Characterization of suramin binding sites on the human group IIA secreted phospholipase A(2) by site-directed mutagenesis and molecular dynamics simulation

Suramin is a polysulphonated naphthylurea with inhibitory activity against the human secreted group IIA phospholipase A(2) (hsPLA2GIIA), and we have investigated suramin binding to recombinant hsPLA2GIIA using site-directed mutagenesis and molecular dynamics (MD) simulations. The changes in suramin binding affinity of 13 cationic residue mutants of the hsPLA2GIIA was strongly correlated with alterations in the inhibition of membrane damaging activity of the protein. Suramin binding to hsPLA2GIIA was also studied by MD simulations, which demonstrated that altered intermolecular potential energy of the suramin/mutant complexes was a reliable indicator of affinity change. Although residues in the C-terminal region play a major role in the stabilization of the hsPLA2GIIA/suramin complex, attractive and repulsive hydrophobic and electrostatic interactions with residues throughout the protein together with the adoption of a bent suramin conformation, all contribute to the stability of the complex. Analysis of the h5PLA2GIIA/suramin interactions allows the prediction of the properties of suramin analogues with improved binding and higher affinities which may be candidates for novel phospholipase A(2) inhibitors. (C) 2012 Elsevier Inc. All rights reserved.

Regulation of SRC-3 localization and dynamics by phosphorylation during ERα-dependent transcriptional activation

Transcription is controlled by promoter-selective transcriptional factors (TFs), which bind to cis-regulatory enhancers elements, termed hormone response elements (HREs), in a specific subset of genes. Regulation by these factors involves either the recruitment of coactivators or corepressors and direct interaction with the basal transcriptional machinery (1). Hormone-activated nuclear receptors (NRs) are well characterized transcriptional factors (2) that bind to the promoters of their target genes and recruit primary and secondary coactivator proteins which possess many enzymatic activities required for gene expression (1,3,4). In the present study, using single-cell high-resolution fluorescent microscopy and high throughput microscopy (HTM) coupled to computational imaging analysis, we investigated transcriptional regulation controlled by the estrogen receptor alpha (ERalpha), in terms of large scale chromatin remodeling and interaction with the associated coactivator SRC-3 (Steroid Receptor Coactivator-3), a member of p160 family (28) primary coactivators. ERalpha is a steroid-dependent transcriptional factor (16) that belongs to the NRs superfamily (2,3) and, in response to the hormone 17-ß estradiol (E2), regulates transcription of distinct target genes involved in development, puberty, and homeostasis (8,16). ERalpha spends most of its lifetime in the nucleus and undergoes a rapid (within minutes) intranuclear redistribution following the addition of either agonist or antagonist (17,18,19). We designed a HeLa cell line (PRL-HeLa), engineered with a chromosomeintegrated reporter gene array (PRL-array) containing multicopy hormone response-binding elements for ERalpha that are derived from the physiological enhancer/promoter region of the prolactin gene. Following GFP-ER transfection of PRL-HeLa cells, we were able to observe in situ ligand dependent (i) recruitment to the array of the receptor and associated coregulators, (ii) chromatin remodeling, and (iii) direct transcriptional readout of the reporter gene. Addition of E2 causes a visible opening (decondensation) of the PRL-array, colocalization of RNA Polymerase II, and transcriptional readout of the reporter gene, detected by mRNA FISH. On the contrary, when cells were treated with an ERalpha antagonist (Tamoxifen or ICI), a dramatic condensation of the PRL-array was observed, displacement of RNA Polymerase II, and complete decreasing in the transcriptional FISH signal. All p160 family coactivators (28) colocalize with ERalpha at the PRL-array. Steroid Receptor Coactivator-3 (SRC-3/AIB1/ACTR/pCIP/RAC3/TRAM1) is a p160 family member and a known oncogenic protein (4,34). SRC-3 is regulated by a variety of posttranslational modifications, including methylation, phosphorylation, acetylation, ubiquitination and sumoylation (4,35). These events have been shown to be important for its interaction with other coactivator proteins and NRs and for its oncogenic potential (37,39). A number of extracellular signaling molecules, like steroid hormones, growth factors and cytokines, induce SRC-3 phosphorylation (40). These actions are mediated by a wide range of kinases, including extracellular-regulated kinase 1 and 2 (ERK1-2), c-Jun N-terminal kinase, p38 MAPK, and IkB kinases (IKKs) (41,42,43). Here, we report SRC-3 to be a nucleocytoplasmic shuttling protein, whose cellular localization is regulated by phosphorylation and interaction with ERalpha. Using a combination of high throughput and fluorescence microscopy, we show that both chemical inhibition (with U0126) and siRNA downregulation of the MAP/ERK1/2 kinase (MEK1/2) pathway induce a cytoplasmic shift in SRC-3 localization, whereas stimulation by EGF signaling enhances its nuclear localization by inducing phosphorylation at T24, S857, and S860, known partecipants in the regulation of SRC-3 activity (39). Accordingly, the cytoplasmic localization of a non-phosphorylatable SRC-3 mutant further supports these results. In the presence of ERalpha, U0126 also dramatically reduces: hormone-dependent colocalization of ERalpha and SRC-3 in the nucleus; formation of ER-SRC-3 coimmunoprecipitation complex in cell lysates; localization of SRC-3 at the ER-targeted prolactin promoter array (PRL-array) and transcriptional activity. Finally, we show that SRC-3 can also function as a cotransporter, facilitating the nuclear-cytoplasmic shuttling of estrogen receptor. While a wealth of studies have revealed the molecular functions of NRs and coregulators, there is a paucity of data on how these functions are spatiotemporally organized in the cellular context. Technically and conceptually, our findings have a new impact upon evaluating gene transcriptional control and mechanisms of action of gene regulators.

Studies of OPA1 pathogenic mechanisms in Dominant Optic Atrophy and novel protein function in mitochondrial DNA stability

The mitochondrion is an essential cytoplasmic organelle that provides most of the energy necessary for eukaryotic cell physiology. Mitochondrial structure and functions are maintained by proteins of both mitochondrial and nuclear origin. These organelles are organized in an extended network that dynamically fuses and divides. Mitochondrial morphology results from the equilibrium between fusion and fission processes, controlled by a family of “mitochondria-shaping” proteins. It is becoming clear that defects in mitochondrial dynamics can impair mitochondrial respiration, morphology and motility, leading to apoptotic cell death in vitro and more or less severe neurodegenerative disorders in vivo in humans. Mutations in OPA1, a nuclear encoded mitochondrial protein, cause autosomal Dominant Optic Atrophy (DOA), a heterogeneous blinding disease characterized by retinal ganglion cell degeneration leading to optic neuropathy (Delettre et al., 2000; Alexander et al., 2000). OPA1 is a mitochondrial dynamin-related guanosine triphosphatase (GTPase) protein involved in mitochondrial network dynamics, cytochrome c storage and apoptosis. This protein is anchored or associated on the inner mitochondrial membrane facing the intermembrane space. Eight OPA1 isoforms resulting from alternative splicing combinations of exon 4, 4b and 5b have been described (Delettre et al., 2001). These variants greatly vary among diverse organs and the presence of specific isoforms has been associated with various mitochondrial functions. The different spliced exons encode domains included in the amino-terminal region and contribute to determine OPA1 functions (Olichon et al., 2006). It has been shown that exon 4, that is conserved throughout evolution, confers functions to OPA1 involved in maintenance of the mitochondrial membrane potential and in the fusion of the network. Conversely, exon 4b and exon 5b, which are vertebrate specific, are involved in regulation of cytochrome c release from mitochondria, and activation of apoptosis, a process restricted to vertebrates (Olichon et al., 2007). While Mgm1p has been identified thanks to its role in mtDNA maintenance, it is only recently that OPA1 has been linked to mtDNA stability. Missense mutations in OPA1 cause accumulation of multiple deletions in skeletal muscle. The syndrome associated to these mutations (DOA-1 plus) is complex, consisting of a combination of dominant optic atrophy, progressive external ophtalmoplegia, peripheral neuropathy, ataxia and deafness (Amati- Bonneau et al., 2008; Hudson et al., 2008). OPA1 is the fifth gene associated with mtDNA “breakage syndrome” together with ANT1, PolG1-2 and TYMP (Spinazzola et al., 2009). In this thesis we show for the first time that specific OPA1 isoforms associated to exon 4b are important for mtDNA stability, by anchoring the nucleoids to the inner mitochondrial membrane. Our results clearly demonstrate that OPA1 isoforms including exon 4b are intimately associated to the maintenance of the mitochondrial genome, as their silencing leads to mtDNA depletion. The mechanism leading to mtDNA loss is associated with replication inhibition in cells where exon 4b containing isoforms were down-regulated. Furthermore silencing of exon 4b associated isoforms is responsible for alteration in mtDNA-nucleoids distribution in the mitochondrial network. In this study it was evidenced that OPA1 exon 4b isoform is cleaved to provide a 10kd peptide embedded in the inner membrane by a second transmembrane domain, that seems to be crucial for mitochondrial genome maintenance and does correspond to the second transmembrane domain of the yeasts orthologue encoded by MGM1 or Msp1, which is also mandatory for this process (Diot et al., 2009; Herlan et al., 2003). Furthermore in this thesis we show that the NT-OPA1-exon 4b peptide co-immuno-precipitates with mtDNA and specifically interacts with two major components of the mitochondrial nucleoids: the polymerase gamma and Tfam. Thus, from these experiments the conclusion is that NT-OPA1- exon 4b peptide contributes to the nucleoid anchoring in the inner mitochondrial membrane, a process that is required for the initiation of mtDNA replication and for the distribution of nucleoids along the network. These data provide new crucial insights in understanding the mechanism involved in maintenance of mtDNA integrity, because they clearly demonstrate that, besides genes implicated in mtDNA replications (i.e. polymerase gamma, Tfam, twinkle and genes involved in the nucleotide pool metabolism), OPA1 and mitochondrial membrane dynamics play also an important role. Noticeably, the effect on mtDNA is different depending on the specific OPA1 isoforms down-regulated, suggesting the involvement of two different combined mechanisms. Over two hundred OPA1 mutations, spread throughout the coding region of the gene, have been described to date, including substitutions, deletions or insertions. Some mutations are predicted to generate a truncated protein inducing haploinsufficiency, whereas the missense nucleotide substitutions result in aminoacidic changes which affect conserved positions of the OPA1 protein. So far, the functional consequences of OPA1 mutations in cells from DOA patients are poorly understood. Phosphorus MR spectroscopy in patients with the c.2708delTTAG deletion revealed a defect in oxidative phosphorylation in muscles (Lodi et al., 2004). An energetic impairment has been also show in fibroblasts with the severe OPA1 R445H mutation (Amati-Bonneau et al., 2005). It has been previously reported by our group that OPA1 mutations leading to haploinsufficiency are associated in fibroblasts to an oxidative phosphorylation dysfunction, mainly involving the respiratory complex I (Zanna et al., 2008). In this study we have evaluated the energetic efficiency of a panel of skin fibroblasts derived from DOA patients, five fibroblast cell lines with OPA1 mutations causing haploinsufficiency (DOA-H) and two cell lines bearing mis-sense aminoacidic substitutions (DOA-AA), and compared with control fibroblasts. Although both types of DOA fibroblasts maintained a similar ATP content when incubated in a glucose-free medium, i.e. when forced to utilize the oxidative phosphorylation only to produce ATP, the mitochondrial ATP synthesis through complex I, measured in digitonin-permeabilized cells, was significantly reduced in cells with OPA1 haploinsufficiency only, whereas it was similar to controls in cells with the missense substitutions. Furthermore, evaluation of the mitochondrial membrane potential (DYm) in the two fibroblast lines DOA-AA and in two DOA-H fibroblasts, namely those bearing the c.2819-2A>C mutation and the c.2708delTTAG microdeletion, revealed an anomalous depolarizing response to oligomycin in DOA-H cell lines only. This finding clearly supports the hypothesis that these mutations cause a significant alteration in the respiratory chain function, which can be unmasked only when the operation of the ATP synthase is prevented. Noticeably, oligomycin-induced depolarization in these cells was almost completely prevented by preincubation with cyclosporin A, a well known inhibitor of the permeability transition pore (PTP). This results is very important because it suggests for the first time that the voltage threshold for PTP opening is altered in DOA-H fibroblasts. Although this issue has not yet been addressed in the present study, several are the mechanisms that have been proposed to lead to PTP deregulation, including in particular increased reactive oxygen species production and alteration of Ca2+ homeostasis, whose role in DOA fibroblasts PTP opening is currently under investigation. Identification of the mechanisms leading to altered threshold for PTP regulation will help our understanding of the pathophysiology of DOA, but also provide a strategy for therapeutic intervention.

Statistical mechanics of protein complexed and condensed DNA

In this thesis I treat various biophysical questions arising in the context of complexed / ”protein-packed” DNA and DNA in confined geometries (like in viruses or toroidal DNA condensates). Using diverse theoretical methods I consider the statistical mechanics as well as the dynamics of DNA under these conditions. In the first part of the thesis (chapter 2) I derive for the first time the single molecule ”equation of state”, i.e. the force-extension relation of a looped DNA (Eq. 2.94) by using the path integral formalism. Generalizing these results I show that the presence of elastic substructures like loops or deflections caused by anchoring boundary conditions (e.g. at the AFM tip or the mica substrate) gives rise to a significant renormalization of the apparent persistence length as extracted from single molecule experiments (Eqs. 2.39 and 2.98). As I show the experimentally observed apparent persistence length reduction by a factor of 10 or more is naturally explained by this theory. In chapter 3 I theoretically consider the thermal motion of nucleosomes along a DNA template. After an extensive analysis of available experimental data and theoretical modelling of two possible mechanisms I conclude that the ”corkscrew-motion” mechanism most consistently explains this biologically important process. In chapter 4 I demonstrate that DNA-spools (architectures in which DNA circumferentially winds on a cylindrical surface, or onto itself) show a remarkable ”kinetic inertness” that protects them from tension-induced disruption on experimentally and biologically relevant timescales (cf. Fig. 4.1 and Eq. 4.18). I show that the underlying model establishes a connection between the seemingly unrelated and previously unexplained force peaks in single molecule nucleosome and DNA-toroid stretching experiments. Finally in chapter 5 I show that toroidally confined DNA (found in viruses, DNAcondensates or sperm chromatin) undergoes a transition to a twisted, highly entangled state provided that the aspect ratio of the underlying torus crosses a certain critical value (cf. Eq. 5.6 and the phase diagram in Fig. 5.4). The presented mechanism could rationalize several experimental mysteries, ranging from entangled and supercoiled toroids released from virus capsids to the unexpectedly short cholesteric pitch in the (toroidaly wound) sperm chromatin. I propose that the ”topological encapsulation” resulting from our model may have some practical implications for the gene-therapeutic DNA delivery process.

Computational strategies to include protein flexibility in Ligand Docking and Virtual Screening

The dynamic character of proteins strongly influences biomolecular recognition mechanisms. With the development of the main models of ligand recognition (lock-and-key, induced fit, conformational selection theories), the role of protein plasticity has become increasingly relevant. In particular, major structural changes concerning large deviations of protein backbones, and slight movements such as side chain rotations are now carefully considered in drug discovery and development. It is of great interest to identify multiple protein conformations as preliminary step in a screening campaign. Protein flexibility has been widely investigated, in terms of both local and global motions, in two diverse biological systems. On one side, Replica Exchange Molecular Dynamics has been exploited as enhanced sampling method to collect multiple conformations of Lactate Dehydrogenase A (LDHA), an emerging anticancer target. The aim of this project was the development of an Ensemble-based Virtual Screening protocol, in order to find novel potent inhibitors. On the other side, a preliminary study concerning the local flexibility of Opioid Receptors has been carried out through ALiBERO approach, an iterative method based on Elastic Network-Normal Mode Analysis and Monte Carlo sampling. Comparison of the Virtual Screening performances by using single or multiple conformations confirmed that the inclusion of protein flexibility in screening protocols has a positive effect on the probability to early recognize novel or known active compounds.

Transcriptional dynamics of the human DMD locus

The human DMD locus encodes dystrophin protein. Absence or reduced levels of dystrophin (DMD or BMD phenotype, respectively) lead to progressive muscle wasting. Little is known about the complex coordination of dystrophin expression and its transcriptional regulation is a field of intense interest. In this work we found that DMD locus harbours multiple long non coding RNAs which orchestrate and control transcription of muscle dystrophin mRNA isoforms. These lncRNAs are tissue-specific and highly expressed during myogenesis, suggesting a possible role in tissue-specific expression of DMD gene isoforms. Their forced ectopic expression in human muscle and neuronal cells leads to a specific and negative regulation of endogenous dystrophin full lenght isoforms. An intriguing aspect regarding the transcription of the DMD locus is the gene size (2.4Mb). The mechanism that ensures the complete synthesis of the primary transcript and the coordinated splicing of 79 exons is still completely unknown. By ChIP-on-chip analyses, we discovered novel regions never been involved before in the transcription regulation of the DMD locus. Specifically, we observed enrichments for Pol II, P-Ser2, P-Ser5, Ac-H3 and 2Me-H3K4 in an intronic region of 3Kb (approximately 21Kb) downstream of the end of DMD exon 52 and in a region of 4Kb spanning the DMD exon 62. Interestingly, this latter region and the TSS of Dp71 are strongly marked by 3Me-H3K36, an histone modification associated with the regulation of splicing process. Furthermore, we also observed strong presence of open chromatin marks (Ac-H3 and 2Me-H3K4) around intron 34 and the exon 45 without presence of RNA pol II. We speculate that these two regions may exert an enhancer-like function on Dp427m promoter, although further investigations are necessary. Finally, we investigated the nuclear-cytoplasmic compartmentalization of the muscular dystrophin mRNA and, specifically, we verified whether the exon skipping therapy could influence its cellular distribution.

Local dynamics and bending mechanics of mesostructured materials

Die vorliegende Arbeit behandelt die Anwendung der Rasterkraftmikroskopie auf die Untersuchung mesostrukturierter Materialien. Mesostrukturierte Materialien setzen sich aus einzelnen mesoskopen Bausteinen zusammen. Diese Untereinheiten bestimmen im Wesentlichen ihr charakteristisches Verhalten auf äußere mechanische oder elektrische Reize, weshalb diesen Materialien eine besondere Rolle in der Natur sowie im täglichen Leben zukommt. Ein genaues Verständnis der Selbstorganisation dieser Materialien und der Wechselwirkung der einzelnen Bausteine untereinander ist daher von essentieller Bedeutung zur Entwicklung neuer Synthesestrategien sowie zur Optimierung ihrer Materialeigenschaften. Die Charakterisierung dieser mesostrukturierten Materialien erfolgt üblicherweise mittels makroskopischer Analysemethoden wie der dielektrischen Breitbandspektroskopie, Thermogravimetrie sowie in Biegungsexperimenten. In dieser Arbeit wird gezeigt, wie sich diese Analysemethoden mit der Rasterkraftmikroskopie verbinden lassen, um mesostrukturierte Materialien zu untersuchen. Die Rasterkraftmikroskopie bietet die Möglichkeit, die Oberfläche eines Materials abzubilden und zusätzlich dazu seine quantitativen Eigenschaften, wie die mechanische Biegefestigkeit oder die dielektrische Relaxation, zu bestimmen. Die Übertragung makroskopischer Analyseverfahren auf den Nano- bzw. Mikrometermaßstab mittels der Rasterkraftmikroskopie erlaubt die Charakterisierung von räumlich sehr begrenzten Proben bzw. von Proben, die nur in einer sehr kleinen Menge (<10 mg) vorliegen. Darüberhinaus umfasst das Auflösungsvermögen eines Rasterkraftmikroskops, welche durch die Größe seines Federbalkens (50 µm) sowie seines Spitzenradius (5 nm) definiert ist, genau den Längenskalenbereich, der einzelne Atome mit der makroskopischen Welt verbindet, nämlich die Mesoskala. In dieser Arbeit werden Polymerfilme, kolloidale Nanofasern sowie Biomineralien ausführlicher untersucht.rnIm ersten Projekt werden mittels Rasterkraftmikroskopie dielektrische Spektren von mischbaren Polymerfilmen aufgenommen und mit ihrer lokalen Oberflächenstruktur korreliert. Im zweiten Projekt wird die Rasterkraftmikroskopie eingesetzt, um Biegeexperimente an kolloidalen Nanofasern durchzuführen und so ihre Brucheigenschaften genauer zu untersuchen. Im letzten Projekt findet diese Methode Anwendung bei der Charakterisierung der Biegeeigenschaften von Biomineralien. Des Weiteren erfolgt eine Analyse der organischen Zusammensetzung dieser Biomineralien. Alle diese Projekte demonstrieren die vielseitige Einsetzbarkeit der Rasterkraftmikroskopie zur Charakterisierung mesostrukturierter Materialien. Die Korrelation ihrer mechanischen und dielektrischen Eigenschaften mit ihrer topographischen Beschaffenheit erlaubt ein tieferes Verständnis der mesoskopischen Materialien und ihres Verhaltens auf die Einwirkung äußerer Stimuli.rn

Micropreparative isoelectric focusing protein separation in a suspended drop

IEF protein binary separations were performed in a 12-μL drop suspended between two palladium electrodes, using pH gradients created by electrolysis of simple buffers at low voltages (1.5-5 V). The dynamics of pH gradient formation and protein separation were investigated by computer simulation and experimentally via digital video microscope imaging in the presence and absence of pH indicator solution. Albumin, ferritin, myoglobin, and cytochrome c were used as model proteins. A drop containing 2.4 μg of each protein was applied, electrophoresed, and allowed to evaporate until it splits to produce two fractions that were recovered by rinsing the electrodes with a few microliters of buffer. Analysis by gel electrophoresis revealed that anode and cathode fractions were depleted from high pI and low pI proteins, respectively, whereas proteins with intermediate pI values were recovered in both fractions. Comparable data were obtained with diluted bovine serum that was fortified with myoglobin and cytochrome c.

Purification and Preliminary Characterization of the TM0727 Protein from Thermotoga maritima

The TM0727 gene of Thermotoga maritima is responsible for encoding what has been reported to be a modulator of DNA gyrase (pmbA). Although the function of pmbA is still unknown, it is believedto be involved in cell division, carbon storage regulation, and the synthesis of the antibiotic peptide microcin B17. It is suggested that it serves together with tldD, a known zinc dependent protease, tomodulate DNA gyrase. TM0727 is believed to be a zinc dependent protease that binds zinc in the central active site of the molecule, located between two equivalent monomeric units. However, thecrystal structure determined by Wilson et al. (2005) did not contain zinc. It therefore remains to be seen if TM0727 requires zinc for activity, or regulation, and if the protein is indeed a protease. To begin studying this protein, the gene was expressed in BL21(DE3) pLysS cells and the induction time was optimized. Using affinity and ion exchange chromatography, the protein has been successfully purified. The purification procedure can be replicated to obtain sufficient protein for characterization. Purification results show that the protein loses stability after 24 hours and remains stable under an imidazole-free lysis workup. Preliminary characterization of TM0727 has focused on understanding the protein’s structuralproperties through tryptophan fluorescence anisotropy measurements. The four tryptophan residues located within the TM0727 dimer fluoresce at different maximum wavelengths and with differentintensities upon excitation with 295nm light. These emission properties are highly sensitive to the environment (solvent, surrounding residues) of each tryptophan residue. The low number oftryptophans allows for a specific monitoring of the protein’s structure as it denatures. As more denaturant is added to the protein, its tryptophan environments have clearly altered. This is indicative of unfolding and increased solvent exposure of the protein. This unfolding has been confirmed with the addition of a fluorescent quencher. Additionally, fluorescence anisotropy measurements have been carried out on the protein to gain a preliminary understanding of the rotational dynamics of the tryptophan residues. These experiments excite the tryptophan residues within the sample using a polarized light source. Polarized emission is then detected, the degree of which depends on the rotational dynamics and local environment of the tryptophan residues. The protein was denatured and the changes in emission were recorded to detect these structural changes. Results have shown a large change in quaternary structure, consistent with a dimer to monomer transition, occurs at 1.5M Guandidine HCl. There has also been an examination of the crystal structure for the location of a potential active site. The inner cavity of the protein was inspected visually to locate a potential location for a catalytic triad, specifically the amino acids found in the active sites of serine, cyteine, and aspartateproteases. It was found that a potential aspartic protease active site may be located between the Asparate286 and Aspartate287 residues. Further investigation is warranted to test this remotepossibility.

A novel computer simulation method for simulating the multiscale transduction dynamics of signal proteins

Signal proteins are able to adapt their response to a change in the environment, governing in this way a broad variety of important cellular processes in living systems. While conventional molecular-dynamics (MD) techniques can be used to explore the early signaling pathway of these protein systems at atomistic resolution, the high computational costs limit their usefulness for the elucidation of the multiscale transduction dynamics of most signaling processes, occurring on experimental timescales. To cope with the problem, we present in this paper a novel multiscale-modeling method, based on a combination of the kinetic Monte-Carlo- and MD-technique, and demonstrate its suitability for investigating the signaling behavior of the photoswitch light-oxygen-voltage-2-Jα domain from Avena Sativa (AsLOV2-Jα) and an AsLOV2-Jα-regulated photoactivable Rac1-GTPase (PA-Rac1), recently employed to control the motility of cancer cells through light stimulus. More specifically, we show that their signaling pathways begin with a residual re-arrangement and subsequent H-bond formation of amino acids near to the flavin-mononucleotide chromophore, causing a coupling between β-strands and subsequent detachment of a peripheral α-helix from the AsLOV2-domain. In the case of the PA-Rac1 system we find that this latter process induces the release of the AsLOV2-inhibitor from the switchII-activation site of the GTPase, enabling signal activation through effector-protein binding. These applications demonstrate that our approach reliably reproduces the signaling pathways of complex signal proteins, ranging from nanoseconds up to seconds at affordable computational costs.

EXPLORING CHANGES IN DETRITAL FLOCCULENT LAYER DYNAMICS DUE TO SHIFTS IN MACROPHYTE COMMUNITIES IN THE NORTHERN EVERGLADES

A shift in plant communities of the Water Conservation Areas (WCAs) within the Everglades has been linked to changes in hydrology and high levels of nutrient loading from surrounding agicultural areas. This has resulted in the encroachment of dense cattail stands (Typha domingensis) into areas that had previously been a ridge and slough landscape populated primarily by native sawgrass (Cladium jamaicense). In order to study ecological management solutions in this area, WCA-2A was broken into study plots; several of which became open water areas through the application of herbicide and burning regimens. The open water areas allowed for Chara spp (a submersed algal species) to replace Typha domingensis as the dominant macrophyte. This study investigated the polymer and ionic profiles of Chara spp, Typha domingensis and Cladium jamaicense and their contributions to detrital flocculent (floc) in the study plots where they are the dominant macrophytes. Floc is not only an important food source for aquatic species; it also supports many algal, fungal and bacterial communities. Data gathered in this study indicated that the floc sample from a phosphorus enriched open water study plot (EO1) where Chara spp was the dominant macrophyte may contain cell wall polymers from sources other than Chara spp (most likely Typha domingensis), while the chemical and polymeric profile of the floc of the study plot where Typha domingensis is the dominant macrophyte (EC1) suggests that the floc layer has contributions from algal sources as well as Typha domingensis. Additionally, monoclonal antibodies to Arabinoglalactan protein (AGP) and (1,4)-β-D galactan were identified as possible biomarkers for distinguishing algal dominated floc layers from layers dominated by emergent vegetation. Calcium labeling could be a useful tool for this as well because of the high amount of Ca2+ associated with Chara spp cell walls. When looking into the soluble phosphorus content of the macrophytes and paired floc samples of WCA-2A, it was found that Chara spp may be contributing a greater amount of Ca-bound phosphorus to floc layers where it is the dominant macrophyte when compared to floc layers from study plots dominated by emergent macrophytes. Floc layers also appear to be acting as a nutrient sink for soluble phosphorus. The findings of this study support the overall hypothesis that the shift from native emergent macrophyte communities to submersed macrophyte communities in study sites of the northern Everglades is affecting the polymeric/chemical profile and ionic content of detrital floc layers. The effects of this shift may contribute to changes in complex flocculent community dynamics.

Recruitment of EB1, a master regulator of microtubule dynamics, to the surface of the Theileria annulata schizont

The apicomplexan parasite Theileria annulata transforms infected host cells, inducing uncontrolled proliferation and clonal expansion of the parasitized cell population. Shortly after sporozoite entry into the target cell, the surrounding host cell membrane is dissolved and an array of host cell microtubules (MTs) surrounds the parasite, which develops into the transforming schizont. The latter does not egress to invade and transform other cells. Instead, it remains tethered to host cell MTs and, during mitosis and cytokinesis, engages the cell's astral and central spindle MTs to secure its distribution between the two daughter cells. The molecular mechanism by which the schizont recruits and stabilizes host cell MTs is not known. MT minus ends are mostly anchored in the MT organizing center, while the plus ends explore the cellular space, switching constantly between phases of growth and shrinkage (called dynamic instability). Assuming the plus ends of growing MTs provide the first point of contact with the parasite, we focused on the complex protein machinery associated with these structures. We now report how the schizont recruits end-binding protein 1 (EB1), a central component of the MT plus end protein interaction network and key regulator of host cell MT dynamics. Using a range of in vitro experiments, we demonstrate that T. annulata p104, a polymorphic antigen expressed on the schizont surface, functions as a genuine EB1-binding protein and can recruit EB1 in the absence of any other parasite proteins. Binding strictly depends on a consensus SxIP motif located in a highly disordered C-terminal region of p104. We further show that parasite interaction with host cell EB1 is cell cycle regulated. This is the first description of a pathogen-encoded protein to interact with EB1 via a bona-fide SxIP motif. Our findings provide important new insight into the mode of interaction between Theileria and the host cell cytoskeleton.

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.

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.

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.