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.

Attractant and repellent signaling conformers of sensory rhodopsin-transducer complexes.

Attractant and repellent signaling conformers of the dual-signaling phototaxis receptor sensory rhodopsin I and its transducer subunit (SRI-HtrI) have recently been distinguished experimentally by the opposite connection of their retinylidene protonated Schiff bases to the outwardly located periplasmic side and inwardly located cytoplasmic side. Here we show that the pK(a) of the outwardly located Asp76 counterion in the outwardly connected conformer is lowered by approximately 1.5 units from that of the inwardly connected conformer. The pK(a) difference enables quantitative determination of the relative amounts of the two conformers in wild-type cells and behavioral mutants prior to photoexcitation, comparison of their absorption spectra, and determination of their relative signaling efficiency. We have shown that the one-photon excitation of the SRI-HtrI attractant conformer causes a Schiff base connectivity switch from inwardly connected to outwardly connected states in the attractant signaling photoreaction. Conversely, a second near-UV photon drives the complex back to the inwardly connected conformer in the repellent signaling photoreaction. The results suggest a model of the color-discriminating dual-signaling mechanism in which phototaxis responses (his-kinase modulation) result from the photointerconversion of the two oppositely connected SRI-HtrI conformers by one-photon and two-photon activation. Furthermore, we find that the related repellent phototaxis SRII-HtrII receptor complex has an outwardly connected retinylidene Schiff base like the repellent signaling forms of the SRI-HtrI complex, indicating the general applicability of macro conformational changes, which can be detected by the connectivity switch, to phototaxis signaling by sensory rhodopsin-transducer complexes.

Evolutionarily conserved role of calcineurin in phosphodegron-dependent degradation of phosphodiesterase 4D.

Calcineurin is a widely expressed and highly conserved Ser/Thr phosphatase. Calcineurin is inhibited by the immunosuppressant drug cyclosporine A (CsA) or tacrolimus (FK506). The critical role of CsA/FK506 as an immunosuppressant following transplantation surgery provides a strong incentive to understand the phosphatase calcineurin. Here we uncover a novel regulatory pathway for cyclic AMP (cAMP) signaling by the phosphatase calcineurin which is also evolutionarily conserved in Caenorhabditis elegans. We found that calcineurin binds directly to and inhibits the proteosomal degradation of cAMP-hydrolyzing phosphodiesterase 4D (PDE4D). We show that ubiquitin conjugation and proteosomal degradation of PDE4D are controlled by a cullin 1-containing E(3) ubiquitin ligase complex upon dual phosphorylation by casein kinase 1 (CK1) and glycogen synthase kinase 3beta (GSK3beta) in a phosphodegron motif. Our findings identify a novel signaling process governing G-protein-coupled cAMP signal transduction-opposing actions of the phosphatase calcineurin and the CK1/GSK3beta protein kinases on the phosphodegron-dependent degradation of PDE4D. This novel signaling system also provides unique functional insights into the complications elicited by CsA in transplant patients.

A Schiff base connectivity switch in sensory rhodopsin signaling.

Sensory rhodopsin I (SRI) in Halobacterium salinarum acts as a receptor for single-quantum attractant and two-quantum repellent phototaxis, transmitting light stimuli via its bound transducer HtrI. Signal-inverting mutations in the SRI-HtrI complex reverse the single-quantum response from attractant to repellent. Fast intramolecular charge movements reported here reveal that the unphotolyzed SRI-HtrI complex exists in two conformational states, which differ by their connection of the retinylidene Schiff base in the SRI photoactive site to inner or outer half-channels. In single-quantum photochemical reactions, the conformer with the Schiff base connected to the cytoplasmic (CP) half-channel generates an attractant signal, whereas the conformer with the Schiff base connected to the extracellular (EC) half-channel generates a repellent signal. In the wild-type complex the conformer equilibrium is poised strongly in favor of that with CP-accessible Schiff base. Signal-inverting mutations shift the equilibrium in favor of the EC-accessible Schiff base form, and suppressor mutations shift the equilibrium back toward the CP-accessible Schiff base form, restoring the wild-type phenotype. Our data show that the sign of the behavioral response directly correlates with the state of the connectivity switch, not with the direction of proton movements or changes in acceptor pK(a). These findings identify a shared fundamental process in the mechanisms of transport and signaling by the rhodopsin family. Furthermore, the effects of mutations in the HtrI subunit of the complex on SRI Schiff base connectivity indicate that the two proteins are tightly coupled to form a single unit that undergoes a concerted conformational transition.

Localized diacylglycerol-dependent stimulation of Ras and Rap1 during phagocytosis.

We describe a role for diacylglycerol in the activation of Ras and Rap1 at the phagosomal membrane. During phagocytosis, Ras density was similar on the surface and invaginating areas of the membrane, but activation was detectable only in the latter and in sealed phagosomes. Ras activation was associated with the recruitment of RasGRP3, a diacylglycerol-dependent Ras/Rap1 exchange factor. Recruitment to phagosomes of RasGRP3, which contains a C1 domain, parallels and appears to be due to the formation of diacylglycerol. Accordingly, Ras and Rap1 activation was precluded by antagonists of phospholipase C and of diacylglycerol binding. Ras is dispensable for phagocytosis but controls activation of extracellular signal-regulated kinase, which is partially impeded by diacylglycerol inhibitors. By contrast, cross-activation of complement receptors by stimulation of Fcgamma receptors requires Rap1 and involves diacylglycerol. We suggest a role for diacylglycerol-dependent exchange factors in the activation of Ras and Rap1, which govern distinct processes induced by Fcgamma receptor-mediated phagocytosis to enhance the innate immune response.

Ligand-signaled upregulation of Enterococcus faecalis ace transcription, a mechanism for modulating host-E. faecalis interaction.

Enterococcus faecalis, the third most frequent cause of bacterial endocarditis, appears to be equipped with diverse surface-associated proteins showing structural-fold similarity to the immunoglobulin-fold family of staphylococcal adhesins. Among the putative E. faecalis surface proteins, the previously characterized adhesin Ace, which shows specific binding to collagen and laminin, was detectable in surface protein preparations only after growth at 46 degrees C, mirroring the finding that adherence was observed in 46 degrees C, but not 37 degrees C, grown E. faecalis cultures. To elucidate the influence of different growth and host parameters on ace expression, we investigated ace expression using E. faecalis OG1RF grown in routine laboratory media (brain heart infusion) and found that ace mRNA levels were low in all growth phases. However, quantitative reverse transcription-PCR showed 18-fold-higher ace mRNA amounts in cells grown in the presence of collagen type IV compared to the controls. Similarly, a marked increase was observed when cells were either grown in the presence of collagen type I or serum but not in the presence of fibrinogen or bovine serum albumin. The production of Ace after growth in the presence of collagen type IV was demonstrated by immunofluorescence microscopy, mirroring the increased ace mRNA levels. Furthermore, increased Ace expression correlated with increased collagen and laminin adhesion. Collagen-induced Ace expression was also seen in three of three other E. faecalis strains of diverse origins tested, and thus it appears to be a common phenomenon. The observation of host matrix signal-induced adherence of E. faecalis may have important implications on our understanding of this opportunistic pathogen.

Role of the cytoplasmic domain in Anabaena sensory rhodopsin photocycling: vectoriality of Schiff base deprotonation.

Light-induced electric signals in intact E. coli cells generated by heterologously expressed full-length and C-terminally truncated versions of Anabaena sensory rhodopsin (ASR) demonstrate that the charge movements within the membrane-embedded part of the molecule are stringently controlled by the cytoplasmic domain. In particular, truncation inverts the direction of proton movement during Schiff base deprotonation from outward to cytoplasmic. Truncation also alters faster charge movements that occur before Schiff base deprotonation. Asp(217) as previously shown by FTIR serves as a proton acceptor in the truncated ASR but not in the full-length version, and its mutation to Asn restores the natural outward direction of proton movement. Introduction of a potential negative charge (Ser(86) to Asp) on the cytoplasmic side favors a cytoplasmic direction of proton release from the Schiff base. In contrast, mutation of the counterion Asp(75) to Glu reverses the photocurrent to the outward direction in the truncated pigment, and in both truncated and full-length versions accelerates Schiff base deprotonation more than 10-fold. The communication between the cytoplasmic domain and the membrane-embedded photoactive site of ASR demonstrated here is likely to derive from the receptor's use of a cytoplasmic protein for signal transduction, as has been suggested previously from binding studies.

Crystal structure of the Anabaena sensory rhodopsin transducer.

We present crystal structures of the Anabaena sensory rhodopsin transducer (ASRT), a soluble cytoplasmic protein that interacts with the first structurally characterized eubacterial retinylidene photoreceptor Anabaena sensory rhodopsin (ASR). Four crystal structures of ASRT from three different spacegroups were obtained, in all of which ASRT is present as a planar (C4) tetramer, consistent with our characterization of ASRT as a tetramer in solution. The ASRT tetramer is tightly packed, with large interfaces where the well-structured beta-sandwich portion of the monomers provides the bulk of the tetramer-forming interactions, and forms a flat, stable surface on one side of the tetramer (the beta-face). Only one of our four different ASRT crystals reveals a C-terminal alpha-helix in the otherwise all-beta protein, together with a large loop from each monomer on the opposite face of the tetramer (the alpha-face), which is flexible and largely disordered in the other three crystal forms. Gel-filtration chromatography demonstrated that ASRT forms stable tetramers in solution and isothermal microcalorimetry showed that the ASRT tetramer binds to ASR with a stoichiometry of one ASRT tetramer per one ASR photoreceptor with a K(d) of 8 microM in the highest affinity measurements. Possible mechanisms for the interaction of this transducer tetramer with the ASR photoreceptor via its flexible alpha-face to mediate transduction of the light signal are discussed.

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.

Morphoproteomic analysis reveals an overexpressed and constitutively activated phospholipase D1-mTORC2 pathway in endometrial carcinoma.

The mammalian target of rapamycin (MTOR) assembles into two distinct complexes: mTOR complex 1 (mTORC1) is predominantly cytoplasmic and highly responsive to rapamycin, whereas mTOR complex 2 (mTORC2) is both cytoplasmic and nuclear, and relatively resistant to rapamycin. mTORC1 and mTORC2 phosphorylatively regulate their respective downstream effectors p70S6K/4EBP1, and Akt. The resulting activated mTOR pathways stimulate protein synthesis, cellular proliferation, and cell survival. Moreover, phospholipase D (PLD) and its product, phosphatidic acid (PA) have been implicated as one of the upstream activators of mTOR signaling. In this study, we investigated the activation status as well as the subcellular distribution of mTOR, and its upstream regulators and downstream effectors in endometrial carcinomas (ECa) and non-neoplastic endometrial control tissue. Our data show that the mTORC2 activity is selectively elevated in endometrial cancers as evidenced by a predominant nuclear localization of the activated form of mTOR (p-mTOR at Ser2448) in malignant epithelium, accompanied by overexpression of nuclear p-Akt (Ser473), as well as overexpression of vascular endothelial growth factor (VEGF)-A isoform, the latter a resultant of target gene activation by mTORC2 signaling via hypoxia-inducible factor (HIF)-2alpha. In addition, expression of PLD1, one of the two major isoforms of PLD in human, is increased in tumor epithelium. In summary, we demonstrate that the PLD1/PA-mTORC2 signal pathway is overactivated in endometrial carcinomas. This suggests that the rapamycin-insensitive mTORC2 pathway plays a major role in endometrial tumorigenesis and that therapies designed to target the phospholipase D pathway and components of the mTORC2 pathway should be efficacious against ECa.

Dopamine-stimulated dephosphorylation of connexin 36 mediates AII amacrine cell uncoupling.

Gap junction proteins form the substrate for electrical coupling between neurons. These electrical synapses are widespread in the CNS and serve a variety of important functions. In the retina, connexin 36 (Cx36) gap junctions couple AII amacrine cells and are a requisite component of the high-sensitivity rod photoreceptor pathway. AII amacrine cell coupling strength is dynamically regulated by background light intensity, and uncoupling is thought to be mediated by dopamine signaling via D(1)-like receptors. One proposed mechanism for this uncoupling involves dopamine-stimulated phosphorylation of Cx36 at regulatory sites, mediated by protein kinase A. Here we provide evidence against this hypothesis and demonstrate a direct relationship between Cx36 phosphorylation and AII amacrine cell coupling strength. Dopamine receptor-driven uncoupling of the AII network results from protein kinase A activation of protein phosphatase 2A and subsequent dephosphorylation of Cx36. Protein phosphatase 1 activity negatively regulates this pathway. We also find that Cx36 gap junctions can exist in widely different phosphorylation states within a single neuron, implying that coupling is controlled at the level of individual gap junctions by locally assembled signaling complexes. This kind of synapse-by-synapse plasticity allows for precise control of neuronal coupling, as well as cell-type-specific responses dependent on the identity of the signaling complexes assembled.

Thoracic aortic disease in tuberous sclerosis complex: molecular pathogenesis and potential therapies in Tsc2+/- mice.

Tuberous sclerosis complex (TSC) is a genetic disorder with pleiotropic manifestations caused by heterozygous mutations in either TSC1 or TSC2. One of the less investigated complications of TSC is the formation of aneurysms of the descending aorta, which are characterized on pathologic examination by smooth muscle cell (SMC) proliferation in the aortic media. SMCs were explanted from Tsc2(+/-) mice to investigate the pathogenesis of aortic aneurysms caused by TSC2 mutations. Tsc2(+/-) SMCs demonstrated increased phosphorylation of mammalian target of rapamycin (mTOR), S6 and p70S6K and increased proliferation rates compared with wild-type (WT) SMCs. Tsc2(+/-) SMCs also had reduced expression of SMC contractile proteins compared with WT SMCs. An inhibitor of mTOR signaling, rapamycin, decreased SMC proliferation and increased contractile protein expression in the Tsc2(+/-) SMCs to levels similar to WT SMCs. Exposure to alpha-elastin fragments also decreased proliferation of Tsc2(+/-) SMCs and increased levels of p27(kip1), but failed to increase expression of contractile proteins. In response to artery injury using a carotid artery ligation model, Tsc2(+/-) mice significantly increased neointima formation compared with the control mice, and the neointima formation was inhibited by treatment with rapamycin. These results demonstrate that Tsc2 haploinsufficiency in SMCs increases proliferation and decreases contractile protein expression and suggest that the increased proliferative potential of the mutant cells may be suppressed in vivo by interaction with elastin. These findings provide insights into the molecular pathogenesis of aortic disease in TSC patients and identify a potential therapeutic target for treatment of this complication of the disease.

Altered transmission of HOX and apoptotic SNPs identify a potential common pathway for clubfoot.

Clubfoot is a common birth defect that affects 135,000 newborns each year worldwide. It is characterized by equinus deformity of one or both feet and hypoplastic calf muscles. Despite numerous study approaches, the cause(s) remains poorly understood although a multifactorial etiology is generally accepted. We considered the HOXA and HOXD gene clusters and insulin-like growth factor binding protein 3 (IGFBP3) as candidate genes because of their important roles in limb and muscle morphogenesis. Twenty SNPs from the HOXA and HOXD gene clusters and 12 SNPs in IGFBP3 were genotyped in a sample composed of non-Hispanic white and Hispanic multiplex and simplex families (discovery samples) and a second sample of non-Hispanic white simplex trios (validation sample). Four SNPs (rs6668, rs2428431, rs3801776, and rs3779456) in the HOXA cluster demonstrated altered transmission in the discovery sample, but only rs3801776, located in the HOXA basal promoter region, showed altered transmission in both the discovery and validation samples (P = 0.004 and 0.028). Interestingly, HOXA9 is expressed in muscle during development. An SNP in IGFBP3, rs13223993, also showed altered transmission (P = 0.003) in the discovery sample. Gene-gene interactions were identified between variants in HOXA, HOXD, and IGFBP3 and with previously associated SNPs in mitochondrial-mediated apoptotic genes. The most significant interactions were found between CASP3 SNPS and variants in HOXA, HOXD, and IGFBP3. These results suggest a biologic model for clubfoot in which perturbation of HOX and apoptotic genes together affect muscle and limb development, which may cause the downstream failure of limb rotation into a plantar grade position.

Sequencing of TGF-beta pathway genes in familial cases of intracranial aneurysm.

BACKGROUND AND PURPOSE: Familial aggregation of intracranial aneurysms (IA) strongly suggests a genetic contribution to pathogenesis. However, genetic risk factors have yet to be defined. For families affected by aortic aneurysms, specific gene variants have been identified, many affecting the receptors to transforming growth factor-beta (TGF-beta). In recent work, we found that aortic and intracranial aneurysms may share a common genetic basis in some families. We hypothesized, therefore, that mutations in TGF-beta receptors might also play a role in IA pathogenesis. METHODS: To identify genetic variants in TGF-beta and its receptors, TGFB1, TGFBR1, TGFBR2, ACVR1, TGFBR3, and ENG were directly sequenced in 44 unrelated patients with familial IA. Novel variants were confirmed by restriction digestion analyses, and allele frequencies were analyzed in cases versus individuals without known intracranial disease. Similarly, allele frequencies of a subset of known SNPs in each gene were also analyzed for association with IA. RESULTS: No mutations were found in TGFB1, TGFBR1, TGFBR2, or ACVR1. Novel variants identified in ENG (p.A60E) and TGFBR3 (p.W112R) were not detected in at least 892 reference chromosomes. ENG p.A60E showed significant association with familial IA in case-control studies (P=0.0080). No association with IA could be found for any of the known polymorphisms tested. CONCLUSIONS: Mutations in TGF-beta receptor genes are not a major cause of IA. However, we identified rare variants in ENG and TGFBR3 that may be important for IA pathogenesis in a subset of families.

Morphoproteomic evidence of constitutively activated and overexpressed mTOR pathway in cervical squamous carcinoma and high grade squamous intraepithelial lesions.

Human papilloma virus (HPV) infection of the uterine cervix is linked to the pathogenesis of cervical cancer. Preclinical in vitro and in vivo studies using HPV-containing human cervical carcinoma cell lines have shown that the mammalian target of rapamycin (mTOR) inhibitor, rapamycin, and epidermal growth factor receptor (EGFR)-tyrosine kinase inhibitor, erlotinib, can induce growth delay of xenografts. Activation of Akt and mTOR are also observed in cervical squamous cell carcinoma and, the expression of phosphorylated mTOR was reported to serve as a marker to predict response to chemotherapy and survival of cervical cancer patients. Therefore, we investigated: a) the expression level of EGFR in cervical squamous cell carcinoma (SCC) and high-grade squamous intraepithelial lesions (HSIL) versus non-neoplastic cervical squamous epithelium; b) the state of activation of the mTOR pathway in these same tissues; and c) any impact of these signal transduction molecules on cell cycle. Formalin-fixed paraffin-embedded tissue microarray blocks containing 20 samples each of normal cervix, HSIL and invasive SCC, derived from a total of 60 cases of cervical biopsies and cervical conizations were examined. Immunohistochemistry was utilized to detect the following antigens: EGFR; mTOR pathway markers, phosphorylated (p)-mTOR (Ser2448) and p-p70S6K (Thr389); and cell cycle associated proteins, Ki-67 and S phase kinase-associated protein (Skp)2. Protein compartmentalization and expression were quantified in regard to proportion (0-100%) and intensity (0-3+). Mitotic index (MI) was also assessed. An expression index (EI) for pmTOR, p-p70S6K and EGFR, respectively was calculated by taking the product of intensity score and proportion of positively staining cells. We found that plasmalemmal EGFR expression was limited to the basal/parabasal cells (2-3+, EI = 67) in normal cervical epithelium (NL), but was diffusely positive in all HSIL (EI = 237) and SCC (EI 226). The pattern of cytoplasmic p-mTOR and nuclear p-p70S6K expression was similar to that of EGFR; all showed a significantly increased EI in HSIL/SCC versus NL (p<0.02). Nuclear translocation of p-mTOR was observed in all SCC lesions (EI = 202) and was significantly increased versus both HSIL (EI = 89) and NL (EI = 54) with p<0.015 and p<0.0001, respectively. Concomitant increases in MI and proportion of nuclear Ki-67 and Skp2 expression were noted in HSIL and SCC. In conclusion, morphoproteomic analysis reveals constitutive activation and overexpression of the mTOR pathway in HSIL and SCC as evidenced by: increased nuclear translocation of pmTOR and p-p70S6K, phosphorylated at putative sites of activation, Ser2448 and Thr389, respectively; correlative overexpression of the upstream signal transducer, EGFR, and increases in cell cycle correlates, Skp2 and mitotic indices. These results suggest that the mTOR pathway plays a key role in cervical carcinogenesis and targeted therapies may be developed for SCC as well as its precursor lesion, HSIL.