ASU Scholarship Showcase

In spite of the recent interest and advances in linear controllability of complex networks, controlling nonlinear network dynamics remains an outstanding problem. Here we develop an experimentally feasible control framework for nonlinear dynamical networks that exhibit multistability. The control objective is to apply parameter perturbation to drive the system from one attractor to another, assuming that the former is undesired and the latter is desired. To make our framework practically meaningful, we consider restricted parameter perturbation by imposing two constraints: it must be experimentally realizable and applied only temporarily. We introduce the concept of attractor network, which allows us to formulate a quantifiable controllability framework for nonlinear dynamical networks: a network is more controllable if the attractor network is more strongly connected. We test our control framework using examples from various models of experimental gene regulatory networks and demonstrate the beneficial role of noise in facilitating control.

Nocturnal cooling of urban areas governs the evolution of thermal state and many thermal-driven environmental issues in cities, especially those suffer strong urban heat island (UHI) effect. Advances in the fundamental understanding of the underlying physics of nighttime UHI involve disentangling complex contributing effects and remains an open challenge. In this study, we develop new numerical algorithms to characterize the thermodynamics of urban nocturnal cooling based on solving the energy balance equations for both the landscape surface and the overlying atmosphere. Further, a scaling law is proposed to relate the UHI intensity to a range of governing mechanisms, including the vertical and horizontal transport of heat in the surface layer, the urban-rural breeze, and the possible urban expansion. The accuracy of proposed methods is evaluated against in-situ urban measurements collected in cities with different geographic and climatic conditions. It is found that the vertical and horizontal contributors modulate the nocturnal UHI at distinct elevation in the atmospheric boundary layer.

Network reconstruction is a fundamental problem for understanding many complex systems with unknown interaction structures. In many complex systems, there are indirect interactions between two individuals without immediate connection but with common neighbors. Despite recent advances in network reconstruction, we continue to lack an approach for reconstructing complex networks with indirect interactions. Here we introduce a two-step strategy to resolve the reconstruction problem, where in the first step, we recover both direct and indirect interactions by employing the Lasso to solve a sparse signal reconstruction problem, and in the second step, we use matrix transformation and optimization to distinguish between direct and indirect interactions. The network structure corresponding to direct interactions can be fully uncovered. We exploit the public goods game occurring on complex networks as a paradigm for characterizing indirect interactions and test our reconstruction approach. We find that high reconstruction accuracy can be achieved for both homogeneous and heterogeneous networks, and a number of empirical networks in spite of insufficient data measurement contaminated by noise. Although a general framework for reconstructing complex networks with arbitrary types of indirect interactions is yet lacking, our approach opens new routes to separate direct and indirect interactions in a representative complex system.

Background: The cytokine MIF (Macrophage Migration Inhibitory Factor) has diverse physiological roles and is present at elevated concentrations in numerous disease states. However, its molecular heterogeneity has not been previously investigated in biological samples. Mass Spectrometric Immunoassay (MSIA) may help elucidate MIF post-translational modifications existing in vivo and provide additional clarity regarding its relationship to diverse pathologies.

Results: In this work, we have developed and validated a fully quantitative MSIA assay for MIF, and used it in the discovery and quantification of different proteoforms of MIF in serum samples, including cysteinylated and glycated MIF. The MSIA assay had a linear range of 1.56-50 ng/mL, and exhibited good precision, linearity, and recovery characteristics. The new assay was applied to a small cohort of human serum samples, and benchmarked against an MIF ELISA assay.

Conclusions: The quantitative MIF MSIA assay provides a sensitive, precise and high throughput method to delineate and quantify MIF proteoforms in biological samples.

Our ability to uncover complex network structure and dynamics from data is fundamental to understanding and controlling collective dynamics in complex systems. Despite recent progress in this area, reconstructing networks with stochastic dynamical processes from limited time series remains to be an outstanding problem. Here we develop a framework based on compressed sensing to reconstruct complex networks on which stochastic spreading dynamics take place. We apply the methodology to a large number of model and real networks, finding that a full reconstruction of inhomogeneous interactions can be achieved from small amounts of polarized (binary) data, a virtue of compressed sensing. Further, we demonstrate that a hidden source that triggers the spreading process but is externally inaccessible can be ascertained and located with high confidence in the absence of direct routes of propagation from it. Our approach thus establishes a paradigm for tracing and controlling epidemic invasion and information diffusion in complex networked systems.

Background: Cystatin C (CysC) is an endogenous cysteine protease inhibitor that can be used to assess the progression of kidney function. Recent studies demonstrate that CysC is a more specific indicator of glomerular filtration rate (GFR) than creatinine. CysC in plasma exists in multiple proteoforms. The goal of this study was to clarify the association of native CysC, CysC missing N-terminal Serine (CysC des-S), and CysC without three N-terminal residues (CysC des-SSP) with diabetic chronic kidney disease (CKD).

Results: Using mass spectrometric immunoassay, the plasma concentrations of native CysC and the two CysC truncation proteoforms were examined in 111 individuals from three groups: 33 non-diabetic controls, 34 participants with type 2 diabetes (DM) and without CKD and 44 participants with diabetic CKD. Native CysC concentrations were 1.4 fold greater in CKD compared to DM group (p = 0.02) and 1.5 fold greater in CKD compared to the control group (p = 0.001). CysC des-S concentrations were 1.55 fold greater in CKD compared to the DM group (p = 0.002) and 1.9 fold greater in CKD compared to the control group (p = 0.0002). CysC des-SSP concentrations were 1.8 fold greater in CKD compared to the DM group (p = 0.008) and 1.52 fold greater in CKD compared to the control group (p = 0.002). In addition, the concentrations of CysC proteoforms were greater in the setting of albuminuria. The truncated CysC proteoform concentrations were associated with estimated GFR independent of native CysC concentrations.

Conclusion: Our findings demonstrate a greater amount of CysC proteoforms in diabetic CKD. We therefore suggest assessing the role of cystatin C proteoforms in the progression of CKD.

Two classes of scaling behaviours, namely the super-linear scaling of links or activities, and the sub-linear scaling of area, diversity, or time elapsed with respect to size have been found to prevail in the growth of complex networked systems. Despite some pioneering modelling approaches proposed for specific systems, whether there exists some general mechanisms that account for the origins of such scaling behaviours in different contexts, especially in socioeconomic systems, remains an open question. We address this problem by introducing a geometric network model without free parameter, finding that both super-linear and sub-linear scaling behaviours can be simultaneously reproduced and that the scaling exponents are exclusively determined by the dimension of the Euclidean space in which the network is embedded. We implement some realistic extensions to the basic model to offer more accurate predictions for cities of various scaling behaviours and the Zipf distribution reported in the literature and observed in our empirical studies. All of the empirical results can be precisely recovered by our model with analytical predictions of all major properties. By virtue of these general findings concerning scaling behaviour, our models with simple mechanisms gain new insights into the evolution and development of complex networked systems.

Insulin-like growth factor 1 (IGF1) is an important biomarker for the management of growth hormone disorders. Recently there has been rising interest in deploying mass spectrometric (MS) methods of detection for measuring IGF1. However, widespread clinical adoption of any MS-based IGF1 assay will require increased throughput and speed to justify the costs of analyses, and robust industrial platforms that are reproducible across laboratories. Presented here is an MS-based quantitative IGF1 assay with performance rating of >1,000 samples/day, and a capability of quantifying IGF1 point mutations and posttranslational modifications. The throughput of the IGF1 mass spectrometric immunoassay (MSIA) benefited from a simplified sample preparation step, IGF1 immunocapture in a tip format, and high-throughput MALDI-TOF MS analysis. The Limit of Detection and Limit of Quantification of the resulting assay were 1.5 μg/L and 5 μg/L, respectively, with intra- and inter-assay precision CVs of less than 10%, and good linearity and recovery characteristics. The IGF1 MSIA was benchmarked against commercially available IGF1 ELISA via Bland-Altman method comparison test, resulting in a slight positive bias of 16%. The IGF1 MSIA was employed in an optimized parallel workflow utilizing two pipetting robots and MALDI-TOF-MS instruments synced into one-hour phases of sample preparation, extraction and MSIA pipette tip elution, MS data collection, and data processing. Using this workflow, high-throughput IGF1 quantification of 1,054 human samples was achieved in approximately 9 hours. This rate of assaying is a significant improvement over existing MS-based IGF1 assays, and is on par with that of the enzyme-based immunoassays. Furthermore, a mutation was detected in ∼1% of the samples (SNP: rs17884626, creating an A→T substitution at position 67 of the IGF1), demonstrating the capability of IGF1 MSIA to detect point mutations and posttranslational modifications.

This study examines the spatial and temporal patterns of the surface urban heat island (SUHI) intensity in the Phoenix metropolitan area and the relationship with land use land cover (LULC) change between 2000 and 2014. The objective is to identify specific regions in Phoenix that have been increasingly heated and cooled to further understand how LULC change influences the SUHI intensity. The data employed include MODerate-resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) 8-day composite June imagery, and classified LULC maps generated using 2000 and 2014 Landsat imagery. Results show that the regions that experienced the most significant LST changes during the study period are primarily on the outskirts of the Phoenix metropolitan area for both daytime and nighttime. The conversion to urban, residential, and impervious surfaces from all other LULC types has been identified as the primary cause of the UHI effect in Phoenix. Vegetation cover has been shown to significantly lower LST for both daytime and nighttime due to its strong cooling effect by producing more latent heat flux and less sensible heat flux. We suggest that urban planners, decision-makers, and city managers formulate new policies and regulations that encourage residential, commercial, and industrial developers to include more vegetation when planning new construction.