Stability and analysis of multi-agent network systems with state-dependent switching typologies have been a fundamental and longstanding challenge in control, social sciences, and many other related fields. These already complex systems become further complicated once one accounts for asymmetry or heterogeneity of the underlying agents/dynamics. Despite extensive progress in analysis of conventional networked decision systems where the network evolution and state dynamics are driven by independent or weakly coupled processes, most of the existing results fail to address multi-agent systems where the network and state dynamics are highly coupled and evolve based on status of heterogeneous agents. Motivated by numerous applications of such dynamics in social sciences, in this paper we provide a new direction toward analysis of dynamic networks of heterogeneous agents under complex time-varying environments. As a result we show how Lyapunov stability and convergence of several challenging problems from opinion dynamics can be established using a simple application of our framework. Moreover, we introduce a new class of asymmetric opinion dynamics, namely nearest neighbor dynamics, and show how our approach can be used to analyze their behavior. In particular, we extend our results to game-theoretic settings and provide new insights toward analysis of complex networked multi-agent systems using exciting field of sequential optimization.This function can also be written in a compact form as f (y, λ) = y T Ly − tr(L) 2 , where L := diag(λ1) − λ, and tr(·) denotes the trace function. Intuitively, the block variable λ is meant to capture the communication network G t , and the block variable y captures the opinion states. Note that if we restrict λ ij s to binary variables in {0, 1}, then λ simply represents the adjacency matrix of a network of n agents where λ ij = 1 if there is a directed edge (i, j) from node i to node j, and λ ij = 0 otherwise. Moreover, for such a binary block variable λ, the matrix L is precisely the Laplacian matrix of the communication network associated with λ. Although we still need to assume that λ ∈ {0, 1} n×n , to avoid complication of handling integral variables, for now we allow λ ij s to vary continuously in the interval [0, 1]. As we shall see soon the integrality of network variables will be automatically achieved during iterations of the BCD method. Now let us consider the BCD method applied to the objective function (5) with block variables x and λ. For a generic time t, let us fix the state variable to y = x(t). Minimizing (5) with respect to the network variable λ ∈ Λ = [0, 1] n×n we obtain, λ t := arg min λ∈[0,1] n 2 i,j