Multiscale feature analysis of salivary gland branching morphogenesis

Pattern formation in developing tissues involves dynamic spatio-temporal changes in cellular organization and subsequent evolution of functional adult structures. Branching morphogenesis is a developmental mechanism by which patterns are generated in many developing organs, which is controlled by underlying molecular pathways. Understanding the relationship between molecular signaling, cellular behavior and resulting morphological change requires quantification and categorization of the cellular behavior. In this study, tissue-level and cellular changes in developing salivary gland in response to disruption of ROCK-mediated signaling by are modeled by building cell-graphs to compute mathematical features capturing structural properties at multiple scales. These features were used to generate multiscale cell-graph signatures of untreated and ROCK signaling disrupted salivary gland organ explants. From confocal images of mouse submandibular salivary gland organ explants in which epithelial and mesenchymal nuclei were marked, a multiscale feature set capturing global structural properties, local structural properties, spectral, and morphological properties of the tissues was derived. Six feature selection algorithms and multiway modeling of the data was performed to identify distinct subsets of cell graph features that can uniquely classify and differentiate between different cell populations. Multiscale cell-graph analysis was most effective in classification of the tissue state. Cellular and tissue organization, as defined by a multiscale subset of cell-graph features, are both quantitatively distinct in epithelial and mesenchymal cell types both in the presence and absence of ROCK inhibitors. Whereas tensor analysis demonstrate that epithelial tissue was affected the most by inhibition of ROCK signaling, significant multiscale changes in mesenchymal tissue organization were identified with this analysis that were not identified in previous biological studies. We here show how to define and calculate a multiscale feature set as an effective computational approach to identify and quantify changes at multiple biological scales and to distinguish between different states in developing tissues.     

Genetic Modification and Recombination of Salivary Gland Organ Cultures

Branching morphogenesis occurs during the development of many organs, and the embryonic mouse submandibular gland (SMG) is a classical model for the study of branching morphogenesis. In the developing SMG, this process involves iterative steps of epithelial bud and duct formation, to ultimately give rise to a complex branched network of acini and ducts, which serve to produce and modify/transport the saliva, respectively, into the oral cavity^1-3. The epithelial-associated basement membrane and aspects of the mesenchymal compartment, including the mesenchyme cells, growth factors and the extracellular matrix, produced by these cells, are critical to the branching mechanism, although how the cellular and molecular events are coordinated remains poorly understood ⁴. The study of the molecular mechanisms driving epithelial morphogenesis advances our understanding of developmental mechanisms and provides insight into possible regenerative medicine approaches. Such studies have been hampered due to the lack of effective methods for genetic manipulation of the salivary epithelium. Currently, adenoviral transduction represents the most effective method for targeting epithelial cells in adult glands in vivo⁵. However, in embryonic explants, dense mesenchyme and the basement membrane surrounding the epithelial cells impedes viral access to the epithelial cells. If the mesenchyme is removed, the epithelium can be transfected using adenoviruses, and epithelial rudiments can resume branching morphogenesis in the presence of Matrigel or laminin-111^6,7. Mesenchyme-free epithelial rudiment growth also requires additional supplementation with soluble growth factors and does not fully recapitulate branching morphogenesis as it occurs in intact glands⁸. Here we describe a technique which facilitates adenoviral transduction of epithelial cells and culture of the transfected epithelium with associated mesenchyme. Following microdissection of the embryonic SMGs, removal of the mesenchyme, and viral infection of the epithelium with a GFP-containing adenovirus, we show that the epithelium spontaneously recombines with uninfected mesenchyme, recapitulating intact SMG glandular structure and branching morphogenesis. The genetically modified epithelial cell population can be easily monitored using standard fluorescence microscopy methods, if fluorescently-tagged adenoviral constructs are used. The tissue recombination method described here is currently the most effective and accessible method for transfection of epithelial cells with a wild-type or mutant vector within a complex 3D tissue construct that does not require generation of transgenic animals.     

LIM kinase regulation of cytoskeletal dynamics is required for salivary gland branching morphogenesis

Coordinated actin microfilament and microtubule dynamics is required for salivary gland development, although the mechanisms by which they contribute to branching morphogenesis are not defined. Because LIM kinase (LIMK) regulates both actin and microtubule organization, we investigated the role of LIMK signaling in mouse embryonic submandibular salivary glands using ex vivo organ cultures. Both LIMK 1 and 2 were necessary for branching morphogenesis and functioned to promote epithelial early- and late-stage cleft progression through regulation of both microfilaments and microtubules. LIMK-dependent regulation of these cytoskeletal systems was required to control focal adhesion protein–dependent fibronectin assembly and integrin β1 activation, involving the LIMK effectors cofilin and TPPP/p25, for assembly of the actin- and tubulin-based cytoskeletal systems, respectively. We demonstrate that LIMK regulates the early stages of cleft formation—cleft initiation, stabilization, and progression—via establishment of actin stability. Further, we reveal a novel role for the microtubule assembly factor p25 in regulating stabilization and elongation of late-stage progressing clefts. This study demonstrates the existence of multiple actin- and microtubule-dependent stabilization steps that are controlled by LIMK and are required in cleft progression during branching morphogenesis.     

An ensemble micro neural network approach for elucidating interactions between zinc finger proteins and their target DNA

Background

The ability to engineer zinc finger proteins binding to a DNA sequence of choice is essential for targeted genome editing to be possible. Experimental techniques and molecular docking have been successful in predicting protein-DNA interactions, however, they are highly time and resource intensive. Here, we present a novel algorithm designed for high throughput prediction of optimal zinc finger protein for 9 bp DNA sequences of choice. In accordance with the principles of information theory, a subset identified by using K-means clustering was used as a representative for the space of all possible 9 bp DNA sequences. The modeling and simulation results assuming synergistic mode of binding obtained from this subset were used to train an ensemble micro neural network. Synergistic mode of binding is the closest to the DNA-protein binding seen in nature, and gives much higher quality predictions, while the time and resources increase exponentially in the trade off. Our algorithm is inspired from an ensemble machine learning approach, and incorporates the predictions made by 100 parallel neural networks, each with a different hidden layer architecture designed to pick up different features from the training dataset to predict optimal zinc finger proteins for any 9 bp target DNA.

Results

The model gave an accuracy of an average 83% sequence identity for the testing dataset. The BLAST e-value are well within the statistical confidence interval of E-05 for 100% of the testing samples. The geometric mean and median value for the BLAST e-values were found to be 1.70E-12 and 7.00E-12 respectively. For final validation of approach, we compared our predictions against optimal ZFPs reported in literature for a set of experimentally studied DNA sequences. The accuracy, as measured by the average string identity between our predictions and the optimal zinc finger protein reported in literature for a 9 bp DNA target was found to be as high as 81% for DNA targets with a consensus sequence GCNGNNGCN reported in literature. Moreover, the average string identity of our predictions for a catalogue of over 100 9 bp DNA for which the optimal zinc finger protein has been reported in literature was found to be 71%.

Conclusions

Validation with experimental data shows that our tool is capable of domain adaptation and thus scales well to datasets other than the training set with high accuracy. As synergistic binding comes the closest to the ideal mode of binding, our algorithm predicts biologically relevant results in sync with the experimental data present in the literature. While there have been disjointed attempts to approach this problem synergistically reported in literature, there is no work covering the whole sample space. Our algorithm allows designing zinc finger proteins for DNA targets of the user’s choice, opening up new frontiers in the field of targeted genome editing. This algorithm is also available as an easy to use web server, ZifNN, at http://web.iitd.ac.in/~sundar/ZifNN/.

     

Cell-Based Multi-Parametric Model of Cleft Progression during Submandibular Salivary Gland Branching Morphogenesis

Cleft formation during submandibular salivary gland branching morphogenesis is the critical step initiating the growth and development of the complex adult organ. Previous experimental studies indicated requirements for several epithelial cellular processes, such as proliferation, migration, cell-cell adhesion, cell-extracellular matrix (matrix) adhesion, and cellular contraction in cleft formation; however, the relative contribution of each of these processes is not fully understood since it is not possible to experimentally manipulate each factor independently. We present here a comprehensive analysis of several cellular parameters regulating cleft progression during branching morphogenesis in the epithelial tissue of an early embryonic salivary gland at a local scale using an on lattice Monte-Carlo simulation model, the Glazier-Graner-Hogeweg model. We utilized measurements from time-lapse images of mouse submandibular gland organ explants to construct a temporally and spatially relevant cell-based 2D model. Our model simulates the effect of cellular proliferation, actomyosin contractility, cell-cell and cell-matrix adhesions on cleft progression, and it was used to test specific hypotheses regarding the function of these parameters in branching morphogenesis. We use innovative features capturing several aspects of cleft morphology and quantitatively analyze clefts formed during functional modification of the cellular parameters. Our simulations predict that a low epithelial mitosis rate and moderate level of actomyosin contractility in the cleft cells promote cleft progression. Raising or lowering levels of contractility and mitosis rate resulted in non-progressive clefts. We also show that lowered cell-cell adhesion in the cleft region and increased cleft cell-matrix adhesions are required for cleft progression. Using a classifier-based analysis, the relative importance of these four contributing cellular factors for effective cleft progression was determined as follows: cleft cell contractility, cleft region cell-cell adhesion strength, epithelial cell mitosis rate, and cell-matrix adhesion strength.     

Prediction of DNA-binding specificity in zinc finger proteins

Zinc finger proteins interact via their individual fingers to three base pair subsites on the target DNA. The four key residue positions -1, 2, 3 and 6 on the alpha-helix of the zinc fingers have hydrogen bond interactions with the DNA. Mutating these key residues enables generation of a plethora of combinatorial possibilities that can bind to any DNA stretch of interest. Exploiting the binding specificity and affinity of the interaction between the zinc fingers and the respective DNA can help to generate engineered zinc fingers for therapeutic purposes involving genome targeting. Exploring the structure-function relationships of the existing zinc finger-DNA complexes can aid in predicting the probable zinc fingers that could bind to any target DNA. Computational tools ease the prediction of such engineered zinc fingers by effectively utilizing information from the available experimental data. A study of literature reveals many approaches for predicting DNA-binding specificity in zinc finger proteins. However, an alternative approach that looks into the physico-chemical properties of these complexes would do away with the difficulties of designing unbiased zinc fingers with the desired affinity and specificity. We present a physico-chemical approach that exploits the relative strengths of hydrogen bonding between the target DNA and all combinatorially possible zinc fingers to select the most optimum zinc finger protein candidate.     

Exploiting the recognition code for elucidating the mechanism of zinc finger protein-DNA interactions

Background

Engineering zinc finger protein motifs for specific binding to double-stranded DNA is critical for targeted genome editing. Most existing tools for predicting DNA-binding specificity in zinc fingers are trained on data obtained from naturally occurring proteins, thereby skewing the predictions. Moreover, these mostly neglect the cooperativity exhibited by zinc fingers.

Methods

Here, we present an ab-initio method that is based on mutation of the key α-helical residues of individual fingers of the parent template for Zif-268 and its consensus sequence (PDB ID: 1AAY). In an attempt to elucidate the mechanism of zinc finger protein-DNA interactions, we evaluated and compared three approaches, differing in the amino acid mutations introduced in the Zif-268 parent template, and the mode of binding they try to mimic, i.e., modular and synergistic mode of binding.

Results

Comparative evaluation of the three strategies reveals that the synergistic mode of binding appears to mimic the ideal mechanism of DNA-zinc finger protein binding. Analysis of the predictions made by all three strategies indicate strong dependence of zinc finger binding specificity on the amino acid propensity and the position of a 3-bp DNA sub-site in the target DNA sequence. Moreover, the binding affinity of the individual zinc fingers was found to increase in the order Finger 1 < Finger 2 < Finger 3, thus confirming the cooperative effect.

Conclusions

Our analysis offers novel insights into the prediction of ZFPs for target DNA sequences and the approaches have been made available as an easy to use web server at http://web.iitd.ac.in/~sundar/zifpredict_ihbe