Skip to main content

Integration of the Gene Ontology into an object-oriented architecture

Abstract

Background

To standardize gene product descriptions, a formal vocabulary defined as the Gene Ontology (GO) has been developed. GO terms have been categorized into biological processes, molecular functions, and cellular components. However, there is no single representation that integrates all the terms into one cohesive model. Furthermore, GO definitions have little information explaining the underlying architecture that forms these terms, such as the dynamic and static events occurring in a process. In contrast, object-oriented models have been developed to show dynamic and static events. A portion of the TGF-beta signaling pathway, which is involved in numerous cellular events including cancer, differentiation and development, was used to demonstrate the feasibility of integrating the Gene Ontology into an object-oriented model.

Results

Using object-oriented models we have captured the static and dynamic events that occur during a representative GO process, "transforming growth factor-beta (TGF-beta) receptor complex assembly" (GO:0007181).

Conclusion

We demonstrate that the utility of GO terms can be enhanced by object-oriented technology, and that the GO terms can be integrated into an object-oriented model by serving as a basis for the generation of object functions and attributes.

Background

Complexity combined with an imprecise terminology has hindered the understanding of biology. A formal and structured vocabulary is now being developed to address this imprecise biology terminology. This vocabulary or Gene Ontology (GO) is being developed by the Gene Ontology Consortium (GOC) [1] to standardize the descriptions of gene products. Ontologies define the basic terms and relations comprising the vocabulary of a topic area, as well as the rules for combining terms and relations to define extensions to the vocabulary [2]. Despite these efforts, the mechanism of representing these terms lacks a unifying architecture that can be applied to the annotation of a gene product. However, computer science has developed a well-defined process and methodology for the development of software models. Adapting this process and methodology can orchestrate the assembly of biological models with integrated gene ontologies. In doing so, a standardized terminology and object-oriented model is created that can facilitate communication between biologists and computer scientists.

The Gene Ontology project is a collaborative effort that addresses the need for a controlled vocabulary that provides a consistent description of gene products in different databases [1]. The GO collaborators are developing three structured, controlled vocabularies that describe gene products, which have been classified into molecular function, biological process, and cellular component domains. GO terms are organized in structures called directed acyclic graphs (DAGs), which differ from hierarchies in that a 'child' (more specialized term) can have many 'parents' (less specialized terms). As part of these graphs, each component is given a GOid (unique identifier), and is associated with a GO definition. Collectively, these agreed upon terms are being developed to help explain various aspects of biology. When applied to a gene, that gene is annotated with a concise description of its molecular function, cellular location and associated biological processes. However, the GOC never intended to represent gene products or correlate ontological terms with these gene products [1]. To address this need, a Gene Ontology Annotation database [3] has been created to associate the GO terms with their gene product counterparts. With sustained effort, the descriptions of these gene products will ultimately be established. Still, much of the current bioinformatics work regarding GO has focused on constructing databases [4–7], applying it to other research areas [8–22], and building tools to mine the GO database. (For a description of some of these tools see [23].)

In addition, there has been an ongoing discussion regarding the depth of information obtained from the Gene Ontology [24]. It has been noted that there remains a need for a unifying architecture that integrates all three GO domains as part of a gene product's annotation. Furthermore, to enhance the Gene Ontology and facilitate its use as a cross-disciplinary tool, several additional issues need to be addressed. First, relationships between the biological processes, molecular functions and cellular components are not readily apparent [25–28]. Second, GO terms lack details. For instance, when one looks at molecular function there is no indication of what is inputted or outputted. Finally, existing tools such as GO-DEV [29] only contain software used for tool development and information retrieval, not software modeled directly after the three domains of the Gene Ontology. However, these issues can be resolved by integrating the Gene Ontology into an object-oriented system.

On a conceptual level, the Gene Ontology has features that support an object-oriented architecture. Consequently, the Gene Ontology can be applied and mapped to the fundamental concepts that form the object-oriented paradigm (i.e. class, object, inheritance, composition, polymorphism, and encapsulation) (Table 1). Furthermore, in an object-oriented sense, biological process terms are equivalent to high-level concepts. However, GO biological process terms do not contain descriptive information about the dynamics or static interactions defined by the terms. By translating a biological process into an object-oriented model the dynamic and static events occurring within a process can be represented. Building a static and dynamic model of a biological process requires defining the components of the process as well as the functions and attributes contained within these components. These components are biological entities (bioentities) that may include individual gene products, whose processes, functions and cellular components are captured in the Gene Ontology, or other higher-level entities such as gene product complexes.

Table 1 The use of object-oriented concepts in the integration of the Gene Ontology into an object-oriented model. Object-oriented terms, their definitions, and corresponding mechanisms of incorporating GO terms into an object-oriented model are shown. A specific example from the manuscript is also given. GO, Gene Ontology; DAG, directed acyclic graph; OOM, object-oriented model

The functions of gene products are the jobs or abilities that it has. In the GO terminology these are described in the molecular function domain. These are analogous to the operations that an object can perform in an object-oriented paradigm. Attributes, which define key properties of a component that when changed may alter the function of that component, may be defined by the cellular component and molecular function sections. For example, the cellular component domain can specify the place in a cell where a gene product is located. When there are multiple cellular components associated with a gene product, however, there is currently no mechanism to designate which cellular component represents the appropriate location.

The unified modeling language has been used to capture various aspects of biology [30–32]. These examples highlight the utility of the unified modeling language as a tool for biological data integration, and indicate that it can be applied to construct large, complex biological models. Therefore, to demonstrate the feasibility of integrating the Gene Ontology into an object-oriented model we have created unified modeling language (UML) representations of a GO biological process, "transforming growth factor beta (TGF-beta) receptor complex assembly" (GO:0007181).

The TGF-beta receptor pathway is involved in numerous cellular events including apoptosis, tumor development, differentiation, and development. These processes stem from the binding of TGF-beta to its cellular receptors. Briefly, dimerized TGF-beta 1 binds to TGF-beta receptor II (TβRII) and then TGF-beta receptor I (TβRI) complexes [33], causing their tetramerization (two type II receptors and two type I receptors) [34–36]. Constitutively activated type II receptor phosphorylates and activates type I receptor. Type I receptor propagates the signal by phosphorylating Smad 2, which is presented by the Smad Anchor for Receptor Activation (SARA) [37]. Phosphorylation of Smad destabilizes the Smad interaction with SARA, releasing it [38]. On TGF-beta stimulation, Smad 2 forms heterotrimeric complexes with Smad 4 and accumulates in the nucleus, binds DNA and remains for several hours [39–42]. Dephosphorylation allows Smad 2 to dissociate from Smad 4 and to be exported to the cytoplasm [43, 44]. If the receptors are no longer active, then the Smads accumulate over time in the cytoplasm [44]. Alternatively, activated Smad 2 is ubiquitinated in the nucleus and undergoes proteasome-mediated degradation [45].

To create a unified model using the Gene Ontology we have taken the biological process term, "transforming growth factor beta (TGF-beta) receptor complex assembly" (GO:0007181), and used object-oriented models to define its dynamic and static architecture. We also show that one can augment the biological process domain terms by using the ontological terms and gene products associated with this process, and integrating them into an object-oriented model. Furthermore, we show that the molecular function, and cellular component domains can serve as a basis for the generation of object functions and attributes to create a standardized, comprehensive, and integrated model encompassing all the Gene Ontology domains.

Results

Converting GO directed acyclic graphs to object-oriented diagrams

The current DAG structure in which the Gene Ontology is represented is not readily amenable to transformation into software code. However, the architecture of directed acyclic graphs mimics that of an object-oriented class diagram. GO terms are presented in a parent-child hierarchy connected by 'is a' (generalizations) and 'part of' (composition) relationships. Read from top to bottom, the GO terms proceed from more specific to less specific. Directed acyclic graphs also allow the properties of multiple parent nodes to be inherited by child nodes, a form of multiple inheritance in object-oriented modeling. In figure 1A, cellular components related to the TGF-beta receptor complex are shown. One can create a UML diagram to mimic these relationships as shown in figure 1B. Since not all cellular components involved in the TGF-beta receptor complex assembly process are present in the current Gene Ontology, additional gene products based on literature searches were added to the object-oriented diagram (Figure 1B, shaded boxes). Relationships are captured as an object-oriented system through containment, composition and inheritance. Cellular components were decomposed into objects and connected via generalizations, which illustrate inheritance. Because these terms inherit all the attributes of their parents, only GO terms at the terminal nodes need be characterized in an object-oriented model. Relationships described as 'part of' were also extrapolated into object-oriented terms as composition.

Figure 1
figure 1

The GO terms associated with the process, TGF-beta receptor complex assembly (GO: GO:0007181). Lines with solid diamonds () at the end indicate composition. These are read from the diamond end, for example, as 'cell' (GO:0005623) contains a 'membrane' (GO:0016020). Lines with open triangles () represent generalizations. These are read from the triangle end, as a 'membrane' (GO:0016020) is a general type of 'plasma membrane' (GO:0005886). A) Directed acyclic graph for the cellular component GO terms associated with TGF-beta receptor complex assembly (GO:0007181). B) Object-oriented representation of the DAG described in Figure 1A (white). Additional cellular components not represented by the current Gene Ontology, but essential to the TGF-beta receptor complex assembly process are shown in the gray boxes (i.e. TGF-beta, TGF-beta receptors, SMAD2)

The functions of gene products were also decomposed into object functions. The creation of object functions involved the transition from gene product functions to standardized GO molecular functions, and then to standardized, fully parameterized object functions. By applying formal ontological terms from the molecular function domain to gene products, object functions can be created with a consistent vocabulary. In table 2 we show the relationships between the function of a gene product defined in our model, and the GO molecular function term most closely corresponding to that cellular function. Here, we first compared ontology terms from the molecular function domain to those ascribed to individual gene products. Due to the incompleteness of the Gene Ontology, some gene product functions were extrapolated from the current literature, and then comparable GO molecular function terms were assigned to the gene products. Next, these molecular function terms were converted to object functions through reverse engineering. We identified the parameters that would normally be input into and output from a cellular reaction. In this way we defined the input and output parameters necessary for an object function. The object function itself was given the GOid that corresponds with its closest matching molecular function as defined by the GOid's definition. Together, object functions were created that are fully parameterized with inputs and outputs and that contain a standardized GO notation.

Table 2 The gene product functions described herein are listed with their associated GO molecular functions and parameters. These gene product functions are mapped to corresponding Gene Ontology molecular functions. These GO functions are integrated into an object-oriented model by amending them with input and output parameters, thereby creating object functions.

We conclude that it is feasible to create standardized functions for objects based on the current literature and an approved ontology. Together, ontological terms can be integrated into an object-oriented model paralleling the relationships, capturing the inherited aspects of the GO terminology, and providing a compact architecture while maintaining a standardized notation.

Sequence diagram generation

The GO biological process term, TGF-beta receptor complex assembly (GO:0007181), contains both static and dynamic features. The events of the TGF-beta receptor complex assembly (GO:0007181) process include TGF-beta binding (GO:0050431) to its receptors and SMAD binding (GO:0046332) and activation (GO:0042301). To capture the dynamic nature of these actions as an object-oriented software system, sequence diagrams were created. The events leading to Smad 2 activation are reflected chronologically in a high-level sequence diagram in Figure 2. The creation of the sequence diagram first entails identifying gene products and their functions by literature searches. Simple or complex bioentities are modeled as objects, which are represented by rectangles with vertical lifelines. Ontology terms taken from the molecular function domain that best corresponded to these functions were incorporated as object functions, which represent the functions of these gene products. These functions are implemented by the methods contained within the objects. Furthermore, these methods allow an object to communicate and interact with other objects, thus capturing cellular activities. To capture interactions between objects, one object can call a method of another object by connecting object lifelines in the sequence diagram (Figure 2). This invocation of a function of one object by another is described as one object sending a message to another object. Alternatively, a message may be passed from an object to itself as in the case of self-checks or autoactivation signals. In this way, real world processes may be captured using an object-oriented approach. For instance, to capture the formation of the TGF-beta and TGF-beta RII complex a GOid that closely corresponds to this ability is chosen as the method name. In this way the method can be cross-referenced to a GO term. Specifically, the method 'GO:0046982 (in Dimerized_TGF-beta, in Dimerized_RII)' references via the GOid, GO:0046982, "protein heterodimerization activity", and shows that a homodimer of TGF-beta and a homodimer of TGF-beta RII are needed to form the complex. Here, each dimer is thought of as a single entity, so the combination of these two entities is best represented as heterodimerization. A value of TGF-beta-TβRII_Complex is returned upon completion of the method as indicated by the return arrow. In contrast, the function call "GO:0042803 (in: SMAD2, in SMAD2)," references a self-call. The GOid can be cross-referenced to "protein homodimerization activity", which requires two SMAD2 components to generate the SMAD2 homodimer, but the message is passed only within the SMAD2 object. Furthermore, a message need not accept any parameters, as in the "translate()" function, which only returns a Boolean value indicating whether the action has occurred. Additional events such as TGF-beta RI activation, and Smad homodimerization, binding and activation are also reflected in figure 2. Together, this diagram demonstrates that the sequence of events occurring in the biological process, TGF-beta receptor complex assembly (GO:0007181), can be represented using the Gene Ontology, and can be integrated as part of the dynamics of an object-oriented software system.

Figure 2
figure 2

An example of the sequence diagram showing the TGF-beta receptor complex assembly (GO:0007181). The binding of TGF-beta to its receptor (GO:0050431), receptor heterotetramerization (RI and RII homodimers, heterodimerizing)(GO:0046982), translation (GO:0043037), transferase activity (GO:0016740), and Smad 2 binding (GO:0046332) and activation (GO:0042301) are shown.

Activity diagram generation

Biological processes are created from a series of complex events. While there may be one main event scenario that most frequently leads to a specific outcome often, alternative scenarios that lead to a process conclusion exist. This is exemplified by the sequence of events found in the TGF-beta receptor complex assembly (GO:0007181). For instance, TGF-beta may initially bind to TGF-beta RII or TGF-beta RIII. To capture these alternative events as part of the dynamic architecture, an activity diagram was created to reflect the initial stages of TGF-beta signaling (Figure 3). Unlike the sequence diagram, which captures main scenario events, the action sequence or flow of the activity diagram can portray alternative outcomes. Taking the example above, if TGF-beta binds to the type III receptor then an alternative flow of events occurs for a time that then returns to the main flow of events. Other possible divergences that were modeled included whether to internalize the TGF-beta receptors via clathrin-dependent or lipid raft-dependent mechanisms. These pathways lead to either complex degradation or signal promotion. Because complex degradation is not specified in our use case, for simplicity, this event is routed to the final state. However, the main success scenario, signal promotion, continues until SMAD2 is released and TGF-beta complex assembly is finished. Together, the dynamic events occurring during the biological process, TGF-beta receptor complex assembly (GO:0007181) are captured.

Figure 3
figure 3

An example of an activity diagram showing the main and alternative flow of events occurring during TGF-beta receptor complex assembly (GO:0007181).

Class diagram generation

The major components of a biological system are bioentities with functions and interactions. Likewise, the center of an object-oriented software system is objects. Complex bioentities formed from multiple gene products along with their relationships, are contained within the biological system encompassing the biological process term, TGF-beta receptor complex assembly (GO:0007181). To represent the components that execute the process, we captured these components as bioentities with functions, and their interactions. The events of the TGF-beta receptor complex assembly (GO:0007181) process include TGF-beta binding (GO:0050431) to its receptors, and SMAD binding (GO:0046332) and activation (GO:0042301). To capture this static architecture, class diagrams were generated that model the bioentities, operations, and interrelationships that occur between TGF-beta, its receptors, and Smad 2. Similarly to figure 1, figure 4 captures the major components of the initial phases of TGF-beta signaling as objects with their associations, using an object-oriented representation. However, unlike figure 1, this object-oriented representation of the components of the main receptor complex is enhanced by the addition of attributes and functions. These objects were given attributes that describe important characteristics that if changed, might alter the function of a component. The functions of the objects, which parallel gene product functions, were generated from the sequence diagrams and were represented using Gene Ontology terms. These functions or operations are a declaration of the methods that an object may use. Together, the models generated using the described object-oriented methodology yield a software system representation of a biological process, TGF-beta receptor complex assembly, capturing both static and dynamic relationships annotated with Gene Ontology terms.

Figure 4
figure 4

An example of a class diagram showing the interactions between the components of the TGF-beta receptor complex (GO:0007181) (grayed). Cellular components containing these gene products are also shown. Due to space constraints functions are named using their corresponding GO ids. The class diagram was also constructed with a descriptive function name [see Additional file 8]. The GOid reference list is: GO:0042803, protein homodimerization activity; GO:0005160, TGF beta receptor binding; GO:0046982, protein heterodimerization activity, GO:0050431, TGF beta binding, GO:0005524, ATP binding; GO:0016740, transferase activity, GO:0046332, Smad binding; GO:0042301, phosphate binding; GO:0003677, DNA binding.

In addition, the UML notation provides a mechanism to specify inheritance that may be used to indicate an object that is the foundation for other objects. For instance, a TGF-beta receptor object might be a generalization of the TGF-beta receptor I object (data not shown). These specific objects inherit the properties of the receptor object. In addition, binary associations containing cardinalities may indicate the number of objects interacting with another. For instance, TGF-beta can interact with one to many receptors, while a receptor can only interact with one TGF-beta at a time (Fig. 4). Cellular compartments where these gene products can be found are also shown. Here, guard conditions are added to distinguish conditions under which each gene product might be found in a particular cellular compartment. In this way, a spatial representation of the TGF-beta receptor complex components is also achieved. These class diagrams demonstrate that the static structure of a biological system can be represented as an object-oriented model with integrated Gene Ontology terms. Collectively, the models generated using the described object-oriented methodology yield a software system representation of a biological system, capturing both static and dynamic relationships annotated with integrated Gene Ontology terms.

Discussion

We have utilized the Gene Ontology to construct an object-oriented representation of the initial steps of TGF-beta signaling, and the gene products contained therein. In doing so, we have provided a standardized framework for the integration of Gene Ontology terms into gene product descriptions. By capturing all of the relevant GO terms in one model, the disjointed GO vocabulary is assembled into a cohesive structure. This cohesive structure encompasses the fundamental concepts of the object-oriented paradigm.

We proposed a solution to three unaddressed issues within the current Gene Ontology. First, while the Gene Ontology has helped to formalize the vocabulary that describes biological systems, it lacks a specific integration method. Currently, when applied to gene products, Gene Ontology terms are only categorically listed. Second, the Gene Ontology domains, biological process, molecular function and cellular component lack coherence. In particular, no association exists between domains. Finally, the current Gene Ontology defines GO terms, but gives no indication of what is necessary to accomplish a particular function, or process. To resolve these problems we defined an object-oriented methodology and architecture that provides a unifying framework to integrate all Gene Ontology domains.

The central dogma of the object-oriented paradigm revolves around several key aspects. Specifically, an object-oriented framework should accommodate the class, object, inheritance, composition, encapsulation and polymorphism concepts. As shown in table 1, gene products and other bioentities can be decomposed into objects, which are created based on template classes. These objects utilize inheritance to acquire the attributes and properties of more general objects. Complex classes can also be disassembled into subclasses using composition. Encapsulation allows the simplification of the model without sacrificing functionality. For instance, we do not need to know specific details regarding how a gene product is translated, just that a process that is encapsulated by the function 'translate()' can create a protein. However, if we wished to delve deeper into the mechanics of the translation process the layered architecture of the object-oriented system would allow us to do so. It is also worth noting that the modular nature of the object-oriented system closely resembles the recently discovered modular structure of biological networks [46–48]. This resemblance further indicates that biological systems can be easily modeled as object-oriented systems. Finally, polymorphism allows one to describe shared functions among different gene products. In this way, a function that may be shared broadly with other gene products can be uniquely specified for a particular gene product.

By applying object-oriented methodologies and concepts the various domains of the Gene Ontology can be coordinated into one model. Currently, the mechanisms in the biological process domain are veiled. There is no indication as to what gene products form the biological process, or what molecular functions are necessary to accomplish the process. Furthermore, the outcome of a specific process is not obvious. As in our example, a process such as TGF-beta receptor complex assembly (GO:0007181) does not give any indication of the components, dynamics or outcomes that occur during this process. However, by incorporating GO terms as attributes and functions we can discern relationships between the three domains. Likewise, the cellular components domain does not provide temporal or spatial clues when applied to gene products. For instance, GO terms 'extracellular' and 'intracellular' may both be associated with a particular gene product. However, the distinction between when a gene product is extracellular and when it is intracellular is not apparent. By applying object-oriented principles we can set extracellular and intracellular to Boolean values, and we can specify which location is the current (true) location of a gene product.

In addition, by using object-oriented principles a GO molecular function term can be augmented with parameters and outcomes. For example, the function "GO:0046982: protein heterodimerization activity" has different input and output parameters depending on the particular protein that contains the function. This type of polymorphic behavior, where one function can be performed in multiple ways is not supported by the Gene Ontology. For example, protein A may heterodimerize with protein B, whereas protein C heterodimerizes with protein D. From the Gene Ontology it is not readily apparent as to what is being inputted into the dimerization function. However, by applying an object-oriented architecture to function "GO:0046982: protein heterodimerization activity" we get "GO:0046982 (in: Protein A, in: Protein B): Protein AB". This format is an improvement to the unparameterized GO term in that the function can be cross-referenced to protein heterodimerization activity via its GO term, and we also see that for protein A to heterodimerize we need both protein A and protein B. In addition, we now observe that a new entity called protein AB is created from this function. By capturing the above details in an object-oriented model the GO term becomes far more useful for both biologists and computer scientists. Using an object-oriented approach the Gene Ontology domains are integrated into one cohesive model.

Integration of the Gene Ontology terms into an object-oriented representation offers several additional benefits. The object-oriented model provides additional levels of detail not found in the Gene Ontology. One of the strengths of object-oriented technology is the ability to capture the dynamics of a system. For example, sequence diagrams can chronologically order events in a biological process. Activity diagrams afford one the opportunity to envision different scenarios that might be occurring in a process. This additional level of detail significantly increases the depth of information that can be applied to the description of a biological process. State-transition diagrams also contribute to the realization of the full dynamics of a process by allowing the visualization of gene product states within a process. Furthermore, UML models can be translated into code, facilitating the creation of simulations.

The standardization of biological system modeling and integration is growing rapidly. A widely accepted example of the drive toward standardization is the Systems Biology Markup Language (SBML) [49], which has been adopted by more than 70 software tools [50]. The Gene Ontology is another example. However, each of the technologies, the Gene Ontology, the object-oriented approach, and SBML, has strengths and weaknesses. The Gene Ontology provides a standardized vocabulary but contains disconnected domains with no details regarding terms. SBML was developed to communicate biological models, with an emphasis on mathematical modeling of biological systems, but does not specify how to construct these models. Object-oriented technologies, on the other hand, provide a well-defined process for model creation and visualization, but have not been standardized for biology. However, the Gene Ontology, object-oriented paradigm, and SBML can form a new synergism when jointly applied to a common biological system model. These technologies are steps toward a unified approach to biological information integration, and studying biological phenomena at the systems level. Together, this unified approach will make biological system integration and analysis consistent, manageable and controllable, which is essential in handling complex systems, as demonstrated by decades of software industry experience.

While the described object-oriented approach can significantly enhance the annotation of gene products using the Gene Ontology, several challenges will need to be addressed. Specifically, object-orientation was not specifically designed for use in biological systems. Therefore, its use in capturing biological systems is not well defined. Furthermore, the Gene Ontology is still expanding and undergoing revisions. Consequently, in the near future it will still be necessary to do literature searches to define all the gene ontologies associated with a gene product. However, automated extraction of information for UML model generation and software implementation for simulations is under development, but is beyond the scope of this paper.

Future systems may also be implemented as software libraries in object-oriented programming languages (C++ and Java) for computer scientists to construct software for various applications and can be distributed as part of the GO-DEV toolkit for Gene Ontology development [29]. In addition, reformatting gene products with Gene Ontology terms will require the cooperation of multiple groups of biologists and computer scientists. However, we must take into consideration that a primary issue with this approach is the lack of people with cross-disciplinary skills able to comprehend both the biology and the computer science. Nonetheless, our own experience has shown that with supervision one biologist without a formal computer science background can learn to model a biological system using UML in a matter of months. Furthermore, automation of some of the annotation process will significantly reduce the human effort, but not eliminate the need for human annotators. Additional standards for automation will also need to be developed to thoroughly specify the process of object-oriented biological system integration. Despite these challenges the ultimate goal of creating a library of UML objects or modules integrated with Gene Ontology attributes and functions is worthwhile. Through this endeavor, biological processes could be assembled from these libraries for the development of simulation tools that will increase the productivity of biologists through increased insight into disease pathways and mechanisms.

Conclusion

Here, we have demonstrated that Gene Ontology terms can be integrated into an object-oriented model. Furthermore, the object-oriented technology and methodologies used for this integration should improve the usability of these terms, and increase the depth of information that they contain. This work also serves as a framework for reverse-engineering biological gene products as objects in an object-oriented system. Together, this should facilitate additional collaborations between biologists and computer scientists.

Methods

UML representations of the TGF-beta receptor complex assembly process were created following a software engineering process consisting of phases of requirement-gathering, analysis and design. UML models were generated using Microsoft Visio Pro. AmiGO [51] was used to determine Gene Ontology links for TGF-beta gene products.

Requirement-gathering phase

Information collection

To define the requirements and collect the information necessary for the generation of the models, two approaches were necessary. First, annotations of the TGF-beta signaling pathways were conducted during an extensive literature review. Second, gene ontologies and Uniprot entries were searched to assign Gene Ontology terms to gene products. The attributes and the interactions of the TGF-beta signaling components were captured using class-responsibility collaboration (CRC) cards as described previously [52] [see Additional files 1, 2, 3, 4].

Use case development

Based on the gathered information, best-case and alternative scenarios were developed within a so-called "use case" to describe the TGF-beta receptor complex assembly process [see Additional file 5]. The use case also serves to define the boundary and scope of the TGF-beta model. For demonstration purposes the boundary of the system was limited to the steps TGF-beta receptor complex assembly. Therefore, alternative events such as receptor ubiquitination and degradation, as well as the specifics of SMAD 2 mobility were not captured in the dynamic models (i.e. sequence diagram).

Analysis phase

Conceptual model generation

To provide an overview of the system and its interrelationships a conceptual based on the information defined in the requirement-gathering phase was generated [see Additional file 6]. This conceptual model integrated biological information, and represented TGF-beta and the cellular components involved in the complex assembly and their relationships in UML notation. By applying object-oriented analysis, the TGF-beta receptor complex assembly was decomposed into objects and component relationships were realized. However, information regarding component properties is hidden through encapsulation. This conceptual model defines the organization of the biological system and provides an overview of the components and their relationships.

Design phase

State diagram generation

The dynamics of the system can also be captured using state diagrams, which can be used to describe the transitions and different states that a cellular component can exist [see Additional file 7]. In addition, multiple concurrent states can be illustrated using this UML notation.

Sequence, activity, class diagram generation

Sequence, activity and class diagrams have been used as an example to demonstrate the feasibility of generating an object-oriented representation of the biological process described by the GO term TGF-beta receptor complex assembly (GO:0007181), with Gene Ontology terms applied to generate these diagrams. Objects representing corresponding gene products are created, and their essential attributes are captured. Interactions among objects are also identified. For each interaction, a corresponding method is generated. This method is matched to a Gene Ontology term. The nature of the interaction determines the method parameters. The sequence of events is captured, and used to generate sequence diagrams. Scenarios are also generated for object interactions, and used to generate activity diagrams. The information captured in the sequence diagram and activity diagrams are used, along with the gene products attributes, to generate class diagrams.

References

  1. Gene Ontology Consortium[http://www.geneontology.org/]

  2. Lambrix P, Habbouche M, Perez M: Evaluation of ontology development tools for bioinformatics. Bioinformatics 2003, 19(12):1564–1571. 10.1093/bioinformatics/btg194

    Article  CAS  PubMed  Google Scholar 

  3. Gene Ontology Annotation[http://www.ebi.ac.uk/goa]

  4. Camon E, Magrane M, Barrell D, Lee V, Dimmer E, Maslen J, Binns D, Harte N, Lopez R, Apweiler R: The Gene Ontology Annotation (GOA) Database: sharing knowledge in Uniprot with Gene Ontology. Nucl Acids Res 2004, 32(90001):D262–266. 10.1093/nar/gkh021

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Dwight SS, Harris MA, Dolinski K, Ball CA, Binkley G, Christie KR, Fisk DG, Issel-Tarver L, Schroeder M, Sherlock G, Sethuraman A, Weng S, Botstein D, Cherry JM: Saccharomyces Genome Database (SGD) provides secondary gene annotation using the Gene Ontology (GO). Nucleic Acids Res 2002, 30(1):69–72. 10.1093/nar/30.1.69

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Harris MA, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, Eilbeck K, Lewis S, Marshall B, Mungall C, Richter J, Rubin GM, Blake JA, Bult C, Dolan M, Drabkin H, Eppig JT, Hill DP, Ni L, Ringwald M, Balakrishnan R, Cherry JM, Christie KR, Costanzo MC, Dwight SS, Engel S, Fisk DG, Hirschman JE, Hong EL, Nash RS, Sethuraman A, Theesfeld CL, Botstein D, Dolinski K, Feierbach B, Berardini T, Mundodi S, Rhee SY, Apweiler R, Barrell D, Camon E, Dimmer E, Lee V, Chisholm R, Gaudet P, Kibbe W, Kishore R, Schwarz EM, Sternberg P, Gwinn M, Hannick L, Wortman J, Berriman M, Wood V, de la Cruz N, Tonellato P, Jaiswal P, Seigfried T, White R: The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res 2004, (32 Database):D258–261.

    Google Scholar 

  7. Lu P, Szafron D, Greiner R, Wishart DS, Fyshe A, Pearcy B, Poulin B, Eisner R, Ngo D, Lamb N: PA-GOSUB: a searchable database of model organism protein sequences with their predicted Gene Ontology molecular function and subcellular localization. Nucleic Acids Res 2005, (33 Database):D147–153.

    Google Scholar 

  8. Adryan B, Schuh R: Gene-Ontology-based clustering of gene expression data. Bioinformatics 2004, 20(16):2851–2852. 10.1093/bioinformatics/bth289

    Article  CAS  PubMed  Google Scholar 

  9. Ahn WS, Kim KW, Bae SM, Yoon JH, Lee JM, Namkoong SE, Kim JH, Kim CK, Lee YJ, Kim YW: Targeted cellular process profiling approach for uterine leiomyoma using cDNA microarray, proteomics and gene ontology analysis. Int J Exp Pathol 2003, 84(6):267–279. 10.1111/j.0959-9673.2003.00362.x

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Arciero C, Somiari SB, Shriver CD, Brzeski H, Jordan R, Hu H, Ellsworth DL, Somiari RI: Functional relationship and gene ontology classification of breast cancer biomarkers. Int J Biol Markers 2003, 18(4):241–272.

    CAS  PubMed  Google Scholar 

  11. Badea L: Functional discrimination of gene expression patterns in terms of the gene ontology. Pac Symp Biocomput 2003, 565–576.

    Google Scholar 

  12. Chou KC, Cai YD: A new hybrid approach to predict subcellular localization of proteins by incorporating gene ontology. Biochem Biophys Res Commun 2003, 311(3):743–747. 10.1016/j.bbrc.2003.10.062

    Article  CAS  PubMed  Google Scholar 

  13. Deng M, Tu Z, Sun F, Chen T: Mapping Gene Ontology to proteins based on protein-protein interaction data. Bioinformatics 2004, 20(6):895–902. 10.1093/bioinformatics/btg500

    Article  CAS  PubMed  Google Scholar 

  14. Feng W, Wang G, Zeeberg BR, Guo K, Fojo AT, Kane DW, Reinhold WC, Lababidi S, Weinstein JN, Wang MD: Development of gene ontology tool for biological interpretation of genomic and proteomic data. AMIA Annu Symp Proc 2003, 839.

    Google Scholar 

  15. Jensen LJ, Gupta R, Staerfeldt HH, Brunak S: Prediction of human protein function according to Gene Ontology categories. Bioinformatics 2003, 19(5):635–642. 10.1093/bioinformatics/btg036

    Article  CAS  PubMed  Google Scholar 

  16. Lagreid A, Hvidsten TR, Midelfart H, Komorowski J, Sandvik AK: Predicting gene ontology biological process from temporal gene expression patterns. Genome Res 2003, 13(5):965–979. 10.1101/gr.1144503

    Article  PubMed  Google Scholar 

  17. Li S, Becich MJ, Gilbertson J: Microarray data mining using gene ontology. Medinfo 2004, 11: 778–782.

    Google Scholar 

  18. Lu X, Zhai C, Gopalakrishnan V, Buchanan BG: Automatic annotation of protein motif function with Gene Ontology terms. BMC Bioinformatics 2004, 5(1):122. 10.1186/1471-2105-5-122

    Article  PubMed Central  PubMed  Google Scholar 

  19. Masseroli M, Martucci D, Pinciroli F: Towards biological knowledge mining by statistical analysis of gene ontology annotations. Medinfo 2004, 2004(CD):1745.

    Google Scholar 

  20. Pinto FR, Cowart LA, Hannun YA, Rohrer B, Almeida JS: Local correlation of expression profiles with gene annotations – proof of concept for a general conciliatory method. Bioinformatics 2005, 21: 1037–1045. 10.1093/bioinformatics/bti074

    Article  CAS  PubMed  Google Scholar 

  21. Schug J, Diskin S, Mazzarelli J, Brunk BP, Stoeckert CJ Jr: Predicting gene ontology functions from ProDom and CDD protein domains. Genome Res 2002, 12(4):648–655. 10.1101/gr.222902

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Vinayagam A, Konig R, Moormann J, Schubert F, Eils R, Glatting KH, Suhai S: Applying Support Vector Machines for Gene Ontology based gene function prediction. BMC Bioinformatics 2004, 5(1):116. 10.1186/1471-2105-5-116

    Article  PubMed Central  PubMed  Google Scholar 

  23. Gene Ontology Tools[http://www.geneontology.org/GO.tools.shtml]

  24. Ashburner M, Mungall CJ, Lewis SE: Ontologies for biologists: a community model for the annotation of genomic data. Cold Spring Harb Symp Quant Biol 2003, 68: 227–235. 10.1101/sqb.2003.68.227

    Article  CAS  PubMed  Google Scholar 

  25. Zhang S, Bodenreider O: Comparing Associative Relationships among Equivalent Concepts Across Ontologies. Medinfo 2004, 11: 459–466.

    Google Scholar 

  26. Smith B, Williams J, Schulze-Kremer S: The ontology of the gene ontology. AMIA Annu Symp Proc 2003, 609–613.

    Google Scholar 

  27. Ogren PV, Cohen KB, Acquaah-Mensah GK, Eberlein J, Hunter L: The compositional structure of Gene Ontology terms. Pac Symp Biocomput 2004, 214–225.

    Google Scholar 

  28. Smith B, Kumar A: Controlled vocabularies in bioinformatics: a case study in the gene ontology. DDT: BIOSILICO 2004, 2(6):246–252. 10.1016/S1741-8364(04)02424-2

    CAS  Google Scholar 

  29. GO-DEV[http://www.godatabase.org/dev/index.html]]

  30. Taylor CF, Paton NW, Garwood KL, Kirby PD, Stead DA, Yin Z, Deutsch EW, Selway L, Walker J, Riba-Garcia I, Mohammed S, Deery MJ, Howard JA, Dunkley T, Aebersold R, Kell DB, Lilley KS, Roepstorff P, Yates JR 3rd, Brass A, Brown AJ, Cash P, Gaskell SJ, Hubbard SJ, Oliver SG: A systematic approach to modeling, capturing, and disseminating proteomics experimental data. Nat Biotechnol 2003, 21(3):247–254. 10.1038/nbt0303-247

    Article  CAS  PubMed  Google Scholar 

  31. Spellman PT, Miller M, Stewart J, Troup C, Sarkans U, Chervitz S, Bernhart D, Sherlock G, Ball C, Lepage M, Swiatek M, Marks WL, Goncalves J, Markel S, Iordan D, Shojatalab M, Pizarro A, White J, Hubley R, Deutsch E, Senger M, Aronow BJ, Robinson A, Bassett D, Stoeckert CJ Jr, Brazma A: Design and implementation of microarray gene expression markup language (MAGE-ML). Genome Biol 2002, 3(9):RESEARCH0046. 10.1186/gb-2002-3-9-research0046

    Article  PubMed Central  PubMed  Google Scholar 

  32. Shegogue D, Zheng WJ: Object-oriented biological system integration: a SARS coronavirus example. Bioinformatics 2005, 21: 2502–9. 10.1093/bioinformatics/bti344

    Article  CAS  PubMed  Google Scholar 

  33. Rodriguez C, Chen F, Weinberg RA, Lodish HF: Cooperative binding of transforming growth factor (TGF)-beta 2 to the types I and II TGF-beta receptors. J Biol Chem 1995, 270(27):15919–15922. 10.1074/jbc.270.27.15919

    Article  CAS  PubMed  Google Scholar 

  34. Brown CB, Boyer AS, Runyan RB, Barnett JV: Requirement of type III TGF-beta receptor for endocardial cell transformation in the heart. Science 1999, 283(5410):2080–2082. 10.1126/science.283.5410.2080

    Article  CAS  PubMed  Google Scholar 

  35. Massague J: TGF-beta signal transduction. Annu Rev Biochem 1998, 67: 753–791. 10.1146/annurev.biochem.67.1.753

    Article  CAS  PubMed  Google Scholar 

  36. Yamashita H, ten Dijke P, Franzen P, Miyazono K, Heldin CH: Formation of hetero-oligomeric complexes of type I and type II receptors for transforming growth factor-beta. J Biol Chem 1994, 269(31):20172–20178.

    CAS  PubMed  Google Scholar 

  37. Tsukazaki T, Chiang TA, Davison AF, Attisano L, Wrana JL: SARA, a FYVE domain protein that recruits Smad2 to the TGFbeta receptor. Cell 1998, 95(6):779–791. 10.1016/S0092-8674(00)81701-8

    Article  CAS  PubMed  Google Scholar 

  38. Xu L, Chen YG, Massague J: The nuclear import function of Smad2 is masked by SARA and unmasked by TGFbeta-dependent phosphorylation. Nat Cell Biol 2000, 2(8):559–562. 10.1038/35019649

    Article  CAS  PubMed  Google Scholar 

  39. Inman GJ, Hill CS: Stoichiometry of active smad-transcription factor complexes on DNA. J Biol Chem 2002, 277(52):51008–51016. 10.1074/jbc.M208532200

    Article  CAS  PubMed  Google Scholar 

  40. Dennler S, Itoh S, Vivien D, ten Dijke P, Huet S, Gauthier JM: Direct binding of Smad3 and Smad4 to critical TGF beta-inducible elements in the promoter of human plasminogen activator inhibitor-type 1 gene. Embo J 1998, 17(11):3091–3100. 10.1093/emboj/17.11.3091

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Yingling JM, Datto MB, Wong C, Frederick JP, Liberati NT, Wang XF: Tumor suppressor Smad4 is a transforming growth factor beta-inducible DNA binding protein. Mol Cell Biol 1997, 17(12):7019–7028.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Zawel L, Dai JL, Buckhaults P, Zhou S, Kinzler KW, Vogelstein B, Kern SE: Human Smad3 and Smad4 are sequence-specific transcription activators. Mol Cell 1998, 1(4):611–617. 10.1016/S1097-2765(00)80061-1

    Article  CAS  PubMed  Google Scholar 

  43. Xu L, Kang Y, Col S, Massague J: Smad2 nucleocytoplasmic shuttling by nucleoporins CAN/Nup214 and Nup153 feeds TGFbeta signaling complexes in the cytoplasm and nucleus. Mol Cell 2002, 10(2):271–282. 10.1016/S1097-2765(02)00586-5

    Article  CAS  PubMed  Google Scholar 

  44. Inman GJ, Nicolas FJ, Hill CS: Nucleocytoplasmic shuttling of Smads 2, 3, and 4 permits sensing of TGF-beta receptor activity. Mol Cell 2002, 10(2):283–294. 10.1016/S1097-2765(02)00585-3

    Article  CAS  PubMed  Google Scholar 

  45. Lo RS, Massague J: Ubiquitin-dependent degradation of TGF-beta-activated smad2. Nat Cell Biol 1999, 1(8):472–478. 10.1038/70258

    Article  CAS  PubMed  Google Scholar 

  46. Papin JA, Reed JL, Palsson BO: Hierarchical thinking in network biology: the unbiased modularization of biochemical networks. Trends Biochem Sci 2004, 29(12):641–647. 10.1016/j.tibs.2004.10.001

    Article  CAS  PubMed  Google Scholar 

  47. Bolouri H, Davidson EH: Modeling transcriptional regulatory networks. Bioessays 2002, 24(12):1118–1129. 10.1002/bies.10189

    Article  CAS  PubMed  Google Scholar 

  48. Wolf DM, Arkin AP: Motifs, modules and games in bacteria. Curr Opin Microbiol 2003, 6(2):125–134. 10.1016/S1369-5274(03)00033-X

    Article  CAS  PubMed  Google Scholar 

  49. Hucka M, Finney A, Sauro HM, Bolouri H, Doyle JC, Kitano H, Arkin AP, Bornstein BJ, Bray D, Cornish-Bowden A, Cuellar AA, Dronov S, Gilles ED, Ginkel M, Gor V, Goryanin II, Hedley WJ, Hodgman TC, Hofmeyr JH, Hunter PJ, Juty NS, Kasberger JL, Kremling A, Kummer U, Le Novere N, Loew LM, Lucio D, Mendes P, Minch E, Mjolsness ED, Nakayama Y, Nelson MR, Nielsen PF, Sakurada T, Schaff JC, Shapiro BE, Shimizu TS, Spence HD, Stelling J, Takahashi K, Tomita M, Wagner J, Wang J: The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models. Bioinformatics 2003, 19(4):524–531. 10.1093/bioinformatics/btg015

    Article  CAS  PubMed  Google Scholar 

  50. Finney A, Hucka M: Systems biology markup language: Level 2 and beyond. Biochem Soc Trans 2003, 31(Pt 6):1472–1473.

    Article  CAS  PubMed  Google Scholar 

  51. AmiGO[http://godatabase.org/]

  52. Shegogue D, Zheng WJ: Capturing biological information with class-responsibility-collaboration cards. Bioinformatics 2005, 21: 415. 10.1093/bioinformatics/bti005

    Article  CAS  PubMed  Google Scholar 

  53. Graham I: Basic Concepts. In Object-oriented Methods, Principles & Practice. Third edition. Harlow, England: Addison-Wesley; 2001:1–37.

    Google Scholar 

Download references

Acknowledgements

Daniel Shegogue is supported by NLM training grant 5-T15-LM007438-02. W. Jim Zheng is partly supported by a grant (DE-FG02-01ER63121) from the Department of Energy.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to W Jim Zheng.

Additional information

Authors' contributions

DS drafted the manuscript, constructed the models and participated in the design of the study. WJZ was the principal investigator, conceived of the project and guided its development. All authors read and approved the final manuscript.

Electronic supplementary material

12859_2004_438_MOESM1_ESM.pdf

Additional File 1: Class-responsibility-collaboration card for TGF-beta 1. Attributes, collaborators and responsibilities of the specified protein are given. The attributes section allows the ordered listing of information not easily captured by the UML notation. The collaborator section lists the cellular components that interact with TGF-beta 1. The responsibilities section specifies the consequence of TGF-beta 1 interacting with its collaborator. This card allows TGF-beta 1 to be decomposed into an object containing attributes, operations and interactions. (PDF 78 KB)

12859_2004_438_MOESM2_ESM.pdf

Additional File 2: Class-responsibility-collaboration card for TGF-beta receptor II. Attributes, collaborators and responsibilities of the specified protein are given. The attributes section allows the ordered listing of information not easily captured by the UML notation. The collaborator section lists the cellular components that interact with TGF-beta receptor II. The responsibilities section specifies the consequence of TGF-beta receptor II interacting with its collaborator. This card allows TGF-beta receptor II to be decomposed into an object containing attributes, operations and interactions. (PDF 79 KB)

12859_2004_438_MOESM3_ESM.pdf

Additional File 3: Class-responsibility-collaboration card for TGF-beta receptor I. Attributes, collaborators and responsibilities of the specified protein are given. The attributes section allows the ordered listing of information not easily captured by the UML notation. The collaborator section lists the cellular components that interact with TGF-beta receptor I. The responsibilities section specifies the consequence of TGF-beta receptor I interacting with its collaborator. This card allows TGF-beta receptor I to be decomposed into an object containing attributes, operations and interactions. (PDF 80 KB)

12859_2004_438_MOESM4_ESM.pdf

Additional File 4: Class-responsibility-collaboration card for SMAD2. Attributes, collaborators and responsibilities of the specified protein are given. The attributes section allows the ordered listing of information not easily captured by the UML notation. The collaborator section lists the cellular components that interact with SMAD2. The responsibilities section specifies the consequence of SMAD2 interacting with its collaborator. This card allows SMAD2 to be decomposed into an object containing attributes, operations and interactions. (PDF 87 KB)

12859_2004_438_MOESM5_ESM.pdf

Additional File 5: Use case describing the events leading to TGF-beta receptor complex assembly. This use case defines the boundaries of the system model. Here, the main success and alternative scenarios leading to the assembly of the TGF-beta receptor complex are described. (PDF 66 KB)

12859_2004_438_MOESM6_ESM.pdf

Additional File 6: Conceptual diagram of the TGF-beta receptor complex, and the proteins that associate with this complex. Gene products have been decomposed into objects. Object attributes and operations are hidden to reduce complexity. Gene products that comprise the receptor complex are shown in blue with their associated relationships in green. As the boundaries of the use case do not include them, other associated proteins that comprise the TGF-beta signaling pathway are not emphasized in further diagrams. (PDF 122 KB)

12859_2004_438_MOESM7_ESM.pdf

Additional File 7: State diagrams for the TGF-beta receptor complex components. State diagrams describe the possible states in which a bioentity can exist. Concurrent states are separated by vertical dashed lines. A) State diagram for TGF-beta. B) State diagram for TGF-beta receptor II. C) State diagram for TGF-beta receptor I. D) State diagram for SMAD2. (PDF 80 KB)

12859_2004_438_MOESM8_ESM.pdf

Additional File 8: An example of a class diagram with expanded function names showing the interactions between the components of the TGF-beta receptor complex (GO:0007181) (grayed). Cellular components containing these gene products are also shown. GO function terms are shown with their corresponding GO ids. This format demonstrates a more user-friendly interface for reading the diagrams, whereas figure 4 is more suitable as a computer readable format. (PDF 194 KB)

Authors’ original submitted files for images

Rights and permissions

Reprints and permissions

About this article

Cite this article

Shegogue, D., Zheng, W.J. Integration of the Gene Ontology into an object-oriented architecture. BMC Bioinformatics 6, 113 (2005). https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2105-6-113

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2105-6-113

Keywords