Using digital electronic design flow to create a Genetic Design Automation tool.

Synthetic bio-systems become increasingly more complex and their development is lengthy and expensive. In the same way, in microelectronics, the design process of very complex circuits has benefited from many years of experience. It is now partly automated through Electronic Design Automation tools. Both areas present analogies that can be used to create a Genetic Design Automation tool inspired from EDA tools used in digital electronics. This tool would allow moving away from a totally manual design of bio-systems to assisted conception. This ambitious project is presented in this paper, with a deep focus on the tool that automatically generates models of bio-systems directly usable in electronic simulators.

Ligand modulation of glial activation: cell permeable, small molecule inhibitors of serine-threonine protein kinases can block induction of interleukin 1 beta and nitric oxide synthase II.

Activated glia (astrocytes and microglia) and their associated neuroinflammatory sequelae have been linked to the disease progression of several neurodegenerative disorders, including Alzheimer's disease. We found that the experimental anti-inflammatory drug K252a, an inhibitor of calmodulin regulated protein kinases (CaMKs), can block induction of both the oxidative stress related enzyme iNOS and the proinflammatory cytokine IL-1 beta in primary cortical glial cultures and the microglial BV-2 cell line. We also found that the profile of CaMKIV and CaMKII isoforms in primary cortical glial cultures and BV-2 cells is distinct from that found in neurons. Knowledge of cellular mechanisms and high throughput screens of a pharmacologically focused chemical library allowed the discovery of novel pyridazine-based compounds that are cell permeable ligand modulators of gene regulating protein kinases involved in the induction of iNOS and IL-1 beta in activated glia. Pyridazine-based compounds are attractive for the development of new therapeutics due to the retention of the remarkable pharmacological properties of K252a and related indolocarbazole alkaloids, and presence of enhanced functional selectivity in a comparatively simple structure amenable to diverse synthetic chemistries.

Mechanism of calcium oscillations in migrating human astrocytoma cells.

Numerous studies show that intracellular calcium controls the migration rate of different mobile cell types. We studied migrating astrocytoma cells from two human cell lines, U-87MG and A172, in order to clarify the mechanisms by which calcium potentially influences cell migration. Using the wound-healing model to assay migration, we showed that four distinct components of migration could be distinguished: (i) a Ca(2+)/serum-dependent process; (ii) a Ca(2+)-dependent/serum-independent process; (iii) a Ca(2+)/serum-independent process; (iv) a Ca(2+)-independent/serum-dependent process. In U-87MG cells which lack a Ca(2+)-dependent/serum-independent component, we found that intracellular Ca(2+) oscillations are involved in Ca(2+)-dependent migration. Removing extracellular Ca(2+) greatly decreased the frequency of migration-associated Ca(2+) oscillations. Furthermore, non-selective inhibition of Ca(2+) channels by heavy metals such as Cd(2+) or La(3+) almost completely abolished changes in intracellular Ca(2+) observed during migration, indicating an essential role for Ca(2+) channels in the generation of these Ca(2+) oscillations. However, specific blockers of voltage-gated Ca(2+) channels, including nitrendipine, omega-conotoxin GVIA, omega-conotoxin MVIIC or low concentrations of Ni(2+) were without effect on Ca(2+) oscillations. We examined the role of internal Ca(2+) stores, showing that thapsigargin-sensitive Ca(2+) stores and InsP(3) receptors are involved in Ca(2+) oscillations, unlike ryanodine-sensitive Ca(2+) stores. Detailed analysis of the spatio-temporal aspect of the Ca(2+) oscillations revealed the existence of Ca(2+) waves initiated at the leading cell edge which propagate throughout the cell. Previously, we have shown that the frequency of Ca(2+) oscillations was reduced in the presence of inhibitory antibodies directed against beta3 integrin subunits. A simple model of a Ca(2+) oscillator is proposed, which may explain how the generation of Ca(2+) oscillations is linked to cell migration.

Analysis of the kinase-related protein gene found at human chromosome 3q21 in a multi-gene cluster: organization, expression, alternative splicing, and polymorphic marker.

We report the amino acid sequence, genomic organization, tissue expression, and alternative splice patterns for the human kinase related protein (KRP) gene, as well as the discovery of a new CA repeat sequence polymorphic marker in an upstream intron of the myosin light chain kinase (MLCK) gene. The KRP/MLCK genetic locus is a prototype for a recently discovered paradigm in which an independently regulated gene for a non-enzymic protein is embedded within a larger gene for a signal transduction enzyme, and both classes of proteins are involved in the regulation of the same cellular structure. The MLCK/KRP gene cluster has been found only in higher vertebrates and is localized to human chromosome 3q21. The determination of the human KRP amino acid sequence through cDNA sequence analysis and its comparison to the exon/intron organization of the human KRP gene revealed an alternative splice pattern at the start of KRP exon 2, resulting in the insertion of a single glutamic acid in the middle of the protein. Examination of tissue distribution using Northern blot analysis revealed that the human expression pattern is more similar to the well-characterized chicken KRP gene expression pattern than to rodent or rabbit. Unexpected differences of the human gene from other species is the apparent expression of the human gene products in adult cardiac muscle, an observation that was pursued further by the production of a site-directed antiserum and immunohistochemistry analysis. The results reported here provide insight into the conserved and variable features of this late evolving genetic paradigm, raise new questions about the molecular aspects of cardiac muscle regulation, and provide tools needed for future clinical studies. The comparative analysis of the MLCK/KRP locus, combined with the recent discovery of a similar genomic relationship among other signal transduction proteins, suggest a diverse distribution of this theme among signal transduction systems in higher vertebrate genomes and indicate the utility of comparative genomics in revealing late evolving genetic paradigms.

Apocalmodulin binds to the myosin light chain kinase calmodulin target site.

The interaction of a 20-residue-long peptide derived from the calmodulin-binding domain of the smooth muscle myosin light chain kinase with calcium-free calmodulin (apocalmodulin) was studied using a combination of isothermal titration calorimetry and differential scanning calorimetry. We showed that: (i) a significant binding between apocalmodulin and the target peptide (RS20) exists in the absence of salt (Ka = 10(6) M-1), (ii) the peptide interacts with the C-terminal lobe of calmodulin and adopts a partly helical conformation, and (iii) the presence of salt weakens the affinity of the peptide for apocalmodulin, emphasizing the importance of electrostatic interactions in the complex. Based on these results and taking into account the work of Bayley et al. (Bayley, P. M., Findlay, W.A., and Martin, S. R. (1996) Protein Sci. 5, 1215-1228), we suggest a physiological role for apocalmodulin.

Thermodynamic analysis of calcium and magnesium binding to calmodulin.

To elucidate some aspects still debated concerning the interaction of Ca2+ and Mg2+ with CaM, the thermodynamic binding parameters of Ca2+-CaM and Mg2+-CaM complexes were characterized by flow dialysis and isothermal microcalorimetry under different experimental conditions. In particular, the enthalpy and entropy changes associated with Ca2+ and Mg2+ binding to their sites were determined, allowing a better understanding of the mechanism underlying cation-CaM interactions. Ca2+-CaM interaction follows an enthalpy-entropy compensation relationship, suggesting that CaM explores a subspace of isoenergetical conformations which is modified by Ca2+ binding. This Ca2+-induced change in CaM dynamics is proposed to play a key role in CaM function, i.e. in its interaction with and/or activation of target proteins. Furthermore, data show that Mg2+ does not act as a direct competitor for Ca2+ binding on the four main Ca2+ binding sites, but rather as an allosteric effector. This implies that the four main Mg2+ binding sites are distinct from the EF-hand Ca2+ binding sites. Finally, Ca2+ is shown to interact with auxiliary binding sites on CaM. These weak affinity sites were thermodynamically characterized. The results presented here challenge the current accepted view of CaM ion binding.

New protein kinase and protein phosphatase families mediate signal transduction in bacterial catabolite repression.

Carbon catabolite repression (CCR) is the prototype of a signal transduction mechanism. In enteric bacteria, cAMP was considered to be the second messenger in CCR by playing a role reminiscent of its actions in eukaryotic cells. However, recent results suggest that CCR in Escherichia coli is mediated mainly by an inducer exclusion mechanism. In many Gram-positive bacteria, CCR is triggered by fructose-1,6-bisphosphate, which activates HPr kinase, presumed to be one of the most ancient serine protein kinases. We here report cloning of the Bacillus subtilis hprK and hprP genes and characterization of the encoded HPr kinase and P-Ser-HPr phosphatase. P-Ser-HPr phosphatase forms a new family of phosphatases together with bacterial phosphoglycolate phosphatase, yeast glycerol-3-phosphatase, and 2-deoxyglucose-6-phosphate phosphatase whereas HPr kinase represents a new family of protein kinases on its own. It does not contain the domain structure typical for eukaryotic protein kinases. Although up to now the HPr modifying/demodifying enzymes were thought to exist only in Gram-positive bacteria, a sequence comparison revealed that they also are present in several Gram-negative pathogenic bacteria.

Calcium signalling in Bacillus subtilis.

Few systematic studies have been devoted to investigating the role of Ca2+ as an intracellular messenger in prokaryotes. Here we report an investigation on the potential involvement of Ca2+ in signalling in Bacillus subtilis, a Gram-positive bacterium. Using aequorin, it is shown that B. subtilis cells tightly regulate intracellular Ca2+ levels. This homeostasis can be changed by an external stimulus such as hydrogen peroxide, pointing to a relationship between oxidative stress and Ca2+ signalling. Also, B. subtilis growth appears to be intimately linked to the presence of Ca2+, as normal growth can be immediately restored by adding Ca2+ to an almost non-growing culture in EGTA containing Luria broth medium. Addition of Fe2+ or Mn2+ also restores growth, but with 5-6 h delay, whereas Mg2+ did not have any effect. In addition, the expression of alkyl hydroperoxide reductase C (AhpC), which is strongly enhanced in bacteria grown in the presence of EGTA, also appears to be regulated by Ca2+. Finally, using 45Ca2+ overlay on membrane electrotransferred two-dimensional gels of B. subtilis, four putative Ca2+ binding proteins were found, including AhpC. Our results provide strong evidence for a regulatory role for Ca2+ in bacterial cells.

Isolation and characterization of the lacA gene encoding beta-galactosidase in Bacillus subtilis and a regulator gene, lacR.

We have isolated transposon insertions in the lacA gene encoding an endogenous beta-galactosidase of Bacillus subtilis. Upstream of the putative operon containing lacA is a negative regulator, lacR, which encodes a product related to a family of regulators that includes the lactose repressor, lacI, of Escherichia coli. New strains with insertions in the lacA gene should be of use in studies using lacZ fusions in B. subtilis.

The Bacillus subtilis crh gene encodes a HPr-like protein involved in carbon catabolite repression.

Carbon catabolite repression (CCR) of several Bacillus subtilis catabolic genes is mediated by ATP-dependent phosphorylation of histidine-containing protein (HPr), a phosphocarrier protein of the phosphoenolpyruvate (PEP): sugar phosphotransferase system. In this study, we report the discovery of a new B. subtilis gene encoding a HPr-like protein, Crh (for catabolite repression HPr), composed of 85 amino acids. Crh exhibits 45% sequence identity with HPr, but the active site His-15 of HPr is replaced with a glutamine in Crh. Crh is therefore not phosphorylated by PEP and enzyme I, but is phosphorylated by ATP and the HPr kinase in the presence of fructose-1,6-bisphosphate. We determined Ser-46 as the site of phosphorylation in Crh by carrying out mass spectrometry with peptides obtained by tryptic digestion or CNBr cleavage. In a B. subtilis ptsH1 mutant strain, synthesis of beta-xylosidase, inositol dehydrogenase, and levanase was only partially relieved from CCR. Additional disruption of the crh gene caused almost complete relief from CCR. In a ptsH1 crh1 mutant, producing HPr and Crh in which Ser-46 is replaced with a nonphosphorylatable alanyl residue, expression of beta-xylosidase was also completely relieved from glucose repression. These results suggest that CCR of certain catabolic operons requires, in addition to CcpA, ATP-dependent phosphorylation of Crh, and HPr at Ser-46.

Conformation and thermal denaturation of apocalmodulin: role of electrostatic mutations.

Scanning microcalorimetry and circular dichroism were used to study conformational state and heat denaturation of Ca2+-free synthetic calmodulin (SynCaM) and three charge reversal mutants. We produced evidence for the major role of the electrostatic potential in the stability and flexibility of SynCaM. The substitution of 118DEE120 by 118KKK120 (SynCaM12A) does not influence the flexibility of the protein; the replacement of 82EEE84 by 82KKK84 (SynCaM8) decreases its level, while the combination of these two mutations in SynCaM18A significantly increases the flexibility. The heat denaturation of apoSynCaM and its mutants is well approximated by two two-state transitions with the lower-temperature transition corresponding to C-terminal lobe melting and the higher-temperature one to N-terminal lobe melting. The difference in transition temperatures for the two lobes decreases in SynCaM8 and increases in SynCaM18A, suggesting a modification in the influence of one lobe to the other. The electrostatic mutations change the parameters of thermal denaturation of SynCaM lobes in a similar way as pH conditions affect thermal transition parameters of multidomain proteins, leading to a linear temperature dependence of transition enthalpy. One domain of the N-terminal lobe in apoSynCaM18A is unfolded in the native state. Near-UV CD spectra point out the invariability of the local structure of aromatic residues upon mutations, although the secondary structure undergoes striking transformations. Cacodylate ions strongly and specifically alter the helical content of SynCaM. Our data unambiguously demonstrate that the two lobes are not independent, and interactions between the lobes are mediated by the electrostatic potential of the molecule.

Integrated mapping and sequencing of a 115 kb DNA fragment from Bacillus subtilis: sequence analysis of a 21 kb segment containing the sigL locus.

A sequence strategy which combines a low redundancy shotgun approach and directed sequencing has been elaborated. Essentially, the sequences, as well as the size of the fragments utilized for a low coverage shotgun approach, were exploited for the construction of a physical map of the region to be sequenced. The latter considerably simplified the subsequent directed sequencing steps. We report the physical mapping of a 115 kb segment which covers nearly 100 kb of the hisA-cysB region of the Bacillus subtilis chromosome and contains previously sequenced genes sigL and sacB. Sequencing and analysis of a 21305 bp segment, which includes the sigL locus, revealed 21 ORFs, apparently belonging to at least seven transcription units. This segment has a G + C content greater than 47%, compared to 43% characteristic of the flanking regions, and mainly consists of genes whose products seem to be involved in the synthesis of an exopolysaccharide. These observations leave open the possibility that the analysed fragment has been acquired through horizontal transfer.

Rational choice of molecular dynamics simulation parameters through the use of the three-dimensional autocorrelation method: application to calmodulin flexibility study.

We examined the effects of several adjustable parameters for use in molecular dynamics simulations of proteins using both standard criteria (radius of gyration, root mean square deviation from starting coordinates, molecular mechanics energy) and a new description of protein conformations by 3-D autocorrelation vectors (3-D ACV). We chose calmodulin (CaM) as a protein model and analysed 23 simulations using different combinations of the four molecular dynamics parameters studied, such as the dielectric constant (epsilon), the heating phase time (H), the thermal bath coupling time (zeta T) and the time step size (delta t). The correctness of the various trajectories generated with different parameter sets was evaluated through geometric analysis and use of a knowledge-based profile method. It is shown that 3-D ACV combined with multivariate statistical analysis provides a convenient way to describe and compare molecular dynamics simulations and constitutes a valuable complementary tool to standard methods. Using these methods, comparison of the various simulations performed on CaM indicated that the best in vacuo parameter set was epsilon = 1 x r, H = 15 ps, zeta T = 0.1 ps and delta t = 1 fs in fairly good agreement with previous less extensive comparisons of molecular dynamics trajectories.

Heme-CO binding to tryptophan-containing calmodulin mutants.

The binding of heme-CO to genetically engineered calmodulin containing a single tryptophan residue has been studied. A tryptophan residue was integrated at one of five positions: 26 or 62 of the N-terminal, 81 in the central helix, or 99 or 135 of the C-terminal. As for the wild type, the mutant calmodulins bind four molecules of heme-CO with an average affinity of 1 microM. (i) Homotropic effect. The quenching of the tryptophan fluorescence by energy transfer to the hemes indicates that there is no preference between the N- or C-terminal pockets for heme binding. The quenching is less than expected for a binomial distribution of four sites. This could indicate a lower energy transfer rate due to a specific orientation factor. The weak quenching as a function of the number of hemes bound may also reveal a cooperativity in the heme binding; the data can be simulated assuming two pairs of sites, where each pocket shows a cooperative binding for two hemes. (ii) Heterotropic effect. As observed for the wild type, addition of melittin does not displace the hemes from the mutant calmodulins; the affinity of heme-CO for the calmodulin.melittin complex is higher than that for calmodulin alone. The affinity of heme-CO for native calmodulin is also higher in the presence of trifluoperazine.

Cloning, expression, and properties of the microtubule-stabilizing protein STOP.

Nerve cells contain abundant subpopulations of cold-stable microtubules. We have previously isolated a calmodulin-regulated brain protein, STOP (stable tubule-only polypeptide), which reconstitutes microtubule cold stability when added to cold-labile microtubules in vitro. We have now cloned cDNA encoding STOP. We find that STOP is a 100.5-kDa protein with no homology to known proteins. The primary structure of STOP includes two distinct domains of repeated motifs. The central region of STOP contains 5 tandem repeats of 46 amino acids, 4 with 98% homology to the consensus sequence. The STOP C terminus contains 28 imperfect repeats of an 11-amino acid motif. STOP also contains a putative SH3-binding motif close to its N terminus. In vitro translated STOP binds to both microtubules and Ca2+-calmodulin. When STOP cDNA is expressed in cells that lack cold-stable microtubules, STOP associates with microtubules at 37 degrees C, and stabilizes microtubule networks, inducing cold stability, nocodazole resistance, and tubulin detyrosination on microtubules in transfected cells. We conclude that STOP must play an important role in the generation of microtubule cold stability and in the control of microtubule dynamics in brain.

Analysis of the ion binding sites of calmodulin by electrospray ionization mass spectrometry.

The binding of Ca2+ and Mg2+ to four calmodulins (SynCaM 1, SynCaM 8, SynCaM 12A, and SynCaM 18A) has been studied by ESI-MS. The mass spectra were recorded by dissolving the apoproteins in methanol/water (20/80, v/v) containing 1 mM CaCl2 or 1 mM MgCl2 and the pH adjusted to 6.0 with ammonia. The carrier solvent was methanol/water (20/80, v/v). In the case of Ca2+ complexation, ESI-MS reveals the presence of three kinds of sites: the first of high affinity corresponding to those determined using flow and equilibrium dialysis techniques and two others with lower affinities. These results clearly confirm the conclusion of Milos et al. [Milos, M., Comte, M., Schaer, J. J., & Cox, J. A. (1989) J. Inorg. Biochem. 36, 11-25] that there should exist between four and six auxiliary sites for Ca2+. Concerning the complexation of magnesium, the four proteins are able to bind two Mg2+ almost certainly on auxiliary cationic sites.