Application of docking methodologies to modeled proteins

Protein docking is essential for structural characterization of protein interactions. Besides providing the structure of protein complexes, modeling of proteins and their complexes is important for understanding the fundamental principles and specific aspects of protein interactions. The accuracy of protein modeling, in general, is still less than that of the experimental approaches. Thus, it is important to investigate the applicability of docking techniques to modeled proteins. We present new comprehensive benchmark sets of protein models for the development and validation of protein docking, as well as a systematic assessment of free and template-based docking techniques on these sets. As opposed to previous studies, the benchmark sets reflect the real case modeling/docking scenario where the accuracy of the models is assessed by the modeling procedure, without reference to the native structure (which would be unknown in practical applications). We also expanded the analysis to include docking of protein pairs where proteins have different structural accuracy. The results show that, in general, the template-based docking is less sensitive to the structural inaccuracies of the models than the free docking. The near-native docking poses generated by the template-based approach, typically, also have higher ranks than those produces by the free docking (although the free docking is indispensable in modeling the multiplicity of protein interactions in a crowded cellular environment). The results show that docking techniques are applicable to protein models in a broad range of modeling accuracy. The study provides clear guidelines for practical applications of docking to protein models.     

Seasonal variation in leaf traits between congeneric savanna and forest trees in Central Brazil: implications for forest expansion into savanna

The ecology of forest and savanna trees species will largely determine the structure and dynamics of the forest–savanna boundaries, but little is known about the constraints to leaf trait variation imposed by selective forces and evolutionary history during the process of savanna invasion by forest species. We compared seasonal patterns in leaf traits related to leaf structure, carbon assimilation, water, and nutrient relations in 10 congeneric species pairs, each containing one savanna species and one forest species. All individuals were growing in dystrophic oxisols in a fire-protected savanna of Central Brazil. We tested the hypothesis that forest species would be more constrained by seasonal drought and nutrient-poor soils than their savanna congeners. We also hypothesized that habitat, rather than phylogeny, would explain more of the interspecific variance in leaf traits of the studied species. We found that throughout the year forest trees had higher specific leaf area (SLA) but lower integrated water use efficiency than savanna trees. Forest and savanna species maintained similar values of predawn and midday leaf water potential along the year. Lower values were measured in the dry season. However, this was achieved by a stronger regulation of stomatal conductance and of CO₂ assimilation on an area basis (A _area) in forest trees, particularly toward the end of the dry season. Relative to savanna trees, forest trees maintained similar (P, K, Ca, and Mg) or slightly higher (N) leaf nutrient concentrations. For the majority of traits, more variance was explained by phylogeny, than by habitat of origin, with the exception of SLA, leaf N concentration, and A _area, which were apparently subjected to different selective pressures in the savanna and forest environments. In conclusion, water shortage during extended droughts would be more limiting for forest trees than nutrient-poor soils.     

Muscle Logic: New Knowledge Resource for Anatomy Enables Comprehensive Searches of the Literature on the Feeding Muscles of Mammals

Background

In recent years large bibliographic databases have made much of the published literature of biology available for searches. However, the capabilities of the search engines integrated into these databases for text-based bibliographic searches are limited. To enable searches that deliver the results expected by comparative anatomists, an underlying logical structure known as an ontology is required.

Development and Testing of the Ontology

Here we present the Mammalian Feeding Muscle Ontology (MFMO), a multi-species ontology focused on anatomical structures that participate in feeding and other oral/pharyngeal behaviors. A unique feature of the MFMO is that a simple, computable, definition of each muscle, which includes its attachments and innervation, is true across mammals. This construction mirrors the logical foundation of comparative anatomy and permits searches using language familiar to biologists. Further, it provides a template for muscles that will be useful in extending any anatomy ontology. The MFMO is developed to support the Feeding Experiments End-User Database Project (FEED, https://feedexp.org/), a publicly-available, online repository for physiological data collected from in vivo studies of feeding (e.g., mastication, biting, swallowing) in mammals. Currently the MFMO is integrated into FEED and also into two literature-specific implementations of Textpresso, a text-mining system that facilitates powerful searches of a corpus of scientific publications. We evaluate the MFMO by asking questions that test the ability of the ontology to return appropriate answers (competency questions). We compare the results of queries of the MFMO to results from similar searches in PubMed and Google Scholar.

Results and Significance

Our tests demonstrate that the MFMO is competent to answer queries formed in the common language of comparative anatomy, but PubMed and Google Scholar are not. Overall, our results show that by incorporating anatomical ontologies into searches, an expanded and anatomically comprehensive set of results can be obtained. The broader scientific and publishing communities should consider taking up the challenge of semantically enabled search capabilities.     

Exocyst SEC3 and Phosphoinositides Define Sites of Exocytosis in Pollen Tube Initiation and Growth

Polarized exocytosis is critical for pollen tube growth, but its localization and function are still under debate. The exocyst vesicle-tethering complex functions in polarized exocytosis. Here, we show that a sec3a exocyst subunit null mutant cannot be transmitted through the male gametophyte due to a defect in pollen tube growth. The green fluorescent protein (GFP)-SEC3a fusion protein is functional and accumulates at or proximal to the pollen tube tip plasma membrane. Partial complementation of sec3a resulted in the development of pollen with multiple tips, indicating that SEC3 is required to determine the site of pollen germination pore formation. Time-lapse imaging demonstrated that SEC3a and SEC8 were highly dynamic and that SEC3a localization on the apical plasma membrane predicts the direction of growth. At the tip, polar SEC3a domains coincided with cell wall deposition. Labeling of GFP-SEC3a-expressing pollen with the endocytic marker FM4-64 revealed the presence of subdomains on the apical membrane characterized by extensive exocytosis. In steady-state growing tobacco (Nicotiana tabacum) pollen tubes, SEC3a displayed amino-terminal Pleckstrin homology-like domain (SEC3a-N)-dependent subapical membrane localization. In agreement, SEC3a-N interacted with phosphoinositides in vitro and colocalized with a phosphatidylinositol 4,5-bisphosphate (PIP₂) marker in pollen tubes. Correspondingly, molecular dynamics simulations indicated that SEC3a-N associates with the membrane by interacting with PIP₂. However, the interaction with PIP₂ is not required for polar localization and the function of SEC3a in Arabidopsis (Arabidopsis thaliana). Taken together, our findings indicate that SEC3a is a critical determinant of polar exocytosis during tip growth and suggest differential regulation of the exocytotic machinery depending on pollen tube growth modes.Pollen tube growth provides a unique model system for studying the role of exocytosis in cell morphogenesis. Pollen tubes are characterized by a highly rapid polarized unidirectional tip growth. Given the relative simplicity of their structure, fast growth rates, haploid genome content, and ability to grow under in vitro culture conditions, pollen tubes provide an extremely attractive system for studying cell morphogenesis. Furthermore, the growth characteristics of pollen tubes resemble those of root hairs, moss protonema, and fungal hyphae and to some extent can be paralleled to neurite growth (Chebli and Geitmann, 2007; Cheung and Wu, 2008; Guan et al., 2013; Hepler and Winship, 2015).It is well established that oscillating polarized exocytosis is fundamental for pollen tube development and determines growth rate (Bove et al., 2008; McKenna et al., 2009; Chebli et al., 2013). Exocytosis is required for the delivery of membrane and cell wall components to the growing tip. Yet, the exact location where exocytosis takes place is under debate. Ultrastructural studies showing the accumulation of vesicles at the tip suggested that exocytosis takes place at the tip (Lancelle et al., 1987; Lancelle and Hepler, 1992; Derksen et al., 1995), which was further supported by studies on the dynamics of cell wall thickness (Rojas et al., 2011), secretion of pectin methyl esterase (PME) and PME inhibitor, and staining of pectin by propidium iodide (PI; Röckel et al., 2008; Rounds et al., 2014). Conversely, based on colabeling with FM1-43 and FM4-64, it was concluded that exocytosis takes place in a subapical collar located in the transition zone between the tip and the shank, as well as at the shank, but not at the tip (Bove et al., 2008; Zonia and Munnik, 2008). In agreement, the pollen tube-specific syntaxin GFP-SYP124 was observed in the inverted cone, 10 to 25 μm away from the tip (Silva et al., 2010), and fluorescence recovery after photobleaching experiments with FM dyes also have indicated that exocytosis takes place at the subapical region (Bove et al., 2008; Moscatelli et al., 2012; Idilli et al., 2013). Yet, based on pollen tube reorientation experiments in a microfluidics device, it was concluded that growth takes place at the tip rather than at a subapical collar located in the transition zone between the apex and the shank (Sanati Nezhad et al., 2014). The tip-based growth is in agreement with exocytosis taking place at the tip. Presumably, part of the disagreement regarding the site of exocytosis resulted from the lack of intracellular markers for exocytosis (Cheung and Wu, 2008; Hepler and Winship, 2015), and as a result, the relationship between the FM dye-labeled inverted cone and exocytotic events during pollen tube growth is not fully understood.In many cell types, the process of secretory vesicles tethering and docking prior to fusion with the plasma membrane is initially mediated by an evolutionarily conserved tethering complex known as the exocyst. The exocyst is a heterooligomeric protein complex composed of eight subunits, SEC3, SEC5, SEC6, SEC8, SEC10, SEC15, EXO70, and EXO84 (TerBush et al., 1996; Guo et al., 1999). Studies originally based on budding yeast (Saccharomyces cerevisiae) have shown that the exocyst functions as an effector of Rab and Rho small GTPases that specifies the sites of vesicle docking and fusion at the plasma membrane in both space and time (Guo et al., 2001; Zhang et al., 2001). Support for the function of the exocyst in vesicle tethering was demonstrated recently by ectopic Sec3p-dependent vesicle recruitment to the mitochondria (Luo et al., 2014).Land plants contain all subunits of the exocyst complex, which were shown to form the functional complex (Elias et al., 2003; Cole et al., 2005; Synek et al., 2006; Hála et al., 2008). Studies in Arabidopsis (Arabidopsis thaliana) and maize (Zea mays) have implicated the exocyst in the regulation of pollen tube and root hair growth, seed coat deposition, response to pathogens, cytokinesis, and meristem and stigma function (Cole et al., 2005; Synek et al., 2006; Hála et al., 2008; Fendrych et al., 2010; Kulich et al., 2010; Pecenková et al., 2011; Safavian and Goring, 2013; Wu et al., 2013; Safavian et al., 2015; Zhang et al., 2016). The growth arrest of pollen tubes in sec8, sec6, sec15a, and sec5a/sec5b single and double mutants (Cole et al., 2005; Hála et al., 2008) or following treatment with the EXO70 inhibitor ENDOSIDIN2 (Zhang et al., 2016), and of root hairs in maize root hairless1 (rth1) SEC3 mutant (Wen et al., 2005), the inhibition of seed coat deposition in the sec8 and exo70A1 mutants (Kulich et al., 2010), and stigmatic papillae function in exo70A1 mutant plants (Safavian and Goring, 2013; Safavian et al., 2015) have implicated the exocyst in polarized exocytosis in plants. Given their function, it was likely that exocyst subunits could be used as markers for polarized exocytosis. Furthermore, it could also be hypothesized that, by studying the mechanisms that underlie the association of the exocyst complex with the plasma membrane, it should be possible to identify mechanisms underlying the regulation of polarized exocytosis (Guan et al., 2013). Moreover, since the interaction of exocytotic vesicles with the exocyst is transient and marks the site(s) of active exocytosis in the membrane, fluorescently labeled exocyst subunits could be used as markers for exocytosis while avoiding potential imaging artifacts stemming from pollen tube tips densely populated with vesicles.We have shown previously that the ROP effector ICR1 can interact with SEC3a and that ROPs can recruit SEC3a-ICR1 complexes to the plasma membrane (Lavy et al., 2007). However, ICR1 is not expressed in pollen tubes, suggesting that SEC3a membrane binding in these cells is likely dependent on other factors. In yeast, the interaction of Sec3p and Exo70p subunits with the plasma membrane is critical for exocyst function (He and Guo, 2009). It has been shown that the membrane binding of both Sec3p and Exo70p is facilitated by their interaction with phosphatidylinositol 4,5-bisphosphate (PIP₂; He et al., 2007; Zhang et al., 2008). The yeast Exo70p interacts with PIP₂ via a number of positively charged residues distributed along the protein, with the highest number located at the C-terminal end (Pleskot et al., 2015). It has been suggested that yeast Sec3p interacts with PIP₂ through N-terminal basic residues (Zhang et al., 2008). These data were further corroborated by x-ray crystallography studies, which showed that the yeast Sec3p N-terminal region forms a Pleckstrin homology (PH) domain fold (Baek et al., 2010; Yamashita et al., 2010), a PIP₂ interaction motif (Lemmon, 2008).The localization of the exocyst subunits has been addressed in several studies. In Arabidopsis root hairs and root epidermis cells, SEC3a-GFP was observed in puncta distributed throughout the cell (Zhang et al., 2013). Studies on the Arabidopsis EXO70 subunits EXO70E2, EXO70A1, and EXO70B1 revealed them to be localized in distinct compartments that were termed exocyst-positive organelles (Wang et al., 2010). The exocyst-positive organelles, visualized mostly by ectopic expression, were shown to be cytoplasmic double membrane organelles that can fuse with the plasma membrane and secrete their contents to the apoplast in an exosome-like manner. It is not yet known whether other exocyst subunits also are localized to the same organelles and what might be the biological function of this putative compartment (Wang et al., 2010; Lin et al., 2015). In differentiating xylem cells, two coiled-coil proteins termed VESICLE TETHERING1 and VESICLE TETHERING2 recruit EXO70A1-positive puncta to microtubules via the GOLGI COMPLEX2 protein (Oda et al., 2015). Importantly, the functionality of the XFP fusion proteins used for the localization studies described above was not tested, and in most cases, the fusion proteins were overexpressed. Therefore, the functional localization of the exocyst is still unclear.Here, we studied the function and subcellular localization of the Arabidopsis exocyst SEC3a subunit using a combination of genetics, cell biology, biochemistry, and structural modeling approaches. Our results show that SEC3a is essential for the determination of pollen tube tip germination site and growth. Partial complementation of sec3a resulted in the formation of pollen with multiple pollen tube tips. In Arabidopsis growing pollen tubes, SEC3a localization is dynamic, and it accumulates in domains of polarized secretion, at or close to the tip plasma membrane (PM). Labeling of GFP-SEC3-expressing pollen with FM4-64 revealed the spatial correlation between polarized exocytosis and endocytic recycling. Furthermore, the association of SEC3a with PM at the tip marks the direction of tube elongation and positively correlates with the deposition of PI-labeled pectins and specific anti-esterified pectin antibodies in the cell wall. In tobacco (Nicotiana tabacum), the mechanisms underlying SEC3a interaction with the PM and its subcellular distribution depend on pollen tube growth mode and involve the interaction with PIP2 through the N-terminal PH domain. Collectively, our results highlight the function of SEC3a as a polarity determinant that links between polarized exocytosis and cell morphogenesis. The correlation between exocyst function and distribution in pollen tubes provides an explanation for some of the current discrepancies regarding the localization of exocytosis.     

Digital image analysis of AgNORs in cervical smears of women with premalignant and malignant lesions of the uterine cervix

We performed a hospital-based, unmatched case-control study to investigate the association between progressive stages of cervical neoplasia and digital analysis of cell proliferation by silver stained nucleolus organizer region associated proteins (AgNORs). We measured cell proliferation levels in the cervical epithelial cells of 10 women with low grade squamous intraepithelial lesions (LG-SIL), eight with high grade squamous intraepithelial lesions (HG-SIL), 11 with cervical cancer (CC) and eight with no cervical lesions (controls) using the AgNORs technique. Cell proliferation was measured by digital image analysis (DIA). DIA revealed increased total areas of AgNORs in HG-SIL and CC compared to LG-SIL and control patients. AgNORs with a kidney or cluster shape exhibited greater areas than those with a spherical or long shape. We propose a cut-off of 118 pixels to differentiate benign (control and LG-SIL) from malignant (HG-SIL and CC) lesions. DIA of AgNORs is a simple and inexpensive method for studying proliferation. The increased total area of AgNORs in malignant lesions provides information regarding cell behavior and may be related to cervical carcinogenesis; however, further validation studies are required to establish its usefulness in cytological analysis.     

Automated discovery of structural signatures of protein fold and function

There are constraints on a protein sequence/structure for it to adopt a particular fold. These constraints could be either a local signature involving particular sequences or arrangements of secondary structure or a global signature involving features along the entire chain. To search systematically for protein fold signatures, we have explored the use of Inductive Logic Programming (ILP). ILP is a machine learning technique which derives rules from observation and encoded principles. The derived rules are readily interpreted in terms of concepts used by experts. For 20 populated folds in SCOP, 59 rules were found automatically. The accuracy of these rules, which is defined as the number of true positive plus true negative over the total number of examples, is 74% (cross-validated value). Further analysis was carried out for 23 signatures covering 30% or more positive examples of a particular fold. The work showed that signatures of protein folds exist, about half of rules discovered automatically coincide with the level of fold in the SCOP classification. Other signatures correspond to homologous family and may be the consequence of a functional requirement. Examination of the rules shows that many correspond to established principles published in specific literature. However, in general, the list of signatures is not part of standard biological databases of protein patterns. We find that the length of the loops makes an important contribution to the signatures, suggesting that this is an important determinant of the identity of protein folds. With the expansion in the number of determined protein structures, stimulated by structural genomics initiatives, there will be an increased need for automated methods to extract principles of protein folding from coordinates.     

An imaging system for standardized quantitative analysis of <Emphasis Type="Italic">C. elegans</Emphasis> behavior

Background  

The nematode Caenorhabditis elegans is widely used for the genetic analysis of neuronal cell biology, development, and behavior. Because traditional methods for evaluating behavioral phenotypes are qualitative and imprecise, there is a need for tools that allow quantitation and standardization of C. elegans behavioral assays.     

Gene expression markers for Caenorhabditis elegans vulval cells

The analysis of cell fate patterning during the vulval development of Caenorhabditis elegans has relied mostly on the direct observation of cell divisions and cell movements (cell lineage analysis). However, reconstruction of the developing vulva from EM serial sections has suggested seven different cell types (vulA, vulB1, vulB2, vulC, vulD, vulE, and vulF), many of which cannot be distinguished based on such observations. Here we report the vulval expression of seven genes, egl-17, cdh-3, ceh-2, zmp-1, B0034.1, T04B2.6 and F47B8.6 based on gfp, cfp and yfp (green fluorescent protein and color variants) reporter fusions. Each gene expresses in a specific subset of vulval cells, and is therefore useful as a marker for vulval cell fates. Together, expressions of markers distinguish six cell types, and reveal a strict temporal control of gene expression in the developing vulva.