Citation: Thomas Schubert, Gernot Längst. Studying epigenetic interactions using MicroScale Thermophoresis (MST)[J]. AIMS Biophysics, 2015, 2(3): 370-380. doi: 10.3934/biophy.2015.3.370
Epigenetics, the inheritable change of gene expression programs without altering the underlying DNA sequence, is most commonly associated with mechanisms that change the packaging of DNA and the landscape of chromatin associated proteins and RNAs. Regulation is correlated with changes in DNA and histone modifications and the specific recognition of these modifications by effector molecules.
Modifications to DNA and histones are dynamically set and removed by chromatin-modifying enzymes in a highly regulated manner. Currently there are four different DNA modifications and a plethory of different histone modifications [1,2,3,4,5]. The DNA and histone code is set by so called writer molecules (methylases, kinases, acetylases and others) that modify histones and DNA. These chromatin marks are specifically recognized by so called reader molecules that were shown to alter chromatin structure by non-covalent interactions with and between nucleosomes. The type and combination of sequestered protein/RNA readers determine the function of the underlying chromatin domain. In a regulated fashion the activity of such chromatin domains is changed by the recruitment of eraser molecules that remove and change local DNA and histone marks (deacetylases, demethylases and others) [6]. Metastable epigenetic regulation relies on a complex re gulatory network specifically binding to modified chromatin by many weak but combinatorial interactions. An understanding of the regulatory processes requires detailed knowledge about the specific binding parameters, with MicroScale Thermophoresis being a novel and versatile method to accelerate their identification.
The epigenetic network is build up by the combinatorial interaction between the recognition domains in Proteins, RNAs and DNAs with their associated modifications. Methods based on high throughput sequencing or mass spectrometry, present a correlative picture of epigenetic coding but fail to provide a dynamic and functional description of the underlying mechanisms. The quantitative description of interactions is absolutely essential to understand the regulatory basis of these networks. On the bi-molecular level binding parameters, such as binding affinities, stoichiometries, thermodynamics and kinetics, have to be unraveled, requiring easy and fast methodologies to study the interactions in a higher throughput.
Currently, researchers in the field of epigenetics use mainly semi-quantitative methods such as electrophoretic mobility shift assays (EMSA) and classical pull down assays. Binding kinetics and binding affinities are measured by surface based Surface Plasmon Resonance (SPR) and labelfree calorimetric approaches that measure the enthalpy of the binding event (Isothermal titration calorimetry, ITC). Fluorescence correlation spectroscopy (FCS) determines binding dependent changes in diffusion times to characterize molecular interactions. FCS measures changes in the translational diffusion, whereas fluorescence anisotropy (FA) is based on measurement of changes in rotational diffusion. MicroScale Thermophoresis (MST) is a rather new technology combining different aspects of the methods mentioned above. The technology provides significant experimental advantages, is easy in setup and requires low amounts of material, as such having the potential to develop into the first choice method to study regulatory/interaction networks in epigenetics.
MicroScale Thermophoresis is based on a physical phenomenon called thermophoresis, which describes the movement of molecules through temperature gradients. This directed movement is dependent on three molecular parameters: i) size; ii) charge; iii) hydration shell. Upon bi-molecular interaction of ligand and target molecule, at least one of these parameters will be altered, resulting in a change of the thermophoretic mobility of the complex, which is used to access binding parameters of the molecular interaction [7,8,9,10,11,12,13].
The changes in thermophoretic movement of one molecule is monitored with an optical system, detecting either the intrinsic fluorescence of tryptophanes within proteins, or the fluorescence signal of a fluorophore, attached to one of the interaction partners. The immobilization-free MicroScale Thermophoresis method records the directed movement of the fluorescent molecules in a µm-sized temperature gradient. The binding parameters are determined by the changing thermophoretic properties of the sample with increasing concentrations of the potential binding partner. The technical setup of the MST device (Monolith NT.115, NanoTemper Technologies GmbH, Munich, Germany) is represented in Figure 1A. The optics focuses in the center of thin glass capillaries, holding 4-6 µl of sample. The microscopic temperature gradient with a diameter of ~50 µm and a temperature difference ΔT of 2-6 °C, is established by an infrared laser, also focusing in the observation window. Upon activation of the laser, thermophoretic depletion or accumulation of molecules can be observed in the region of elevated temperature, which is quantified by the Soret coefficient: ST: chot/ccold = exp(−ST ΔT).
Figure 1B represents a typical MST experiment. The initial fluorescence in the sample is recorded and serves as a base line in the absence of the temperature gradient. After 5 seconds, the IR-laser is activated, which leads to the establishment of the temperature gradient. As a result of this, an initial steep drop of the fluorescence signal—the so-called Temperature- or T-Jump—reflecting the temperature dependence of the fluorophore quantum yield is observed. A slower thermophoresis-driven depletion of fluorophores follows the T-Jump. Once the IR-laser is deactivated, a reverse T-Jump and concomitant back diffusion of fluorescent molecules can be observed.
As mentioned above, the thermophoresis of a molecule is highly sensitive with respect to changes in size, charge and hydration shell. Hence, the binding of a ligand to the molecule will alter the thermophoretic behaviour, which can be used to study equilibrium constants, such as the dissociation constant Kd. To get this information, a serial dilution of the ligand is prepared, mixed with a constant concentration of labelled target molecule, loaded into capillaries and analyzed in the instrument by subsequent scanning of each capillary. The results of a typical binding experiment are illustrated in Figure 1C.
The strengths of the MicroScale Thermophoresis in speed, quality and costs are summarized in Figure 2. The MST system allows a fast and flexible setup and optimization of the assays. Very helpful to obtain high quality data-in short time-with internal quality controls, allowing to detect sticking and aggregation/precipitation effects in real time. Upon detection of these effects, technical conditions and buffers can be easily optimized to ensure ideal binding conditions. Besides the rapid measurements (Kd in 15 min), the immobilization free MicroScale Thermophoresis offers free choice of buffers, with the possibility to measure in lysates and sera. Due to the fact, that the thermophoresis of molecules depends on the charge and hydration shell, besides size, there are no limitations in the size of the measured interaction partners. Additional strengths of MST are the wide dynamic range, from pM to mM, to reveal binding affinities and the low sample consumption required to allow measurements.
The quantification of transcription factor binding affinities is an ideal application of MST. Besides others [14,15,16,17,18], the study of Pham and colleagues [19] on the hematopoietic master transcription factor PU.1 nicely demonstrates the accuracy of MST in the characterization of TF binding site affinities. Binding site selection of the transcription factor PU.1 was analyzed during human macrophage differentiation by ChIP-Seq. Bioinformatic analysis of the ChIP-Seq data, allowed to identify bound Pu.1 sites and to predict binding site motifs with high Pu.1 occupancy in the cells. A score list with binding motifs exhibiting best binding was derived and these sites were used for binding site affinity measurements by MST. Studying the binding affinity of 75 individual DNA sequence motifs revealed an inverse correlation between the ChIP-Seq score list and the Kd values determined in vitro. This demonstrates that motif log-odds scores represent a good measure for motif binding site affinity, implying that binding site selection in vivo is in part driven by the affinity of the transcription factor to its site [19]. This study identified three major categories of consensus sites for PU.1 binding. (i) non-bound sites within inaccessible chromatin domains exhibiting low binding affinity, (ii) PU.1 bound, autonomous binding sites with high affinity which are DNase I inaccessible, (iii) PU.1 bound, autonomous binding sites with intermediate / low affinity which are DNase I accessible. Here binding is likely stabilized by cooperativeness with neighboring transcription factor binding sites [19].
Besides the binding affinity of an interaction, its specificity and selectivity can be determined using the MicroScale Thermophoresis. Using the MST, the study of Schubert and colleagues demonstrated, that the effector protein Df31 (Decondensation factor 31), which plays a role as RNP complex in the decondensation of Drosophila chromatin, is able to bind specifically to Cy5 labelled ssRNA but not to ssDNA [20]. Furthermore, the study identified specific RNAs (snoRNAs), enriched at open, accessible chromatin, which do preferentially interact with Df31. MST competition assays using differentially labelled RNAs revealed the specificity of the interaction, which were later validated by Df31 pull down assays and RT PCR. In addition, MicroScale Thermophoresis revealed an additional interaction surface in Df31, targeting the histones. Df31-EGFP interacts with the core histones H3 and H4, but not with H2A and H2B.
Histone interactions allow the specific binding of Df31 to nucleosomes, as depicted in Figure 3. Mononucleosomes were reconstituted on a 601 positioning sequence carrying a Cy5 label on the DNA and incubated with increasing amounts of Df31 to analyse the binding affinity (Figure 3).
Histone tails are known to be modified post-translationally, with these modifications serving as epigenetic marks. Modifications such as methylation, acetylation, phosphorylation or ubiquitinylation are set and recognized by a large number of proteins-domains that have been named “reader” domains. Lysine methylation at positions 4 and 36 of histone 3 do correlate with transcriptional activity, whereas methylation of Lysine at positions 9 and 27 on the same histone tail correlate with transcriptional repression. The acetylation of lysine residues in histones is known to correlate with accessible chromatin and active transcription [4,21,22,23]. The histone code hypothesis suggests a modification landscape that enables the recruitment of specific effector proteins that in turn determine the chromatin structure and the activity of the underlying DNA. As indicated by several studies, the MST is perfectly suited to characterize protein-histone and protein-histone tail interactions [24,25,26,27].
Alpatov and colleagues studied the role of the fragile X mental retardation protein (FMRP) in the fragile X syndrome [28]. FMRP was identified as a chromatin binding protein that uses its tandem tudor domain for binding. MST experiments using MLA nucleosomes (methyl-lysine analogs) showed that the protein has a preference towards K3Kc79me (Kd = 135 ± 28 nM) and H3Kc27me1 (Kd = 102 ± 11 nM) modified nucleosomes. Unmodified H3 nucleosomes were bound with an affinity of about 1 mM. The study describes an important role of FMRP in DNA damage response (DDR), which is mediated in a chromatin-dependent manner. FMRP regulates genomic stability at the chromatin interface and potentially impacts on gametogenesis and some other developmental aspects of the fragile X syndrome [28].
The study of Josling and colleagues using the malaria parasite Plasmodium falciparum, demonstrates how the parasite specific Bromodomain protein PfBDP1 is involved in the transcriptional regulation of invasion-related genes [29]. Knockdown of the protein results in a drastic inhibition of growth and caused defects in erythrocyte invasion. The study identified 186 downregulated genes affecting different time points of the parasites live cycle. Interestingly, invasion related genes were predominantly found to be downregulated, strengthening the hypothesis that PfBDP1 is involved in the transcriptional regulation these genes. ChIP experiments revealed that the protein binds to 5xintergenic regions of invasion-related genes in schizonts, but not in the ring stages of the parasite. Label-free MST analyses demonstrated that the bromodomain protein PfBDP1 binds specifically to H3K9Ac and H3K14Ac peptides with a Kd = 110.79 ± 10.31 µM (R2 = 0.995) and Kd = 126.17 ± 14.39 µM (R2 = 0.993) respectively. These reported affinities are nicely in line with previously published affinities for bromodomain-histone tail interactions. Furthermore, BDP1 bound to unmodified H3 and H3K16Ac peptides much reduced affinity and did not bind to H3K9me2 and H3K4me2 peptides.
As demonstrated in different studies, the MST can be used in the development of small compound drugs binding to epigenetic proteins [30,31]. Due to the specific role of PfBDP1 in the parasite it presents an interesting drug target. Two known bromodomain inhibitors were tested for binding to the protein in our laboratory at the 2bind GmbH, Regensburg, Germany. The compound PFl-1, which is a selective BET (bromodomain containing protein) inhibitor for BRD4 with a described IC50 of 0.22 mM was shown not to bind to the PfBDP1. However, the compound I-BET151 (GSK1210151A), which is a selective inhibitor for BRD2 (IC50 = 0.5 mM), BRD3 (IC50 = 0.25 mM) and BRD4 (IC50 = 0.79 mM), was identified to bind to PfBDP1 in vitro with an affinity of 260 ± 53 nM (R2 = 0.978) (Figure 4).
The study of Tan and colleagues on the proliferation cell nuclear antigen (PCNA), which takes part in various cellular processes, demonstrates that the MST is a useful tool to characterize the interaction between an epigenetic effector protein and a small compound [30]. Initially, an in silico screen of a compound library using the crystal structure of PCNA revealed a set of 200 compounds that could potentially inhibit PCNA. Two compounds (PCNA-I1 and 2) showed an IC50 value in the low mM range in tumor cell proliferation assays. A subsequent similarity search in the ZINC chemical data base based on the PCNA-I1 compound identified in total 8 additional candidates (PCNA-I3-10). MicroScale Thermophoresis was used to directly determine the binding affinity between recombinant, fluorescently labelled PCNA and the compound PCNA-I1. The Kd of this interaction was determined to be 0.41 ± 0.17 mM. Rabbit IgG was not bound by the compound. Surface plasmon resonance experiments strengthened the MST data, as the affinity was calculated to be 0.17 mM. Native gel electrophoresis indicated that the compound was binding to the trimeric form of PCNA, potentially stabilizing this state. Additionally the compound was shown to release PCNA from chromatin in a dose and time dependent manner. Growth assays using different tumor types and non-transformed cells indicated that the growth of tumor cells is affected at lower compound concentrations than the non-transformed cells. A cell cycle arrest was identified to be the major reason of the growth deficiency upon treatment with the PCNA-inhibitor 1 [30].
Beginning 2014 more than 500 publications using the MicroScale Thermophoresis (MST) could be found in literature, characterizing different kinds of molecular interactions in terms of binding affinity (pM to mM range), stoichiometry and thermodynamics. The MST represents a versatile and precise technology with a broad application range from screening assays, binding assays, competition assays, assays with multiple binding partners, buffer optimization assays to stability assays. The high information content of these assays is widely used in drug development, antibody characterization and aptamer development. This review summarized recent studies demonstrating the suitability of this technology to characterize epigenetic networks. The fast and flexible MST assay setup is strengthened by integrated quality controls detecting sticking effects and aggregation/precipitation effects allowing for fast optimization of reaction conditions to ensure high quality data. Freedom of reaction buffer enables the user to work under previously challenging conditions (sera, lysates). The technology is perfectly suited to study protein-protein, protein-nucleic acids, which are characteristic for interaction within chromatin environment. In addition, the method allows for fast quantification of protein-small compound interactions, making it the perfect tool for the characterization of drugs targeting epigenetic players.
The authors declare no conflict of interest.
[1] |
Rothbart SB, Strahl BD (2014) Interpreting the language of histone and DNA modifications. Biochim Biophys Acta 1839: 627-643. doi: 10.1016/j.bbagrm.2014.03.001
![]() |
[2] |
Baylin SB, Jones PA (2011) A decade of exploring the cancer epigenome-biological and translational implications. Nat Rev Cancer 11: 726-734. doi: 10.1038/nrc3130
![]() |
[3] |
Wu H, Zhang Y (2011) Mechanisms and functions of Tet protein-mediated 5-methylcytosine oxidation. Genes Dev 25: 2436-2452. doi: 10.1101/gad.179184.111
![]() |
[4] |
Kouzarides T (2007) Chromatin modifications and their function. Cell 128: 693-705. doi: 10.1016/j.cell.2007.02.005
![]() |
[5] |
Tan M, Luo H, Lee S, et al. (2011) Identification of 67 histone marks and histone lysine crotonylation as a new type of histone modification. Cell 146: 1016-1028. doi: 10.1016/j.cell.2011.08.008
![]() |
[6] |
Zhang G, Pradhan S (2014) Mammalian epigenetic mechanisms. IUBMB Life 66: 240-256. doi: 10.1002/iub.1264
![]() |
[7] |
Baaske P, Wienken CJ, Reineck P, et al. (2010) Optical thermophoresis for quantifying the buffer dependence of aptamer binding. Angew Chem Int Ed Engl 49: 2238-2241. doi: 10.1002/anie.200903998
![]() |
[8] |
Seidel SA, Dijkman PM, Lea WA, et al. (2013) Microscale thermophoresis quantifies biomolecular interactions under previously challenging conditions. Methods 59: 301-315. doi: 10.1016/j.ymeth.2012.12.005
![]() |
[9] |
Seidel SA, Wienken CJ, Geissler S, et al. (2012) Label-free microscale thermophoresis discriminates sites and affinity of protein-ligand binding. Angew Chem Int Ed Engl 51: 10656-10659. doi: 10.1002/anie.201204268
![]() |
[10] |
Wienken CJ, Baaske P, Rothbauer U, et al. (2010) Protein-binding assays in biological liquids using microscale thermophoresis. Nat Commun 1: 100. doi: 10.1038/ncomms1093
![]() |
[11] |
Jerabek-Willemsen M, André T, Wanner R, et al. (2014) MicroScale Thermophoresis: Interaction analysis and beyond. J Mol Struct 1077: 101-113. doi: 10.1016/j.molstruc.2014.03.009
![]() |
[12] |
Jerabek-Willemsen M, Wienken CJ, Braun D, et al. (2011) Molecular interaction studies using microscale thermophoresis. Assay Drug Dev Technol 9: 342-53. doi: 10.1089/adt.2011.0380
![]() |
[13] | Duhr S, Baaske P (2013) Microscale Thermophoresis: Immobilization-free Binding Assays. GIT Laboratory J. Available from: http://www.laboratory-journal.com/applications/biological-methods/microscale-thermophoresis-immobilization-free-binding-assays. |
[14] |
Diermeier SD, Németh A, Rehli M, et al. (2013) Chromatin-Specific Regulation of Mammalian rDNA Transcription by Clustered TTF-I Binding Sites. PLoS Genet 9: e1003786. doi: 10.1371/journal.pgen.1003786
![]() |
[15] |
Timofeeva OA, Chasovskikh S, Lonskaya I, et al. (2012) Mechanisms of unphosphorylated STAT3 transcription factor binding to DNA. J Biol Chem 287: 14192-14200. doi: 10.1074/jbc.M111.323899
![]() |
[16] |
Gjerstorff MF, Relster MM, Greve KBV, et al. (2014) SSX2 is a novel DNA-binding protein that antagonizes polycomb group body formation and gene repression. Nucleic Acids Res 42: 11433-11446 doi: 10.1093/nar/gku852
![]() |
[17] |
Posse V, Hoberg E, Dierckx A, et al. (2014) The amino terminal extension of mammalian mitochondrial RNA polymerase ensures promoter specific transcription initiation. Nucleic Acids Res 42: 3638-3647. doi: 10.1093/nar/gkt1397
![]() |
[18] |
Silvers R, Saxena K, Kudlinzki D, et al. (2012) Recombinant expression and purification of human TATA binding protein using a chimeric fusion. Protein Expr Purif 85: 142-147. doi: 10.1016/j.pep.2012.07.006
![]() |
[19] |
Pham TH, Minderjahn J, Schmidl C, et al. (2013) Mechanisms of in vivo binding site selection of the hematopoietic master transcription factor PU.1. Nucleic Acids Res 41: 6391-6402. doi: 10.1093/nar/gkt355
![]() |
[20] |
Schubert T, Pusch MC, Diermeier S, et al. (2012) Df31 protein and snoRNAs maintain accessible higher-order structures of chromatin. Mol Cell 48: 434-444. doi: 10.1016/j.molcel.2012.08.021
![]() |
[21] |
Wang Y, Fischle W, Cheung W, et al. (2004) Beyond the double helix: writing and reading the histone code. Novartis Found Symp 259: 3-17; discussion 17-21, 163-169. doi: 10.1002/0470862637.ch2
![]() |
[22] |
Jenuwein T, Allis CD (2001) Translating the histone code. Science 293: 1074-1080. doi: 10.1126/science.1063127
![]() |
[23] |
Taverna SD, Li H, Ruthenburg AJ, et al. (2007) How chromatin-binding modules interpret histone modifications: lessons from professional pocket pickers. Nat Struct Mol Biol 14: 1025-1040. doi: 10.1038/nsmb1338
![]() |
[24] | Southall SM, Cronin NB, Wilson JR (2013) A novel route to product specificity in the Suv4-20 family of histone H4K20 methyltransferases. Nucleic Acids Res 42: 661-671. |
[25] | Guo R, Zheng L, Park J, et al. (2014) BS69/ZMYND11 reads and connects histone h3.3 lysine 36 trimethylation-decorated chromatin to regulated pre-mrna processing. Mol Cell 56: 298-310. |
[26] |
Greer E, Beese-Sims S, Brookes E, et al. (2014) A Histone Methylation Network Regulates Transgenerational Epigenetic Memory in C. elegans. Cell Reports 7: 113-126. doi: 10.1016/j.celrep.2014.02.044
![]() |
[27] |
Allam R, Scherbaum CR, Darisipudi MN, et al. (2012) Histones from Dying Renal Cells Aggravate Kidney Injury via TLR2 and TLR4. J Am Soc Nephrol 23: 1375-1388. doi: 10.1681/ASN.2011111077
![]() |
[28] |
Alpatov R, Lesch B, Nakamoto-Kinoshita M, et al. (2014) A Chromatin-Dependent Role of the Fragile X Mental Retardation Protein FMRP in the DNA Damage Response. Cell 157: 869-881. doi: 10.1016/j.cell.2014.03.040
![]() |
[29] |
Josling GA, Petter M, Oehring SC, et al. (2015) A plasmodium falciparum bromodomain protein regulates invasion gene expression. Cell Host Microbe 17: 741-751. doi: 10.1016/j.chom.2015.05.009
![]() |
[30] |
Tan Z, Wortman M, Dillehay KL, et al. (2012) Small molecule targeting of PCNA chromatin association inhibits tumor cell growth. Mol Pharmacol 81: 811-819. doi: 10.1124/mol.112.077735
![]() |
[31] |
Zheng YC, Duan YC, Ma JL, et al. (2013) Triazole-dithiocarbamate based selective lysine specific demethylase 1 (LSD1) inactivators inhibit gastric cancer cell growth, invasion, and migration. J Med Chem 56: 8543-8560. doi: 10.1021/jm401002r
![]() |
1. | Vladimir B. Teif, Andrey G. Cherstvy, Chromatin and epigenetics: current biophysical views, 2016, 3, 2377-9098, 88, 10.3934/biophy.2016.1.88 | |
2. | Philine Kirsch, Alwin M. Hartman, Anna K. H. Hirsch, Martin Empting, Concepts and Core Principles of Fragment-Based Drug Design, 2019, 24, 1420-3049, 4309, 10.3390/molecules24234309 | |
3. | Anuroop Venkateswaran Venkatasubramani, Katy McLaughlin, Giovanny Rodriguez Blanco, Vladimir Larionov, Alexander Kagansky, Pilot RNAi screening using mammalian cell-based system identifies novel putative silencing factors including Kat5/Tip60, 2015, 2, 2377-9098, 570, 10.3934/biophy.2015.4.570 | |
4. | Aleksandra Kosianova, Vladlena Tiasto, Margarita Yatsunskaya, Yuri Khotimchenko, Alexander Kagansky, Natural molecules as modulators of epigenetic silencing in human cells for cancer care and aging, 2020, 65, 2587-5779, 10.21638/spbu03.2020.405 | |
5. | Yan Yang, Yingai Yin, Sai Wang, Yiyang Dong, Simultaneous determination of zearalenone and ochratoxin A based on microscale thermophoresis assay with a bifunctional aptamer, 2021, 1155, 00032670, 338345, 10.1016/j.aca.2021.338345 | |
6. | Ivan Corbeski, Velten Horn, Ramon A. van der Valk, Ulric B. le Paige, Remus T. Dame, Hugo van Ingen, 2018, Chapter 11, 978-1-4939-8674-3, 177, 10.1007/978-1-4939-8675-0_11 | |
7. | Ghassan Al Hamoui Dit Banni, Rouba Nasreddine, Syntia Fayad, Cyril Colas, Axel Marchal, Reine Nehmé, Investigation of lipase-ligand interactions in porcine pancreatic extracts by microscale thermophoresis, 2021, 1618-2642, 10.1007/s00216-021-03314-7 | |
8. | Sarah Kreuz, Stefan‐Sebastian David, Lorena Viridiana Cortes Medina, Wolfgang Fischle, 2022, 9783527348640, 213, 10.1002/9783527830503.ch11 | |
9. | Austė Kanapeckaitė, Neringa Burokienė, Asta Mažeikienė, Graeme S. Cottrell, Darius Widera, Biophysics is reshaping our perception of the epigenome: from DNA-level to high-throughput studies, 2021, 1, 26670747, 100028, 10.1016/j.bpr.2021.100028 | |
10. | D. R. Martin, A. T. Mutombwera, A. M. Madiehe, M. O. Onani, M. Meyer, R. Cloete, Molecular modeling and simulation studies of SELEX-derived high-affinity DNA aptamers to the Ebola virus nucleoprotein, 2024, 0739-1102, 1, 10.1080/07391102.2024.2302922 | |
11. | Ivan Corbeski, Velten Horn, Ramon A. van der Valk, Ulric le Paige, Remus T. Dame, Hugo van Ingen, 2024, Chapter 17, 978-1-0716-3929-0, 357, 10.1007/978-1-0716-3930-6_17 |