The functional relationship between yawning and vigilance

BACKGROUND: Although yawning is a ubiquitous and phylogenetically old phenomenon, its origin and purpose remain unclear. The study aimed at testing the widely held hypothesis that yawning is triggered by drowsiness and brings about a reversal or suspension of the process of falling asleep. METHODS: Subjects complaining of excessive sleepiness were spontaneously yawning while trying to stay awake in a quiet and darkened room. Changes in their electroencephalogram (EEG) and heart rate variability (HRV) associated with yawning were compared to changes associated with isolated voluntary body movements. Special care was taken to remove eye blink- and movement-artefacts from the recorded signals. RESULTS: Yawns were preceded and followed by a significantly greater delta activity in EEG than movements (p< or =0.008). After yawning, alpha rhythms were attenuated, decelerated, and shifted towards central brain regions (p< or =0.01), whereas after movements, they were attenuated and accelerated (p<0.02). A significant transient increase of HRV occurred after the onset of yawning and movements, which was followed by a significant slow decrease peaking 17s after onset (p<0.0001). No difference in HRV changes was found between yawns and movements. CONCLUSIONS: Yawning occurred during periods with increased drowsiness and sleep pressure, but was not followed by a measurable increase of the arousal level of the brain. It was neither triggered nor followed by a specific autonomic activation. Our results therefore confirm that yawns occur due to sleepiness, but do not provide evidence for an arousing effect of yawning.

Indocyanine green plasma disappearance rate during the anhepatic phase of orthotopic liver transplantation

Non-invasive pulse spectrophotometry to measure indocyanine green (ICG) elimination correlates well with the conventional invasive ICG clearance test. Nevertheless, the precision of this method remains unclear for any application, including small-for-size liver remnants. We therefore measured ICG plasma disappearance rate (PDR) during the anhepatic phase of orthotopic liver transplantation using pulse spectrophotometry. Measurements were done in 24 patients. The median PDR after exclusion of two outliers and two patients with inconstant signal was 1.55%/min (95% confidence interval [CI]=0.8-2.2). No correlation with patient age, gender, body mass, blood loss, administration of fresh frozen plasma, norepinephrine dose, postoperative albumin (serum), or difference in pre and post transplant body weight was detected. In conclusion, we found an ICG-PDR different from zero in the anhepatic phase, an overestimation that may arise in particular from a redistribution into the interstitial space. If ICG pulse spectrophotometry is used to measure functional hepatic reserve, the verified average difference from zero (1.55%/min) determined in our study needs to be taken into account.

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.

Method of Detecting System Function by Measuring Frequency Response

Methods of rapidly measuring an impedance spectrum of an energy storage device in-situ over a limited number of logarithmically distributed frequencies are described. An energy storage device is excited with a known input signal, and aresponse is measured to ascertain the impedance spectrum. An excitation signal is a limited time duration sum-of-sines consisting of a select number offrequencies. In one embodiment, magnitude and phase of each frequency ofinterest within the sum-of-sines is identified when the selected frequencies and sample rate are logarithmic integer steps greater than two. This technique requires a measurement with a duration of one period of the lowest frequency. In another embodiment, where selected frequencies are distributed in octave steps, the impedance spectrum can be determined using a captured time record that is reduced to a half-period of the lowest frequency.

T1(Gd) gives comparable information as Delta T1 relaxation rate in dGEMRIC evaluation of cartilage repair tissue

OBJECTIVES: To evaluate the relationship between T1 after intravenous contrast administration (T1Gd) and Delta relaxation rate (DeltaR1) = (1/T1(Gd) - 1/T1o) in the delayed Gadolinium-Enhanced MRI of cartilage (dGEMRIC) evaluation of cartilage repair tissue. MATERIALS AND METHODS: Thirty single MR examinations from 30 patients after matrix-associated autologous chondrocyte transplantations of the knee joint with different postoperative intervals were examined using an 8-channel knee-coil at 3T. T1 mapping using a 3D GRE sequence with a 35/10 degrees flip angle excitation pulse combination was performed before and after contrast administration (dGEMRIC technique). T1 postcontrast (T1(Gd)) and the DeltaR1 (relative index of pre- and postcontrast R1 value) were calculated for repair tissue and the weight-bearing normal appearing control cartilage. For evaluation of the different postoperative intervals, MR exams were subdivided into 3 groups (up to 12 months, 12-24 months, more than 24 months). For statistical analysis Spearman correlation coefficients were calculated. RESULTS: The mean value for T1 postcontrast was 427 +/- 159 ms, for DeltaR1 1.85 +/- 1.0; in reference cartilage 636 +/- 181 ms for T1 postcontrast and 0.83 +/- 0.5 for DeltaR1.The correlation coefficients were highly significant between T1 (Gd) and DeltaR1 for repair tissue (0.969) as well as normal reference cartilage (0.928) in total, and for the reparative cartilage in the early, middle postoperative, and late postoperative interval after surgery (R values: -0.986, -0.970, and -0.978, respectively). Using either T1(Gd) or DeltaR1, the 2 metrics resulted in similar conclusions regarding the time course of change of repair tissue and control tissue, namely that highly significant (P > 0.01) differences between cartilage repair tissue and reference cartilage were found for all follow-up groups. Additionally, for both metrics highly significant differences (P < 0.01) between early follow up and the 2 later postoperative groups for cartilage repair tissue were found. No statistical differences were found between the 2 later follow-up groups of reparative cartilage either for T1 (Gd) or DeltaR1. CONCLUSION: The high correlation between T1 (Gd) and DeltaR1 and the comparable conclusions reached utilizing metric implies that T1 mapping before intravenous administration of MR contrast agent is not necessary for the evaluation of repair tissue. This will help to reduce costs, inconvenience for the patients, simplifies the examination procedure, and makes dGEMRIC more attractive for follow-up of patients after cartilage repair surgeries.

Unanticipated structural and functional properties of delta-subunit-containing GABAA receptors

GABA(A) receptors mediate inhibitory neurotransmission in the mammalian brain via synaptic and extrasynaptic receptors. The delta (delta)-subunit-containing receptors are expressed exclusively extra-synaptically and mediate tonic inhibition. In the present study, we were interested in determining the architecture of receptors containing the delta-subunit. To investigate this, we predefined the subunit arrangement by concatenation. We prepared five dual and three triple concatenated subunit constructs. These concatenated dual and triple constructs were used to predefine nine different GABA(A) receptor pentamers. These pentamers composed of alpha(1)-, beta(3)-, and delta-subunits were expressed in Xenopus oocytes and maximal currents elicited in response to 1 mm GABA were determined in the presence and absence of THDOC (3alpha, 21-dihydroxy-5alpha-pregnane-20-one). beta(3)-alpha(1)-delta/alpha(1)-beta(3) and beta(3)-alpha(1)-delta/beta(3)-alpha(1) resulted in the expression of large currents in response to GABA. Interestingly, the presence of the neurosteroid THDOC uncovered alpha(1)-beta(3)-alpha(1)/beta(3)-delta receptors, additionally. The functional receptors were characterized in detail using the agonist GABA, THDOC, Zn(2+), and ethanol and their properties were compared with those of non-concatenated alpha(1)beta(3) and alpha(1)beta(3)delta receptors. Each concatenated receptor isoform displayed a specific set of properties, but none of them responded to 30 mm ethanol. We conclude from the investigated receptors that delta can assume multiple positions in the receptor pentamer. The GABA dose-response properties of alpha(1)-beta(3)-alpha(1)/beta(3)-delta and beta(3)-alpha(1)-delta/alpha(1)-beta(3) match most closely the properties of non-concatenated alpha(1)beta(3)delta receptors. Furthermore, we show that the delta-subunit can contribute to the formation of an agonist site in alpha(1)-beta(3)-alpha(1)/beta(3)-delta receptors.

Dose–effect relationship in routine outpatient psychotherapy: Does treatment duration matter?

Objective: There is an ongoing debate concerning how outcome variables change during the course of psychotherapy. We compared the dose–effect model, which posits diminishing effects of additional sessions in later treatment phases, against a model that assumes a linear and steady treatment progress through termination. Method: Session-by-session outcome data of 6,375 outpatients were analyzed, and participants were categorized according to treatment length. Linear and log-linear (i.e., negatively accelerating) latent growth curve models (LGCMs) were estimated and compared for different treatment length categories. Results: When comparing the fit of the various models, the log-linear LGCMs assuming negatively accelerating treatment progress consistently outperformed the linear models irre- spective of treatment duration. The rate of change was found to be inversely related to the length of treatment. Conclusion: As proposed by the dose–effect model, the expected course of improvement in psychotherapy appears to follow a negatively accelerated pattern of change, irrespective of the duration of the treatment. However, our results also suggest that the rate of change is not constant across various treatment lengths. As proposed by the “good enough level” model, longer treatments are associated with less rapid rates of change.

Molecular evolution of primate immunodeficiency viruses and hepatitis delta virus

Primate immunodeficiency viruses, or lentiviruses (HIV-1, HIV-2, and SIV), and hepatitis delta virus (HDV) are RNA viruses characterized by rapid evolution. Infection by primate immunodeficiency viruses usually results in the development of acquired immunodeficiency syndrome (AIDS) in humans and AIDS-like illnesses in Asian macaques. Similarly, hepatitis delta virus infection causes hepatitis and liver cancer in humans. These viruses are heterogeneous within an infected patient and among individuals. Substitution rates in the virus genomes are high and vary in different lineages and among sites. Methods of phylogenetic analysis were applied to study the evolution of primate lentiviruses and the hepatitis delta virus. The following results have been obtained: (1) The substitution rate varies among sites of primate lentivirus genes according to the two parameter gamma distribution, with the shape parameter $\alpha$ being close to 1. (2) Primate immunodeficiency viruses fall into species-specific lineages. Therefore, viral transmissions across primate species are not as frequent as suggested by previous authors. (3) Primate lentiviruses have acquired or lost their pathogenicity several times in the course of evolution. (4) Evidence was provided for multiple infections of a North American patient by distinct HIV-1 strains of the B subtype. (5) Computer simulations indicate that the probability of committing an error in testing HIV transmission depends on the number of virus sequences and their length, the divergence times among sequences, and the model of nucleotide substitution. (6) For future investigations of HIV-1 transmissions, using longer virus sequences and avoiding the use of distant outgroups is recommended. (7) Hepatitis delta virus strains are usually related according to the geographic region of isolation. (8) Evolution of HDV is characterized by the rate of synonymous substitution being lower than the nonsynonymous substitution rate and the rate of evolution of the noncoding region. (9) There is a strong preference for G and C nucleotides at the third codon positions of the HDV coding region. ^

Functional studies on DNA double-strand break repair proteins (MmRad51, Ku80) in mice

RecA in Escherichia coli and it's homologue, ScRad51 in Saccharomyces cerevisiae, play important roles in recombinational repair. ScRad51 homologues have been discovered in a wide range of organisms including Schizosaccharomyces pombe, lily, chicken, mouse and human. To date there is no direct evidence to describe that mouse Rad51(MmRad51) is involved in DNA double-strand break repair. In order to elucidate the role of MmRad51 in vivo, it was mutated by the embryonic stem (ES) cell/gene targeting technology in mice. The mutant embryos arrested in development shortly after implantation. There was a decrease in cell proliferation followed by programmed cell death, and trophectoderm-derived cells were sensitive to $\gamma$-radiation. Severe chromosome loss was observed in most mitotically dividing cells. The mutant embryos lived longer and developed further in a p53 mutant background; however, double-mutant embryonic fibroblasts failed to proliferate in tissue culture, reflecting the embryos limited life span. Based on these data, MmRad51 repairs DNA damage induced by $\gamma$-radiation, is needed to maintain euplody, and plays an important role in proliferating cells.^ Ku is a heterodimer of 70 and 80 kDs subunit, which binds to DNA ends and other altered DNA structures such as hairpins, nicks, and gaps. In addition, Ku is required for DNA-PK activity through a direct association. Although the biochemical properties of Ku and DNA-PKcs have been characterized in cells, their physiological functions are not clear. In order to understand the function of Ku in vivo, we generated mice homozygous for a mutation of the Ku80 gene. Ku80-deficient mice, like scid mice, showed severe immunodeficiency due to a impairment of V(D)J recombination. Mutant mice were semiviable and runted, cells derived from mutant embryos displayed hypersensitivity to $\gamma$-radiation, a decreased growth rate, a slow entry into S phase, altered colony size distributions, and a short life span. Based on these results, mutant cells and mice appeared to prematurely age. ^

Coherence and phase locking in the scalp EEG and between LORETA model sources, and microstates as putative mechanisms of brain temporo-spatial functional organization

Brain electric mechanisms of temporary, functional binding between brain regions are studied using computation of scalp EEG coherence and phase locking, sensitive to time differences of few milliseconds. However, such results if computed from scalp data are ambiguous since electric sources are spatially oriented. Non-ambiguous results can be obtained using calculated time series of strength of intracerebral model sources. This is illustrated applying LORETA modeling to EEG during resting and meditation. During meditation, time series of LORETA model sources revealed a tendency to decreased left-right intracerebral coherence in the delta band, and to increased anterior-posterior intracerebral coherence in the theta band. An alternate conceptualization of functional binding is based on the observation that brain electric activity is discontinuous, i.e., that it occurs in chunks of up to about 100 ms duration that are detectable as quasi-stable scalp field configurations of brain electric activity, called microstates. Their functional significance is illustrated in spontaneous and event-related paradigms, where microstates associated with imagery- versus abstract-type mentation, or while reading positive versus negative emotion words showed clearly different regions of cortical activation in LORETA tomography. These data support the concept that complete brain functions of higher order such as a momentary thought might be incorporated in temporal chunks of processing in the range of tens to about 100 ms as quasi-stable brain states; during these time windows, subprocesses would be accepted as members of the ongoing chunk of processing.

Rapid method to express and purify human membrane protein using the Xenopus oocyte system for functional and low-resolution structural analysis.

Progress toward elucidating the 3D structures of eukaryotic membrane proteins has been hampered by the lack of appropriate expression systems. Recent work using the Xenopus oocyte as a novel expression system for structural analysis demonstrates the capability of providing not only the significant amount of protein yields required for structural work but also the expression of eukaryotic membrane proteins in a more native and functional conformation. There is a long history using the oocyte expression system as an efficient tool for membrane transporter and channel expression in direct functional analysis, but improvements in robotic injection systems and protein yield optimization allow the rapid scalability of expressed proteins to be purified and characterized in physiologically relevant structural states. Traditional overexpression systems (yeast, bacteria, and insect cells) by comparison require chaotropic conditions over several steps for extraction, solubilization, and purification. By contrast, overexpressing within the oocyte system for subsequent negative-staining transmission electron microscopy studies provides a single system that can functionally assess and purify eukaryotic membrane proteins in fewer steps maintaining the physiological properties of the membrane protein.

Intrinsic characteristics of the proton pump in the luminal membrane of a tight urinary epithelium. The relation between transport rate and delta mu H.

A number of tight urinary epithelia, as exemplified by the turtle bladder, acidify the luminal solution by active transport of H+ across the luminal cell membrane. The rate of active H+ transport (JH) decreases as the electrochemical potential difference for H+ [delta mu H = mu H(lumen) - mu H(serosa)] across the epithelium is increased. The luminal cell membrane has a low permeability for H+ equivalents and a high electrical resistance compared with the basolateral cell membrane. Changes in JH thus reflect changes in active H+ transport across the luminal membrane. To examine the control of JH by delta mu H in the turtle bladder, transepithelial electrical potential differences (delta psi) were imposed at constant acid-base conditions or the luminal pH was varied at delta psi = 0 and constant serosal PCO2 and pH. When the luminal compartment was acidified from pH 7 to 4 or was made electrically positive, JH decreased as a linear function of delta mu H as previously described. When the luminal compartment was made alkaline from pH 7 to 9 or was made electrically negative, JH reached a maximal value, which was the same whether the delta mu H was imposed as a delta pH or a delta psi. The nonlinear JH vs. delta mu H relation does not result from changes in the number of pumps in the luminal membrane or from changes in the intracellular pH, but is a characteristic of the H+ pumps themselves. We propose a general scheme, which, because of its structural features, can account for the nonlinearity of the JH vs. delta mu H relations and, more specifically, for the kinetic equivalence of the effects of the chemical and electrical components of delta mu H. According to this model, the pump complex consists of two components: a catalytic unit at the cytoplasmic side of the luminal membrane, which mediates the ATP-driven H+ translocation, and a transmembrane channel, which mediates the transfer of H+ from the catalytic unit to the luminal solution. These two components may be linked through a buffer compartment for H+ (an antechamber).

Bayesian Adaptive Designs for Early Phase Clinical Trials

My dissertation focuses mainly on Bayesian adaptive designs for phase I and phase II clinical trials. It includes three specific topics: (1) proposing a novel two-dimensional dose-finding algorithm for biological agents, (2) developing Bayesian adaptive screening designs to provide more efficient and ethical clinical trials, and (3) incorporating missing late-onset responses to make an early stopping decision. Treating patients with novel biological agents is becoming a leading trend in oncology. Unlike cytotoxic agents, for which toxicity and efficacy monotonically increase with dose, biological agents may exhibit non-monotonic patterns in their dose-response relationships. Using a trial with two biological agents as an example, we propose a phase I/II trial design to identify the biologically optimal dose combination (BODC), which is defined as the dose combination of the two agents with the highest efficacy and tolerable toxicity. A change-point model is used to reflect the fact that the dose-toxicity surface of the combinational agents may plateau at higher dose levels, and a flexible logistic model is proposed to accommodate the possible non-monotonic pattern for the dose-efficacy relationship. During the trial, we continuously update the posterior estimates of toxicity and efficacy and assign patients to the most appropriate dose combination. We propose a novel dose-finding algorithm to encourage sufficient exploration of untried dose combinations in the two-dimensional space. Extensive simulation studies show that the proposed design has desirable operating characteristics in identifying the BODC under various patterns of dose-toxicity and dose-efficacy relationships. Trials of combination therapies for the treatment of cancer are playing an increasingly important role in the battle against this disease. To more efficiently handle the large number of combination therapies that must be tested, we propose a novel Bayesian phase II adaptive screening design to simultaneously select among possible treatment combinations involving multiple agents. Our design is based on formulating the selection procedure as a Bayesian hypothesis testing problem in which the superiority of each treatment combination is equated to a single hypothesis. During the trial conduct, we use the current values of the posterior probabilities of all hypotheses to adaptively allocate patients to treatment combinations. Simulation studies show that the proposed design substantially outperforms the conventional multi-arm balanced factorial trial design. The proposed design yields a significantly higher probability for selecting the best treatment while at the same time allocating substantially more patients to efficacious treatments. The proposed design is most appropriate for the trials combining multiple agents and screening out the efficacious combination to be further investigated. The proposed Bayesian adaptive phase II screening design substantially outperformed the conventional complete factorial design. Our design allocates more patients to better treatments while at the same time providing higher power to identify the best treatment at the end of the trial. Phase II trial studies usually are single-arm trials which are conducted to test the efficacy of experimental agents and decide whether agents are promising to be sent to phase III trials. Interim monitoring is employed to stop the trial early for futility to avoid assigning unacceptable number of patients to inferior treatments. We propose a Bayesian single-arm phase II design with continuous monitoring for estimating the response rate of the experimental drug. To address the issue of late-onset responses, we use a piece-wise exponential model to estimate the hazard function of time to response data and handle the missing responses using the multiple imputation approach. We evaluate the operating characteristics of the proposed method through extensive simulation studies. We show that the proposed method reduces the total length of the trial duration and yields desirable operating characteristics for different physician-specified lower bounds of response rate with different true response rates.