key: cord-0164402-6x56ys4w authors: Chirom, Keilash; Malik, Md. Zubbair; Somvanshi, Pallavi; Singh, R. K. Brojen title: Network medicine in ovarian cancer: Topological properties to drug discovery date: 2021-08-15 journal: nan DOI: nan sha: a0c6fed216462a1f16f02566a93f2f37da4e4942 doc_id: 164402 cord_uid: 6x56ys4w The investigation of topological properties of ovarian cancer network (OCN) and the roles of hubs involved in it by digging the network at various levels of organization are important to understand how OCN is organized to understand disease states. The OCN constructed from the experimentally verified genes exhibits fractal nature in the topological properties of the network and deeply rooted communities. Also, the network properties at all levels of organization obey one parameter scaling law which lacks centrality lethality rule. We then showed that $langle krangle$ can be taken as a scaling parameter, where, power law exponent can be estimated from the ratio of network diameters, $lambda=1+rholeft[frac{d_c}{d_k}right]$. The betweenness centrality shows two distinct behaviors one shown by high degree hubs and the other by segregated low degree nodes. The $C_B$ distribution follows power law behavior with the exponent connected to exponents of distributions of high and low degree nodes by, $thetasimleft[1+frac{1}{2}left(frac{1}{epsilon_s}+frac{1}{epsilon_h}right)(lambda-1)right]$. Absence of rich-club formation leads to the missing of a number of attractors in the network causing formation of weakly tied diverse functional modules to keep optimal network efficiency. The hubs knockout experiment shows that provincial hubs take major responsibility to keep network integrity and organization. The identified key regulators are found to be provincial and connector hubs. Further, two key regulators, EPCAM and CD44 are found to be maximally over expressed at various cancer stages (II-IV). They are also positively correlated with immune infiltrates (CD4+ T cells). Finally, few potential drugs are identified related to the key regulators. Cancer is generally believed to be genetic disease [1] and interferes a number of interacting spatio-temporal signaling gene regulatory networks (SGRN), which are responsible for molecular basis communications, cellular and other processes in the whole complex cellular networks [2] . Since the SGRN are composed of basic motif circuits [3] , the overall signal processing by them is the integration of the signals by these circuits [4] which may reflect fractal or scale-free properties in the network organization with hierarchical features [5] given by self-affine process [7] of any property F (x) of a scale factor α: F (αx) F (x) = Γ(α) ∼ α D , such that, Γ(α 1 α 2 ...α N ) = Γ 1 (α 1 )Γ 2 (α 2 )...Γ N (α N ) : Γ k (α k ) ∼ α D k k with 0 ≤ α, α 1 , α 2 , ..., α N ≤ 1, where, D, D k are self-similarity dimension. Hence, the static and dynamical changes in the SGRN due to changes in signal processing when the cellular system move from normal to disease state are reflected in the topological properties of the network [8] . Because the changes in the signal patterns in the network is driven by genetic and epigenetic perturbations which alter significant rewiring of the these signaling networks [10, 11] . This perturbation may enhance or diminish the functional activity of a gene or group of genes causing significant changes in the signaling patterns leading to variations in the cellular phenotype [12] which may exhibit in the topological properties of the network [13] . On the other hand, Kauffman's cancer attractor theory [14, 15] developed on the basis of his works and earlier works of Max Delbruck [16] , and Jacob and Monod [17] , is another important theory to explain the disruption in signaling mechanism at cancer state [15] . The theory claims that the interacting genes in each SGRN spatio-temporally regulate among themselves to reach an equilibrium state or attractor which corresponds to balancing of their gene expression patterns [14] . Hence, the state of a complex network having N interacting genes with their expression levels x 1 , x 2 , ..., x N at any instant of time ′ t ′ can be defined as an N -dimensional vector Λ(t) = [x 1 , x 2 , ..., x N ] −1 in the state space [15] . In this state space, there could be a large number of attractors, where, each attractor corresponds to an equilibrium state S * surrounded by non-stationary unstable states [14, 18] and encodes specific genetic activities [15] . The attractors could be of two types, one may correspond to stable network corresponding to normal mature cell type, and the other may correspond to cancer attractor separated by epigenetic barrier which can be lowered by mutation load. Hence, in the state corresponding to cancer attractor, the property of the signal processing mechanism in the network will be quite different and may be measured in the topological properties of the network. It has been proposed that signaling entropy [19] , which can be calculated from local signaling entropy of the network [13] , could be one potential parameter to characterize the variations in the signaling patterns in the network. However, capturing these altered signaling patterns with network topological properties and their quantification are still debatable and not fully understood yet. Ovarian cancer is known as silent killer because early diagnosis of this cancer type is quite difficult and generally more than 70% cases are not able to be detected until at significantly advanced stage [20, 21] . One reason could be inability to recognize the symptoms of this cancer because the data indicates that the origin of this cancer is generally from fallopian tube and not from ovaries [22] . Hence, this cancer type is one of the most inheritable cancers, the fifth cause of cancer in women and the fourth cause of women death in developed countries [23, 24] . Further, since genetic instability is the hallmark of cancer disease that causes genetic heterogeneity triggering neoplastic progression [25] , it becomes quite difficult to which key gene/genes would be targeted for possible cure of cancer diseases [26] . However, identification of key genes of ovarian cancer should be done systematically from various sensitive experimental and theoretical perspectives because those identified key genes are the hopes of possible drug targets in systematic drug design that may lead to cure and prevent from this dreaded disease. There have been few techniques to identify the key regulators, namely, network theory which was used to identify the five novel key genes (AKT1, KRAS, EPCAM, CD44 and MCAM) in ovarian cancer network [29] , new biomarker (CREB1 gene) that was proposed by analyzing various databases with experiments [30] , and GWAS (genome-wide association studies) based identification of biomarkers specially for African women [31] etc. In a particular network, hubs (high degree nodes), specially key regulators, are the main agents which are involved in many signaling processes, in maintaining topological properties, stability and self-organization of the network [5, 32] . For instance, it has been found that a group of important hub genes involved in prostate cancer are significantly reduced their gene expressions and hence quite disturb in their regulatory mechanisms and signal processing [33] , regulatory mechanism of a group of genes in few pathways (eg glucose metabolism) in Lung cancer are found to be disturbed [34] , significant dysregulation of signaling hubs of lipid-regulated genes was found in Hepatocellular carcinoma [35] etc. If the cancer cause dysregulation or disturbance in the signaling hubs in the complex gene regulatory network, then there should be significant changes in the topological properties of the network. The question is how to capture the changes in the information patterns in the genes and associated networks. In the complex protein-protein interaction (PPI), one possible technique to estimate the information contained in the network is to measure entropy of the network from which information contained in the network can be estimated [36, 37] . The entropy of a network G(V, N ) with N possible states and adjacency matrix [U ] N ×N , which is the amount of uncertainty one can observe in accessing the states of the network [37] , can be measured by, [38] . The maximum value of H of U , which have N 2 elements {x i ; i = 1, 2, ..., N 2 } and x i = 0, 1(binary), is given by, H max = log 2 [N 2 ] from which one can estimate the information contain in the network G by, I = H max − H [40] . There are various ways to calculate entropy of a network, namely, graph entropy [41] : define probability of a node n to have inner links e as p i (n) = e |N (n)| and probability of this node n to have outer links by p o (n) = 1 − p i (n), then H(G) = n∈V h(n), h(n) = −p i (n)lnp i (n) − p o (n)lnp o (n), such that, I = H max − H(G). Next method is to calculate signaling entropy [13] for PPI network with each node represented by gene expression data, where, the Pearson correlation value between any pair of nodes (i, j) C ij can be related to weight w ij assigned to the edge by w ij = 1 2 [1 + C ij ] from which one can obtain probability distribution p ij = w ij j∈N (i) w ij to calculate entropy and On the other hand, simple network centrality, namely, betweenness centrality (BC) metric can be used as a significantly sensitive parameter to measure the capability of a node to spread information in a complex network [60] . The reason could be due to the fact that nodes having high centrality values can reach the other nodes in the network on short paths [61] enhancing the speed of information processing, and vice versa. This BC can then be correlated with the connectivity (degree) k in the network as a power law C B (k) ∼ k ǫ [62] . From this relation, one can estimate maximum BC as C max B ∼ k ǫ γ−1 max , where, maximum connectivity is a function of network size N by, k max ∼ N network follows scale free then hubs try to control and manipulate the whole network such that their removal cause network breakdown [102] , whereas, if the network obeys hierarchical features, hubs are busy in correlating nodes in its own module/community as well as with the communities such that their removal should not cause network breakdown [5, 119] . These hubs are heterogeneous in nature and can be classified into various types depending on their participation in the network and roles [123] . Hence, for a particular disease, it is quite important to categorize the identified key/driver regulating genes/proteins to understand their roles and potential in controlling the disease. Then, the key regulators should be studied systematically in terms of their binding affinity to various available drugs and medicinal compounds for possible drug discovery. In this work, we address structural behavior of the OC network to understand the information processing and organization of the network along with the nature of the hubs with their classification using network theoretical approach. The behavior of the constituent communities at various levels of organization in order to understand how OC network works at fundamental level. Then the changes in the gene expression levels of the key regulators of the OC and at their stages of the cancer. Then, the key regulators are allowed to interact with the FDA approved available drugs to identify drug which could be used to prevent/cure OC disease. We now present our analysis of the ovarian cancer (OC) data we have mined for study. Our main focus is three fold, first, we calculated topological properties of the complex network we have constructed using experimentally verified list of genes/proteins (Methods) to capture the changes in the calculated topological properties at cancer state and their roles in the network organization. Second, the roles of the identified novel key regulators in the cancer state, specially, in signal transduction from information theoretical approach. Third, search for novel natural compounds (drug (FDA approved) or food nutrients) with respect to the identified key regulators. Theory of fractal properties of the ovarian cancer network We considered experimentally verified pre-processed 4818 genes with 16,320 interactions collected from the patients' data mined from the ovarian cancer dedicated and related database (see the workflow and Methods for detail information of mining and pre-processing the data). The complex OC network is analyzed using Grivan and Newman community finding method [42] , and found that the network has five levels of hierarchical organization of communities/sub-communities till triangular motif level (Fig. 2) , which are the basic functional unit of network organization [3] and quite abundant in real networks [43, 44] , and follows interesting systems level hierarchical organization in the network [29] as also supported by the topological results of the network and its communities. Let us represent the hierarchical scale free OC network by a graph of size N G(V, E), where, V = {n i ; i = 1, 2, ..., N, n i ∈ V } is the set of nodes, and E = {e ij ; i, j = 1, 2, ..., N, e ij ∈ E} is the set of edges. The topological properties of the network, namely, probability of degree distribution P (k), clustering coefficient C(k), connectivity C N (k), centrality measurement parameters, betweenness C B (k), closeness C C (k) and eigenvector C E (k) centralities follow power law properties, where, the power law exponents have their own values. We also found that the data of each topological distribution function of the primary network and the constituent communities collapse into a single curve after interpolation of the data following the MacKinnon and Kramer's numerical procedure which indicates that each topological distribution function (U (k) → P/C/C N /C B /C C /C E ) from primary network to communities/sub-communities down to motif level follow one parameter scaling law [45] given by, Here, η was calculated by fitting with the scaled data and found to be approximately proportional to the minimum path length of the corresponding network/community/sub-community, and we found that H → constant = b. This shows that U (k) ∼ bk λ , such that F (k) = U b ∼ k λ . If we take b as the normalization constant such that b = where, N j is the size of the jth community/sub-community with j = 0 corresponds to primary network N 0 → N , then, 0 ≤ F (k) ≤ 1 and F behaves just like probability distribution function, such that, λ = {−γ, −α, −β, ǫ, κ, δ} corresponding to P (k), C(k), C N (k), C B (k), C C (k), C E (k) respectively. Now ifT is an operator which coarse grain the system with a scale factor ω, then it is clearly seen that:T F (x) = F (ωx) = (ωx) λ = ω λ F (x) [7] . If one consider the property F (k) of the primary OC network and the same scaled property F ′ (k ′ ) of any one of the communities/sub-communities of the network, then one parameter scaling law (2) provide the renormalization relation F (k)dk = F ′ (k ′ )dk ′ : k → k ′ : k = ωk ′ which gives us the following renormalization relation, This means that the topological properties of the OC network remains the same if the network is scale up or down by a scale factor ω which carries the fractal properties in the network organization [7] indicating self-organizing behavior in the OC network organization [18] . Theorem 1: The communities which organize complex hierarchical OC network exhibit self-similar or fractal topological properties. Proof: Since the OC network G(V, N ) is hierarchical scale free topology, the network can be decomposed into M communities g i (V i , E i ); i = 1, 2, ..., M ; ∀V i ∈ V, ∀E i ∈ E with random size and structures (one can use various community finding methods, namely, Girvan and Newman method [42] , Louvain method [47] , constant potts model method [48] and other related methods [49] ). The topological properties of these communities F i follow one parameter scaling behavior (2) with the primary OC network and obeys renormalized scaling law (3) . Generally, these communities can be dependent on different scaling parameters ω i , such that, Then, one can write: have the following scaling law, This scaling law indicates that the topological properties of the communities of OC network have fractal properties with the primary network as well as among the communities also following the scaling law in ′ k ′ given by equation (4) . The scaling law still holds true for the organization of smaller communities at various levels of organization till motif level. Numerical proof of the theorem can be done by the analysis of the topological properties OC network from primary network down to the motif level (Fig. 2) . The behavior of P (k) as the function of k for different networks at various levels of organization (primary, communities, sub-communities down to motif level as shown in Fig. 2 b) are approximately similar as shown by the nearly parallel fitted lines to the data points, where, the calculated values of γs for all networks are approximately the same (Fig. 2 a. ). This clearly indicates that all the P (k) data are collapsed to a single curve obeying the one parameter scaling behavior given by equation (2) . Similar one parameter scaling nature can be found in the behaviors of C(k) and C N (k) (Fig. 2 a middle and lower panels). The one parameter scaling property can also be observed in the behaviors of centrality measurement parameters (C B (k), C C (k) and C E (k)) where we discussed the behaviors of high degree nodes in the respective networks ( Fig. 2 c all panels). The reason is that signal processing in each respective network is mainly regulated by high degree hubs [5, 32] . However, there is significant roles of the low degree nodes in regulating the network behavior which we will be discussing also in this work. Then the data of the constituent networks fall in a single curve obeying one parameter scaling law given by (2) . Lastly, we calculated network connectance or density [50] given by, Ω = 2EN N (N −1) ∼k N , where, E N is the total number of edges in the network of size N andk is the mean degree of the network. If the network is dense then lim N →∞ Ω → constant, otherwise, the network is sparse, where, lim N →∞ Ω → 0. We calculated Ω for all communities at various levels of organization (Fig. 2 b) and found that all the data points collapsed to a single curve which obey power law behavior, Ω(N ) ∼ N −1 with good goodness of fit. The behavior can be seen in two ways, first, Ω rapidly decreases as N increases till N ∼ 100 then becomes nearly steady state which approximately converges to small value, but does not go to zero. This indicates that the complex OC network is weakly compact as N increases. Second, since the data points beautifully falls in a single curve which obey power law behavior, the wiring/rewiring of the edges among the nodes of the communities/sub-communities is done proportionately to keep the self-similarity in the topological properties of the network. Theorem 2:Hierarchical scale free OC network does not obey centrality lethality rule. Proof: Communities at a certain level of organization are loosely connected to minimize wiring energy cost to maximize information processing among them [51] , but within each community nodes are connected densely with large within-community clustering coefficients and significantly small path lengths [52] to allow each individual community function independently [53] . Hence, we assume that the topological properties of the communities are nearly independent of each other and we consider the network G(V, E) has M communities at a certain level of organization, such that, Now, if we consider degree k as random variable which is integrable in the set G(V, N ), and since the the topological , one can assume that this function behaves like probability distribution function. Taking distribution function of equation (5) and applying theorem of addition of probability [54] , we have, Now, if most important hub/hubs in the OC network is is/are removed, the most probable changes the perturbation can make in the network is that the community/communities g c+1 , g c+2 , ..., g c+n in which the nodes present it/them will be affected the most causing rewiring of the nodes in the network [49] , and in the worst scenario the community or some of the communities may collapse loosing all corresponding topological properties [55] . Hence, even though the affected communities become breakdown due to the perturbation, then, corresponding topological properties will be vanished: c+n j=c+1 F (k j ) → 0, the other remaining (M − n) communities will be survived because of their tight interaction keeping themselves small worldness property [52, 53] , and they will reorganize the OC network with the scaling transformation of the random variable k by k → k ′ : k ′ = θk and k i by The above equation indicates that removing hub/hubs below a critical number of hubs does not cause the collapse of the OC network even though few communities are broken down due to the perturbation, but the network reorganizes with modified topological properties. Hence, emergence of tightly bound nearly independent communities/subcommunities are the main reason in keeping complex network integration, stability, adaptation and management of external and internal perturbations and most importantly self-organization of the network, where, the activities of the leading hubs are well regulated by the emergent functional communities/sub-communities. Proof from the OC network data analysis is as follows. We carry out hubs knock-out experiment by systematically removing the first 10-400 leading hubs from the primary OC network, then studied the changes in the topological properties of the resulting networks and compared with those of primary OC network ( Fig. 3 and Fig. 4 ). We found the following observations. The systematic removal of the hubs from the primary OC network upto certain critical hub number N T C , the fractal properties of the resulting networks strictly kept preserved with change in the power of the topological distribution function: However, if the number of hubs removed is larger than this critical value, N C N T C then the k ′ of the resulting networks become random such that the F become randomly varied. For example, the exponent of power law behavior of P (k) become γ 3 and γ show random variation ( Fig. 3 left panel and inset figure) . For the case N C N T C and depending on the values γ one can classify as : (a) hierarchical nature for γ 2ln2; (b) scale free nature for the regime 2ln2 γ 3, and (c) random for γ 3. For clustering coefficient after N C 80 it is difficult to see the behavior because the calculated clustering coefficient values for those resulting networks from which large number of hubs are generally small and nearly overlapped and fitting to those data becomes quite meaningless ( Fig. 3 middle panel) . The scenario is distinctly visible in the behavior of neighborhood connectivity, where, the power law exponent β → 0 for N C N T C otherwise finite even though there is quite fluctuations in the values (Fig. 3 panel 4 and inset). In the case of centrality measurement functions also the exponent of the power law distributions are either goes to zero (ǫ, κ → 0) or become random (for δ) for N C N T C ( Fig. 4 and insets). The goodness of the fit to each data measured by r 2 value is also quite reasonably good (0.97 r 2 0.7). From these results we can conclude the following behavior of F , The scale free network which is obtained here is categorized based on the value of γ (2 γ3), however, the rewired network still have strongly modular structure in the network, where, C(k) and C N (k) still exhibit power law behaviors. This modular structure becomes vanishing as γ → 3, and then as γ 3 the topological properties become randomly distributed. Proposition: Self-similarity process in OC network is given by the following scaling law: where, d c and d k are network diameters in cth community, and ρ is proportionality constant. Proof: The average diameter, which is the shortest distance between two nodes in a network, can be taken as network metric and has the form in hierarchical network [8, 56, 57] where, N is the network size from which scaling law of cth community of the network is given by, where, D c is the fractal dimension of cth community which can be related to the diameter of the cth community of the network. l c and l 0 are shortest length between two nodes in cth community and characteristics length respectively. Here, k c is taken as the scaling parameter to understand scaling law in complex OC network. Now network mass of the cth community can be obtained as [9] , where, scaling of k to k c with scaling factor ω[ k c ] is given by renormalization theory given below [9] , k → k c : Renormalization theory allows to express the distribution function F (k) in terms of the number of nodes having 'k' links n(k) in primary network with the number of nodes having k c in cth community Putting the power law distribution of F (k) ∼ k −λ and scaling law in ω and equation (13) we have, The fractal dimensions D c and D k given in equation (12) can be expressed as: Now, by using scaling in ω and equation (14) , one can reach the following relation, The equations (12) and (15), we could relate the fractal dimension to the diameter of the network and mean degree k , which is dependent on various network topological properties, can be taken as a scaling parameter to understand fundamental laws of network properties and their organization. Further, the distribution function F can be obtained as follows, The above equation indicates that if one knows the mean value of degree of a network and network diameter one can able to know the fractal properties of a network. Significant role of lower degree hubs at cancer state in OC network One way to understand the role of hubs in regulating a complex network and its constituent communities/subcommunities is to study the network behavior as a function of degree (k) because k relates directly with various important topological distribution functions [5, 32, 50] . One of the most important functions of the hubs involved in a complex network is the capability of triggering information flow among the nodes directly or importantly influencing the information flow among the nodes which are not connected directly, and can be characterized by behavior of centrality measures [58] . The behavior of betweenness centrality C B as a function of k is two-fold in OC network ( Fig. 4 left panel) . First, the primary network and the networks in which significant number of hubs are not removed (number of highest degree hubs removed 150) C B follows two behaviors, (a) low degree nodes follow power law with larger exponent ǫ: C Fig. 4 left panel inset). This indicates that lower degree nodes can able to communicate others in efficient manner with shorter paths which may be due to the possibility of these small degree nodes may be segregated with smaller paths [59, 60] . On the other hand, high degree nodes have longer paths as compared to small degree nodes, which could be due to sparsely distribution of high degree nodes buried inside the network and proper link needed to have by each of them might be reduced significantly at cancer state. Now, one can have two sets of nodes of the whole network, It has been shown that for large complex distribution of betweenness centrality follows power law: , hence, the joint probability distribution of betweenness centrality for the whole network can be written as, This equation (17) shows the connection among the power law exponents θ, ǫ s , ǫ h and γ. Now, as the number of high degree hubs are removed ǫ s → ǫ h , and at a certain critical number of hubs removed N c , we found that ǫ decreases such that lim Nc→large ǫ → 0 which indicates that the path length becomes constant indicating information flow becomes constant. Similar behavior is found in the behavior of closeness centrality (CC) which is the measure of distance between a node to all other nodes in the network segregated with close shortest paths from it [64] . The nature of closeness centrality of high degree nodes as a function of degree exhibits power law C [h] C (k) ∼ k δ h , but the exponent κ h has comparatively smaller value as compared to that of lower degree nodes, C [s] C (k) ∼ k δs , δ h δ s (Fig. 4 right panel) . Since CC indicates the measure of geographical distance to spread the information from a node in the network [60, 64] , the rate of information transfer coordinated by high degree hubs is significantly reduced as compared to those of low degree nodes at OC state. One of the key reasons could be the direct drastic changes in the expression of some of the driver genes (high degree hubs) in the genetic network at cancer state, where, rest of the existing nodes (mainly low degree nodes) in the network are indirectly influenced [29, 65] . However, as the low degree nodes generally tie with high degree hubs [66] , the low degree nodes might know the changes in the information processing of the high degree hubs in the network and probably might act accordingly to keep efficiency of information processing in the network. Further, as the number of high degree hubs removed is increased δ s → δ h = δ and also lim Nc→large δ → 0 indicating information spread becomes constant. The nature of eigenvector centrality (EC) is also quite similar to that of betweenness centrality and obeys power law behavior, where, the exponent of the power law function of the low degree nodes of the OC network, C [s] E (k) ∼ k κs , is larger than that of high degree nodes, C [h] E (k) ∼ k κ h , i.e. κ h κ s (Fig. 4 middle panel and inset) . The importance of EC is that this network metric not only reveals the influence of a node in the neighboring nodes but also the score of the connected nodes [67] . Hence, the significant decrease in the value of κ h of the high degree hubs indicates that the signal propagation and influence of the hubs at various layers of neighbors of nodes starting from each one of them connecting till the end of the network [68] is quite reduced at cancer state. Whereas, the low degree nodes, which received change in influence of the OC network, maintain faster rate where κ s is large as compared to κ s . As seen before in other centrality measures we have: κ s → κ h = κ as lim Nc→large κ → 0, such that, C E (k) → constant. at large N c limit. The behavior of neighborhood connectivity C N (k) (Fig. 3 right panel) obeyed power law nature with negative exponent which is the signature of disassortativity nature [69, 70] . This disassortativity nature, which is again supported by the rich-club coefficient φ and normalized rich-club coefficient φ N results (Fig. 5) , reveals that there is no evidence of rich-club formation of the high degree hubs [71] . Since the normalized rich-club coefficient φ N 1 for the primary network, the high degree hubs are sparsely distributed in the OC network topology and this value of φ N is still found to be smaller for those networks after few hubs are knockout N c 10 indicating the correlation among the high degree hubs become decreased significantly as more number of highest rank hubs are removed from the network [72] . This weakening in the correlation among the rich degree hubs in the complex OC network reduces the number of attractors of the nodes in the network and loosen the connectivity among the low degree and peripheral nodes in the network causing drastic functional diversity [73] . Hence, the dominance of controlling of the network by these rich-degree hubs such signal processing, functional structuring and many other regulating mechanisms will be less probable in OC network. Even though the hubs having k ≥ 30 have the tendency to tie tightly among themselves (increasing features of φ and φ N in Fig. 5) , their controlling dominance is weak (φ N ). Since the OC network is hierarchical topology, the loosely bound low degree nodes and moderate degree hubs to the high degree hubs might have segregated at cancer state and organized forming diverse loosely connected communities/subcommunities [71] to maintain optimal efficiency of the system's properties. Hence, the signaling behavior of high degree hubs could be different from that of lower degree nodes (ǫ s ǫ h , κ s κ h and δ s δ h in Fig. 4 ) in order to keep network organization optimal even at cancer state. Now, in order to identify the type of the hubs which take prominent responsibility for OC network regulation we calculated participation ratio P i of the nodes with respect to Z i score (see Methods) showing the different types of hubs for different number of leading hubs removed (Fig. 6 ). The analysis of the evolution of the different types of hubs in (P i , Z i )-plane as a function of N c show that provincial hubs are more prominent as compared to other types of hubs, namely, kinless and connector hubs. For this we defined the probability of finding each type of hubs P H as the ratio of the number of each type of hub N u ; ∀u ∈ {k, q, p} (N k → number of kinless hubs, N q → number of connector hubs and N p → number of provincial hubs) to the size of the network N as given below, From the Fig. 7 we could see that probability of finding kinless hubs is quite low as compared to other types of the hubs (connector and provincial hubs) and could able to see when one removed N c = 25, 30 and surprisingly there is no signature of it for N c 30 and N c 25. However, probability of finding of connector and provincial hubs could be seen starting from primary network and in the regime 0 ≤ N c ≤ 50 connector hubs abundance dominate that of provincial hubs in the network. This indicates that even though first leading 50 hubs from the OC network are removed or their regulatory functions are suppressed due to signal manipulation at cancer state, these connector and provincial hubs take over main responsibilities to keep network integrity, stability, regulation and preserve network topological properties. As the OC severity increases (cancer state moves to stage I, II, III etc.) the number of genes expressed/suppressed start increasing [76] most probably causing in increase in the number of manipulated high degree hubs in the OC network [77] . In such situation cancer driven hub manipulation, both connector and provincial hubs take major responsibility in network regulation of silencing of number of hubs N c ≤ 100. If the number of manipulated hubs is N c 100 only provincial hubs take main responsibility of network integrity and regulation. The key regulators of OC network [29] : CD44, EPCAM, KRAS and MCAM are found to be provincial hubs, and AKT1 falls in connector hub (Fig. 5b lowest panel) . All the key regulators do no change their hub identity in the knockout experiment of N C = 10, except CD44, where, it changes from provincial to connector as around N C = 40 hubs are knockout from the OC network. This indicates that provincial and connector hubs are quite important in understanding OC network regulation in maintaining network stability, properties and self-organization when the network is perturbed externally and internally by cancer signal. Patterns of mRNA expressions of the key regulators show diverse activities Now, we studied the activities of the key genes [29] from the mRNA expressions of the five key regulators using GEPIA (Gene Expression Profiling Interactive Analysis) online server (see Methods) where p-value P 0.01. The mRNA expression levels of these key regulators show diverse activities at OC state as compared to that of normal expressions (Fig. 6a) . We found that EPCAM was over expressed at ovarian cancer state, whereas, MCAM was down regulated but CD44 was significantly up regulated as compared to that of normal state. However, AKT1 and KRAS were slightly up regulated as compared to that of normal state. This is quite strongly correlated to the observed survival analysis of these key regulators [29] , where, some regulators (AKT1, CD44, MCAM) act as cancer suppressor, whereas, some of them (EPCAM and KRAS) work in the cancer progression. Hence, expressions AKT1, CD44, EPCAM and KRAS can be correlated to increasing or decreasing the risk of disease progression, and hence, they can be potential prognostics markers or drug targets in OC carcinogenesis. Then, the mRNA expression levels of the key regulators are calculated and clearly show that at stage-II all the key genes unregulated maximally, and then slowly down regulated as the cancer stage approaches from II → IV (Fig. 6b) . Moreover, another important aspect of these results is that the overall mRNA expression level EPCAM at any OC state and it's expression level at any stage of the OC is maximum, whereas, those of KRAS are minimum. Further, the expression levels of EPCAM at any stages (II-IV) of OC remain equally over expressed. Across all stages (I-III), EPCAM was the most frequently expressed biomarker [80] . EPCAM could therefore be an ideal tumor antigen candidate for detecting cancer cells. Overexpression of EPCAM was connected to advanced stage of disease and poor overall survival in some tumor types, suggesting EPCAM as a possible prognostic marker. AKT activation is prevalent in high-grade, late-stage serous ovarian carcinomas (6-9), according to immunohistochemical investigations, and may thus play a role in promoting tumor development [81] [82] [83] [84] [85] . The TIMER web tool [86] analysis showed that the expression of AKT1, CD44, EPCAM, KRAS and MCAM genes was significantly correlated with one or more tumor-infiltrating immune subsets (Fig. 7) . We found that the immune infiltrates, namely, CD4+ T-cells show positive correlations with the expressions of AKT1, CD44, MCAM and KRAS; macrophages have positive correlation with KRAS and MCAM; and neutrophils with the expression of KRAS. This indicates that the key regulators have the most significant correlation with immune infiltrates. Important aspects of KRAS is that it's expression has positive correlation with maximal infiltrates (Purity, CD4+ T cell, macrophage, neutrophil, dendritic cell). Hence, it is quite important to notice that these key regulators involve in immune infiltrates and regulations connecting to OC with immune system [93] . The identified novel key regulators of OC network using network theoretical approach are shown to be potential drug targets to prevent/cure this disease because of their high connectivity with the rest of the nodes (both high and low degree) in the network and the ability to rapidly transfer information in the network. Key regulator genes could be the most appropriate target in complex biological network for drug detection because of their high connection and capacity to rapidly transfer information [94] . Hence, we use these five regulators for possible identification of available list of the FDA approved drugs. Now a total of 219 drug-gene interactions for these key genes were identified (Supplementary Table) . We found that out of these available list of drugs, Uprosertib, Progesterone, Solitomab and Regorafenib drugs showed significant strong interactions with the key genes AKT1, CD44, EPCAM and KRAS respectively having the highest score (Fig. 8A) . Many of the drugs targeting these key regulators can be used either individually or in combination. Uprosertib is a reversible pan-AKT inhibitor that interacts with adenosine triphosphate [95] [96] [97] . A recent study demonstrated that ovulation, the release of hormones (e.g. progesterone and estrogen), inflammatory factors and growth factors among others, promoted the migration of intrauterine-injected malignant cells towards the ovarian stromal compartment to form ovarian tumors [98] . Solitomab is an EPCAM antibody that is being developed for the treatment of patients with multiple EPCAM-expressing solid tumors [99, 100] . Because of its capacity to engage resting, polyclonal, CD8+, and CD4+ T lymphocytes for very potent, redirected destruction of target tumor cells, solitomab has exhibited exceptional anticancer activity in preclinical ovarian tumor xenograft models [100] . Regorafenib is an anti-angiogenic oral multikinase inhibitor that has been approved for use in a variety of solid cancers [101] . We then did molecular docking analysis of the filtered drug compounds with the AKT1, CD44, EPCAM, KRAS and MCAM proteins as shown in Fig. 8B . Based on their binding affinities and visual inspection of docked complexes for their ability to shape H-bonds and interactions key genes and ligands (Fig. 8B) . Binding affinities of AKT1, CD44, EPCAM and KRAS with drug Uprosertib (-7.723 kcal/mol), Progesterone (-3.715 kcal/mol), Solitomab (-6.877 kcal/mol) and Regorafenib (-3.934) respectively, showed greatest affinity of binding among these. The key protein (AKT1, CD44, EPCAM and KRAS) is bind with drug with hydrogen bond and other non-polar bond with residues of key proteins (Fig. 3B ). Identification of driver genes in OC and understanding their roles are still open questions and critical study in this field is necessary to prevent this disease and also finding a possible drug target/s to cure this disease. The study can be done from diverse directions of experimental techniques, theoretical methods, data analysis etc with interdisciplinary perspective. We mined experimentally verified genes, which involved in OC, from various dedicated databases and analyzed that data from network theoretical approach to understand how OC network is organized, what is the role of the hubs, what could be the mechanism of signal manipulation in OC etc. The OC network carries fractal properties which is reflected in the topological properties of the network as power law distributions. If one goes deeper inside the network through communities/subcommunities till basic building block motif level (OC network has five hierarchical levels of organization) one can observe self-similar power law distributions which could be the origin of fractal nature in the OC network. These communities/subcommunities can be think of as coarse grained functional modules which can be scale down to basic functional units whose functional properties are self-similar because the topological properties of these functional modules at all levels of organization obey one parameter scaling law where all the curves collapse to a single curve by interpolating the curves. We proved it from both from data analysis and general fractal theory approach [46] and breakdown of any of the functional module caused by removing by knocking out few leading hubs from the network do not cause OC network breakdown. This property reveals that OC network does not follow centrality lethality rule [102] and self-organized. Understanding the role/s of the hubs are important because driver genes are generally one of these hubs and they are meant for potential drug target to enable to eradicate OC disease. The hubs knockout experiment we performed clearly shows us few important features in the OC network. First, removal of significant number of hubs from the OC network allows the transition from hierarchical network to scale free topology. Continued removal of hubs bring changes from the scale free topology to random feature which is the signature of network breakdown. This indicates the significant role of high degree hubs in shaping the OC network topology as well as keeping network stability. Second, one can observe the decrease in the rate of information processing of the leading hubs which could be due to the manipulation of activity of hub nodes as compared to the activity of lower degree hubs. This manipulation of activities of high degree hubs at cancer state leads to decrease in the number of 'attractors' in the network structure causing segregation of lower degree hubs with short paths to become more active. However, the nodes in the OC network try to keep the network properties, stability and organization after rewiring among them. There are different types of hubs in OC network but surprisingly it is not the high rank hubs but moderate rank 'provincial' and 'connector' hubs which are the keepers and safeguarders of the OC network at cancer state when the cancer signal spreads all over the entire network. There is no signature of formation of 'rich-clubs' by the high rank hubs in the OC network causing weakening the ties of the lower degree hubs in the network periphery with the hubs affecting in the functional diversity of the network. Hence, there is quite complex regulatory mechanisms in network organization when the system is at cancer state. Further, when the OC moves to severity condition i.e. stage I, II, III etc which could be modeled by allowing the network to destabilize more high degree hubs [76] the network regulation becomes more complicated. Hence, we need to study further on the cancer data at various severity stages so that one can detect OC at early stage for prevention. Another important aspect is to study the proper analysis of these key genes and clinical trials for possible drug discovery to eradicate this cancer type. We then studied the gene expression patterns of these key regulators at OC state and their behaviors at different stages of the OC. The analysis indicates that EPCAM has maximum over expression overall at various stages of the cancer state which indicates this gene behaves as cancer enhancer in the OC progression [103] . It is also found that CD44 is also highly expressed as compared to normal state and could also be another important target gene [104] . Similarly, other remaining regulators are also important for analysis in terms of drug targets and other analysis. These regulators are also positively correlated with immune infiltrates such as CD4+ T cells etc indicating the correlation of OC with immune system. Although CD44 expression has been observed during the progression of ovarian carcinoma, the questions of whether high CD44 expression indicates a relatively favourable prognosis [106] or aggressive tumor behaviour and a poor prognosis [106] and whether CD44 has any prognostic significance [107] [108] [109] have remained unresolved. To better understand the significance of CD44 expression in epithelial ovarian cancer, researchers used immunohistochemistry to examine a large number of samples and examined the prognostic value of CD44 using multivariate and univariate analyses [110] . The finding points CD44 playing a role in establishing treatment resistance by stimulating the PI3K/Akt signalling pathway [112] . In preclinical models, an RNA-based bispecific CD44-EPCAM aptamer was created and tested. By combining single CD44 and EPCAM aptamers with a double-stranded RNA adaptor, it can block CD44 and EPCAM at the same time. In vitro and in vivo investigations, the bispecific CD44-EPCAM aptamer was much more effective than either single CD44 or EPCAM aptamer at suppressing cell proliferation or inducing apoptosis in ovarian cancer cells [111] . The findings show that the KRAS-variant is a new genetic marker for OC risk, and that it could explain genetic risk in a number of hereditary ovarian cancer and breast syndrome families, according to the researchers [113] . In an in vitro situation similar to the patient's tumour, proliferation is aided by MER-mediated ERK and AKT signalling, which enhances phosphorylated nuclear ER. When MAPK activity is increased, inhibiting both the Mutated, Estrogen Receptor, and ERK pathways is required to stop proliferation, which may explain the use of aromatase inhibitors rather than tamoxifen [114] . When there are significant genomic co-alterations, such as KRAS or PI3K, it has been suggested that personalised combination therapy may be a key to establishing tumur control [117] , especially when certain gene product pathways, such as KRAS or PI3K, are activated [118] . It could be a useful indicator for predicting ovarian cancer metastasis, recurrence, and drug resistance, as well as a diagnostic, therapeutic, and prognostic tool. However, it's surprising that the stronger key regulators expression is, the earlier the clinical stage is, implying that it's a protooncogene. However, in the context of prognosis, it also functions as a tumour suppressor gene. This issue will be the subject of our future research. Further, network theoretical approach was used to identify potential drugs from the list of FDA approved drugs which have high affinity to bind with the key proteins and high score in molecular docking analysis. The identified drugs could be useful for the prevention of this disease and need to be done critical study on it both from theoretical and experimental perspectives. We propose a systematic investigation on the identified key/driver genes and need to do various clinical trials and experiments. The methods and algorithms we have used to study OC network are summarized in the flowchart in the Fig. 1 . We chose six highly cited ovarian cancer databases, namely, COSMIC database, Gene cards database, Ovarian kaleidoscope database, Dragon database of ovarian cancer, Curated ovarian database and OC-Gene database, and mined 4818 experimentally verified genes which are involved in OC. We then curated the identified list of the genes by using CGI code written in Perl and removed gene duplication and involved aliases in the list. From this method we could able to retrieve 660 genes out of 2000. This curated list of the genes is then cleaned manually as well as using Cytoscape 3.7.1 plugin by mapping to UniProt (January 2016), and we arrived at final list of 600 genes. Now following the principle of one gene one protein idea we constructed OC primary network using the curated 600 genes using GeneMANIA App [115] by uploading the 600 genes in Cytoscape 3.7.1 [116] . The primary network so constructed is composed of 4818 with 16320 physical interactions among the involved nodes. The knockout experiments of the high degree hubs were carried out by removing the first 10-400 leading hubs along with their associated edges from the primary OC network systematically and analyzed the topological properties of the resulting networks. Consider a complex network of size N defined by a graph G(V, E), where, V = {n i ; i = 1, 2, ..., N ; ∀n i ∈ V } and E = {e ij ; i, j = 1, 2, ..., N ; ∀e ij ∈ E} are sets of nodes and edges respectively. Since the PPI we constructed is undirected and unweighted network, G(V, E) can be constructed by defining an adjacency matrix [A] N ×N with matrix elements, A ij = 1 if there is an edge between 'i'th and 'j'th nodes in the network, otherwise A ij = 1, and A is symmetric matrix [50] . We used Cytoscape3.6.0 with modules Network analyzer and CytoNCA [38, 78, 79] and Igraph for calculating the topological properties of the OC network. We employed Louvain community detection algorithm [47] for digging communities/subcommunities at various levels of organization embedded in OC network. The topological properties of a complex network which are calculated for analysis are the following: Degree (k): This network metric is the total number of links a node has in a network. If k i is the degree of its node of the network, then it can be obtained by, This network metric characterize how much a node can influence other remaining nodes in the network [50] . Larger the value of k i more influential the its node is and vice versa. Probability of degree distribution, P (k): It indicates the measure of the probability of finding nodes of degree k in the network and can be calculated as follows, where, n k is the total number of nodes with degree k in the network. This network metric can able to characterize topological features of different types of networks, such as, scale free and hierarchical network follow power law distribution with k, whereas, small world or random network obey Poisson distribution with k [5, 32, 119] . Clustering coefficient (local) C i (k): It is a network metric that characterize the ability of a node to have local internal links with its nearest neighbors, which can be measured by by the ratio of the number of linked pairs of neighbors of a node (L i ) to total number of possible pairs of neighbors of that node [50] . It is given by, Higher the value of C i local nodes are densely connected. However, higher degree nodes generally have lower value of clustering coefficient [32] . For hierarchical network C i can be conjectured by C i (k) ∼ k α with α ≤ 1, whereas, for scale free, random and small world network C i (k) → constant [5, 32, 119] . One can also calculate average clustering coefficient of the network by using, Neighborhood connectivity C N (k): This metric characterize the average ties of a node with the nearest neighbours of degree k in a network defined as follows [69] , where, P (s|k) is conditional probability of the links of a node with k connections to another node having s connections. If the C N (k) exhibit power law behavior C N (k) ∼ k ±β with β ≤ 0.5, then the network falls to have in hierarchical features [69] , where, +β behavior is the signature of assortativity property of the network, whereas, −β indicates disassortativity nature of the network. Networks with assortativity nature is the evidence of rich − club formation of high degree hubs in the network [71] . Rich-club coefficient φ(k): This metric is the measure of how well high degree hubs establish connections among themselves to control the network regulation, and can be defined as the ratio of the number of links connected among the rich-club member hubs (U k ) to the maximum possible links they (M k ) can have [71] , as given below, Now the normalized rich-club coefficient can then be written as follows [71] , where, φ rand (k) is the calculated rich-club coefficient of the generated random networks with similar size and degree sequence. If the network have signature of rich-club then φ norm (k) 1, otherwise (φ norm (k) 1), there is rich-club formation. Betweenness centrality C B : It is the measure of the influence of a node has to propagate an information in the other nodes of the network [50, 120, 122] which can be quantified as follows, where, d ij (v) denotes the number of geodesic paths from node i to node j passing through node v and d ij is the all possible number of geodesic paths from node i to node j. Hence, normalized betweenness centrality C B is given by, where, M denotes the number of node pairs, excluding v Higher the value of C B faster the information processing throughout the network and vice versa. Closeness centrality C C : It can be regarded as the measure of how far an information can be propagated from a certain node to the rest of the other nodes in the network [50, 120, 121] . This metric of any node i can be calculated from the harmonic mean of the all possible geodesic paths d ij linking to the node i, Hence, larger the value of C C farther the information can be reached, and vice versa [50] . Eigenvector centrality C E : This network metric is the measure of strength of influence of a node to the rest of the nodes in the network having diverse influential capabilities [67] . The eigenvector centrality of a node i in a network can be obtained, where, nn(i) represents the nearest neighbors of node i in the network with eigenvalue λ and eigenvector v i in the eigenvalue equation, Av i = λv i where, A is the network adjacency matrix. C E values of all nodes in a network are always positive C E 0, and it is considered to be the extension of degree centrality [50] . Techniques to classify hubs Classification of different types of hubs in a network can be done by the following network metrics. Z-score parameter (Z i ): This parameter is generally known as within module Z-score parameter which measures how strongly a particular node in a module is linked to the rest of the nodes in the module [123] . The Z-score of a node i having k i links in certain module "c" can be defined as, where, k c and σ kc are the expectation value and standard deviation of degrees k of all nodes in the module respectively. This Z−score can able to distinguish hubs from non-hub nodes depending on its value, namely, modular − hubs if Z i ≥ 2.5 (nodes having large number of links) and non − hubs if Z i ≤ 2.5 (lower degree nodes) [123] . Participation co-efficient (P i ): The local and global roles of the identified modular hubs by Z-score can be further classified depending on the participation of the hubs locally within their own modules and globally by establishing links to other modules other than its own module. This participation co-efficient of node i in a network having N M modules, which have degree k i , can be defined as [123] , where, k is denotes the degree of node i in c module. If P i → 1 then the node i is uniformly linked to all other modules, whereas, if P i → 0 the node links to the nodes within its own module [123] . Generally, nodes which have P i 0.7 has larger participation in other modules. Now, different types of hubs can be classified depending on the values of Z-score and P i in (P i , Z i )-plane as follows: 1. kinless hubs: Hubs in the range Z i ≥ 2.5 & P i > 0. 70, 2. connector hubs: in the area Z i ≥ 2.5 & 0.30 < P i ≤ 0.75, and 3. provincial hubs: in the area Z i ≥ 2.5 & P i ≤ 0.30 in the (P i , Z i ) space [123] . The GEPIA [124] for Cancer Genomics was used to investigate the expression pattern of the key genes. Tumor or normal differential expression analysis, profiling according to cancer subtypes or pathological stages, patient survival analysis, similar gene finding, correlation analysis, and dimensionality reduction analysis are all customizable functions in GEPIA [124] . Correlation of gene expression with tumor-infiltrating immune cells: Considering the importance of immune dysregulation in ovarian cancer, we explored the correlation between AKT1, CD44, EPCAM, KRAS and MCAM mRNA expression and tumor-infiltrating immune cells. The web tool TIMER (https://cistrome.shinyapps.io/timer/) [86] was used. Six tumor-infiltrating immune subsets, neutrophils, macrophages, CD8+ T cells, CD4+ T-cells and including B cells were analysed. The DGIdb [87] database was screened out for discovering therapeutic targets for key genes. The DGIdb was used to produce a similar drug-gene interaction network. DGIdb is a database that can anticipate the physical and functional interactions between key genes and drugs in a query. The interactions in the DGIdb database come from five different sources: automatic high-throughput lab experiment data, text mining, co-expression interaction data, genomic context interaction prediction, and previous knowledge from other databases. The X-ray structures of some key genes of ovarian cancer like AKT1, CD44, EPCAM, KRAS and MCAM are available on protein data bank (PDB) [88] corresponding to PDB ids 9S9W, 1POZ, 4MZV, 4EPT and 6LYN, respectively. Prior to docking studies, the proteins were processed using Schrodinger Protein Preparation Wizard during which missing hydrogens were added, hydrogen bonds were optimized and water molecules were deleted [89] . Nutraceuticals are bioactive phytochemicals that deliver health benefits and are relatively safe to be used for prevention and treatment against disease. A large number of drugs are FDA approved pharmaceutical drugs. The chemical structures of high score ligands from gene-drug interaction of key genes were obtained from PubChem [90] . The chemical structures of ligands were prepared with LigPrep module Schrodinger which generated diverse, accurate and energy minimized conformations using OPLS-2005 force field [89] . Docking studies using Glide: Glide uses a range of filters and performs a thorough search of the conformational, orientational and positional space for ligand in the binding site of receptor. The receptor was represented by creating a cubic grid of size 30x30x30Å centred on the active site residues for each protein by Schrodinger's Receptor Grid Generation program [91] . Further the ligands were screened against each receptor using high throughput virtual screening (HTVS) approach of Glide [91] . The high ranking ligands were then subjected to stringent screening via extra precision (XP) [91] method to eliminate false positives and acquire precise binding modes of compounds. Maestro was used to compute interactions between the docked protein-ligand complexes. 2D diagrams were generated depicting hydrogen bonds and hydrophobically interacting residues. Then PyMOL [92] was used to visualize and study the resultant receptor-ligand complexes. FIG. 8: Identification of potential drug compounds which have optimal binding affinity with the key regulators: (A) Interaction network of key regulators with FDA-approved drugs. The green circle represents drugs, blue color diamonds represents key regulators and red color drug shows highest score. (b) Molecular docking of key regulators and drugs with highest score. Uprosertib, Progesterone, Solitomab and Regorafenib drugs showed significant strong interactions with the key genes AKT1, CD44, EPCAM and KRAS respectively having the highest score. KC and RKBS are financially supported by UPE-II, under sanction no. 101 and DBT-COE. MZM financially supported by DHR The dynamic control of signal transduction networks in cancer cells The evolution of hierarchical gene regulatory networks Network motifs: simple building blocks of complex networks SnapShot: network motifs Network biology: Understanding the cells's functional organization Complex networks: Structure and dynamics The Multifractal Model of Asset Returns. Discussion papers of the Cowles Foundation for Economics. Yale University: Cowles Foundation 114-1166 Complex networks: Structure and dynamics Self-similarity of complex networks Rewiring makes the difference Regulation patterns in signaling networks of cancer Navigating cancer network attractors for tumor-specific therapy Increased signaling entropy in cancer requires the scale-free property of protein interaction networks Metabolic stability and epigenesis in randomly constructed genetic nets Cancer attractors: a systems view of tumors from a gene network dynamics and developmental perspective Unités biologiques douées de continuité génétique Teleonomic mechanisms in cellular metabolism, growth, and differentiation The origins of order Entropy rate of diffusion processes on complex networks Screening for ovarian cancer Surgical prevention strategies in ovarian cancer The changing view of high-grade serous ovarian cancer Epithelial ovarian cancer risk: A review of the current genetic landscape The hallmarks of cancer Cancer as an evolutionary and ecological process Exploring novel key regulators in breast cancer network Assessment of the key regulatory genes and their Enterology for Turner Syndrome employing network approach Methodology of predicting novel key regulators in ovarian cancer network: A network theoretical approach Identification of Novel Biomarkers and Candidate Drug in Ovarian Cancer Identification of novel epithelial ovarian cancer loci in women of African ancestry Statistical mechanics of complex networks Deregulation of an imprinted gene network in prostate cancer The metabolic landscape of lung cancer: new insights in a disturbed glucose metabolism Dysregulated signaling hubs of liver lipid metabolism reveal hepatocellular carcinoma pathogenesis On measuring the complexity of networks: Kolmogorov complexity versus entropy Review of graph and network complexity from an algorithmic information perspective A mathematical theory of communication Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks Three approaches to the quantitative definition of information Decomposing protein interactome networks by graph entropy Community structure in social and biological networks Structure and function of the feed-forward loop network motif Network motifs: theory and experimental approaches One-Parameter Scaling of Localization Length and Conductance in Disordered Systems The Multifractal Model of Asset Returns Cowles Foundation discussion paper no. 1164 Fast unfolding of communities in large networks Community detection in networks with positive and negative links Networks: An introduction A small world of weak ties provides optimal global integration of self-similar modules in functional brain networks The human connectome: a complex network Complexity and self-organization Geometric fractal growth model for scale-free networks The diameter of random sparse graphs A first look at graph theory A measure of betweenness centrality based on random walks Betweenness centrality in large complex networks Universal behavior of load distribution in scale-free networks Large-scale topological and dynamical properties of the Internet Centrality in social networks conceptual clarification A compendium of mutational cancer driver genes Eigenvector-like measures of centrality for asymmetric relations Justification and application of eigenvector centrality Dynamical and Correlation Properties of the Internet Assortative Mixing in Networks Detecting rich-club ordering in complex networks Rich-club organization of the human connectome Rich club organization supports a diverse set of functional network configurations Kinless hubs are potential target genes in prostate cancer network Identification of key regulators in Prostate cancer from gene expression datasets of patients Gene expression patterns in ovarian carcinomas The biology of ovarian cancer: new opportunities for translation Topological analysis and interactive visualization of biological networks and protein structures CytoNCA: A cytoscape plugin for centrality analysis and evaluation of protein interaction networks Epithelial cell adhesion molecule expression (CD326) in cancer: a short review Frequent activation of AKT2 and induction of apoptosis by inhibition of phosphoinositide-3-OH kinase/Akt pathway in human ovarian cancer AKT and mTOR phosphorylation is frequently detected in ovarian cancer and can be targeted to disrupt ovarian tumor cell growth A specific role for AKT3 in the genesis of ovarian cancer through modulation of G(2)-M phase transition Frequent loss of PTEN expression is linked to elevated phosphorylated Akt levels, but not associated with p27 and cyclin D1 expression, in primary epithelial ovarian carcinomas Akt1 and akt3 exert opposing roles in the regulation of vascular tumor growth TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells DGIdb 3.0: a redesign and expansion of the drug-gene interaction database The Protein Data Bank Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments PubChem substance and compound databases Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening Ligand docking and binding site analysis with PyMOL and Autodock/Vina Ovarian cancer creates a suppressive microenvironment to escape immune elimination Brain disease network analysis to elucidate the neurological manifestations of COVID-19 Characterization of a chemical affinity probe targeting Akt kinases Discovery of novel AKT inhibitors with enhanced anti-tumor effects in combination with the MEK inhibitor Phase I dose-escalation trial of the oral AKT inhibitor uprosertib in combination with the oral MEK1/MEK2 inhibitor trametinib in patients with solid tumors Epigenetic modification restores functional PR expression in endometrial cancer cells MT110: a novel bispecific single-chain antibody construct with high efficacy in eradicating established tumors Solitomab, an epithelial cell adhesion molecule/CD3 bispecific antibody (BiTE), is highly active against primary chemotherapy-resistant ovarian cancer cell lines in vitro and fresh tumor cells ex vivo Phase I study of regorafenib and sildenafil in advanced solid tumors Lethality and centrality in protein networks Simultaneous targeting of CD44 and EpCAM with a bispecific aptamer effectively inhibits intraperitoneal ovarian cancer growth Up-regulation of CD44 in the development of metastasis, recurrence and drug resistance of ovarian cancer Genomic and transcriptomic profiling expands precision cancer medicine: The WINTHER trial Expression of the hyaluronan receptor, CD44s, in epithelial ovarian cancer is an independent predictor of survival Expression of CD44 in effusions of patients diagnosed with serous ovarian carcinoma: diagnostic and prognostic implications CD44 expression in benign, premalignant, and malignant ovarian neoplasms: relation to tumour development and progression CD44 variant expression is a common feature of epithelial ovarian cancer: lack of association with standard prognostic factors CD44 expression indicates favorable prognosis in epithelial ovarian cancer Simultaneous targeting of CD44 and EpCAM with a bispecific aptamer effectively inhibits intraperitoneal ovarian cancer growth Akt signaling is sustained by a CD44 splice isoform-mediated positive feedback loop A KRAS-variant in ovarian cancer acts as a genetic marker of cancer risk Estrogen Receptor-Positive Low-Grade Serous Ovarian Cancer: Unraveling an Exceptional Response Mystery GeneMANIA: a real-time multiple association network integration algorithm for predicting gene function Cytoscape: a software environment for integrated models of biomolecular interaction networks Molecular profiling of cancer patients enables personalized combination therapy: The I-PREDICT study Prognostic implications of RAS alterations in diverse malignancies and impact of targeted therapies Hierarchical organization in complex networks A Graph-theoretic perspective on centrality Communication Patterns in Task-Oriented Groups A Set of Measures of Centrality Based on Betweenness Functional cartography of complex metabolic networks GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses Correlation of the OC key regulators with various cancer-infiltrating immune subsets: Plots of the gene expression levels of the OC key regulators (AKT1, CD44, EPCAM, KRAS and MCAM) as a function of infiltration levels for various cancer-infiltrating immune subsets via the TIMER database Author Contributions RKBS conceived the model. KC, MZM and RKBS did the numerical experiment. KC, MZM and RKBS prepared the figures of the numerical results, analyzed and interpreted the simulation results. RKBS, KC and MZM wrote the manuscript. KC, MZM, PS and RKBS involved in the study and reviewed the manuscript. Competing financial interests: The authors declare no competing financial interests.