Mathematical modeling for intracellular transport and binding of HIV-1 Gag proteins

Mathematical modeling and quantitative analysis of HIV-1 Gag trafficking and polymerization

Gag, as the major structural protein of HIV-1, is necessary for the assembly of the HIV-1 sphere shell. An in-depth understanding of its trafficking and polymerization is important for gaining further insights into the mechanisms of HIV-1 replication and the design of antiviral drugs. We developed a mathematical model to simulate two biophysical processes, specifically Gag monomer and dimer transport in the cytoplasm and the polymerization of monomers to form a hexamer underneath the plasma membrane. Using experimental data, an optimization approach was utilized to identify the model parameters, and the identifiability and sensitivity of these parameters were then analyzed. Using our model, we analyzed the weight of the pathways involved in the polymerization reactions and concluded that the predominant pathways for the formation of a hexamer might be the polymerization of two monomers to form a dimer, the polymerization of a dimer and a monomer to form a trimer, and the polymerization of two trimers to form a hexamer. We then deduced that the dimer and trimer intermediates might be crucial in hexamer formation. We also explored four theoretical combined methods for Gag suppression, and hypothesized that the N-terminal glycine residue of the MA domain of Gag might be a promising drug target. This work serves as a guide for future theoretical and experimental efforts aiming to understand HIV-1 Gag trafficking and polymerization, and might help accelerate the efficiency of anti-AIDS drug design.

The human immunodeficiency virus (HIV-1) is a retrovirus that causes acquired immunodeficiency syndrome (AIDS), an infectious disease with high annual mortality. Gag protein is the major structural protein of HIV-1 and can self-assemble into the HIV-1 sphere shell. Therefore, an in-depth understanding of Gag protein trafficking and polymerization is important for gaining further insights into the mechanisms of HIV-1 replication and the design of antiviral drugs. Through mathematical modeling, optimization and quantitative analysis, we hypothesized the budding and release time of virus-like particles and revealed that the dimer and trimer intermediates might be crucial in hexamer formation. We also concluded that the predominant pathways in hexamer formation might be the polymerization of two monomers to form a dimer, the polymerization of a dimer and a monomer to form a trimer, and the polymerization of two trimers to form a hexamer. Our analysis also suggested that the N-terminal glycine residue of the MA domain of Gag might be a promising drug target. These results serve as a guide for further theoretical and experimental efforts aiming to understand HIV-1 Gag trafficking and polymerization and could aid anti-AIDS drug design.

Crosstalk between pathways enhances the controllability of signalling networks

The control of complex networks is one of the most challenging problems in the fields of biology and engineering. In this study, the authors explored the controllability and control energy of several signalling networks, which consisted of many interconnected pathways, including networks with a bow-tie architecture. On the basis of the theory of structure controllability, they revealed that biological mechanisms, such as cross-pathway interactions, compartmentalisation and so on make the networks easier to fully control. Furthermore, using numerical simulations for two realistic examples, they demonstrated that the control energy of normal networks with crosstalk is lower than in networks without crosstalk. These results indicate that the biological networks are optimally designed to achieve their normal functions from the viewpoint of the control theory. The authors' work provides a comprehensive understanding of the impact of network structures and properties on controllability.

Modeling HIV-1 viral capsid nucleation by dynamical systems

There are two stages generally recognized in the viral capsid assembly: nucleation and elongation. This paper focuses on the nucleation stage and develops mathematical models for HIV-1 viral capsid nucleation based on six-species dynamical systems. The Particle Swarm Optimization (PSO) algorithm is used for parameter fitting to estimate the association and dissociation rates from biological experiment data. Numerical simulations of capsid protein (CA) multimer concentrations demonstrate a good agreement with experimental data. Sensitivity and elasticity analysis of CA multimer concentrations with respect to the association and dissociation rates further reveals the importance of CA trimer-of- dimers in the nucleation stage of viral capsid self- assembly.

Optimal control strategy for abnormal innate immune response

Innate immune response plays an important role in control and clearance of pathogens following viral infection. However, in the majority of virus-infected individuals, the response is insufficient because viruses are known to use different evasion strategies to escape immune response. In this study, we use optimal control theory to investigate how to control the innate immune response. We present an optimal control model based on an ordinary-differential-equation system from a previous study, which investigated the dynamics and regulation of virus-triggered innate immune signaling pathways, and we prove the existence of a solution to the optimal control problem involving antiviral treatment or/and interferon therapy. We conduct numerical experiments to investigate the treatment effects of different control strategies through varying the cost function and control efficiency. The results show that a separate treatment, that is, only inhibiting viral replication (u₁(t)) or enhancing interferon activity (u₂(t)), has more advantages for controlling viral infection than a mixed treatment, that is, controlling both (u₁(t)) and (u₂(t)) simultaneously, including the smallest cost and operability. These findings would provide new insight for developing effective strategies for treatment of viral infectious diseases.

A new asynchronous parallel algorithm for inferring large-scale gene regulatory networks

The reconstruction of gene regulatory networks (GRNs) from high-throughput experimental data has been considered one of the most important issues in systems biology research. With the development of high-throughput technology and the complexity of biological problems, we need to reconstruct GRNs that contain thousands of genes. However, when many existing algorithms are used to handle these large-scale problems, they will encounter two important issues: low accuracy and high computational cost. To overcome these difficulties, the main goal of this study is to design an effective parallel algorithm to infer large-scale GRNs based on high-performance parallel computing environments. In this study, we proposed a novel asynchronous parallel framework to improve the accuracy and lower the time complexity of large-scale GRN inference by combining splitting technology and ordinary differential equation (ODE)-based optimization. The presented algorithm uses the sparsity and modularity of GRNs to split whole large-scale GRNs into many small-scale modular subnetworks. Through the ODE-based optimization of all subnetworks in parallel and their asynchronous communications, we can easily obtain the parameters of the whole network. To test the performance of the proposed approach, we used well-known benchmark datasets from Dialogue for Reverse Engineering Assessments and Methods challenge (DREAM), experimentally determined GRN of Escherichia coli and one published dataset that contains more than 10 thousand genes to compare the proposed approach with several popular algorithms on the same high-performance computing environments in terms of both accuracy and time complexity. The numerical results demonstrate that our parallel algorithm exhibits obvious superiority in inferring large-scale GRNs.

Deciphering deterioration mechanisms of complex diseases based on the construction of dynamic networks and systems analysis

The early diagnosis and investigation of the pathogenic mechanisms of complex diseases are the most challenging problems in the fields of biology and medicine. Network-based systems biology is an important technique for the study of complex diseases. The present study constructed dynamic protein-protein interaction (PPI) networks to identify dynamical network biomarkers (DNBs) and analyze the underlying mechanisms of complex diseases from a systems level. We developed a model-based framework for the construction of a series of time-sequenced networks by integrating high-throughput gene expression data into PPI data. By combining the dynamic networks and molecular modules, we identified significant DNBs for four complex diseases, including influenza caused by either H3N2 or H1N1, acute lung injury and type 2 diabetes mellitus, which can serve as warning signals for disease deterioration. Function and pathway analyses revealed that the identified DNBs were significantly enriched during key events in early disease development. Correlation and information flow analyses revealed that DNBs effectively discriminated between different disease processes and that dysfunctional regulation and disproportional information flow may contribute to the increased disease severity. This study provides a general paradigm for revealing the deterioration mechanisms of complex diseases and offers new insights into their early diagnoses.