Discs in misaligned binary systems

We perform SPH simulations to study precession and changes in alignment between the circumprimary disc and the binary orbit in misaligned binary systems. We find that the precession process can be described by the rigid-disc approximation, where the disc is considered as a rigid body interacting with the binary companion only gravitationally. Precession also causes change in alignment between the rotational axis of the disc and the spin axis of the primary star. This type of alignment is of great important for explaining the origin of spin-orbit misaligned planetary systems. However, we find that the rigid-disc approximation fails to describe changes in alignment between the disc and the binary orbit. This is because the alignment process is a consequence of interactions that involve the fluidity of the disc, such as the tidal interaction and the encounter interaction. Furthermore, simulation results show that there are not only alignment processes, which bring the components towards alignment, but also anti-alignment processes, which tend to misalign the components. The alignment process dominates in systems with misalignment angle near 90 degrees, while the anti-alignment process dominates in systems with the misalignment angle near 0 or 180 degrees. This means that highly misaligned systems will become more aligned but slightly misaligned systems will become more misaligned.Comment: 15 pages, 16 figures, 1 table, accepted for publication in MNRA

Organic geochemistry of the Boltysh impact crater, Ukraine

The Boltysh crater has been know for several decades and was originally drilled in the 1960s - 1980s in a study of economic oil shale deposits. Unfortunately, the cores were not curated and have been lost. However we have recently re-drilled the impact crater and have recovered a near continuous record of ~400m of organic rich sediments deposited in a deep isolated lake which overly the basement rocks spanning a period ~10 Ma. The Boltysh impact crater, centred at 48°54–N and 32°15–E is a complex impact structure formed on the basement rocks of the Ukrainian shield. The age of the impact is 65.17±0.64 Ma [1]. At 24km diameter, the impact is unlikely to have contributed substantially to the worldwide devastation at the end of the Cretaceous. However, the precise age of the Boltysh impact relative to the Chicxulub impact and its location on a stable low lying coastal plain which allowed formation of the postimpact crater lake make it a particularly important locality. After the impact, the crater quickly filled with water, and the crater lake received sediment input from the surrounding land surface for a period >10 Ma [2]. These strata contain a valuable record of Paleogene environmental change in central Europe, and one of very few terrestrial records of the KT event. This preeminent record of the Paleogene of central Europe can help us to answer several related scientific questions. What is the relative age of Boltysh compared with Chicxulub? How long was the hydrothermal system active for after the impact event? How did the devastated area surrounding the crater recover, and how rapid was the recovery? The first sediments to be deposited in the crater lake were a series of relatively thin turbidites, the sediments then become organic rich shales and oil shales. Within the core there is ~400 m of organic rich shales/oil shales spanning a period of ~10 Ma some of which contain macrofossils such as ostracods, fish and plant fossils. Preliminary palynological studies suggest initial sedimentation was slow after the impact followed by more rapid sedimentation through the Late Paleocene. Hydrocarbons extracted from these samples are commonly dominated by terrestrial n-alkanes (Fig 1), Hopanes (including 3-methylhopanes) and steranes are also abundant and indicate the immaturity of the samples. The immaturity of samples is also evident from the abundance of hopenes, sterenes and oleanenes especially in the upper section of the core. In some of the oil shales the hopenes and sterenes are the most abundant hydrocarbons present. There is variation in the distribution of hydrocarbons/biomarkers and palynology throughout the core caused by changing inputs and environmental conditions

A model of estrogen-related gene expression reveals non-linear effects in transcriptional response to tamoxifen

SynthSys is a Centre for Integrative Systems Biology (CISB) funded by BBSRC and EPSRC, reference BB/D019621/1.Background: Estrogen receptors alpha (ER) are implicated in many types of female cancers, and are the common target for anti-cancer therapy using selective estrogen receptor modulators (SERMs, such as tamoxifen). However, cell-type specific and patient-to-patient variability in response to SERMs (from suppression to stimulation of cancer growth), as well as frequent emergence of drug resistance, represents a serious problem. The molecular processes behind mixed effects of SERMs remain poorly understood, and this strongly motivates application of systems approaches. In this work, we aimed to establish a mathematical model of ER-dependent gene expression to explore potential mechanisms underlying the variable actions of SERMs. Results: We developed an equilibrium model of ER binding with 17 beta-estradiol, tamoxifen and DNA, and linked it to a simple ODE model of ER-induced gene expression. The model was parameterised on the broad range of literature available experimental data, and provided a plausible mechanistic explanation for the dual agonism/antagonism action of tamoxifen in the reference cell line used for model calibration. To extend our conclusions to other cell types we ran global sensitivity analysis and explored model behaviour in the wide range of biologically plausible parameter values, including those found in cancer cells. Our findings suggest that transcriptional response to tamoxifen is controlled in a complex non-linear way by several key parameters, including ER expression level, hormone concentration, amount of ER-responsive genes and the capacity of ER-tamoxifen complexes to stimulate transcription (e. g. by recruiting co-regulators of transcription). The model revealed non-monotonic dependence of ER-induced transcriptional response on the expression level of ER, that was confirmed experimentally in four variants of the MCF-7 breast cancer cell line. Conclusions: We established a minimal mechanistic model of ER-dependent gene expression, that predicts complex non-linear effects in transcriptional response to tamoxifen in the broad range of biologically plausible parameter values. Our findings suggest that the outcome of a SERM's action is defined by several key components of cellular micro-environment, that may contribute to cell-type-specific effects of SERMs and justify the need for the development of combinatorial biomarkers for more accurate prediction of the efficacy of SERMs in specific cell types.Publisher PDFPeer reviewe

Optimizing the ensemble for equilibration in broad-histogram Monte Carlo simulations

We present an adaptive algorithm which optimizes the statistical-mechanical ensemble in a generalized broad-histogram Monte Carlo simulation to maximize the system's rate of round trips in total energy. The scaling of the mean round-trip time from the ground state to the maximum entropy state for this local-update method is found to be O([N log N]^2) for both the ferromagnetic and the fully frustrated 2D Ising model with N spins. Our new algorithm thereby substantially outperforms flat-histogram methods such as the Wang-Landau algorithm.Comment: 6 pages, 5 figure

Meteorology of Jupiter's Equatorial Hot Spots and Plumes from Cassini

We present an updated analysis of Jupiter's equatorial meteorology from Cassini observations. For two months preceding the spacecraft's closest approach, the Imaging Science Subsystem (ISS) onboard regularly imaged the atmosphere. We created time-lapse movies from this period in order to analyze the dynamics of equatorial hot spots and their interactions with adjacent latitudes. Hot spots are quasi-stable, rectangular dark areas on visible-wavelength images, with defined eastern edges that sharply contrast with surrounding clouds, but diffuse western edges serving as nebulous boundaries with adjacent equatorial plumes. Hot spots exhibit significant variations in size and shape over timescales of days and weeks. Some of these changes correspond with passing vortex systems from adjacent latitudes interacting with hot spots. Strong anticyclonic gyres present to the south and southeast of the dark areas appear to circulate into hot spots. Impressive, bright white plumes occupy spaces in between hot spots. Compact cirrus-like 'scooter' clouds flow rapidly through the plumes before disappearing within the dark areas. These clouds travel at 150-200 m/s, much faster than the 100 m/s hot spot and plume drift speed. This raises the possibility that the scooter clouds may be more illustrative of the actual jet stream speed at these latitudes. Most previously published zonal wind profiles represent the drift speed of the hot spots at their latitude from pattern matching of the entire longitudinal image strip. If a downward branch of an equatorially-trapped Rossby waves controls the overall appearance of hot spots, however, the westward phase velocity of the wave leads to underestimates of the true jet stream speed.Comment: 33 pages, 11 figures; accepted for publication in Icarus; for supplementary movies, please contact autho

Kinetic modelling of in vitro data of PI3K, mTOR1, PTEN enzymes and on-target inhibitors Rapamycin, BEZ235, and LY294002

The phosphatidylinositide 3-kinases (PI3K) and mammalian target of rapamycin-1 (mTOR1) are two key targets for anti-cancer therapy. Predicting the response of the PI3K/AKT/mTOR1 signalling pathway to targeted therapy is made difficult because of network complexities. Systems biology models can help explore those complexities but the value of such models is dependent on accurate parameterisation. Motivated by a need to increase accuracy in kinetic parameter estimation, and therefore the predictive power of the model, we present a framework to integrate kinetic data from enzyme assays into a unified enzyme kinetic model. We present exemplar kinetic models of PI3K and mTOR1, calibrated on in vitro enzyme data and founded on Michaelis-Menten (MM) approximation. We describe the effects of an allosteric mTOR1 inhibitor (Rapamycin) and ATP-competitive inhibitors (BEZ2235 and LY294002) that show dual inhibition of mTOR1 and PI3K. We also model the kinetics of phosphatase and tensin homolog (PTEN), which modulates sensitivity of the PI3K/AKT/mTOR1 pathway to these drugs. Model validation with independent data sets allows investigation of enzyme function and drug dose dependencies in a wide range of experimental conditions. Modelling of the mTOR1 kinetics showed that Rapamycin has an IC50 independent of ATP concentration and that it is a selective inhibitor of mTOR1 substrates S6K1 and 4EBP1: it retains 40% of mTOR1 activity relative to 4EBP1 phosphorylation and inhibits completely S6K1 activity. For the dual ATP-competitive inhibitors of mTOR1 and PI3K, LY294002 and BEZ235, we derived the dependence of the IC50 on ATP concentration that allows prediction of the IC50 at different ATP concentrations in enzyme and cellular assays. Comparison of the drug effectiveness in enzyme and cellular assays showed that some features of these drugs arise from signalling modulation beyond the on-target action and MM approximation and require a systems-level consideration of the whole PI3K/PTEN/AKT/mTOR1 network in order to understand mechanisms of drug sensitivity and resistance in different cancer cell lines. We suggest that using these models in systems biology investigation of the PI3K/AKT/mTOR1 signalling in cancer cells can bridge the gap between direct drug target action and the therapeutic response to these drugs and their combinations

Customizing the therapeutic response of signaling networks to promote antitumor responses by drug combinations

Drug resistance, de novo and acquired, pervades cellular signaling networks (SNs) from one signaling motif to another as a result of cancer progression and/or drug intervention. This resistance is one of the key determinants of efficacy in targeted anti-cancer drug therapy. Although poorly understood, drug resistance is already being addressed in combination therapy by selecting drug targets where SN sensitivity increases due to combination components or as a result of de novo or acquired mutations. Additionally, successive drug combinations have shown low resistance potential. To promote a rational, systematic development of combination therapies, it is necessary to establish the underlying mechanisms that drive the advantages of combination therapies, and design methods to determine drug targets for combination regimens. Based on a joint systems analysis of cellular SN response and its sensitivity to drug action and oncogenic mutations, we describe an in silico method to analyze the targets of drug combinations. Our method explores mechanisms of sensitizing the SN through a combination of two drugs targeting vertical signaling pathways. We propose a paradigm of SN response customization by one drug to both maximize the effect of another drug in combination and promote a robust therapeutic response against oncogenic mutations. The method was applied to customize the response of the ErbB/PI3K/PTEN/AKT pathway by combination of drugs targeting HER2 receptors and proteins in the down-stream pathway. The results of a computational experiment showed that the modification of the SN response from hyperbolic to smooth sigmoid response by manipulation of two drugs in combination leads to greater robustness in therapeutic response against oncogenic mutations determining cancer heterogeneity. The application of this method in drug combination co-development suggests a combined evaluation of inhibition effects together with the capability of drug combinations to suppress resistance mechanisms before they become clinically manifest

Quantitative analysis of NRF2 pathway reveals key elements of the regulatory circuits underlying antioxidant response and proliferation of ovarian cancer cells

Cells are constantly exposed to Reactive Oxygen Species (ROS) produced both endogenously to meet physiological requirements and from exogenous sources. While endogenous ROS are considered as important signalling molecules, high uncontrollable ROS are detrimental. It is unclear how cells can achieve a balance between maintaining physiological redox homeostasis and robustly activate the antioxidant system to remove exogenous ROS. We have utilised a Systems Biology approach to understand how this robust adaptive system fulfils homeostatic requirements of maintaining steady-state ROS and growth rate, while undergoing rapid readjustment under challenged conditions. Using a panel of human ovarian and normal cell lines, we experimentally quantified and established interrelationships between key elements of ROS homeostasis. The basal levels of NRF2 and KEAP1 were cell line specific and maintained in tight correlation with their growth rates and ROS. Furthermore, perturbation of this balance triggered cell specific kinetics of NRF2 nuclear–cytoplasmic relocalisation and sequestration of exogenous ROS. Our experimental data were employed to parameterise a mathematical model of the NRF2 pathway that elucidated key response mechanisms of redox regulation and showed that the dynamics of NRF2-H2O2 regulation defines a relationship between half-life, total and nuclear NRF2 level and endogenous H2O2 that is cell line specific