Search

Cancers, such as squamous cell carcinoma, frequently invade as multicellular units. However, these invading units can be organised in a variety of ways, ranging from thin discontinuous strands to thick ‘pushing’ collectives. Here we employ an integrated experimental and computational approach to identify the factors that determine the mode of collective cancer cell invasion. We find that matrix proteolysis is linked to the formation of wide strands but has little effect on the maximum extent of invasion. Cell-cell junctions also favour wide strands, but our analysis also reveals a requirement for cell-cell junctions for efficient invasion in response to uniform directional cues. Unexpectedly, the ability to generate wide invasive strands is coupled to the ability to grow effectively when surrounded by extracellular matrix in three-dimensional assays. Combinatorial perturbation of both matrix proteolysis and cell-cell adhesion demonstrates that the most aggressive cancer behaviour, both in terms of invasion and growth, is achieved at high levels of cell-cell adhesion and high levels of proteolysis. Contrary to expectation, cells with canonical mesenchymal traits – no cell-cell junctions and high proteolysis – exhibit reduced growth and lymph node metastasis. Thus, we conclude that the ability of squamous cell carcinoma cells to invade effectively is also linked to their ability to generate space for proliferation in confined contexts. These data provide an explanation for the apparent advantage of retaining cell-cell junctions in squamous cell carcinomas.

PCR-based methods to amplify the 3' untranslated region (3'-UTR) of the heat shock protein 70 (type I) gene (HSP70-I) have previously been used for typing of Leishmania but not with Leishmania (Mundinia) martiniquensis and L. (Mundinia) orientalis, newly identified human pathogens. Here, the 3'-UTRs of HSP70-I of L. martiniquensis, L. orientalis, and 10 other species were sequenced and analyzed. PCR-Restriction Fragment Length Polymorphism (RFLP) analysis targeting the 3'-UTR of HSP70-I was developed. Also, the detection limit of HSP70-I-3'-UTR PCR methods was compared with two other commonly used targets: the 18S small subunit ribosomal RNA (SSU-rRNA) gene and the internal transcribed spacer 1 region of the rRNA (ITS1-rRNA) gene. Results showed that HSP70-I-3'-UTR PCR methods could be used to identify and differentiate between L. martiniquensis (480-2 bp) and L. orientalis (674 bp) and distinguished them from parasites of the subgenus Viannia and of the subgenus Leishmania. PCR-RFLP patterns of the 3'-UTR of HSP70-I fragments digested with BsuRI restriction enzyme successfully differentiated L. martiniquensis, L. orientalis, L. braziliensis, L. guyanensis = L. panamensis, L. mexicana = L. aethiopica = L. tropica, L. amazonensis, L. major, and L. donovani = L. infantum. For the detection limit, the HSP70-I-3'-UTR PCR method could detect the DNA of L. martiniquensis and L. orientalis at the same concentration, 1 pg/μL, at a similar level to the SSU-rRNA PCR. The PCR that amplified ITS1-rRNA was more sensitive (0.01 pg/μL) than that of the HSP70-I-3'-UTR PCR. However, the sizes of both SSU-rRNA and ITS1-rRNA PCR amplicons could not differentiate between L. martiniquensis and L. orientalis. This is the first report of using HSP70-I-3'-UTR PCR based methods to identify the parasites causing leishmaniasis in Thailand. Also, the BsuRI-PCR-RFLP method can be used for differentiating some species within other subgenera.

The higher-order patterning of extra-cellular matrix in normal and pathological tissues has profound consequences on tissue function. Whilst studies have documented both how fibroblasts create and maintain individual matrix fibers and how cell migration is altered by the fibers they interact with, a model unifying these two aspects of tissue organization is lacking. Here we use computational modelling to understand the effect of this interconnectivity between fibroblasts and matrix at the mesoscale level. We created a unique adaptation to the Vicsek flocking model to include feedback from a second layer representing the matrix, and use experimentation to parameterize our model and validate model-driven hypotheses. Our two-layer model demonstrates that feedback between fibroblasts and matrix increases matrix diversity creating higher-order patterns. The model can quantitatively recapitulate matrix patterns of tissues in vivo. Cells follow matrix fibers irrespective of when the matrix fibers were deposited, resulting in feedback with the matrix acting as temporal 'memory' to collective behaviour, which creates diversity in topology. We also establish conditions under which matrix can be remodelled from one pattern to another. Our model elucidates how simple rules defining fibroblast-matrix interactions are sufficient to generate complex tissue patterns.

Amphotericin B is an increasingly important tool in efforts to reduce the global disease burden posed by Leishmania parasites. With few other chemotherapeutic options available for the treatment of leishmaniasis, the potential for emergent resistance to this drug is a considerable threat. Here we characterised four novel amphotericin B-resistant Leishmania mexicana lines. All lines exhibited altered sterol biosynthesis, and hypersensitivity to pentamidine. Whole genome sequencing demonstrated resistance-associated mutation of the sterol biosynthesis gene sterol C5-desaturase in one line. However, in three out of four lines, RNA-seq revealed loss of expression of sterol C24-methyltransferase (SMT) responsible for drug resistance and altered sterol biosynthesis. Additional loss of the miltefosine transporter was associated with one of those lines. SMT is encoded by two tandem gene copies, which we found to have very different expression levels. In all cases, reduced overall expression was associated with loss of the 3' untranslated region of the dominant gene copy, resulting from structural variations at this locus. Local regions of sequence homology, between the gene copies themselves, and also due to the presence of SIDER1 retrotransposon elements that promote multi-gene amplification, correlate to these structural variations. Moreover, in at least one case loss of SMT expression was not associated with loss of virulence in primary macrophages or in vivo. Whilst such repeat sequence-mediated instability is known in Leishmania genomes, its presence associated with resistance to a major antileishmanial drug, with no evidence of associated fitness costs, is a significant concern.

Leishmania parasites are transmitted to vertebrate hosts by female phlebotomine sand flies as they bloodfeed by lacerating the upper capillaries of the dermis with their barbed mouthparts. In the sand fly midgut secreted proteophosphoglycans from Leishmania form a biological plug known as the promastigote secretory gel (PSG), which blocks the gut and facilitates the regurgitation of infective parasites. The interaction between the wound created by the sand fly bite and PSG is not known. Here we nanoinjected a sand fly egested dose of PSG into BALB/c mouse skin that lead to the differential expression of 7,907 transcripts. These transcripts were transiently up-regulated during the first 6 hours post-wound and enriched for pathways involved in inflammation, cell proliferation, fibrosis, epithelial cell differentiation and wound remodelling. We found that PSG significantly accelerated wound healing in vitro and in mice; which was associated with an early up-regulation of transcripts involved in inflammation (IL-1β, IL-6, IL-10, TNFα) and inflammatory cell recruitment (CCL2, CCL3, CCL4, CXCL2), followed 6 days later by enhanced expression of transcripts associated with epithelial cell proliferation, fibroplasia and fibrosis (FGFR2, EGF, EGFR, IGF1). Dermal expression of IGF1 was enhanced following an infected sand fly bite and was acutely responsive to the deposition of PSG but not the inoculation of parasites or sand fly saliva. Antibody blockade of IGF1 ablated the gel's ability to promote wound closure in mouse ears and significantly reduced the virulence of Leishmania mexicana infection delivered by an individual sand fly bite. Dermal macrophages recruited to air-pouches on the backs of mice revealed that IGF1 was pivotal to the PSG's ability to promote macrophage alternative activation and Leishmania infection. Our data demonstrate that through the regurgitation of PSG Leishmania exploit the wound healing response of the host to the vector bite by promoting the action of IGF1 to drive the alternative activation of macrophages.

Accurate protein structure prediction from amino acid sequence is still an unsolved problem. The most reliable methods centre on template based modelling. However, the accuracy of these models entirely depends on the availability of experimentally resolved homologous template structures. In order to generate more accurate models, extensive physics based molecular dynamics (MD) refinement simulations are performed to sample many different conformations to find improved conformational states. In this study, we propose a deep recurrent network model, called DeepTrajectory, that is able to identify these improved conformational states, with high precision, from a variety of different MD based sampling protocols. The proposed model learns the temporal patterns of features computed from MD trajectory data in order to classify whether each recorded simulation snapshot is an improved quality conformational state, decreased quality conformational state or whether there is no perceivable change in state with respect to the starting conformation. The model was trained and tested on 904 trajectories from 42 different protein systems with a cumulative number of more than 1.7 million snapshots. We show that our model outperforms other state of the art machine-learning algorithms that do not consider temporal dependencies. To our knowledge, DeepTrajectory is the first implementation of a time-dependent deep-learning protocol that is re-trainable and able to adapt to any new MD based sampling procedure, thereby demonstrating how a neural network can be used to learn the latter part of the protein folding funnel.

Abstract Phlebotomine sand flies are the subject of much research because of the role of their females as the only proven natural vectors of Leishmania species, the parasitic protozoans that are the causative agents of the neglected tropical disease leishmaniasis. Activity in this field was highlighted by the eighth International Symposium on Phlebotomine Sand flies (ISOPS) held in September 2014, which prompted this review focusing on vector control. Topics reviewed include: Taxonomy and phylogenetics, Vector competence, Genetics, genomics and transcriptomics, Eco-epidemiology, and Vector control. Research on sand flies as leishmaniasis vectors has revealed a diverse array of zoonotic and anthroponotic transmission cycles, mostly in subtropical and tropical regions of Africa, Asia and Latin America, but also in Mediterranean Europe. The challenge is to progress beyond descriptive eco-epidemiology, in order to separate vectors of biomedical importance from the sand fly species that are competent vectors but lack the vectorial capacity to cause much human disease. Transmission modelling is required to identify the vectors that are a public health priority, the ones that must be controlled as part of the integrated control of leishmaniasis. Effective modelling of transmission will require the use of entomological indices more precise than those usually reported in the leishmaniasis literature.

Abstract Background Phlebotomine sand flies transmit the haemoflagellate Leishmania, the causative agent of human leishmaniasis. The Leishmania promastigotes are confined to the gut lumen and are exposed to the gut microbiota within female sand flies. Here we study the colonisation resistance of yeast and bacteria in preventing the establishment of a Leishmania population in sand flies and the ability of Leishmania to provide colonisation resistance towards the insect bacterial pathogen Serratia marcescens that is also pathogenic towards Leishmania. Methods We isolated microorganisms from wild-caught and laboratory-reared female Lutzomyia longipalpis, identified as Pseudozyma sp. Asaia sp. and Ochrobactrum intermedium. We fed the females with a sugar meal containing the microorganisms and then subsequently fed them with a bloodmeal containing Leishmania mexicana and recorded the development of the Leishmania population. Further experiments examined the effect of first colonising the sand fly gut with L. mexicana followed by feeding with, Serratia marcescens, an insect bacterial pathogen. The mortality of the flies due to S. marcescens was recorded in the presence and absence of Leishmania. Results There was a reduction in the number of flies harbouring a Leishmania population that had been pre-fed with Pseudozyma sp. and Asaia sp. or O. intermedium. Experiments in which L. mexicana colonised the sand fly gut prior to being fed an insect bacterial pathogen, Serratia marcescens, showed that the survival of flies with a Leishmania infection was significantly higher compared to flies without Leishmania infection. Conclusions The yeast and bacterial colonisation experiments show that the presence of sand fly gut microorganisms reduce the potential for Leishmania to establish within the sand fly vector. Sand flies infected with Leishmania were able to survive an attack by the bacterial pathogen that would have killed the insect and we concluded that Leishmania may benefit its insect host whilst increasing the potential to establish itself in the sand fly vector. We suggest that the increased ability of the sand fly to withstand a bacterial entomopathogen, due to the presence of the Leishmania, may provide an evolutionary pressure for the maintenance of the Leishmania-vector association.

Understanding changes in mosquito salivary proteins during the time that sporozoite maturation occurs and after blood feeding may give information regarding the roles of salivary proteins during the malarial transmission. Anopheles dissidens (formerly Anopheles barbirostris species A1) is a potential vector of Plasmodium vivax in Thailand. In this study, analyses of the proteomic profiles of female An. dissidens salivary glands during adult development and after blood feeding were carried out using two-dimensional gel electrophoresis coupled with nano-liquid chromatography-mass spectrometry. Results showed at least 17 major salivary gland proteins present from day one to day 21 post emergence at 8 different time points sampled. Although there was variation observed, the patterns of protein expression could be placed into one of four groups. Fifteen protein spots showed significant depletion after blood feeding with the percentages of the amount of depletion ranging from 8.5% to 68.11%. The overall results identified various proteins, including a putative mucin-like protein, an anti-platelet protein, a long form D7 salivary protein, a putative gVAG protein precursor, a D7-related 3.2 protein, gSG7 salivary proteins, and a gSG6 protein. These results allow better understanding of the changes of the salivary proteins during the adult mosquito development. They also provide candidate proteins to investigate any possible link or not between sporozoite maturation, or survival of skin stage sporozoites, and salivary proteins.

Mitotic chromosomes were one of the first cell biological structures to be described, yet their molecular architecture remains poorly understood. We have devised a simple biophysical model of a 300 kb-long nucleosome chain, the size of a budding yeast chromosome, constrained by interactions between binding sites of the chromosomal condensin complex, a key component of interphase and mitotic chromosomes. Comparisons of computational and experimental (4C) interaction maps, and other biophysical features, allow us to predict a mode of condensin action. Stochastic condensin-mediated pairwise interactions along the nucleosome chain generate native-like chromosome features and recapitulate chromosome compaction and individualization during mitotic condensation. Higher order interactions between condensin binding sites explain the data less well. Our results suggest that basic assumptions about chromatin behavior go a long way to explain chromosome architecture and are able to generate a molecular model of what the inside of a chromosome is likely to look like.

Since 1996, there have been several case reports of autochthonous visceral leishmaniasis in Thailand. Here we report a case in a 52-year-old Thai male from northern Thailand, who presented with subacute fever, huge splenomegaly and pancytopenia. Bone marrow aspiration revealed numerous amastigotes within macrophages. Isolation of Leishmania LSCM1 into culture and DNA sequence analysis (ribosomal RNA ITS-1 and large subunit of RNA polymerase II) revealed the parasites to be members of the Leishmania enriettii complex, and apparently identical to L. martiniquensis previously reported from the Caribbean island of Martinique. This is the first report of visceral leishmaniasis caused by L. martiniquensis from the region. Moreover, the majority of parasites previously identified as "L. siamensis" also appear to be L. martiniquensis.

Malaria sporozoites must invade the salivary glands of mosquitoes for maturation before transmission to vertebrate hosts. The duration of the sporogonic cycle within the mosquitoes ranges from 10 to 21 days depending on the parasite species and temperature. During blood feeding salivary gland proteins are injected into the vertebrate host, along with malaria sporozoites in the case of an infected mosquito. To identify salivary gland proteins depleted after blood feeding of female Anopheles campestris-like, a potential malaria vector of Plasmodium vivax in Thailand, two-dimensional gel electrophoresis and nano-liquid chromatography-mass spectrometry techniques were used. Results showed that 19 major proteins were significantly depleted in three to four day-old mosquitoes fed on a first blood meal. For the mosquitoes fed the second blood meal on day 14 after the first blood meal, 14 major proteins were significantly decreased in amount. The significantly depleted proteins in both groups included apyrase, 5'-nucleotidase/apyrase, D7, D7-related 1, short form D7r1, gSG6, anti-platelet protein, serine/threonine-protein kinase rio3, putative sil1, cyclophilin A, hypothetical protein Phum_PHUM512530, AGAP007618-PA, and two non-significant hit proteins. To our knowledge, this study presents for the first time the salivary gland proteins that are involved in the second blood feeding on the day corresponding to the transmission period of the sporozoites to new mammalian hosts. This information serves as a basis for future work concerning the possible role of these proteins in the parasite transmission and the physiological processes that occur during the blood feeding.

Predicting the effects of mutations on the kinetic rate constants of protein-protein interactions is central to both the modeling of complex diseases and the design of effective peptide drug inhibitors. However, while most studies have concentrated on the determination of association rate constants, dissociation rates have received less attention. In this work we take a novel approach by relating the changes in dissociation rates upon mutation to the energetics and architecture of hotspots and hotregions, by performing alanine scans pre- and post-mutation. From these scans, we design a set of descriptors that capture the change in hotspot energy and distribution. The method is benchmarked on 713 kinetically characterized mutations from the SKEMPI database. Our investigations show that, with the use of hotspot descriptors, energies from single-point alanine mutations may be used for the estimation of off-rate mutations to any residue type and also multi-point mutations. A number of machine learning models are built from a combination of molecular and hotspot descriptors, with the best models achieving a Pearson's Correlation Coefficient of 0.79 with experimental off-rates and a Matthew's Correlation Coefficient of 0.6 in the detection of rare stabilizing mutations. Using specialized feature selection models we identify descriptors that are highly specific and, conversely, broadly important to predicting the effects of different classes of mutations, interface regions and complexes. Our results also indicate that the distribution of the critical stability regions across protein-protein interfaces is a function of complex size more strongly than interface area. In addition, mutations at the rim are critical for the stability of small complexes, but consistently harder to characterize. The relationship between hotregion size and the dissociation rate is also investigated and, using hotspot descriptors which model cooperative effects within hotregions, we show how the contribution of hotregions of different sizes, changes under different cooperative effects.

Gauging the systemic effects of non-synonymous single nucleotide polymorphisms (nsSNPs) is an important topic in the pursuit of personalized medicine. However, it is a non-trivial task to understand how a change at the protein structure level eventually affects a cell's behavior. This is because complex information at both the protein and pathway level has to be integrated. Given that the idea of integrating both protein and pathway dynamics to estimate the systemic impact of missense mutations in proteins remains predominantly unexplored, we investigate the practicality of such an approach by formulating mathematical models and comparing them with experimental data to study missense mutations. We present two case studies: (1) interpreting systemic perturbation for mutations within the cell cycle control mechanisms (G2 to mitosis transition) for yeast; (2) phenotypic classification of neuron-related human diseases associated with mutations within the mitogen-activated protein kinase (MAPK) pathway. We show that the application of simplified mathematical models is feasible for understanding the effects of small sequence changes on cellular behavior. Furthermore, we show that the systemic impact of missense mutations can be effectively quantified as a combination of protein stability change and pathway perturbation.

Phlebotomine sand flies are vectors of Leishmania that are acquired by the female sand fly during blood feeding on an infected mammal. Leishmania parasites develop exclusively in the gut lumen during their residence in the insect before transmission to a suitable host during the next blood feed. Female phlebotomine sand flies are blood feeding insects but their life style of visiting plants as well as animals, and the propensity for larvae to feed on detritus including animal faeces means that the insect host and parasite are exposed to a range of microorganisms. Thus, the sand fly microbiota may interact with the developing Leishmania population in the gut. The aim of the study was to investigate and identify the bacterial diversity associated with wild adult female Lutzomyia sand flies from different geographical locations in the New World. The bacterial phylotypes recovered from 16S rRNA gene clone libraries obtained from wild caught adult female Lutzomyia sand flies were estimated from direct band sequencing after denaturing gradient gel electrophoresis of bacterial 16 rRNA gene fragments. These results confirm that the Lutzomyia sand flies contain a limited array of bacterial phylotypes across several divisions. Several potential plant-related bacterial sequences were detected including Erwinia sp. and putative Ralstonia sp. from two sand fly species sampled from 3 geographically separated regions in Brazil. Identification of putative human pathogens also demonstrated the potential for sand flies to act as vectors of bacterial pathogens of medical importance in addition to their role in Leishmania transmission.

The prediction of protein-protein kinetic rate constants provides a fundamental test of our understanding of molecular recognition, and will play an important role in the modeling of complex biological systems. In this paper, a feature selection and regression algorithm is applied to mine a large set of molecular descriptors and construct simple models for association and dissociation rate constants using empirical data. Using separate test data for validation, the predicted rate constants can be combined to calculate binding affinity with accuracy matching that of state of the art empirical free energy functions. The models show that the rate of association is linearly related to the proportion of unbound proteins in the bound conformational ensemble relative to the unbound conformational ensemble, indicating that the binding partners must adopt a geometry near to that of the bound prior to binding. Mirroring the conformational selection and population shift mechanism of protein binding, the models provide a strong separate line of evidence for the preponderance of this mechanism in protein-protein binding, complementing structural and theoretical studies.

The binding of Leishmania promastigotes to the midgut epithelium is regarded as an essential part of the life-cycle in the sand fly vector, enabling the parasites to persist beyond the initial blood meal phase and establish the infection. However, the precise nature of the promastigote stage(s) that mediate binding is not fully understood.To address this issue we have developed an in vitro gut binding assay in which two promastigote populations are labelled with different fluorescent dyes and compete for binding to dissected sand fly midguts. Binding of procyclic, nectomonad, leptomonad and metacyclic promastigotes of Leishmania infantum and L. mexicana to the midguts of blood-fed, female Lutzomyia longipalpis was investigated. The results show that procyclic and metacyclic promastigotes do not bind to the midgut epithelium in significant numbers, whereas nectomonad and leptomonad promastigotes both bind strongly and in similar numbers. The assay was then used to compare the binding of a range of different parasite species (L. infantum, L. mexicana, L. braziliensis, L. major, L. tropica) to guts dissected from various sand flies (Lu. longipalpis, Phlebotomus papatasi, P. sergenti). The results of these comparisons were in many cases in line with expectations, the natural parasite binding most effectively to its natural vector, and no examples were found where a parasite was unable to bind to its natural vector. However, there were interesting exceptions: L. major and L. tropica being able to bind to Lu. longipalpis better than L. infantum; L. braziliensis was able to bind to P. papatasi as well as L. major; and significant binding of L. major to P. sergenti and L. tropica to P. papatasi was observed.The results demonstrate that Leishmania gut binding is strictly stage-dependent, is a property of those forms found in the middle phase of development (nectomonad and leptomonad forms), but is absent in the early blood meal and final stages (procyclic and metacyclic forms). Further they show that although gut binding may be necessary for parasite establishment, in several vector-parasite pairs the specificity of such in vitro binding alone is insufficient to explain overall vector specificity. Other significant barriers to development must exist in certain refractory Leishmania parasite-sand fly vector combinations. A re-appraisal of the specificity of the Leishmania-sand fly relationship is required.

Vascular abnormalities contribute to many diseases such as cancer and diabetic retinopathy. In angiogenesis new blood vessels, headed by a migrating tip cell, sprout from pre-existing vessels in response to signals, e.g., vascular endothelial growth factor (VEGF). Tip cells meet and fuse (anastomosis) to form blood-flow supporting loops. Tip cell selection is achieved by Dll4-Notch mediated lateral inhibition resulting, under normal conditions, in an interleaved arrangement of tip and non-migrating stalk cells. Previously, we showed that the increased VEGF levels found in many diseases can cause the delayed negative feedback of lateral inhibition to produce abnormal oscillations of tip/stalk cell fates. Here we describe the development and implementation of a novel physics-based hierarchical agent model, tightly coupled to in vivo data, to explore the system dynamics as perpetual lateral inhibition combines with tip cell migration and fusion. We explore the tipping point between normal and abnormal sprouting as VEGF increases. A novel filopodia-adhesion driven migration mechanism is presented and validated against in vivo data. Due to the unique feature of ongoing lateral inhibition, 'stabilised' tip/stalk cell patterns show sensitivity to the formation of new cell-cell junctions during fusion: we predict cell fates can reverse. The fusing tip cells become inhibited and neighbouring stalk cells flip fate, recursively providing new tip cells. Junction size emerges as a key factor in establishing a stable tip/stalk pattern. Cell-cell junctions elongate as tip cells migrate, which is shown to provide positive feedback to lateral inhibition, causing it to be more susceptible to pathological oscillations. Importantly, down-regulation of the migratory pathway alone is shown to be sufficient to rescue the sprouting system from oscillation and restore stability. Thus we suggest the use of migration inhibitors as therapeutic agents for vascular normalisation in cancer.

In nature the prevalence of Leishmania infection in whole sand fly populations can be very low (