1
|
Chromosomal Arrangement of Phosphorelay Genes Couples Sporulation and DNA Replication. Cell 2015; 162:328-337. [PMID: 26165942 DOI: 10.1016/j.cell.2015.06.012] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2014] [Revised: 03/27/2015] [Accepted: 05/21/2015] [Indexed: 12/21/2022]
Abstract
Genes encoding proteins in a common regulatory network are frequently located close to one another on the chromosome to facilitate co-regulation or couple gene expression to growth rate. Contrasting with these observations, here, we demonstrate a functional role for the arrangement of Bacillus subtilis sporulation network genes on opposite sides of the chromosome. We show that the arrangement of two sporulation network genes, one located close to the origin and the other close to the terminus, leads to a transient gene dosage imbalance during chromosome replication. This imbalance is detected by the sporulation network to produce cell-cycle coordinated pulses of the sporulation master regulator Spo0A∼P. This pulsed response allows cells to decide between sporulation and continued vegetative growth during each cell cycle spent in starvation. The simplicity of this coordination mechanism suggests that it may be widely applicable in a variety of gene regulatory and stress-response settings. VIDEO ABSTRACT.
Collapse
|
2
|
Stekel DJ, Jenkins DJ. Evolution of resource and energy management in biologically realistic gene regulatory network models. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2012; 751:301-28. [PMID: 22821464 DOI: 10.1007/978-1-4614-3567-9_14] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
Abstract
We describe the use of computational models of evolution of artificial gene regulatory networks to understand the topologies of biological gene regulatory networks. We summarize results from three complementary approaches that explicitly represent biological processes of transcription, translation, metabolism and gene regulation: a fine-grained model that allows detailed molecular interactions, a coarse-grained model that allows rapid evolution of many generations, and a fixed-architecture model that allows for comparison of different hypotheses. In the first two cases, we are able to evolve networks towards the biological fitness objectives of survival and reproduction. A theme that emerges is that the control of cell energy and resources is a major driver of gene network topology and function. This is demonstrated in the fine-grained model with the emergence of biologically realistic mRNA and protein turnover rates that optimize energy usage and cell division time, and the evolution of basic repressor activities; in the fixed architecture model with a negative self-regulating gene evolving major efficiencies in mRNA usage; and in the coarse-grained model by the need for the inclusion of basal gene expression to obtain biologically plausible networks and the emergence of global regulators keeping all cellular systems under negative control. In summary, we demonstrate the value of biologically realistic computer evolution techniques, and the importance of energy and resource management in driving the topology and function of gene regulatory networks.
Collapse
Affiliation(s)
- Dov J Stekel
- School of Biosciences, University of Nottingham, Leicestershire, UK.
| | | |
Collapse
|
3
|
Iber D. Inferring Biological Mechanisms by Data-Based Mathematical Modelling: Compartment-Specific Gene Activation during Sporulation in Bacillus subtilis as a Test Case. Adv Bioinformatics 2012; 2011:124062. [PMID: 22312331 PMCID: PMC3270535 DOI: 10.1155/2011/124062] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2011] [Revised: 10/12/2011] [Accepted: 11/03/2011] [Indexed: 11/27/2022] Open
Abstract
Biological functionality arises from the complex interactions of simple components. Emerging behaviour is difficult to recognize with verbal models alone, and mathematical approaches are important. Even few interacting components can give rise to a wide range of different responses, that is, sustained, transient, oscillatory, switch-like responses, depending on the values of the model parameters. A quantitative comparison of model predictions and experiments is therefore important to distinguish between competing hypotheses and to judge whether a certain regulatory behaviour is at all possible and plausible given the observed type and strengths of interactions and the speed of reactions. Here I will review a detailed model for the transcription factor σ(F), a regulator of cell differentiation during sporulation in Bacillus subtilis. I will focus in particular on the type of conclusions that can be drawn from detailed, carefully validated models of biological signaling networks. For most systems, such detailed experimental information is currently not available, but accumulating biochemical data through technical advances are likely to enable the detailed modelling of an increasing number of pathways. A major challenge will be the linking of such detailed models and their integration into a multiscale framework to enable their analysis in a larger biological context.
Collapse
Affiliation(s)
- Dagmar Iber
- Department for Biosystems Science and Engineering, Switzerland and Swiss Institute of Bioinformatics (SIB), ETH Zurich, Mattenstraße 26, Basel 4058, Switzerland
| |
Collapse
|
4
|
Jenkins DJ, Stekel DJ. A new model for investigating the evolution of transcription control networks. ARTIFICIAL LIFE 2009; 15:259-291. [PMID: 19254178 DOI: 10.1162/artl.2009.stekel.006] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
Biological systems show unbounded capacity for complex behaviors and responses to their environments. This principally arises from their genetic networks. The processes governing transcription, translation, and gene regulation are well understood, as are the mechanisms of network evolution, such as gene duplication and horizontal gene transfer. However, the evolved networks arising from these simple processes are much more difficult to understand, and it is difficult to perform experiments on the evolution of these networks in living organisms because of the timescales involved. We propose a new framework for modeling and investigating the evolution of transcription networks in realistic, varied environments. The model we introduce contains novel, important, and lifelike features that allow the evolution of arbitrarily complex transcription networks. Molecular interactions are not specified; instead they are determined dynamically based on shape, allowing protein function to freely evolve. Transcriptional logic provides a flexible mechanism for defining genetic regulatory activity. Simulations demonstrate a realistic life cycle as an emergent property, and that even in simple environments lifelike and complex regulation mechanisms are evolved, including stable proteins, unstable mRNA, and repressor activity. This study also highlights the importance of using in silico genetics techniques to investigate evolved model robustness.
Collapse
Affiliation(s)
- Dafyd J Jenkins
- Centre for Systems Biology, School of Biosciences, University of Birmingham, Edgbaston, UK.
| | | |
Collapse
|
5
|
Igoshin OA, Brody MS, Price CW, Savageau MA. Distinctive topologies of partner-switching signaling networks correlate with their physiological roles. J Mol Biol 2007; 369:1333-52. [PMID: 17498739 PMCID: PMC2727513 DOI: 10.1016/j.jmb.2007.04.021] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2006] [Revised: 04/09/2007] [Accepted: 04/09/2007] [Indexed: 10/23/2022]
Abstract
Regulatory networks controlling bacterial gene expression often evolve from common origins and share homologous proteins and similar network motifs. However, when functioning in different physiological contexts, these motifs may be re-arranged with different topologies that significantly affect network performance. Here we analyze two related signaling networks in the bacterium Bacillus subtilis in order to assess the consequences of their different topologies, with the aim of formulating design principles applicable to other systems. These two networks control the activities of the general stress response factor sigma(B) and the first sporulation-specific factor sigma(F). Both networks have at their core a "partner-switching" mechanism, in which an anti-sigma factor forms alternate complexes either with the sigma factor, holding it inactive, or with an anti-anti-sigma factor, thereby freeing sigma. However, clear differences in network structure are apparent: the anti-sigma factor for sigma(F) forms a long-lived, "dead-end" complex with its anti-anti-sigma factor and ADP, whereas the genes encoding sigma(B) and its network partners lie in a sigma(B)-controlled operon, resulting in positive and negative feedback loops. We constructed mathematical models of both networks and examined which features were critical for the performance of each design. The sigma(F) model predicts that the self-enhancing formation of the dead-end complex transforms the network into a largely irreversible hysteretic switch; the simulations reported here also demonstrate that hysteresis and slow turn off kinetics are the only two system properties associated with this complex formation. By contrast, the sigma(B) model predicts that the positive and negative feedback loops produce graded, reversible behavior with high regulatory capacity and fast response time. Our models demonstrate how alterations in network design result in different system properties that correlate with regulatory demands. These design principles agree with the known or suspected roles of similar networks in diverse bacteria.
Collapse
Affiliation(s)
- Oleg A. Igoshin
- Department of Biomedical Engineering, University of California, Davis, CA, 95616
| | - Margaret S. Brody
- Department of Food Science, University of California, Davis, CA, 95616
| | - Chester W. Price
- Department of Food Science, University of California, Davis, CA, 95616
| | - Michael A. Savageau
- Department of Biomedical Engineering, University of California, Davis, CA, 95616
- Corresponding author: e-mail: ; phone 1(530) 754-8375; fax: 1(530) 7545739
| |
Collapse
|
6
|
Igoshin OA, Price CW, Savageau MA. Signalling network with a bistable hysteretic switch controls developmental activation of the sigma transcription factor in Bacillus subtilis. Mol Microbiol 2006; 61:165-84. [PMID: 16824103 DOI: 10.1111/j.1365-2958.2006.05212.x] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
Abstract
The sporulation process of the bacterium Bacillus subtilis unfolds by means of separate but co-ordinated programmes of gene expression within two unequal cell compartments, the mother cell and the smaller forespore. sigmaF is the first compartment-specific transcription factor activated during this process, and it is controlled at the post-translational level by a partner-switching mechanism that restricts sigmaF activity to the forespore. The crux of this mechanism lies in the ability of the anti-sigma factor SpoIIAB (AB) to form alternative complexes either with sigmaF, holding it in an inactive form, or with the anti-anti-sigma factor SpoIIAA (AA) and a nucleotide, either ATP or ADP. In the complex with AB and ATP, AA is phosphorylated on a serine residue and released, making AB available to capture sigmaF in an inactive complex. Subsequent activation of sigmaF requires the intervention of the SpoIIE serine phosphatase to dephosphorylate AA, which can then attack the AB-sigmaF complex to induce the release of sigmaF. By incorporating biochemical, biophysical and genetic data from the literature we have constructed an integrative mathematical model of this partner-switching network. The model predicts that the self-enhancing formation of a long-lived complex of AA, AB and ADP transforms the network into an essentially irreversible hysteretic switch, thereby explaining the sharp, robust and irreversible activation of sigmaF in the forespore compartment. The model also clarifies the contributions of the partly redundant mechanisms that ensure correct spatial and temporal activation of sigmaF, reproduces the behaviour of various mutants and makes strong, testable predictions.
Collapse
Affiliation(s)
- Oleg A Igoshin
- Department of Biomedical Engineering, One Shields Avenue, University of California, Davis, CA 95616, USA
| | | | | |
Collapse
|
7
|
Iber D, Clarkson J, Yudkin MD, Campbell ID. The mechanism of cell differentiation in Bacillus subtilis. Nature 2006; 441:371-4. [PMID: 16710423 DOI: 10.1038/nature04666] [Citation(s) in RCA: 46] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2006] [Accepted: 02/21/2006] [Indexed: 11/09/2022]
Abstract
Sporulation in Bacillus subtilis serves as a model for the development of two different cell types from a single cell. Although much information has been accumulated about the mechanisms that initiate the developmental programmes, important questions remain that can be answered only by quantitative analysis. Here we develop, with the help of existing and new experimental results, a mathematical model that reproduces published in vitro experiments and explains how the activation of the key transcription factor is regulated. The model identifies the difference in volume between the two cell types as the primary trigger for determining cell fate. It shows that this effect depends on the allosteric behaviour of a key protein kinase and on a low rate of dephosphorylation by the corresponding phosphatase; both predicted effects are confirmed experimentally.
Collapse
Affiliation(s)
- Dagmar Iber
- Mathematical Institute, Centre for Mathematical Biology, University of Oxford, Oxford OX1 3LB, UK.
| | | | | | | |
Collapse
|
8
|
Yudkin MD, Clarkson J. Differential gene expression in genetically identical sister cells: the initiation of sporulation in Bacillus subtilis. Mol Microbiol 2005; 56:578-89. [PMID: 15819616 DOI: 10.1111/j.1365-2958.2005.04594.x] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
Abstract
Early in sporulation, the cell divides asymmetrically to give two sister compartments, a smaller prespore and a larger mother cell. Differential gene expression in these compartments depends on the regulation of the first sporulation-specific sigma factor, sigma(F), which is activated only in the prespore. Regulation relies on the interactions of four proteins -sigma(F), its antisigma SpoIIAB (which also has protein kinase activity), the anti-antisigma SpoIIAA and the protein phosphatase SpoIIE. Before asymmetric division, and in the mother cell after division, sigma(F) is held in an inactive complex with SpoIIAB and ATP; SpoIIAA is in its phosphorylated form. To disrupt the complex so as to liberate sigma(F) in the prespore, dephosphorylated SpoIIAA is needed, and this is made available by SpoIIE. Thereafter, SpoIIAB and SpoIIE are active simultaneously in the prespore, cycling SpoIIAA through phosphorylated and non-phosphorylated forms. This cycle detains SpoIIAB in a state in which it cannot inhibit sigma(F). Results from biophysical techniques, mathematical simulations and enzyme kinetics have now helped to elucidate the dynamics of the protein-protein interactions involved. An understanding of these dynamics largely accounts for the regulation of sigma(F). We show that the system is tuned to be highly efficient in its use of components and extremely economical in conserving ATP.
Collapse
Affiliation(s)
- Michael D Yudkin
- Microbiology Unit, Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX13QU, UK.
| | | |
Collapse
|
9
|
Shu JC, Clarkson J, Yudkin M. Studies of SpoIIAB mutant proteins elucidate the mechanisms that regulate the developmental transcription factor sigmaF in Bacillus subtilis. Biochem J 2005; 384:169-78. [PMID: 15294015 PMCID: PMC1134100 DOI: 10.1042/bj20040923] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
SigmaF, the first compartment-specific sigma factor of sporulation, is regulated by an anti-sigma factor, SpoIIAB (AB) and its antagonist SpoIIAA (AA). AB can bind to sigmaF in the presence of ATP or to AA in the presence of ADP; in addition, AB can phosphorylate AA. The ability of AB to switch between its two binding partners regulates sigmaF. Early in sporulation, AA activates sigmaF by releasing it from its complex with AB. We have previously proposed a reaction scheme for the phosphorylation of AA by AB which accounts for AA's regulatory role. A crucial feature of this scheme is a conformational change in AB that accompanies its switch in binding partner. In the present study, we have studied three AB mutants, all of which have amino-acid replacements in the nucleotide-binding region; AB-E104K (Glu104-->Lys) and AB-T49K (Thr49-->Lys) fail to activate sigmaF, and AB-R105A (Arg105-->Ala) activates it prematurely. We used techniques of enzymology, surface plasmon resonance and fluorescence spectroscopy to analyse the defects in each mutant. AB-E104K was deficient in binding to AA, AB-T49K was deficient in binding to ADP and AB-R105A bound ADP exceptionally strongly. Although the release of sigmaF from all three mutant proteins was impaired, and all three failed to undergo the wild-type conformational change when switching binding partners, the phenotypes of the mutant cells were best accounted for by the properties of the respective AB species in forming complexes with AA and ADP. The behaviour of the mutants enables us to propose convincing mechanisms for the regulation of sigmaF in wild-type bacteria.
Collapse
Affiliation(s)
- Jwu-Ching Shu
- Microbiology Unit, Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K
| | - Joanna Clarkson
- Microbiology Unit, Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K
| | - Michael D. Yudkin
- Microbiology Unit, Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K
- To whom correspondence should be addressed (email )
| |
Collapse
|
10
|
Perry S, Shearer N, Little R, Dixon R. Mutational analysis of the nucleotide-binding domain of the anti-activator NifL. J Mol Biol 2005; 346:935-49. [PMID: 15701508 DOI: 10.1016/j.jmb.2004.12.033] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2004] [Revised: 12/14/2004] [Accepted: 12/15/2004] [Indexed: 11/19/2022]
Abstract
The NifL regulatory protein controls transcription of nitrogen fixation genes in Azotobacter vinelandii by modulating the activity of the transcriptional activator NifA through direct protein-protein interactions. The ability of NifL to integrate the antagonistic signals of redox and nitrogen status is achieved via the involvement of discrete domains in signalling specific environmental cues. NifL senses the redox status via an FAD co-factor located within the amino-terminal PAS domain and responds to the fixed nitrogen status by interaction with the signal transduction protein GlnK, which binds to the C-terminal GHKL domain of NifL. The GHKL domain binds adenosine nucleotides and is similar to the core catalytic domain of the histidine protein kinases. Binding of ADP to this domain increases the inhibitory activity of NifL and the formation of protein complexes with NifA. This inhibition is antagonised by the binding of 2-oxoglutarate, a key metabolic signal of the carbon status, to the amino-terminal GAF domain of NifA. In this study we have examined the properties of three mutations within conserved residues in the GHKL domain of NifL that impair signal transduction. All three mutations decrease the affinity of NifL for ADP significantly, but the mutant proteins exhibit discrete properties. The N419D mutation prevents inhibition of NifA activity by NifL both in vivo and in vitro. In contrast, the G455A and G480A mutations eliminate the redox response, but the mutant proteins retain some sensitivity to the fixed nitrogen status and the ability to interact with the GlnK signal transduction protein. Our data suggest that the absence of the redox switch in the G455A and G480A mutants is a consequence of their inability to override the allosteric effect of 2-oxoglutarate on NifA activity. Overall, these results demonstrate that the binding of adenosine nucleotides to the GHKL domain of NifL plays an important role in counteracting the response of NifA to 2-oxoglutarate, under conditions that are inappropriate for nitrogen fixation.
Collapse
Affiliation(s)
- Susan Perry
- Department of Molecular Microbiology, John Innes Centre, Norwich NR4 7UK, UK
| | | | | | | |
Collapse
|
11
|
Abstract
Differentiation of vegetative Bacillus subtilis into heat resistant spores is initiated by the activation of the key transcription regulator Spo0A through the phosphorelay. Subsequent events depend on the cell compartment-specific action of a series of RNA polymerase sigma factors. Analysis of genes in the Spo0A regulon has helped delineate the mechanisms of axial chromatin formation and asymmetric division. There have been considerable advances in our understanding of critical controls that act to regulate the phosphorelay and to activate the sigma factors.
Collapse
Affiliation(s)
- Patrick J Piggot
- Department of Microbiology and Immunology, Temple University School of Medicine. 3400N. Broad St., Philadelphia, Pennsylvania 19140, USA
| | | |
Collapse
|
12
|
Clarkson J, Campbell ID, Yudkin MD. Efficient regulation of sigmaF, the first sporulation-specific sigma factor in B.subtilis. J Mol Biol 2004; 342:1187-95. [PMID: 15351644 DOI: 10.1016/j.jmb.2004.07.090] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2004] [Revised: 07/23/2004] [Accepted: 07/26/2004] [Indexed: 10/26/2022]
Abstract
Differential gene expression is established in the prespore and mother-cell compartments of Bacillus subtilis through the successive activation of a series of cell-type-specific sigma factors. Crucial to the success of this process is the control of the first prespore-specific sigma factor, sigmaF. sigmaF is regulated by the proteins SpoIIAB, SpoIIAA and SpoIIE. SpoIIAB forms an inhibitory complex with sigmaF, which can be dissociated by interaction with SpoIIAA. During this interaction SpoIIAA is phosphorylated. SpoIIE is a membrane-bound phosphatase that dephosphorylates SpoIIAA, thereby re-activating it. It is not understood how sigmaF is activated specifically in the prespore but not in the mother cell. Here, we use a recently developed fluorescence spectroscopy technique to follow in real time the formation of sigmaF.SpoIIAB complexes and their dissociation by SpoIIAA. We show that complete activation of sigmaF is induced by a tenfold increase in SpoIIE activity. This result demonstrates that relatively small changes in SpoIIE activity, which could arise from asymmetric septation, can achieve the all-or-nothing response in sigmaF activity required by the cell. For long-term sigmaF activation, we find that sustained SpoIIE activity is required to counteract the activity of SpoIIAB. Even though the continual phosphorylation and dephosphorylation of SpoIIAA by these two enzymes will expend some ATP, the formation of SpoIIAA.SpoIIAB.ADP complexes greatly diminishes the rate of the phosphorylation reaction, and thus minimizes the wastage of energy. These features provide a very efficient system for regulating sigmaF.
Collapse
Affiliation(s)
- Joanna Clarkson
- Microbiology Unit, Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, UK
| | | | | |
Collapse
|
13
|
Clarkson J, Campbell ID, Yudkin MD. Physical evidence for the induced release of the Bacillus subtilis transcription factor, sigma(F), from its inhibitory complex. J Mol Biol 2004; 340:203-9. [PMID: 15201047 DOI: 10.1016/j.jmb.2004.04.061] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2004] [Revised: 04/13/2004] [Accepted: 04/16/2004] [Indexed: 11/22/2022]
Abstract
The release of the transcription factor sigma(F) from its inhibitory complex with SpoIIAB is a key regulatory step in the control of sporulation in Bacillus subtilis as it initiates a pattern of differential gene expression in the mother cell and prespore compartments. The sigma(F).SpoIIAB complex is dissociated by the unphosphorylated form of the protein SpoIIAA, the alternative binding partner of SpoIIAB. Here, we employ fluorescence spectroscopy to examine the mechanism by which SpoIIAA acts on the sigma(F).SpoIIAB complex. We constructed a mutant of sigma(F), sigma(F)-W46L, which displayed a reproducible fluorescence response on binding to SpoIIAB. Using this mutant we were able to quantify the amount of sigma(F) bound to SpoIIAB in real time. The results provide physical evidence for the "induced release" mechanism of sigma(F) activation. We demonstrate that SpoIIAA interacts directly with the sigma(F).SpoIIAB complex, greatly decreasing the affinity of SpoIIAB for sigma(F) and thus causing the release of the latter. We also demonstrate that sigma(F) is released before SpoIIAA is phosphorylated and that release occurs on a similar time scale to the binding of SpoIIAA to SpoIIAB.
Collapse
Affiliation(s)
- Joanna Clarkson
- Microbiology Unit, Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, UK
| | | | | |
Collapse
|