Circadian (∼24 h) clocks are self-sustained endogenous oscillators with which organisms keep track of daily and seasonal time. Circadian clocks frequently rely on interlocked transcriptional-translational feedback loops to generate rhythms that are robust against intrinsic and extrinsic perturbations. To investigate the dynamics and mechanisms of the intracellular feedback loops in circadian clocks, a number of mathematical models have been developed. The majority of the models use Hill functions to describe transcriptional repression in a way that is similar to the Goodwin model. Recently, a new class of models with protein sequestration-based repression has been introduced. Here, the author discusses how this new class of models differs dramatically from those based on Hill-type repression in several fundamental aspects: conditions for rhythm generation, robust network designs and the periods of coupled oscillators. Consistently, these fundamental properties of circadian clocks also differ among Neurospora, Drosophila, and mammals depending on their key transcriptional repression mechanisms (Hill-type repression or protein sequestration). Based on both theoretical and experimental studies, this review highlights the importance of careful modelling of transcriptional repression mechanisms in molecular circadian clocks.

In this study, a model of basal ganglia (BG) is applied to develop a deep brain stimulation controller to reduce Parkinson's tremor. Conventionally, one area in BG is stimulated, with no feedback, to control Parkinson's tremor. In this study, a new architecture is proposed to develop feedback controller as well as to stimulate two areas of BG simultaneously. To this end, two controllers are designed and implemented in globus pallidus internal (GPi) and subthalamic nucleus (STN) in the brain. A proportional controller and a backstepping controller are designed and implemented in GPi and STN, respectively. The proposed controllers deliver suitable stimulatory control signals to GPi and STN based on hand tremor amplitude (as the feedback). When tremor reduces, these controllers decrease the stimulatory energy intensity proportionally. Therefore, additional stimulatory signal is not delivered to the brain. Subsequently, the side effects from the excessive stimulation intensity become much less. Comparing with one area stimulation, the results show that stimulating two areas of BG results in reduction of the level of the stimulation intensity. It is observed that these two controllers are both robust in terms of changing the system parameters.

Signalling pathway analysis is a popular approach that is used to identify significant cancer-related pathways based on differentially expressed genes (DEGs) from biological experiments. The main advantage of signalling pathway analysis lies in the fact that it assesses both the number of DEGs and the propagation of signal perturbation in signalling pathways. However, this method simplifies the interactions between genes by categorising them only as activation (+1) and suppression (−1), which does not encompass the range of interactions in real pathways, where interaction strength between genes may vary. In this study, the authors used newly developed signalling pathway impact analysis (SPIA) methods, SPIA based on Pearson correlation coefficient (PSPIA), and mutual information (MSPIA), to measure the interaction strength between pairs of genes. In analyses of a colorectal cancer dataset, a lung cancer dataset, and a pancreatic cancer dataset, PSPIA and MSPIA identified more candidate cancer-related pathways than were identified by SPIA. Generally, MSPIA performed better than PSPIA.

Combined antiretroviral therapy (cART) suppress HIV-1 viral replication, such that viral load in plasma remains below the limit of detection in standard assays. However, intermittent episodes of transient viremia (blips) occur in a set of HIV-patients. Given that follicular hyperplasia occurs during lymphoid inflammation as a normal response to infection, it is hypothesised that when the diameter of the lymph node follicle (LNF) increases and crosses a critical size, a viral blip occurs due to cryptic viremia. To study this hypothesis, a theoretical analysis of a mathematical model is performed to find the conditions for virus suppression in all compartments and different scenarios of LNF size changes are simulated. According to the analysis, blips with duration of around 30 days arise when the diameter rise rate is between 0.02 and 0.03 days⁻¹. Moreover, the final diameter of the site is directly related to the steady states of the virus load after the occurrence of a blip. When the value of R ₀ is around 2.1, to have a steady-state below the limit of detection after the viral blip, the maximum final diameters should be greater than 0.7 mm so that there is a relative loss of connection between compartments.