Simulating properties of in vitro epithelial cell morphogenesis

How do individual epithelial cells (ECs) organize into multicellular structures? ECs are studied in vitro to help answer that question. Characteristic growth features include stable cyst formation in embedded culture, inverted cyst formation in suspension culture, and lumen formation in overlay culture. Formation of these characteristic structures is believed to be a consequence of an intrinsic program of differentiation and de-differentiation. To help discover how such a program may function, we developed an in silico analogue in which space, events, and time are discretized. Software agents and objects represent cells and components of the environment. "Cells" act independently. The "program" governing their behavior is embedded within each in the form of axioms and an inflexible decisional process. Relationships between the axioms and recognized cell functions are specified. Interactions between "cells" and environment components during simulation give rise to a complex in silico phenotype characterized by context-dependent structures that mimic counterparts observed in four different in vitro culture conditions: a targeted set of in vitro phenotypic attributes was matched by in silico attributes. However, for a particular growth condition, the analogue failed to exhibit behaviors characteristic of functionally polarized ECs. We solved this problem by following an iterative refinement method that improved the first analogue and led to a second: it exhibited characteristic differentiation and growth properties in all simulated growth conditions. It is the first model to simultaneously provide a representation of nonpolarized and structurally polarized cell types, and a mechanism for their interconversion. The second analogue also uses an inflexible axiomatic program. When specific axioms are relaxed, growths strikingly characteristic of cancerous and precancerous lesions are observed. In one case, the simulated cause is aberrant matrix production. Analogue design facilitates gaining deeper insight into such phenomena by making it easy to replace low-resolution components with increasingly detailed and realistic components.     

Synthesis of cyclic ADP-carbocylcic-xylose and its 3"-O-methyl analogue as stable and potent Ca2+ -mobilizing agents

We previously showed that 3"-deoxy-cyclic ADP-carbocyclic-ribose (3"-deoxy-cADPcR, 3) is a stable and highly potent analogue of cyclic ADP-ribose (cADPR, 1), a Ca2+ -mobilizing second messenger. From these results, we newly designed another 3"-modified analogues of cADPcR and identified the N1-"xylo"-type carbocyclic analogue, i.e., cADPcX (4), as one of the most potent cADPR-related compounds reported so far.     

Experimental study of informal rewards in peer production

We test the effects of informal rewards in online peer production. Using a randomized, experimental design, we assigned editing awards or "barnstars" to a subset of the 1% most productive Wikipedia contributors. Comparison with the control group shows that receiving a barnstar increases productivity by 60% and makes contributors six times more likely to receive additional barnstars from other community members, revealing that informal rewards significantly impact individual effort.     

The number of loci that affect milk production traits in dairy cattle

We have used the results of an experiment mapping quantitative trait loci (QTL) affecting milk yield and composition to estimate the total number of QTL affecting these traits. We did this by estimating the number of segregating QTL within a half-sib daughter design using logic similar to that used to estimate the "false discovery rate" (FDR). In a half-sib daughter design with six sire families we estimate that the average sire was heterozygous for approximately 5 QTL per trait. Also, in most cases only one sire was heterozygous for any one QTL; therefore at least 30 QTL were likely to be segregating for these milk production traits in this Holstein population.     

Exploration of the binding proteins of perfluorooctane sulfonate by a T7 phage display screen

Perfluorooctane sulfonate (PFOS) is a pollutant widely found throughout nature and is toxic to animals. We created a PFOS analogue on a polyethylene glycol polyacrylamide copolymer and isolated peptides that preferentially bound the PFOS analogue using a T7 phage display system. Bioinformatic analysis using the FASTAskan program on the RELIC bioinformatics server showed several human proteins that likely bound PFOS. Among them, we confirmed binding between PFOS and a recombinant soluble form of monocyte differentiation antigen CD14 (sCD14) by a surface plasmon biosensor. Furthermore, PFOS inhibited TNF-α production induced by the sCD14 in mouse macrophage RAW264.7 cells.     

Parameter estimation and determinability analysis applied to Drosophila gap gene circuits

Background

Our objective was to discover in silico axioms that are plausible representations of the operating principles realized during characteristic growth of EMT6/Ro mouse mammary tumor spheroids in culture. To reach that objective we engineered and iteratively falsified an agent-based analogue of EMT6 spheroid growth. EMT6 spheroids display consistent and predictable growth characteristics, implying that individual cell behaviors are tightly controlled and regulated. An approach to understanding how individual cell behaviors contribute to system behaviors is to discover a set of principles that enable abstract agents to exhibit closely analogous behaviors using only information available in an agent's immediate environment. We listed key attributes of EMT6 spheroid growth, which became our behavioral targets. Included were the development of a necrotic core surrounded by quiescent and proliferating cells, and growth data at two distinct levels of nutrient.

Results

We then created an analogue made up of quasi-autonomous software agents and an abstract environment in which they could operate. The system was designed so that upon execution it could mimic EMT6 cells forming spheroids in culture. Each agent used an identical set of axiomatic operating principles. In sequence, we used the list of targeted attributes to falsify and revise these axioms, until the analogue exhibited behaviors and attributes that were within prespecified ranges of those targeted, thereby achieving a level of validation.

Conclusion

The finalized analogue required nine axioms. We posit that the validated analogue's operating principles are reasonable representations of those utilized by EMT6/Ro cells during tumor spheroid development.     

Modification of the phylloquinone in the A1 binding site in photosystem I studied using time-resolved FTIR difference spectroscopy and density functional theory

A phylloquinone molecule (2-methyl-3-phytyl-1,4-naphthoquinone) occupies the A1 binding site in photosystem I. Previously, we have obtained A1(-)/A1 FTIR difference spectra using labeled and unlabeled photosystem I particles and proposed assignments for many of the bands in the spectra [Sivakumar, V., Wang, R., and Hastings, G. (2005) Biochemistry 44, 1880-1893]. In particular, we suggested that a negative/positive band at 1654/1495 cm(-1) in A1(-)/A1 FTIR DS is due to a C=O/C-:O mode of the neutral/anionic phylloquinone, respectively. To test this hypothesis, we have obtained A1(-)/A1 FTIR DS for menG mutant PS I particles. In menG mutant PS I, phylloquinone in the A1 binding site is replaced with an analogue in which the methyl group at position 2 of the quinone ring is replaced with a hydrogen atom (2-phytyl-1,4-naphthoquinone). In A1(-)/A1 FTIR DS obtained using menG mutant PS I particles, we find that the 1654/1495 cm(-1) bands are upshifted by approximately 6 cm(-1). To test if such upshifts are likely for C=O/C-:O modes of neutral/anionic phylloquinone, we have used density functional theory to calculate the "anion minus neutral" infrared difference spectra for both phylloquinone and its methyl-less analogue. We have also undertaken calculations in which the C4=O carbonyl group of phylloquinone and its methyl-less analogue are hydrogen bonded (to a water or leucine molecule). We find that, irrespective of the hydrogen bonding state of the C4=O group, the C=O/C-:O modes of neutral/reduced phylloquinone are indeed expected to be upshifted by at least 6 cm(-1) upon replacement of the methyl group at position 2 with hydrogen. The calculations also suggest that certain C=C/C-:C modes of neutral/reduced phylloquinone do not shift upon replacement of the methyl group. On the basis of these calculated results, we suggest which bands in the A1(-)/A1 FTIR DS may be associated with C=C/C-:C modes of neutral/reduced phylloquinone, respectively.     

Design and synthesis of chiral peptidic nucleic acids

Due to the increasing interest in the use ofoligonucleotide analogues as antisense and antigenedrugs, we designed a chiral analogue constituted of apeptidic frame bearing nucleobases in suitablepositions (C-PNA). We recently reported the synthesisof four nonnatural -amino acids with the DNAbases in the lateral chain. In this paper we presentan improved synthesis of the Fmoc monomers and theirpolymerisation to polypeptidic oligonucleotideanalogues using a modification of the standardprotocol for solid phase peptide synthesis.     

Template-directed interference footprinting for RNA based on inosine-specific cleavage

We report here the development of a Template-directed Interference (TDI) footprinting assay for RNA. The TDI nucleotide analogue inosine (I) lacks the exocyclic amine of G and is a suitable probe for the role of this group in RNA structure and function. Using an I-specific cleavage protocol we identified three functionally significant G residues in the Tetrahymena ribozyme. These residues are proximal to the active site of the folded intron and likely contribute to the positioning of substrates at the catalytic core.     

When You Find Yourself in a Hole, Stop Digging: Comments on Earth Systems Engineering

In an earlier "Forum" article in this journal, Brad Allenby outlined his views of a new approach to managing the unintended consequences of human activity, "earth systems engineering." He argues that we must develop the tools, institutions, and moral and ethical systems to allow us to "assume an active management role for most global systems." I believe this to be a significant departure from a core concept of industrial ecology: learning from ecosystems how the natural world operates to be able to more effectively design and manage coupled human-natural systems. Such lessons are more likely to lead away from tightly managed, centralized approaches, and favor approaches with as little intervention as feasible. More important, I believe that we are far less likely to learn how to implement earth systems engineering than simpler approaches, hence less likely to minimize environmental damage.     

Habitat structural complexity mediates food web dynamics in a freshwater macrophyte community

A considerable amount of research has investigated the influence of habitat structure on predator success, yet few studies have explored the implications for community structure and food-web dynamics. The relative importance of macrophyte structure and fish predation on the composition of the macroinvertebrate and periphyton communities in a lowland river was investigated using a multifactorial caging experiment. We hypothesised that: (1) fish predators are less effective in a more structurally complex macrophyte analogue; (2) strong direct and indirect effects of fish predators (e.g. trophic cascades) are less likely to occur in a structurally complex habitat; and (3) the strength of these patterns is influenced by the composition of the prevailing community assemblage. We measured the abundance and composition of the macroinvertebrate and periphyton communities associated with three different-shaped macrophyte analogues, under different fish predator treatments and at different times. Macrophyte analogue architecture had strong, consistent effects on both the macroinvertebrate and periphyton communities; both were most abundant and diverse on the most structurally complex plant analogue. In contrast, the fish predators affected only a subset of the macroinvertebrate community and there was a suggestion of minor indirect effects on periphyton community composition. Contrary to expectations, the fish predators had their strongest effects in the most structurally complex macrophyte analogue. We conclude that in this system, macrophyte shape strongly regulates the associated freshwater assemblage, resulting in a diverse community structure less likely to exhibit strong effects of fish predation.     

The camptothecin-resistant topoisomerase I mutant F361S is cross-resistant to antitumor rebeccamycin derivatives. A model for topoisomerase I inhibition by indolocarbazoles.

DNA topoisomerase I is a major cellular target for antitumor indolocarbazole derivatives (IND) such as the antibiotic rebeccamycin and the synthetic analogue NB-506 which is undergoing phase I clinical trials. We have investigated the mechanism of topoisomerase I inhibition by a rebeccamycin analogue, R-3, using the wild-type human topoisomerase I and a well-characterized recombinant enzyme, F361S. The catalytic activity of this mutant remains fully intact, but the enzyme is resistant to inhibition by camptothecin (CPT). Here we show that the mutated enzyme is cross-resistant to the rebeccamycin analogue. Despite their profound structural differences, CPT and R-3 interfere similarly with the activity of the wild-type and mutant topoisomerase I enzymes, and the drug-induced cleavable complexes are equally sensitive to the NaCl concentration. CPT and IND likely recognize identical structural elements of the topoisomerase I-DNA covalent complex; however, differences do exist in terms of sequence-specificity of topoisomerase I-mediated DNA cleavage. For the first time, a molecular model showing that CPT and IND share common steric and electronic features is proposed. The model helps to identify a specific pharmacophore for topoisomerase I inhibitors.     

Efficient algorithms for protein sequence design and the analysis of certain evolutionary fitness landscapes.

Protein sequence design is a natural inverse problem to protein structure prediction: given a target structure in three dimensions, we wish to design an amino acid sequence that is likely fold to it. A model of Sun, Brem, Chan, and Dill casts this problem as an optimization on a space of sequences of hydrophobic (H) and polar (P) monomers; the goal is to find a sequence that achieves a dense hydrophobic core with few solvent-exposed hydrophobic residues. Sun et al. developed a heuristic method to search the space of sequences, without a guarantee of optimality or near-optimality; Hart subsequently raised the computational tractability of constructing an optimal sequence in this model as an open question. Here we resolve this question by providing an efficient algorithm to construct optimal sequences; our algorithm has a polynomial running time, and performs very efficiently in practice. We illustrate the implementation of our method on structures drawn from the Protein Data Bank. We also consider extensions of the model to larger amino acid alphabets, as a way to overcome the limitations of the binary H/P alphabet. We show that for a natural class of arbitrarily large alphabets, it remains possible to design optimal sequences efficiently. Finally, we analyze some of the consequences of this sequence design model for the study of evolutionary fitness landscapes. A given target structure may have many sequences that are optimal in the model of Sun et al.; following a notion raised by the work of J. Maynard Smith, we can ask whether these optimal sequences are "connected" by successive point mutations. We provide a polynomial-time algorithm to decide this connectedness property, relative to a given target structure. We develop the algorithm by first solving an analogous problem expressed in terms of submodular functions, a fundamental object of study in combinatorial optimization.     

Synthesis and bioorthogonal coupling chemistry of a novel cyclopentenone-containing unnatural tyrosine analogue

Herein we report the synthesis of a novel amino acid with orthogonal functionality to the natural amino acid side chains. Tyrosine was O-alkylated with a cyclic 5-membered α,β-unsaturated ketone ring (5). We have established that this amino acid analogue can undergo cycloaddition reactions in aqueous media with in situ generated nitrones. Nitrone formation occurred by micellar catalysis can undergo aqueous 1,3-dipolar cycloaddition reactions with the unnatural Tyr. We also performed a linear free energy analysis of the one pot bioconjugation reaction in water using cyclopentenone as a model for the Tyr analogue and seven different aryl nitrones. We found that the Hammett ρ value was −0.94, suggesting that the reaction occurs in a concerted fashion with a slight positive charge buildup in the transition state. The Hammett ρ value also suggests that the bioconjugation reaction is tolerant of different substituents and thus may be useful for introducing novel functionality into peptides and proteins containing the Tyr analogue 5. The aqueous 1,3-dipolar cycloaddition reactions, that use nitrones to trap the O-alkylated Tyr 5, establish a novel strategy for rapid, water compatible bioconjugation reactions.     

Crystal structure of Trypanosoma cruzi glyceraldehyde-3-phosphate dehydrogenase complexed with an analogue of 1,3-bisphospho-d-glyceric acid.

We report here the first crystal structure of a stable isosteric analogue of 1,3-bisphospho-d-glyceric acid (1,3-BPGA) bound to the catalytic domain of Trypanosoma cruzi glycosomal glyceraldehyde-3-phosphate dehydrogenase (gGAPDH) in which the two phosphoryl moieties interact with Arg249. This complex possibly illustrates a step of the catalytic process by which Arg249 may induce compression of the product formed, allowing its expulsion from the active site. Structural modifications were introduced into this isosteric analogue and the respective inhibitory effects of the resulting diphosphorylated compounds on T. cruzi and Trypanosoma brucei gGAPDHs were investigated by enzymatic inhibition studies, fluorescence spectroscopy, site-directed mutagenesis, and molecular modelling. Despite the high homology between the two trypanomastid gGAPDHs (> 95%), we have identified specific interactions that could be used to design selective irreversible inhibitors against T. cruzi gGAPDH.     

Power in pairs: assessing the statistical value of paired samples in tests for differential expression

Background

When genomics researchers design a high-throughput study to test for differential expression, some biological systems and research questions provide opportunities to use paired samples from subjects, and researchers can plan for a certain proportion of subjects to have paired samples. We consider the effect of this paired samples proportion on the statistical power of the study, using characteristics of both count (RNA-Seq) and continuous (microarray) expression data from a colorectal cancer study.

Results

We demonstrate that a higher proportion of subjects with paired samples yields higher statistical power, for various total numbers of samples, and for various strengths of subject-level confounding factors. In the design scenarios considered, the statistical power in a fully-paired design is substantially (and in many cases several times) greater than in an unpaired design.

Conclusions

For the many biological systems and research questions where paired samples are feasible and relevant, substantial statistical power gains can be achieved at the study design stage when genomics researchers plan on using paired samples from the largest possible proportion of subjects. Any cost savings in a study design with unpaired samples are likely accompanied by underpowered and possibly biased results.

     

Biological effects of a de novo designed myxoma virus peptide analogue: evaluation of cytotoxicity on tumor cells

Background

The Resonant Recognition Model (RRM) is a physico-mathematical model that interprets protein sequence linear information using digital signal processing methods. In this study the RRM concept was employed for structure-function analysis of myxoma virus (MV) proteins and the design of a short bioactive therapeutic peptide with MV-like antitumor/cytotoxic activity.

Methodology/Principal Findings

The analogue RRM-MV was designed by RRM as a linear 18 aa 2.3 kDa peptide. The biological activity of this computationally designed peptide analogue against cancer and normal cell lines was investigated. The cellular cytotoxicity effects were confirmed by confocal immunofluorescence microscopy, by measuring the levels of cytoplasmic lactate dehydrogenase (LDH) and by Prestoblue cell viability assay for up to 72 hours in peptide treated and non-treated cell cultures. Our results revealed that RRM-MV induced a significant dose and time-dependent cytotoxic effect on murine and human cancer cell lines. Yet, when normal murine cell lines were similarly treated with RRM-MV, no cytotoxic effects were observed. Furthermore, the non-bioactive RRM designed peptide RRM-C produced negligible cytotoxic effects on these cancer and normal cell lines when used at similar concentrations. The presence/absence of phosphorylated Akt activity in B16F0 mouse melanoma cells was assessed to indicate the possible apoptosis signalling pathway that could be affected by the peptide treatment. So far, Akt activity did not seem to be significantly affected by RRM-MV as is the case for the original viral protein.

Conclusions/Significance

Our findings indicate the successful application of the RRM concept to design a bioactive peptide analogue (RRM-MV) with cytotoxic effects on tumor cells only. This 2.345 kDa peptide analogue to a 49 kDa viral protein may be suitable to be developed as a potential cancer therapeutic. These results also open a new direction to the rational design of therapeutic agents for future cancer treatment.     

Somatostatin inhibits PDGF-stimulated Ras activation in human neuroblastoma cells.

The main physiological role of somatostatin (SST) is the control of hormone secretion. Recently, SST has been shown to exert antiproliferative effects on some human tumors via both direct and indirect mechanisms. We have previously found that in the human neuroblastoma cell line SY5Y the SST analogue lanreotide (BIM 23014) inhibited serum-stimulated cell proliferation and MAP kinase activity. Here, we examine the effect of SST on PDGF-induced Ras activation. We found that SST suppressed PDGF-induced Ras activation in a pertussis toxin (PTx)-independent and peroxovanadate-dependent manner. Ras-specific GTPase activating protein (GAP) activities were not altered by SST treatment. On the contrary, PDGF-induced PDGF receptor phosphorylation was decreased by SST in a PTx-independent, peroxovanadate-dependent manner, likely accounting for the SST-mediated inhibition of PDGF-induced Ras activation.     

Diminished plasma cGMP during weightlessness

We performed an experiment within the project "RLF" (Russian long-term flight) on a cosmonaut onboard the space station MIR. For creating an analogue to orthostatic stress, we used lower body negative pressure (LBNP) as stimulus. Decrease in central and peripheral baroreceptor load by LBNP can be used as a cardiovascular countermeasure in cosmonauts or for inducing endocrine responses. Altered steady-state plasma concentration values of volume sensitive hormones have been observed inflight as well as postflight. Within this project we measured plasma ANP and cGMP as second messenger. Changes in plasma cGMP concentration are generally considered to be a good indicator of those in ANP activity. However, in our experiments depression of cGMP during space flight was more impressive than ANP decline. We are not aware of previous measurements of plasma cGMP under these conditions, and believe to be the first to report complete suppression of plasma cGMP during long-term stay in space.     

Humanised models of cancer in molecular medicine: the experimental control of disanalogy

This paper explores the epistemology of extrapolation from model organisms to humans in molecular medicine. We take into account two common views on the issue, the homology view and the disanalogy view. In response to both interpretations, we argue that the foundational basis of extrapolations cannot simply be provided by homology and that relevant disanalogies can, thanks to the techniques of molecular biology, be experimentally controlled and exploited to allow useful and reliable extrapolations. The case of "humanised mice" in the context of cancer stem cell research provides evidence of how animal models can be construed to approximate bona fide causal analogue models of human diseases. To supplement this view we show how the epistemology of model organisms needs to take into account the engineering side of molecular medicine. Model organisms are often manipulated to create analogies or remove disanalogies with the target system. We maintain that highlighting this feature is fundamental to explain what warrants extrapolation in the search for the molecular causes of disease.