Dynamics of a minimal model of interlocked positive and negative feedback loops of transcriptional regulation by cAMP-response element binding proteins.

cAMP-response element binding (CREB) proteins are involved in transcriptional regulation in a number of cellular processes (e.g., neural plasticity and circadian rhythms). The CREB family contains activators and repressors that may interact through positive and negative feedback loops. These loops can be generated by auto- and cross-regulation of expression of CREB proteins, via CRE elements in or near their genes. Experiments suggest that such feedback loops may operate in several systems (e.g., Aplysia and rat). To understand the functional implications of such feedback loops, which are interlocked via cross-regulation of transcription, a minimal model with a positive and negative loop was developed and investigated using bifurcation analysis. Bifurcation analysis revealed diverse nonlinear dynamics (e.g., bistability and oscillations). The stability of steady states or oscillations could be changed by time delays in the synthesis of the activator (CREB1) or the repressor (CREB2). Investigation of stochastic fluctuations due to small numbers of molecules of CREB1 and CREB2 revealed a bimodal distribution of CREB molecules in the bistability region. The robustness of the stable HIGH and LOW states of CREB expression to stochastic noise differs, and a critical number of molecules was required to sustain the HIGH state for days or longer. Increasing positive feedback or decreasing negative feedback also increased the lifetime of the HIGH state, and persistence of this state may correlate with long-term memory formation. A critical number of molecules was also required to sustain robust oscillations of CREB expression. If a steady state was near a deterministic Hopf bifurcation point, stochastic resonance could induce oscillations. This comparative analysis of deterministic and stochastic dynamics not only provides insights into the possible dynamics of CREB regulatory motifs, but also demonstrates a framework for understanding other regulatory processes with similar network architecture.

Dynamics of a minimal model of interlocked positive and negative feedback loops of transcriptional regulation by cAMP-response element binding proteins.

cAMP-response element binding (CREB) proteins are involved in transcriptional regulation in a number of cellular processes (e.g., neural plasticity and circadian rhythms). The CREB family contains activators and repressors that may interact through positive and negative feedback loops. These loops can be generated by auto- and cross-regulation of expression of CREB proteins, via CRE elements in or near their genes. Experiments suggest that such feedback loops may operate in several systems (e.g., Aplysia and rat). To understand the functional implications of such feedback loops, which are interlocked via cross-regulation of transcription, a minimal model with a positive and negative loop was developed and investigated using bifurcation analysis. Bifurcation analysis revealed diverse nonlinear dynamics (e.g., bistability and oscillations). The stability of steady states or oscillations could be changed by time delays in the synthesis of the activator (CREB1) or the repressor (CREB2). Investigation of stochastic fluctuations due to small numbers of molecules of CREB1 and CREB2 revealed a bimodal distribution of CREB molecules in the bistability region. The robustness of the stable HIGH and LOW states of CREB expression to stochastic noise differs, and a critical number of molecules was required to sustain the HIGH state for days or longer. Increasing positive feedback or decreasing negative feedback also increased the lifetime of the HIGH state, and persistence of this state may correlate with long-term memory formation. A critical number of molecules was also required to sustain robust oscillations of CREB expression. If a steady state was near a deterministic Hopf bifurcation point, stochastic resonance could induce oscillations. This comparative analysis of deterministic and stochastic dynamics not only provides insights into the possible dynamics of CREB regulatory motifs, but also demonstrates a framework for understanding other regulatory processes with similar network architecture.

Differential dynamic properties of scleroderma fibroblasts in response to perturbation of environmental stimuli.

Diseases are believed to arise from dysregulation of biological systems (pathways) perturbed by environmental triggers. Biological systems as a whole are not just the sum of their components, rather ever-changing, complex and dynamic systems over time in response to internal and external perturbation. In the past, biologists have mainly focused on studying either functions of isolated genes or steady-states of small biological pathways. However, it is systems dynamics that play an essential role in giving rise to cellular function/dysfunction which cause diseases, such as growth, differentiation, division and apoptosis. Biological phenomena of the entire organism are not only determined by steady-state characteristics of the biological systems, but also by intrinsic dynamic properties of biological systems, including stability, transient-response, and controllability, which determine how the systems maintain their functions and performance under a broad range of random internal and external perturbations. As a proof of principle, we examine signal transduction pathways and genetic regulatory pathways as biological systems. We employ widely used state-space equations in systems science to model biological systems, and use expectation-maximization (EM) algorithms and Kalman filter to estimate the parameters in the models. We apply the developed state-space models to human fibroblasts obtained from the autoimmune fibrosing disease, scleroderma, and then perform dynamic analysis of partial TGF-beta pathway in both normal and scleroderma fibroblasts stimulated by silica. We find that TGF-beta pathway under perturbation of silica shows significant differences in dynamic properties between normal and scleroderma fibroblasts. Our findings may open a new avenue in exploring the functions of cells and mechanism operative in disease development.

Heart rate variability in response to pain stimulus in VLBW infants followed longitudinally during NICU stay.

The objective of this longitudinal study, conducted in a neonatal intensive care unit, was to characterize the response to pain of high-risk very low birth weight infants (<1,500 g) from 23 to 38 weeks post-menstrual age (PMA) by measuring heart rate variability (HRV). Heart period data were recorded before, during, and after a heel lanced or wrist venipunctured blood draw for routine clinical evaluation. Pain response to the blood draw procedure and age-related changes of HRV in low-frequency and high-frequency bands were modeled with linear mixed-effects models. HRV in both bands decreased during pain, followed by a recovery to near-baseline levels. Venipuncture and mechanical ventilation were factors that attenuated the HRV response to pain. HRV at the baseline increased with post-menstrual age but the growth rate of high-frequency power was reduced in mechanically ventilated infants. There was some evidence that low-frequency HRV response to pain improved with advancing PMA.

Laser-induced transient grating analysis of dynamics of interaction between sensory rhodopsin II D75N and the HtrII transducer.

The interaction between sensory rhodopsin II (SRII) and its transducer HtrII was studied by the time-resolved laser-induced transient grating method using the D75N mutant of SRII, which exhibits minimal visible light absorption changes during its photocycle, but mediates normal phototaxis responses. Flash-induced transient absorption spectra of transducer-free D75N and D75N joined to 120 amino-acid residues of the N-terminal part of the SRII transducer protein HtrII (DeltaHtrII) showed only one spectrally distinct K-like intermediate in their photocycles, but the transient grating method resolved four intermediates (K(1)-K(4)) distinct in their volumes. D75N bound to HtrII exhibited one additional slower kinetic species, which persists after complete recovery of the initial state as assessed by absorption changes in the UV-visible region. The kinetics indicate a conformationally changed form of the transducer portion (designated Tr*), which persists after the photoreceptor returns to the unphotolyzed state. The largest conformational change in the DeltaHtrII portion was found to cause a DeltaHtrII-dependent increase in volume rising in 8 micros in the K(4) state and a drastic decrease in the diffusion coefficient (D) of K(4) relatively to those of the unphotolyzed state and Tr*. The magnitude of the decrease in D indicates a large structural change, presumably in the solvent-exposed HAMP domain of DeltaHtrII, where rearrangement of interacting molecules in the solvent would substantially change friction between the protein and the solvent.

Intracranial EEG reveals a time- and frequency-specific role for the right inferior frontal gyrus and primary motor cortex in stopping initiated responses.

Inappropriate response tendencies may be stopped via a specific fronto/basal ganglia/primary motor cortical network. We sought to characterize the functional role of two regions in this putative stopping network, the right inferior frontal gyrus (IFG) and the primary motor cortex (M1), using electocorticography from subdural electrodes in four patients while they performed a stop-signal task. On each trial, a motor response was initiated, and on a minority of trials a stop signal instructed the patient to try to stop the response. For each patient, there was a greater right IFG response in the beta frequency band ( approximately 16 Hz) for successful versus unsuccessful stop trials. This finding adds to evidence for a functional network for stopping because changes in beta frequency activity have also been observed in the basal ganglia in association with behavioral stopping. In addition, the right IFG response occurred 100-250 ms after the stop signal, a time range consistent with a putative inhibitory control process rather than with stop-signal processing or feedback regarding success. A downstream target of inhibitory control is M1. In each patient, there was alpha/beta band desynchronization in M1 for stop trials. However, the degree of desynchronization in M1 was less for successfully than unsuccessfully stopped trials. This reduced desynchronization on successful stop trials could relate to increased GABA inhibition in M1. Together with other findings, the results suggest that behavioral stopping is implemented via synchronized activity in the beta frequency band in a right IFG/basal ganglia network, with downstream effects on M1.

Differential dynamic properties of scleroderma fibroblasts in response to perturbation of environmental stimuli.

Diseases are believed to arise from dysregulation of biological systems (pathways) perturbed by environmental triggers. Biological systems as a whole are not just the sum of their components, rather ever-changing, complex and dynamic systems over time in response to internal and external perturbation. In the past, biologists have mainly focused on studying either functions of isolated genes or steady-states of small biological pathways. However, it is systems dynamics that play an essential role in giving rise to cellular function/dysfunction which cause diseases, such as growth, differentiation, division and apoptosis. Biological phenomena of the entire organism are not only determined by steady-state characteristics of the biological systems, but also by intrinsic dynamic properties of biological systems, including stability, transient-response, and controllability, which determine how the systems maintain their functions and performance under a broad range of random internal and external perturbations. As a proof of principle, we examine signal transduction pathways and genetic regulatory pathways as biological systems. We employ widely used state-space equations in systems science to model biological systems, and use expectation-maximization (EM) algorithms and Kalman filter to estimate the parameters in the models. We apply the developed state-space models to human fibroblasts obtained from the autoimmune fibrosing disease, scleroderma, and then perform dynamic analysis of partial TGF-beta pathway in both normal and scleroderma fibroblasts stimulated by silica. We find that TGF-beta pathway under perturbation of silica shows significant differences in dynamic properties between normal and scleroderma fibroblasts. Our findings may open a new avenue in exploring the functions of cells and mechanism operative in disease development.

Motile response patterns and ultrastructural observations of the isolated outer hair cell

The tonotopic organization of the mammalian cochlea is accompanied by structural gradients which include the somatic lengths of outer hair cells (OHCs). These receptors rest upon the vibrating portion of the basilar membrane and have been reported to exhibit motile responses following chemical and electrical stimulation. These movements were examined in detail in this dissertation. It was found that isolated OHCs cultured in vitro respond to chemical depolarization with slow tonic movements, and to electrical waveforms with bi-directional, frequency following movements extending from DC to at least 10 kHz.^ Slow contractions were also elicited following electrical stimulation, bath incubation in carbachol (a cholinergic agonist), and increases in extracellular K+ concentration as little as 50 mM.^ Isolated OHCs display anatomical features which are remarkable when contrasted with those prepared from intact receptor organs. A complex structure located between the cuticular plate and the nuclear membrane was consistently observed and was examined by serial cross-sections which revealed a network of non-membrane bound densities. This corresponded to a granular complex seen at the light microscope level. The complex was composed of dense regions of organelles, striated structures embedded within the core, and a circumferential network of microtubules residing in the peri-nuclear portion of the cell. In cells which had lost their nuclear attachment to the terminal synaptic body, the granular complex could be made to contract without effecting any change in cellular length, implying that the complex may be the driving force behind certain aspects of the motile response.^ Most cells displayed movements which revealed asymmetries analogous to those reported for OHC receptor potentials in vivo. The contraction phase (for longer cells) was shown to have a small time constant (approximately 400 microseconds) and saturated with limited displacements. The expansion phase had time constants as large as 1.3 milliseconds but yielded displacements as much as 60 percent larger than those seen for contractions.^ Additional waveform characteristics seen in the in vivo response could be emulated either by biasing the cell's resting length with either direct current, triggering contractions via large electrical displacements, or incubation with depolarizing compounds.^ Alternatively, short (20-30 um) cells revealed more linear response characteristics to the probe stimulus. Partial saturation was achieved and revealed a DC component which was opposite in polarity to that seen in longer cells. (Abstract shortened with permission of author.) ^

Comparison of compliance and immune response of drug users on the standard and accelerated hepatitis B vaccine schedules - Houston, Texas

Despite the availability of hepatitis B vaccine for over two decades, drug users and other high-risk adult populations have experienced low vaccine coverage. Poor compliance has limited efforts to reduce transmission of hepatitis B infection in this population. Evidence suggests that immunological response in drug users is impaired compared to the general population, both in terms of lower seroprotection rates and antibodies levels.^ The current study investigated the effectiveness of the multi-dose hepatitis B vaccine and compared the effect of the standard and accelerated vaccine schedules in a not-in-treatment, drug-using adult population in the city of Houston, USA.^ A population of drug-users from two communities in Houston, susceptible to hepatitis B, was sampled by outreach workers and referral methodology. Subjects were randomized either to the standard hepatitis vaccine schedule (0, 1-, 6-month) or to an accelerated schedule (0, 1-, 2-month). Antibody levels were detected through laboratory analyses at various time-points. The participants were followed for two years and seroconversion rates were calculated to determine immune response.^ A four percent difference in the overall compliance rate was observed between the standard (73%) and accelerated schedules (77%). Logistic regression analyses showed that drug users living on the streets were twice as likely to not complete all three vaccine doses (p=0.028), and current speedball use was also associated with non-completion (p=0.002). Completion of all three vaccinations in the multivariate analysis was also correlated with older age. Drug users on the accelerated schedule were 26% more likely to achieve completion, although this factor was marginally significant (p=0.085).^ Cumulative adequate protective response was gained by 65% of the HBV susceptible subgroup by 12-months and was identical for both the standard and accelerated schedules. Excess protective response (>=100 mIU/mL) occurred with greater frequency at the later period for the standard schedule (36% at 12-months compared to 14% at six months), while the greater proportion of excess protective response for the accelerated schedule occurred earlier (34% at 6 months compared to 18% at 12-months). Seroconversion at the adequate protective response level of 10 mIU/mL was reached by the accelerated schedule group at a quicker rate (62% vs. 49%), and with a higher mean titer (104.8 vs. 64.3 mIU/mL), when measured at six months. Multivariate analyses indicated a 63% increased risk of non-response for older age and confirmed the existence of an accelerating decline in immune response to vaccination manifesting after 40 years (p=0.001). Injecting more than daily was also highly associated with the risk of non-response (p=0.016).^ The substantial increase in the seroprotection rate at six months may be worth the trade-off against the faster antibody titer decrease and is recommended for enhancing compliance and seroconversion. Utilization of the accelerated schedule with the primary objective of increasing compliance and seroconversion rates during the six months after the first dose may confer early protective immunity and reduce the HBV vulnerability of drug users who continue, or have recently initiated, increased high risk drug use and sexual behaviors.^

Antibody response to Panton-Valentine Leukocidin, a Staphylococcus aureus virulence factor, in pediatric patients from Texas Children's Hospital in Houston, Texas

Staphylococcus aureus is an important human pathogen of global health significance, whose frequency is increasing and whose persistence and versatility allow it to remain established in communities worldwide. An observed significant increase in infections, particularly in children with no predisposing risk factors or medical conditions, led to an investigation into pediatric humoral immune response to Panton-Valentine Leukocidin (PVL) and to other antigens expressed by S. aureus that represent the important classes of virulence activities. Patients who were diagnosed with staphylococcal infections were enrolled (n=60), and serum samples collected at the time of admission were analyzed using ELISA and Western blot to screen for immune response to the panel of recombinant proteins. The dominant circulating immunoglobulin titers in this pediatric population were primarily IgG, were specific, and were directed against LukF and LukS, while suppression of other important virulence factors in the presence of PVL was suggested. Patients with invasive infections (osteomyelitis, pneumonia or myositis) had higher titers against LukF and LukS compared to patients with non-invasive infections (abscesses, cellulitis or lymphadenitis). In patients with osteomyelitis, antibody responses to LukF and LukS were higher than antibody responses to any other virulence factor examined. This description of immune response to selected virulence factors of S. aureus caused by isolates of the USA300 lineage in children is novel. Antibody titers also correlated with markers of inflammation. The significance of these correlations remains to be understood.^

TRANSMISSION DYNAMICS OF ENTERIC BACTERIA IN DAY-CARE CENTERS, HOUSTON, TEXAS

An investigation was undertaken to evaluate the role of fomites in the transmission of diarrhea in day-care centers (DCC) and to elucidate the paths by which enteric organisms spread within this setting.^ During a nine-month period (December 1980-August 1981) extensive culturing of inanimate objects, as well as children and staff was done routinely each month and again repeated during diarrhea outbreaks. Air was sampled from the classrooms and toilets using a Single-Stage Sieve Sampler (Ross Industries, Midland, VA.). Stool samples were collected from both ill and well children and staff in the affected rooms only during outbreaks. Environmental samples were processed for Shigella, salmonella and fecal coliforms while stools were screened for miscellaneous enteropathogens.^ A total of 11 outbreaks occurred in the 5 DCC during the study period. Enteric pathogens were recovered in 7 (64%) of the outbreaks. Multiple pathogens were identified in 3 outbreaks. The most frequently identified pathogen in stools was Giardia lamblia which was recovered in 5 (45%) of the outbreaks. Ten of the 11 (91%) outbreaks occurred in children less than 12 months of age.^ Environmental microbiology studies together with epidemiologic information revealed that enteric organisms were transmitted from person-to-person. On routine sampling, fecal coliforms were most frequently isolated from tap handles and diaper change areas. Contamination with fetal coliforms was wide-spread during diarrhea outbreaks. Fecal coliforms were recovered with significantly greater frequency from hands, toys and other classroom objects during outbreaks than during non-outbreak period. Salmonella typhimurium was recovered from a table top during an outbreak of Salmonellosis. There was no association between the level of enteric microbial contamination in the toilet areas and the occurrence of outbreaks. No evidence was found to indicate that enteric organisms were spread by the airborne route via aerosols.^ Toys, other classroom objects and contaminated hands probably play a major role in the transmission of enteropathogens during day-care center outbreaks. The presence of many enteric agents in the environment undoubtedly explains the polymicrobial etiology of the day-care center associated diarrhea outbreaks. ^

APPLICATION OF THE BRAINSTEM EVOKED RESPONSE (BER) TO NEUROTOXICITY EVALUATION OF THE AUDITORY SYSTEM: ASSESSMENT OF INORGANIC LEAD-INDUCED OTOTOXICITY

An experimental procedure was developed using the Brainstem Evoked Response (BER) electrophysiological technique to assess the effect of neurotoxic substances on the auditory system. The procedure utilizes Sprague-Dawley albino rats who have had dural electrodes implanted in their skulls, allowing neuroelectric evoked potentials to be recorded from their brainstems. Latency and amplitude parameters derived from the evoked potentials help assess the neuroanatomical integrity of the auditory pathway in the brainstem. Moreover, since frequency-specific auditory stimuli are used to evoke the neural responses, additional audiometric information is obtainable. An investigation on non-exposed control animals shows the BER threshold curve obtained by tests at various frequencies very closely approximates that obtained by behavioral audibility tests. Thus, the BER appears to be a valid measure of both functional and neuroanatomical integrity of the afferent auditory neural pathway.^ To determine the usefulness of the BER technique in neurobehavioral toxicology research, a known neurotoxic agent, Pb, was studied. Female Sprague-Dawley rats were dosed for 45 days with low levels of Pb acetate in their drinking water, after which BER recordings were obtained. The Pb dosages were determined from the findings of an earlier pilot study. One group of 6 rats received normal tap water, one group of 7 rats received a solution of 0.1% Pb, and another group of 7 rats received a solution of 0.2% Pb. After 45 days, the three groups exhibited blood Pb levels of 4.5 (+OR-) 0.43 (mu)g/100 ml, 37.8 (+OR-) 4.8 (mu)g/100 ml and 47.3 (+OR-) 2.7 (mu)g/100 ml, respectively.^ The results of the BER recording indicated evoked response waveform latency abnormalities in both the Pb-treated groups when midrange frequency (8 kHz to 32 kHz) stimuli were used. For the most part, waveform amplitudes did not vary significantly from control values. BER recordings obtained after a 30-day recovery period indicated the effects seen in the 0.1% Pb group had disappeared. However, those anomalies exhibited by the 0.2% Pb group either remained or increased in number. This outcome indicates a longer lasting or possibly irreversible effect on the auditory system from the higher dose of Pb. The auditory pathway effect appears to be in the periphery, at the level of the cochlea or the auditory (VIII) nerve. The results of this research indicate the BER technique is a valuable and sensitive indicator of low-level toxic effects on the auditory system.^

Evidence for the role of agonist binding frequency in receptor -mediated activation of adenylate cyclase

$\beta$-adrenergic receptor-mediated activation of adenylate cyclase exhibits an agonist-specific separation between the dose/response curve (characterized by the EC$\sb{50}$) and the dose/binding curve (characterized by the K$\sb{\rm d}$). Cyclase activity can be near-maximal when receptor occupancy is quite low (EC$\sb{50}$ $\ll$ K$\sb{\rm d}$). This separation between the binding and response curves can be explained by the assumption that the rate of cyclase activation is proportional to the concentration of agonist-bound receptors, since the receptor is mobile and can activate more than one cyclase (the Collision Coupling Model of Tolkovsky and Levitzki). Here it is established that agonist binding frequency plays an additional role in adenylate cyclase activation in S49 murine lymphoma cells. Using epinephrine (EC$\sb{50}$ = 10 nM, K$\sb{\rm d}$ = 2 $\mu$M), the rate of cyclase activation decreased by 80% when a small (1.5%) receptor occupancy was restricted (by addition of the antagonist propranolol) to a small number (1.5%) of receptors rather than being proportionally distributed among the cell's entire population of receptors. Thus adenylate cyclase activity is not proportional to receptor occupancy in all circumstances. Collisions between receptor and cyclase pairs apparently occur a number of times in rapid sequence (an encounter); the high binding frequency of epinephrine ensures that discontiguous regions of the cell surface experience some period of agonist-bound receptor activity per small unit time minimizing "wasted" collisions between activated cyclase and bound receptor within an encounter. A contribution of agonist binding frequency to activation is thus possible when: (1) the mean lifetime of the agonist-receptor complex is shorter than the mean encounter time, and (2) the absolute efficiency (intrinsic ability to promote cyclase activation per collision) of the agonist-receptor complex is high. These conclusions are supported by experiments using agonists of different efficiencies and binding frequencies. These results are formalized in the Encounter Coupling Model of adenylate cyclase activation, which takes into explicit account the agonist binding frequency, agonist affinity for the $\beta$-adrenergic receptor, agonist efficiency, encounter frequency and the encounter time between receptor and cyclase. ^