A system for studying evolution of life-like virtual organisms

BACKGROUND: Fitness landscapes, the dependences of fitness on the genotype, are of critical importance for the evolution of living beings. Unfortunately, fitness landscapes that are relevant to the evolution of complex biological functions are very poorly known. As a result, the existing theory of evolution is mostly based on postulated fitness landscapes, which diminishes its usefulness. Attempts to deduce fitness landscapes from models of actual biological processes led, so far, to only limited success. RESULTS: We present a model system for studying the evolution of biological function, which makes it possible to attribute fitness to genotypes in a natural way. The system mimics a very simple cell and takes into account the basic properties of gene regulation and enzyme kinetics. A virtual cell contains only two small molecules, an organic nutrient A and an energy carrier X, and proteins of five types – two transcription factors, two enzymes, and a membrane transporter. The metabolism of the cell consists of importing A from the environment and utilizing it in order to produce X and an unspecified end product. The genome may carry an arbitrary number of genes, each one encoding a protein of one of the five types. Both major mutations that affect whole genes and minor mutations that affect individual characteristics of genes are possible. Fitness is determined by the ability of the cell to maintain homeostasis when its environment changes. The system has been implemented as a computer program, and several numerical experiments have been performed on it. Evolution of the virtual cells usually involves a rapid initial increase of fitness, which eventually slows down, until a fitness plateau is reached. The origin of a wide variety of genetic networks is routinely observed in independent experiments performed under the same conditions. These networks can have different, including very high, levels of complexity and often include large numbers of non-essential genes. CONCLUSION: The described system displays a rich repertoire of biologically sensible behaviors and, thus, can be useful for investigating a number of unresolved issues in evolutionary biology, including evolution of complexity, modularity and redundancy, as well as for studying the general properties of genotype-to-fitness maps. REVIEWERS: This article was reviewed by Drs. Eugene Koonin, Shamil Sunyaev and Arcady Mushegian

The Bacillus anthracis protein MprF is required for synthesis of lysylphosphatidylglycerols and for resistance to cationic antimicrobial peptides

During inhalational anthrax, Bacillus anthracis survives and replicates in alveolar macrophages, followed by rapid invasion into the host's bloodstream, where it multiplies to cause heavy bacteremia. B. anthracis must therefore defend itself from host immune functions encountered during both the intracellular and the extracellular stages of anthrax infection. In both of these niches, cationic antimicrobial peptides are an essential component of the host's innate immune response that targets B. anthracis. However, the genetic determinants of B. anthracis contributing to resistance to these peptides are largely unknown. Here we generated Tn917 transposon mutants in the ΔANR strain (pXO1(−) pXO2(−)) of B. anthracis and screened them for altered protamine susceptibility. A protamine-sensitive mutant identified carried the transposon inserted in the BA1486 gene encoding a putative membrane protein homologous to MprF known in several gram-positive pathogens. A mutant strain with the BAS1375 gene (the orthologue of BA1486) deleted in the Sterne 34F2 strain (pXO1(+) pXO2(−)) of B. anthracis exhibited hypersusceptibility not only to protamine but also to α-helical cathelicidin LL-37 and β-sheet defensin human neutrophil peptide 1 compared to the wild-type Sterne strain. Analysis of membrane lipids using isotopic labeling demonstrated that the BAS1375 deletion mutant is unable to synthesize lysinylated phosphatidylglycerols, and this defect is rescued by genetic complementation. Further, we determined the structures of these lysylphosphatidylglycerols by using various mass spectrometric analyses. These results demonstrate that in B. anthracis a functional MprF is required for the biosynthesis of lysylphosphatidylglycerols, which is critical for resistance to cationic antimicrobial peptides

Crystallization and Preliminary X-ray Diffraction Studies on the DNA-Binding Domain of the Multidrug Transporter Activation Protein (MtaN) from Bacillus subtilis

The N-terminal DNA-binding domain of the multidrug transporter activation protein (MtaN) was crystallized by the hanging-drop vapour-diffusion method using lithium chloride as a precipitant. The crystals are orthorhombic and belong to the space group I212121, with unit-cell parameters a = 49.4, b = 67.8, c = 115.0 Å. Diffraction data have been collected at 100 K to 2.75 Å resolution at a synchrotron-radiation source

Crystal Structure of MtaN, a Global Multidrug Transporter Gene Activator

MtaN (Multidrug Transporter Activation, N terminus) is a constitutive, transcriptionally active 109-residue truncation mutant, which contains only the N-terminal DNA-binding and dimerization domains of MerR family member Mta. The 2.75 Å resolution crystal structure of apo-MtaN reveals a winged helix-turn-helix protein with a protruding 8-turn helix (α5) that is involved in dimerization by the formation of an antiparallel coiled-coil. The hydrophobic core and helices α1 through α4 are structurally homologous to MerR family member BmrR bound to DNA, whereas one wing (Wing 1) is shifted. Differences between the orientation of α5 with respect to the core and the revolution of the antiparallel coiled-coil lead to significantly altered conformations of MtaN and BmrR dimers. These shifts result in a conformation of MtaN that appears to be incompatible with the transcription activation mechanism of BmrR and suggest that additional DNA-induced structural changes are necessary

Monitoring of gene knockouts: genome-wide profiling of conditionally essential genes

Monitoring of gene knockouts is a new microarray-based genetic technique used for genome-wide identification of conditionally essential genes in bacteri

Nucleotide Biosynthesis Is Critical for Growth of Bacteria in Human Blood

Proliferation of bacterial pathogens in blood represents one of the most dangerous stages of infection. Growth in blood serum depends on the ability of a pathogen to adjust metabolism to match the availability of nutrients. Although certain nutrients are scarce in blood and need to be de novo synthesized by proliferating bacteria, it is unclear which metabolic pathways are critical for bacterial growth in blood. In this study, we identified metabolic functions that are essential specifically for bacterial growth in the bloodstream. We used two principally different but complementing techniques to comprehensively identify genes that are required for the growth of Escherichia coli in human serum. A microarray-based and a dye-based mutant screening approach were independently used to screen a library of 3,985 single-gene deletion mutants in all non-essential genes of E. coli (Keio collection). A majority of the mutants identified consistently by both approaches carried a deletion of a gene involved in either the purine or pyrimidine nucleotide biosynthetic pathway and showed a 20- to 1,000-fold drop in viable cell counts as compared to wild-type E. coli after 24 h of growth in human serum. This suggests that the scarcity of nucleotide precursors, but not other nutrients, is the key limitation for bacterial growth in serum. Inactivation of nucleotide biosynthesis genes in another Gram-negative pathogen, Salmonella enterica, and in the Gram-positive pathogen Bacillus anthracis, prevented their growth in human serum. The growth of the mutants could be rescued by genetic complementation or by addition of appropriate nucleotide bases to human serum. Furthermore, the virulence of the B. anthracis purE mutant, defective in purine biosynthesis, was dramatically attenuated in a murine model of bacteremia. Our data indicate that de novo nucleotide biosynthesis represents the single most critical metabolic function for bacterial growth in blood and reveal the corresponding enzymes as putative antibiotic targets for the treatment of bloodstream infections

Mesenchymal Stromal Cells Primed with Paclitaxel Provide a New Approach for Cancer Therapy

BACKGROUND: Mesenchymal stromal cells may represent an ideal candidate to deliver anti-cancer drugs. In a previous study, we demonstrated that exposure of mouse bone marrow derived stromal cells to Doxorubicin led them to acquire anti-proliferative potential towards co-cultured haematopoietic stem cells (HSCs). We thus hypothesized whether freshly isolated human bone marrow Mesenchymal stem cells (hMSCs) and mature murine stromal cells (SR4987 line) primed in vitro with anti-cancer drugs and then localized near cancer cells, could inhibit proliferation. METHODS AND PRINCIPAL FINDINGS: Paclitaxel (PTX) was used to prime culture of hMSCs and SR4987. Incorporation of PTX into hMSCs was studied by using FICT-labelled-PTX and analyzed by FACS and confocal microscopy. Release of PTX in culture medium by PTX primed hMSCs (hMSCsPTX) was investigated by HPLC. Culture of Endothelial cells (ECs) and aorta ring assay were used to test the anti-angiogenic activity of hMSCsPTX and PTX primed SR4987(SR4987PTX), while anti-tumor activity was tested in vitro on the proliferation of different tumor cell lines and in vivo by co-transplanting hMSCsPTX and SR4987PTX with cancer cells in mice. Nevertheless, despite a loss of cells due to chemo-induced apoptosis, both hMSCs and SR4987 were able to rapidly incorporate PTX and could slowly release PTX in the culture medium in a time dependent manner. PTX primed cells acquired a potent anti-tumor and anti-angiogenic activity in vitro that was dose dependent, and demonstrable by using their conditioned medium or by co-culture assay. Finally, hMSCsPTX and SR4987PTX co-injected with human cancer cells (DU145 and U87MG) and mouse melanoma cells (B16) in immunodeficient and in syngenic mice significantly delayed tumor takes and reduced tumor growth. CONCLUSIONS: These data demonstrate, for the first time, that without any genetic manipulation, mesenchymal stromal cells can uptake and subsequently slowly release PTX. This may lead to potential new tools to increase efficacy of cancer therapy

The Roots of Bioinformatics in Theoretical Biology

From the late 1980s onward, the term “bioinformatics” mostly has been used to refer to computational methods for comparative analysis of genome data. However, the term was originally more widely defined as the study of informatic processes in biotic systems. In this essay, I will trace this early history (from a personal point of view) and I will argue that the original meaning of the term is re-emerging

Sensitive and Specific Fluorescent Probes for Functional Analysis of the Three Major Types of Mammalian ABC Transporters

An underlying mechanism for multi drug resistance (MDR) is up-regulation of the transmembrane ATP-binding cassette (ABC) transporter proteins. ABC transporters also determine the general fate and effect of pharmaceutical agents in the body. The three major types of ABC transporters are MDR1 (P-gp, P-glycoprotein, ABCB1), MRP1/2 (ABCC1/2) and BCRP/MXR (ABCG2) proteins. Flow cytometry (FCM) allows determination of the functional expression levels of ABC transporters in live cells, but most dyes used as indicators (rhodamine 123, DiOC2(3), calcein-AM) have limited applicability as they do not detect all three major types of ABC transporters. Dyes with broad coverage (such as doxorubicin, daunorubicin and mitoxantrone) lack sensitivity due to overall dimness and thus may yield a significant percentage of false negative results. We describe two novel fluorescent probes that are substrates for all three common types of ABC transporters and can serve as indicators of MDR in flow cytometry assays using live cells. The probes exhibit fast internalization, favorable uptake/efflux kinetics and high sensitivity of MDR detection, as established by multidrug resistance activity factor (MAF) values and Kolmogorov-Smirnov statistical analysis. Used in combination with general or specific inhibitors of ABC transporters, both dyes readily identify functional efflux and are capable of detecting small levels of efflux as well as defining the type of multidrug resistance. The assay can be applied to the screening of putative modulators of ABC transporters, facilitating rapid, reproducible, specific and relatively simple functional detection of ABC transporter activity, and ready implementation on widely available instruments

A Single Acidic Residue Can Guide Binding Site Selection but Does Not Govern QacR Cationic-Drug Affinity

Structures of the multidrug-binding repressor protein QacR with monovalent and bivalent cationic drugs revealed that the carboxylate side-chains of E90 and E120 were proximal to the positively charged nitrogens of the ligands ethidium, malachite green and rhodamine 6G, and therefore may contribute to drug neutralization and binding affinity. Here, we report structural, biochemical and in vivo effects of substituting these glutamate residues. Unexpectedly, substitutions had little impact on ligand affinity or in vivo induction capabilities. Structures of QacR(E90Q) and QacR(E120Q) with ethidium or malachite green took similar global conformations that differed significantly from all previously described QacR-drug complexes but still prohibited binding to cognate DNA. Strikingly, the QacR(E90Q)-rhodamine 6G complex revealed two mutually exclusive rhodamine 6G binding sites. Despite multiple structural changes, all drug binding was essentially isoenergetic. Thus, these data strongly suggest that rather than contributing significantly to ligand binding affinity, the role of acidic residues lining the QacR multidrug-binding pocket is primarily to attract and guide cationic drugs to the “best available” positions within the pocket that elicit QacR induction