Diversity of neuronal connexins in the mammalian retina

Many neurons in the mammalian retina are electrically coupled by intercellular channels or gap junctions, which are assembled from a family of proteins called connexins. Numerous studies indicate that gap junctions differ in properties such as conductance and tracer permeability. For example, A-type horizontal cell gap junctions are permeable to Lucifer Yellow, but B-type horizontal cell gap junctions are not. This suggests the two cell types express different connexins. My hypothesis is that multiple neuronal connexins are expressed in the mammalian retina in a cell type specific manner. Immunohistochemical techniques and confocal microscopy were used to localize certain connexins within well-defined neuronal circuits. The results of this study can be summarized as follows: AII amacrine cells, which receive direct input from rod bipolar cells, are well-coupled to neighboring AIIs. In addition, AII amacrine cells also form gap junctions with ON cone bipolar cells. This is a complex heterocellular network. In both rabbit and primate retina, connexin36 occurs at dendritic crossings in the AII matrix as well as between AIIs and ON cone bipolar cells. Coupling in the AII network is thought to reduce noise in the rod pathway while AII/bipolar gap junctions are required for the transmission of rod signals to ON ganglion cells. In the outer plexiform layer, connexin36 forms gap junctions between cones and between rods and cones via cone telodendria. Cone to cone coupling is thought to reduce noise and is partly color selective. Rod to cone coupling forms an alternative rod pathway thought to operate at intermediate light intensity. A-type horizontal cells in the rabbit retina are strongly coupled via massive low resistance gap junctions composed from Cx50. Coupling dramatically extends the receptive field of horizontal cells and the modulation of coupling is thought to change the strength of the feedback signal from horizontal cells to cones. Finally, there are other coupled networks, such as B-type horizontal cells and S1/S2 amacrine cells, which do not use either connexin36 or Cx50. These results confirm the hypothesis that multiple neuronal connexins are expressed in the mammalian retina and these connexins are localized to particular retinal circuits. ^

Computer simulation of glioma growth and morphology.

Despite major advances in the study of glioma, the quantitative links between intra-tumor molecular/cellular properties, clinically observable properties such as morphology, and critical tumor behaviors such as growth and invasiveness remain unclear, hampering more effective coupling of tumor physical characteristics with implications for prognosis and therapy. Although molecular biology, histopathology, and radiological imaging are employed in this endeavor, studies are severely challenged by the multitude of different physical scales involved in tumor growth, i.e., from molecular nanoscale to cell microscale and finally to tissue centimeter scale. Consequently, it is often difficult to determine the underlying dynamics across dimensions. New techniques are needed to tackle these issues. Here, we address this multi-scalar problem by employing a novel predictive three-dimensional mathematical and computational model based on first-principle equations (conservation laws of physics) that describe mathematically the diffusion of cell substrates and other processes determining tumor mass growth and invasion. The model uses conserved variables to represent known determinants of glioma behavior, e.g., cell density and oxygen concentration, as well as biological functional relationships and parameters linking phenomena at different scales whose specific forms and values are hypothesized and calculated based on in vitro and in vivo experiments and from histopathology of tissue specimens from human gliomas. This model enables correlation of glioma morphology to tumor growth by quantifying interdependence of tumor mass on the microenvironment (e.g., hypoxia, tissue disruption) and on the cellular phenotypes (e.g., mitosis and apoptosis rates, cell adhesion strength). Once functional relationships between variables and associated parameter values have been informed, e.g., from histopathology or intra-operative analysis, this model can be used for disease diagnosis/prognosis, hypothesis testing, and to guide surgery and therapy. In particular, this tool identifies and quantifies the effects of vascularization and other cell-scale glioma morphological characteristics as predictors of tumor-scale growth and invasion.

Computer simulation of glioma growth and morphology.

Despite major advances in the study of glioma, the quantitative links between intra-tumor molecular/cellular properties, clinically observable properties such as morphology, and critical tumor behaviors such as growth and invasiveness remain unclear, hampering more effective coupling of tumor physical characteristics with implications for prognosis and therapy. Although molecular biology, histopathology, and radiological imaging are employed in this endeavor, studies are severely challenged by the multitude of different physical scales involved in tumor growth, i.e., from molecular nanoscale to cell microscale and finally to tissue centimeter scale. Consequently, it is often difficult to determine the underlying dynamics across dimensions. New techniques are needed to tackle these issues. Here, we address this multi-scalar problem by employing a novel predictive three-dimensional mathematical and computational model based on first-principle equations (conservation laws of physics) that describe mathematically the diffusion of cell substrates and other processes determining tumor mass growth and invasion. The model uses conserved variables to represent known determinants of glioma behavior, e.g., cell density and oxygen concentration, as well as biological functional relationships and parameters linking phenomena at different scales whose specific forms and values are hypothesized and calculated based on in vitro and in vivo experiments and from histopathology of tissue specimens from human gliomas. This model enables correlation of glioma morphology to tumor growth by quantifying interdependence of tumor mass on the microenvironment (e.g., hypoxia, tissue disruption) and on the cellular phenotypes (e.g., mitosis and apoptosis rates, cell adhesion strength). Once functional relationships between variables and associated parameter values have been informed, e.g., from histopathology or intra-operative analysis, this model can be used for disease diagnosis/prognosis, hypothesis testing, and to guide surgery and therapy. In particular, this tool identifies and quantifies the effects of vascularization and other cell-scale glioma morphological characteristics as predictors of tumor-scale growth and invasion.

Cellular dynamic simulator: an event driven molecular simulation environment for cellular physiology.

In this paper, we present the Cellular Dynamic Simulator (CDS) for simulating diffusion and chemical reactions within crowded molecular environments. CDS is based on a novel event driven algorithm specifically designed for precise calculation of the timing of collisions, reactions and other events for each individual molecule in the environment. Generic mesh based compartments allow the creation / importation of very simple or detailed cellular structures that exist in a 3D environment. Multiple levels of compartments and static obstacles can be used to create a dense environment to mimic cellular boundaries and the intracellular space. The CDS algorithm takes into account volume exclusion and molecular crowding that may impact signaling cascades in small sub-cellular compartments such as dendritic spines. With the CDS, we can simulate simple enzyme reactions; aggregation, channel transport, as well as highly complicated chemical reaction networks of both freely diffusing and membrane bound multi-protein complexes. Components of the CDS are generally defined such that the simulator can be applied to a wide range of environments in terms of scale and level of detail. Through an initialization GUI, a simple simulation environment can be created and populated within minutes yet is powerful enough to design complex 3D cellular architecture. The initialization tool allows visual confirmation of the environment construction prior to execution by the simulator. This paper describes the CDS algorithm, design implementation, and provides an overview of the types of features available and the utility of those features are highlighted in demonstrations.

MODELLING β2AR REGULATION

The β2 adrenergic receptor (β2AR) regulates smooth muscle relaxation in the vasculature and airways. Long- and Short-acting β-agonists (LABAs/SABAs) are widely used in treatment of chronic obstructive pulmonary disorder (COPD) and asthma. Despite their widespread clinical use we do not understand well the dominant β2AR regulatory pathways that are stimulated during therapy and bring about tachyphylaxis, which is the loss of drug effects. Thus, an understanding of how the β2AR responds to various β-agonists is crucial to their rational use. Towards that end we have developed deterministic models that explore the mechanism of drug- induced β2AR regulation. These mathematical models can be classified into three classes; (i) Six quantitative models of SABA-induced G protein coupled receptor kinase (GRK)-mediated β2AR regulation; (ii) Three phenomenological models of salmeterol (a LABA)-induced GRK-mediated β2AR regulation; and (iii) One semi-quantitative, unified model of SABA-induced GRK-, protein kinase A (PKA)-, and phosphodiesterase (PDE)-mediated regulation of β2AR signalling. The various models were constrained with all or some of the following experimental data; (i) GRK-mediated β2AR phosphorylation in response to various LABAs/SABAs; (ii) dephosphorylation of the GRK site on the β2AR; (iii) β2AR internalisation; (iv) β2AR recycling; (v) β2AR desensitisation; (vi) β2AR resensitisation; (vii) PKA-mediated β2AR phosphorylation in response to a SABA; and (viii) LABA/SABA induced cAMP profile ± PDE inhibitors. The models of GRK-mediated β2AR regulation show that plasma membrane dephosphorylation and recycling of the phosphorylated β2AR are required to reconcile with the measured dephosphorylation kinetics. We further used a consensus model to predict the consequences of rapid pulsatile agonist stimulation and found that although resensitisation was rapid, the β2AR system retained the memory of prior stimuli and desensitised much more rapidly and strongly in response to subsequent stimuli. This could explain tachyphylaxis of SABAs over repeated use in rescue therapy of asthma patients. The LABA models show that the long action of salmeterol can be explained due to decreased stability of the arrestin/β2AR/salmeterol complex. This could explain long action of β-agonists used in maintenance therapy of asthma patients. Our consensus model of PKA/PDE/GRK-mediated β2AR regulation is being used to identify the dominant β2AR desensitisation pathways under different therapeutic regimens in human airway cells. In summary our models represent a significant advance towards understanding agonist-specific β2AR regulation that will aid in a more rational use of the β2AR agonists in the treatment of asthma.

Cerebellar mechanisms for motor learning: Testing predictions from a large-scale computer simulation

The cerebellum is the major brain structure that contributes to our ability to improve movements through learning and experience. We have combined computer simulations with behavioral and lesion studies to investigate how modification of synaptic strength at two different sites within the cerebellum contributes to a simple form of motor learning—Pavlovian conditioning of the eyelid response. These studies are based on the wealth of knowledge about the intrinsic circuitry and physiology of the cerebellum and the straightforward manner in which this circuitry is engaged during eyelid conditioning. Thus, our simulations are constrained by the well-characterized synaptic organization of the cerebellum and further, the activity of cerebellar inputs during simulated eyelid conditioning is based on existing recording data. These simulations have allowed us to make two important predictions regarding the mechanisms underlying cerebellar function, which we have tested and confirmed with behavioral studies. The first prediction describes the mechanisms by which one of the sites of synaptic modification, the granule to Purkinje cell synapses (gr → Pkj) of the cerebellar cortex, could generate two time-dependent properties of eyelid conditioning—response timing and the ISI function. An empirical test of this prediction using small, electrolytic lesions of the cerebellar cortex revealed the pattern of results predicted by the simulations. The second prediction made by the simulations is that modification of synaptic strength at the other site of plasticity, the mossy fiber to deep nuclei synapses (mf → nuc), is under the control of Purkinje cell activity. The analysis predicts that this property should confer mf → nuc synapses with resistance to extinction. Thus, while extinction processes erase plasticity at the first site, residual plasticity at mf → nuc synapses remains. The residual plasticity at the mf → nuc site confers the cerebellum with the capability for rapid relearning long after the learned behavior has been extinguished. We confirmed this prediction using a lesion technique that reversibly disconnected the cerebellar cortex at various stages during extinction and reacquisition of eyelid responses. The results of these studies represent significant progress toward a complete understanding of how the cerebellum contributes to motor learning. ^

Functional architecture of mammalian horizontal cells

In the rabbit retina, there are two kinds of horizontal cells (HCs). The A-type HC is a large axonless cell which contacts cones exclusively. The B-type HC is an axon bearing cell. While the somatic dendrites of B-type HCs also contact cones, the axon expands into an elaborately branched structure, the axon terminal (AT), which contacts a large number of rods. It is difficult to label the different HCs selectively by immunochemical methods. Therefore, we developed dye injection methods to label each type of HC. Then it was possible, (1) to describe the detailed structure of the AT (2) to identify the glutamate receptors mediating cone input to A and B-type HCs and rod input to ATs and (3) to test the hypothesis that the B-type HCs are coupled via Cx57 gap junctions. ^ To obtain well filled examples of single HCs, it was necessary to block gap junction coupling to stop the spread of Neurobiotin through the network. We used dye coupling in A-type HCs to screen a series of potential gap junction antagonists. One of these compounds, meclofenamic acid (MFA), was potent, water soluble and easily reversible. This compound may be a useful tool to manipulate gap junction coupling. ^ In the presence of MFA, Neurobiotin passed down the axon of B-type HCs to reveal the detailed structure of the AT. We observed that only one AT ending entered each rod spherule invagination. This observation was confirmed by calculation and two dye injections. ^ Glutamate is the neurotransmitter used by both rods and cones. AMPA receptors were colocalized with the dendrites of A and B-type HCs at each cone pedicle. In addition, AMPA receptors were located on the AT ending at each rod spherule. Thus rod and cone input to HCs is mediated by AMPA receptors. ^ A-type and B-type HCs may express different connexins because they have different dye-coupling properties. Recently, we found that connexin50 (Cx50) is expressed by A-type HCs. B-type HCs and B-type ATs are also independently coupled. Cx57 was expressed in the OPL and double label studies showed that Cx 57 was colocalized with the AT matrix but not with the somatic dendrites of B-type HCs. ^ In summary, we have identified a useful gap junction antagonist, MFA. There is one AT ending at each rod spherule, rods inputs to ATs is mediated by AMPA receptors and coupling in the AT matrix is mediated by Cx57. This confirms that HCs with different properties use distinct connexins. The properties of ATs described in this research are consistent. The connections and properties reported here suggest that ATs functions as rod HCs and provide a negative feedback signal to rods. ^

Outcomes of tuberculosis positive non-injection drug users in Harris County from October 1995 through September 2004

Background. The population-based Houston Tuberculosis Initiative (HTI) study has enrolled and gathered demographic, social, behavioral, and disease related data on more than 80% of all reported Mycobacterium Tuberculosis (MTB) cases and 90% of all culture positive patients in Houston/Harris County over a 9 year period (from October 1995-September 2004). During this time period 33% (n=1210) of HTI MTB cases have reported a history of drug use. Of those MTB cases reporting a history of drug use, a majority of them (73.6%), are non-injection drug users (NIDUs). ^ Other than HIV, drug use is the single most important risk factor for progression from latent to infectious tuberculosis (TB). In addition, drug use is associated with increased transmission of active TB, as seen by the increased number of clonally related strains or clusters (see definition on page 30) found in this population. The deregulatory effects of drug use on immune function are well documented. Associations between drug use and increased morbidity have been reported since the late 1970's. However, limited research focused on the immunological consequence of non-injection drug use and its relation to tuberculosis infection among TB patients is available. ^ Methods. TB transmission patterns, symptoms, and prevalence of co-morbidities were a focus of this project. Smoking is known to suppress Nitric Oxide (NO) production and interfere with immune function. In order to limit any possible confounding due to smoking two separate analyses were done. Non-injection drug user smokers (NIDU-S) were compared to non-drug user smokers (NDU-S) and non-injection drug user non-smokers (NIDU-NS) were compared to non-drug user non-smokers (NDU-NS) individually. Specifically proportions, chi-square p-values, and (where appropriate) odds ratios with 95% confidence intervals were calculated to assess characteristics and potential associations of co-morbidities and symptoms of TB among NIDUs HTI TB cases. ^ Results. Significant differences in demographic characteristics and risk factors were found. In addition drug users were found to have a decreased risk for cancer, diabetes mellitus, and chronic pulmonary disease. They were at increased risk of having HIV/AIDS diagnosis, liver disease, and trauma related morbidities. Drug users were more likely to have pulmonary TB disease, and a significantly increased amount of clonally related strains of TB or "clusters" were seen in both smokers and non-smoker drug users when compared to their non-drug user counterparts. Drug users are more likely to belong to print groups (clonally related TB strains with matching spoligotypes) including print one and print three and the Beijing family group, s1. Drug users were found to be no more likely to experience drug resistance to TB therapy and were likely to be cured of disease upon completion of therapy. ^ Conclusion. Drug users demographic and behavioral risk factors put them at an increased risk contracting and spreading TB disease throughout the community. Their increased levels of clustering are evidence of recent transmission and the significance of certain print groups among this population indicate the transmission is from within the social family. For these reasons a focus on this "at risk population" is critical to the success of future public health interventions. Successful completion of directly observed therapy (DOT), the tracking of TB outbreaks and incidence through molecular characterization, and increased diagnostic strategies have led to the stabilization of TB incidence in Houston, Harris County over the past 9 years and proven that the Houston Tuberculosis Initiative has played a critical role in the control and prevention of TB transmission. ^

Ligand-stimulated internalization of a Dictyostelium discoideum G protein-coupled cAMP receptor

The social amoeba, Dictyostelium discoideum, undergoes a remarkable starvation-induced program of development that transforms a population of unicellular amoebae into a fruiting body composed of resistant spores suspended on a stalk. During this development, secreted cAMP drives chemotaxis of the amoebae, leading to their aggregation, and subsequent differentiation and morphogenesis. Four sequentially expressed G protein-coupled receptors (GPCRs) for cAMP play critical roles in this process. The first of these, cAR1, is essential for aggregation as it mediates chemotaxis as well as the propagation of secreted cAMP waves throughout aggregating populations. Ligand-induced internalization has been shown to regulate a variety of GPCRs. However, little was known at the outset of this study about the role of internalization in the regulation of cAR1 function or, for that matter, in developmental systems in general. For this study, cAMP-induced cAR1 internalization was assessed by measuring (1) the reduction of cell surface binding sites for [ 3H]cAMP and (2) the redistribution of YFP-tagged receptors to the cell's interior, cAMP was found to induce little or no loss of ligand binding (LLB) in vegetative cells. However, the ability to induce LLB increased progressively over the initial 6 hrs of development, reaching ∼70% in cells undergoing aggregation. Despite these reductions in surface binding, detectable cAR1-YFP redistribution could be induced by cAMP only after the cells reached the mound stage (10 hrs) and was found to occur naturally by the ensuing slug stage (18 hrs). Site-directed substitution of a cluster of 5 serines in the receptor's cytoplasmic tail that was previously shown to be the principal site of cAMP-induced cAR1 phosphorylation impaired both LLB and receptor redistribution and furthermore resulted in mound-stage developmental arrest, suggesting that phosphorylation of cAR1 is a prerequisite for its internalization and that cAR1 internalization is required for post-aggregative development. To assess the involvement of clathrin mediated endocytosis, Dictyostelium cells lacking the clathrin light chain gene (clc-) or either of two dynamin genes were examined and found to be defective in LLB and, in the case of clc- cells, also cAR1 redistribution and turnover. Furthermore, cAR1 overexpression in clc- cells (like the serine mutant in wild-type cells) promoted developmental arrest in mounds. The mound-arrest phenotype was also recapitulated in a wild-type background by the specific expression of cAR1 in prestalk cells (but not prespore cells), suggesting that development depends critically on internalization and clearance of cAR1 from these cells. Persistent cAR1 expression following aggregation was found to be associated with aberrant expression of prestalk and prespore genes, which may adversely affect development in the prestalk cell lineage. The PI3 kinase-TORC2 signal transduction pathway, known to be important for Dictyostelium chemotaxis and internalization of yeast pheromone receptors, was examined using chemical inhibitors and null cells and found to be necessary for cAR1 internalization. In conclusion, cAR1 was shown to be similar to other GPCRs in that its internalization depends on phosphorylation of cytoplasmic domain serines, utilizes clathrin and dynamin, and involves the TORC2 complex. In addition, the findings presented here that cAR1 internalization is both developmentally regulated and required for normal development represent a novel regulatory paradigm that might pertain to other GPCRs known to play important roles in the development of humans and other metazoans. ^

A simulation study of the standard design, the rolling six design, the CRM, and the modified CRM in Phase I clinical trials

The Phase I clinical trial is considered the "first in human" study in medical research to examine the toxicity of a new agent. It determines the maximum tolerable dose (MTD) of a new agent, i.e., the highest dose in which toxicity is still acceptable. Several phase I clinical trial designs have been proposed in the past 30 years. The well known standard method, so called the 3+3 design, is widely accepted by clinicians since it is the easiest to implement and it does not need a statistical calculation. Continual reassessment method (CRM), a design uses Bayesian method, has been rising in popularity in the last two decades. Several variants of the CRM design have also been suggested in numerous statistical literatures. Rolling six is a new method introduced in pediatric oncology in 2008, which claims to shorten the trial duration as compared to the 3+3 design. The goal of the present research was to simulate clinical trials and compare these phase I clinical trial designs. Patient population was created by discrete event simulation (DES) method. The characteristics of the patients were generated by several distributions with the parameters derived from a historical phase I clinical trial data review. Patients were then selected and enrolled in clinical trials, each of which uses the 3+3 design, the rolling six, or the CRM design. Five scenarios of dose-toxicity relationship were used to compare the performance of the phase I clinical trial designs. One thousand trials were simulated per phase I clinical trial design per dose-toxicity scenario. The results showed the rolling six design was not superior to the 3+3 design in terms of trial duration. The time to trial completion was comparable between the rolling six and the 3+3 design. However, they both shorten the duration as compared to the two CRM designs. Both CRMs were superior to the 3+3 design and the rolling six in accuracy of MTD estimation. The 3+3 design and rolling six tended to assign more patients to undesired lower dose levels. The toxicities were slightly greater in the CRMs.^

A simulation study of the effects of small center enrollment in multi -center clinical trials

Multi-center clinical trials are very common in the development of new drugs and devices. One concern in such trials, is the effect of individual investigational sites enrolling small numbers of patients on the overall result. Can the presence of small centers cause an ineffective treatment to appear effective when treatment-by-center interaction is not statistically significant?^ In this research, simulations are used to study the effect that centers enrolling few patients may have on the analysis of clinical trial data. A multi-center clinical trial with 20 sites is simulated to investigate the effect of a new treatment in comparison to a placebo treatment. Twelve of these 20 investigational sites are considered small, each enrolling less than four patients per treatment group. Three clinical trials are simulated with sample sizes of 100, 170 and 300. The simulated data is generated with various characteristics, one in which treatment should be considered effective and another where treatment is not effective. Qualitative interactions are also produced within the small sites to further investigate the effect of small centers under various conditions.^ Standard analysis of variance methods and the "sometimes-pool" testing procedure are applied to the simulated data. One model investigates treatment and center effect and treatment-by-center interaction. Another model investigates treatment effect alone. These analyses are used to determine the power to detect treatment-by-center interactions, and the probability of type I error.^ We find it is difficult to detect treatment-by-center interactions when only a few investigational sites enrolling a limited number of patients participate in the interaction. However, we find no increased risk of type I error in these situations. In a pooled analysis, when the treatment is not effective, the probability of finding a significant treatment effect in the absence of significant treatment-by-center interaction is well within standard limits of type I error. ^