· · · Institutional Repository

Background: Humans normally dissipate significant energy during walking, largely at the transitions between steps. The ankle then acts to restore energy during push-off, which may be the reason that ankle impairment nearly always leads to poorer walking economy. The replacement of lost energy is necessary for steady gait, in which mechanical energy is constant on average, external dissipation is negligible, and no net work is performed over a stride. However, dissipation and replacement by muscles might not be necessary if energy were instead captured and reused by an assistive device. Methodology/Principal Findings: We developed a microprocessor-controlled artificial foot that captures some of the energy that is normally dissipated by the leg and ‘‘recycles’’ it as positive ankle work. In tests on subjects walking with an artificially-impaired ankle, a conventional prosthesis reduced ankle push-off work and increased net metabolic energy expenditure by 23% compared to normal walking. Energy recycling restored ankle push-off to normal and reduced the net metabolic energy penalty to 14%. Conclusions/Significance: These results suggest that reduced ankle push-off contributes to the increased metabolic energy expenditure accompanying ankle impairments, and demonstrate that energy recycling can be used to reduce such cost.

Background: The authors present a procedural extension of the popular Implicit Association Test (IAT; [1]) that allows for indirect measurement of attitudes on multiple dimensions (e.g., safe–unsafe; young–old; innovative–conventional, etc.) rather than on a single evaluative dimension only (e.g., good–bad). Methodology/Principal Findings: In two within-subjects studies, attitudes toward three automobile brands were measured on six attribute dimensions. Emphasis was placed on evaluating the methodological appropriateness of the new procedure, providing strong evidence for its reliability, validity, and sensitivity. Conclusions/Significance: This new procedure yields detailed information on the multifaceted nature of brand associations that can add up to a more abstract overall attitude. Just as the IAT, its multi-dimensional extension/application (dubbed md-IAT) is suited for reliably measuring attitudes consumers may not be consciously aware of, able to express, or willing to share with the researcher [2,3].

Binding selectivity and cross-reactivity within one of the largest and most abundant interaction domain families, the PDZ family, has long been enigmatic. The complete human PDZ domain complement (the PDZome) consists of 267 domains and we applied here a Bayesian selectivity model to predict hundreds of human PDZ domain interactions, using target sequences of 22,997 non-redundant proteins. Subsequent analysis of these binding scores shows that PDZs can be divided into two genome-wide clusters that coincide well with the division between canonical class 1 and 2 PDZs. Within the class 1 PDZs we observed binding overlap at unprecedented levels, mediated by two residues at positions 1 and 5 of the second ahelix of the binding pocket. Eight PDZ domains were subsequently selected for experimental binding studies and to verify the basics of our predictions. Overall, the PDZ domain class 1 cross-reactivity identified here implies that auxiliary mechanisms must be in place to overcome this inherent functional overlap and to minimize cross-selectivity within the living cell. Indeed, when we superimpose PDZ domain binding affinities with gene ontologies, network topology data and the domain position within a PDZ superfamily protein, functional overlap is minimized and PDZ domains position optimally in the binding space. We therefore propose that PDZ domain selectivity is achieved through cellular context rather than inherent binding specificity.

Specific interactions of the genome with the nuclear lamina (NL) are thought to assist chromosome folding inside the nucleus and to contribute to the regulation of gene expression. High-resolution mapping has recently identified hundreds of large, sharply defined lamina-associated domains (LADs) in the human genome, and suggested that the insulator protein CTCF may help to demarcate these domains. Here, we report the detailed structure of LADs in Drosophila cells, and investigate the putative roles of five insulator proteins in LAD organization. We found that the Drosophila genome is also organized in discrete LADs, which are about five times smaller than human LADs but contain on average a similar number of genes. Systematic comparison to new and published insulator binding maps shows that only SU(HW) binds preferentially at LAD borders and at specific positions inside LADs, while GAF, CTCF, BEAF-32 and DWG are mostly absent from these regions. By knockdown and overexpression studies we demonstrate that SU(HW) weakens genome – NL interactions through a local antagonistic effect, but we did not obtain evidence that it is essential for border formation. Our results provide insights into the evolution of LAD organization and identify SU(HW) as a fine-tuner of genome – NL interactions.

Background: To support the development of early warning and surveillance systems of emerging zoonoses, we present a general method to prioritize pathogens using a quantitative, stochastic multi-criteria model, parameterized for the Netherlands. Methodology/Principal Findings: A risk score was based on seven criteria, reflecting assessments of the epidemiology and impact of these pathogens on society. Criteria were weighed, based on the preferences of a panel of judges with a background in infectious disease control. Conclusions/Significance: Pathogens with the highest risk for the Netherlands included pathogens in the livestock reservoir with a high actual human disease burden (e.g. Campylobacter spp., Toxoplasma gondii, Coxiella burnetii) or a low current but higher historic burden (e.g. Mycobacterium bovis), rare zoonotic pathogens in domestic animals with severe disease manifestations in humans (e.g. BSE prion, Capnocytophaga canimorsus) as well as arthropod-borne and wildlife associated pathogens which may pose a severe risk in future (e.g. Japanese encephalitis virus and West-Nile virus). These agents are key targets for development of early warning and surveillance.

Tumorigenesis is a multi-step process in which normal cells transform into malignant tumors following the accumulation of genetic mutations that enable them to evade the growth control checkpoints that would normally suppress their growth or result in apoptosis. It is therefore important to identify those combinations of mutations that collaborate in cancer development and progression. DNA copy number alterations (CNAs) are one of the ways in which cancer genes are deregulated in tumor cells. We hypothesized that synergistic interactions between cancer genes might be identified by looking for regions of co-occurring gain and/or loss. To this end we developed a scoring framework to separate truly cooccurring aberrations from passenger mutations and dominant single signals present in the data. The resulting regions of high co-occurrence can be investigated for between-region functional interactions. Analysis of high-resolution DNA copy number data from a panel of 95 hematological tumor cell lines correctly identified co-occurring recombinations at the T-cell receptor and immunoglobulin loci in T- and B-cell malignancies, respectively, showing that we can recover truly cooccurring genomic alterations. In addition, our analysis revealed networks of co-occurring genomic losses and gains that are enriched for cancer genes. These networks are also highly enriched for functional relationships between genes. We further examine sub-networks of these networks, core networks, which contain many known cancer genes. The core network for co-occurring DNA losses we find seems to be independent of the canonical cancer genes within the network. Our findings suggest that large-scale, low-intensity copy number alterations may be an important feature of cancer development or maintenance by affecting gene dosage of a large interconnected network of functionally related genes.

The cell envelope of mycobacteria, a group of Gram positive bacteria, is composed of a plasma membrane and a Gramnegative-like outer membrane containing mycolic acids. In addition, the surface of the mycobacteria is coated with an illcharacterized layer of extractable, non-covalently linked glycans, lipids and proteins, collectively known as the capsule, whose occurrence is a matter of debate. By using plunge freezing cryo-electron microscopy technique, we were able to show that pathogenic mycobacteria produce a thick capsule, only present when the cells were grown under unperturbed conditions and easily removed by mild detergents. This detergent-labile capsule layer contains arabinomannan, a-glucan and oligomannosyl-capped glycolipids. Further immunogenic and proteomic analyses revealed that Mycobacterium marinum capsule contains high amounts of proteins that are secreted via the ESX-1 pathway. Finally, cell infection experiments demonstrated the importance of the capsule for binding to cells and dampening of pro-inflammatory cytokine response. Together, these results show a direct visualization of the mycobacterial capsular layer as a labile structure that contains ESX-1-secreted proteins.

Determination of an accurate glenohumeral-joint rotation center (GH-JRC) from marker data is essential for kinematic and dynamic analysis of shoulder motions. Previous studies have focused on the evaluation of the different functional methods for the estimation of the GH-JRC for healthy subjects. The goal of this paper is to compare two widely used functional methods, namely the instantaneous helical axis (IHA) and symmetrical center of rotation (SCoRE) methods, for estimating the GH-JRC in vivo for patients with implanted shoulder hemiarthroplasty. The motion data of five patients were recorded while performing three different dynamic motions (circumduction, abduction, and forward flexion). The GH-JRC was determined using the CT-images of the subjects (geometric GH-JRC) and was also estimated using the two IHA and SCoRE methods. The rotation centers determined using the IHA and SCoRE methods were on average 1.4760.62 cm and 2.0760.55 cm away from geometric GH-JRC, respectively. The two methods differed significantly (two-tailed p-value from paired t-Test ,0.02, post-hoc power ,0.30). The SCoRE method showed a significant lower (two-tailed p-value from paired t-Test ,0.03, post-hoc power ,0.68) repeatability error calculated between the different trials of each motion and each subject and averaged across all measured subjects (0.6260.10 cm for IHA vs. 0.4360.12 cm for SCoRE). It is concluded that the SCoRE appeared to be a more repeatable method whereas the IHA method resulted in a more accurate estimation of the GH-JRC for patients with endoprostheses.

Background: Articular cartilage has very limited intrinsic regenerative capacity, making cell-based therapy a tempting approach for cartilage repair. Cell tracking can be a major step towards unraveling and improving the repair process of these therapies. We studied superparamagnetic iron oxides (SPIO) for labeling human bone marrow-derived mesenchymal stem cells (hBMSCs) regarding effectivity, cell viability, long term metabolic cell activity, chondrogenic differentiation and hBMSC secretion profile. We additionally examined the capacity of synovial cells to endocytose SPIO from dead, labeled cells, together with the use of magnetic resonance imaging (MRI) for intra-articular visualization and quantification of SPIO labeled cells. Methodology/Prinicipal Findings: Efficacy and various safety aspects of SPIO cell labeling were determined using appropriate assays. Synovial SPIO re-uptake was investigated in vitro by co-labeling cells with SPIO and green fluorescent protein (GFP). MRI experiments were performed on a clinical 3.0T MRI scanner. Two cell-based cartilage repair techniques were mimicked for evaluating MRI traceability of labeled cells: intra-articular cell injection and cell implantation in cartilage defects. Cells were applied ex vivo or in vitro in an intra-articular environment and immediately scanned. SPIO labeling was effective and did not impair any of the studied safety aspects, including hBMSC secretion profile. SPIO from dead, labeled cells could be taken up by synovial cells. Both injected and implanted SPIO-labeled cells could accurately be visualized by MRI in a clinically relevant sized joint model using clinically applied cell doses. Finally, we quantified the amount of labeled cells seeded in cartilage defects using MR-based relaxometry. Conclusions: SPIO labeling appears to be safe without influencing cell behavior. SPIO labeled cells can be visualized in an intra-articular environment and quantified when seeded in cartilage defects.

Psoriasis is characterized by hyperproliferation of keratinocytes and by infiltration of activated Th1 and Th17 cells in the (epi)dermis. By expression microarray, we previously found the GATA3 transcription factor significantly downregulated in lesional psoriatic skin. Since GATA3 serves as a key switch in both epidermal and T helper cell differentiation, we investigated its function in psoriasis. Because psoriatic skin inflammation shares many characteristics of epidermal regeneration during wound healing, we also studied GATA3 expression under such conditions. Psoriatic lesional skin showed decreased GATA3 mRNA and protein expression compared to non-lesional skin. GATA3 expression was also markedly decreased in inflamed skin of mice with a psoriasiform dermatitis induced with imiquimod. Tape-stripping of nonlesional skin of patients with psoriasis, a standardized psoriasis-triggering and skin regeneration-inducing technique, reduced the expression of GATA3. In wounded skin of mice, low GATA3 mRNA and protein expression was detected. Taken together, GATA3 expression is downregulated under regenerative and inflammatory hyperproliferative skin conditions. GATA3 expression could be re-induced by successful narrow-band UVB treatment of both human psoriasis and imiquimodinduced psoriasiform dermatitis in mice. The prototypic Th2 cytokine IL-4 was the only cytokine capable of inducing GATA3 in skin explants from healthy donors. Based on these findings we argue that GATA3 serves as a key regulator in psoriatic inflammation, keratinocyte hyperproliferation and skin barrier dysfunction.

In the study of metabolic networks, optimization techniques are often used to predict flux distributions, and hence, metabolic phenotype. Flux balance analysis in particular has been successful in predicting metabolic phenotypes. However, an inherent limitation of a stoichiometric approach such as flux balance analysis is that it can predict only flux distributions that result in maximal yields. Hence, previous attempts to use FBA to predict metabolic fluxes in Lactobacillus plantarum failed, as this lactic acid bacterium produces lactate, even under glucose-limited chemostat conditions, where FBA predicted mixed acid fermentation as an alternative pathway leading to a higher yield. In this study we tested, however, whether longterm adaptation on an unusual and poor carbon source (for this bacterium) would select for mutants with optimal biomass yields. We have therefore adapted Lactobacillus plantarum to grow well on glycerol as its main growth substrate. After prolonged serial dilutions, the growth yield and corresponding fluxes were compared to in silico predictions. Surprisingly, the organism still produced mainly lactate, which was corroborated by FBA to indeed be optimal. To understand these results, constraint-based elementary flux mode analysis was developed that predicted 3 out of 2669 possible flux modes to be optimal under the experimental conditions. These optimal pathways corresponded very closely to the experimentally observed fluxes and explained lactate formation as the result of competition for oxygen by the other flux modes. Hence, these results provide thorough understanding of adaptive evolution, allowing in silico predictions of the resulting flux states, provided that the selective growth conditions favor yield optimization as the winning strategy.

The genomes of four Lactobacillus reuteri strains isolated from human breast milk and the gastrointestinal tract have been recently sequenced as part of the Human Microbiome Project. Preliminary genome comparisons suggested that these strains belong to two different clades, previously shown to differ with respect to antimicrobial production, biofilm formation, and immunomodulation. To explain possible mechanisms of survival in the host and probiosis, we completed a detailed genomic comparison of two breast milk–derived isolates representative of each group: an established probiotic strain (L. reuteri ATCC 55730) and a strain with promising probiotic features (L. reuteri ATCC PTA 6475). Transcriptomes of L. reuteri strains in different growth phases were monitored using strain-specific microarrays, and compared using a panmetabolic model representing all known metabolic reactions present in these strains. Both strains contained candidate genes involved in the survival and persistence in the gut such as mucus-binding proteins and enzymes scavenging reactive oxygen species. A large operon predicted to encode the synthesis of an exopolysaccharide was identified in strain 55730. Both strains were predicted to produce health-promoting factors, including antimicrobial agents and vitamins (folate, vitamin B12). Additionally, a complete pathway for thiamine biosynthesis was predicted in strain 55730 for the first time in this species. Candidate genes responsible for immunomodulatory properties of each strain were identified by transcriptomic comparisons. The production of bioactive metabolites by human-derived probiotics may be predicted using metabolic modeling and transcriptomics. Such strategies may facilitate selection and optimization of probiotics for health promotion, disease prevention and amelioration.

In recent years increasing evidence appeared that breast cancer may not constitute a single disease at the molecular level, but comprises a heterogeneous set of subtypes. This suggests that instead of building a single monolithic predictor, better predictors might be constructed that solely target samples of a designated subtype, which are believed to represent more homogeneous sets of samples. An unavoidable drawback of developing subtype-specific predictors, however, is that a stratification by subtype drastically reduces the number of samples available for their construction. As numerous studies have indicated sample size to be an important factor in predictor construction, it is therefore questionable whether the potential benefit of subtyping can outweigh the drawback of a severe loss in sample size. Factors like unequal class distributions and differences in the number of samples per subtype, further complicate comparisons. We present a novel experimental protocol that facilitates a comprehensive comparison between subtype-specific predictors and predictors that do not take subtype information into account. Emphasis lies on careful control of sample size as well as class and subtype distributions. The methodology is applied to a large breast cancer compendium involving over 1500 arrays, using a state-of-the-art subtyping scheme. We show that the resulting subtype-specific predictors outperform those that do not take subtype information into account, especially when taking sample size considerations into account.

Many theories of driver behaviour suggest that unconscious or implicit emotions play a functional role in the shaping and control of behaviour. This has not been experimentally tested however. Therefore, in this study the effects of emotive masked images on driver behaviour were examined. While driving a simulator, participants were repeatedly exposed to negative or neutral emotionally laden target images that were sandwich masked by emotionally neutral images. These images were encountered across two different trials each of which consisted of 3–4 minutes of driving on a rural road. The results indicate an effect of the negative target images primarily in reducing the extent of familiarisation occurring between the first and second experimental drives. This is evident in a reduced decrease in heart rate and a reduced increase in high band heart rate variability and actual travelling speed from the first to second drives if the negative target image was presented in the second drive. In addition to these findings there was no clear effect of the target image on subjective ratings of effort or feelings of risk. There was however an effect of gender, with the majority of the effects found in the study being limited to the larger female dataset. These findings suggest that unconscious or implicit emotional stimuli may well influence driver behaviour without explicit awareness.

Mariprofundus ferrooxydans PV-1 has provided the first genome of the recently discovered Zetaproteobacteria subdivision. Genome analysis reveals a complete TCA cycle, the ability to fix CO2, carbon-storage proteins and a sugar phosphotransferase system (PTS). The latter could facilitate the transport of carbohydrates across the cell membrane and possibly aid in stalk formation, a matrix composed of exopolymers and/or exopolysaccharides, which is used to store oxidized iron minerals outside the cell. Two-component signal transduction system genes, including histidine kinases, GGDEF domain genes, and response regulators containing CheY-like receivers, are abundant and widely distributed across the genome. Most of these are located in close proximity to genes required for cell division, phosphate uptake and transport, exopolymer and heavy metal secretion, flagellar biosynthesis and pilus assembly suggesting that these functions are highly regulated. Similar to many other motile, microaerophilic bacteria, genes encoding aerotaxis as well as antioxidant functionality (e.g., superoxide dismutases and peroxidases) are predicted to sense and respond to oxygen gradients, as would be required to maintain cellular redox balance in the specialized habitat where M. ferrooxydans resides. Comparative genomics with other Fe(II) oxidizing bacteria residing in freshwater and marine environments revealed similar content, synteny, and amino acid similarity of coding sequences potentially involved in Fe(II) oxidation, signal transduction and response regulation, oxygen sensation and detoxification, and heavy metal resistance. This study has provided novel insights into the molecular nature of Zetaproteobacteria.

Background Magnetic resonance imaging (MRI), together with histology, is widely used to diagnose and to monitor treatment in oncology. Spatial correspondence between these modalities provides information about the ability of MRI to characterize cancerous tissue. However, registration is complicated by deformations during pathological processing, and differences in scale and information content. Methodology/Principal Findings This study proposes a methodology for establishing an accurate 3D relation between histological sections and high resolution in vivo MRI tumor data. The key features of the methodology are: 1) standardized acquisition and processing, 2) use of an intermediate ex vivo MRI, 3) use of a reference cutting plane, 4) dense histological sampling, 5) elastic registration, and 6) use of complete 3D data sets. Five rat pancreatic tumors imaged by T2*-w MRI were used to evaluate the proposed methodology. The registration accuracy was assessed by root mean squared (RMS) distances between manually annotated landmark points in both modalities. After elastic registration the average RMS distance decreased from 1.4 to 0.7 mm. The intermediate ex vivo MRI and the reference cutting plane shared by all three 3D images (in vivo MRI, ex vivo MRI, and 3D histology data) were found to be crucial for the accurate co-registration between the 3D histological data set and in vivo MRI. The MR intensity in necrotic regions, as manually annotated in 3D histology, was significantly different from other histologically confirmed regions (i.e., viable and hemorrhagic). However, the viable and the hemorrhagic regions showed a large overlap in T2*-w MRI signal intensity. Conclusions The established 3D correspondence between tumor histology and in vivo MRI enables extraction of MRI characteristics for histologically confirmed regions. The proposed methodology allows the creation of a tumor database of spatially registered multi-spectral MR images and multi-stained 3D histology.

After oral exposure, prions are thought to enter Peyer's patches via M cells and accumulate first upon follicular dendritic cells (FDCs) before spreading to the nervous system. How prions are actually initially acquired from the gut lumen is not known. Using high-resolution immunofluorescence and cryo-immunogold electron microscopy, we report the trafficking of the prion protein (PrP) toward Peyer's patches of wild-type and PrP-deficient mice. PrP was transiently detectable at 1 day post feeding (dpf) within large multivesicular LAMP1-positive endosomes of enterocytes in the follicle-associated epithelium (FAE) and at much lower levels within M cells. Subsequently, PrP was detected on vesicles in the late endosomal compartments of macrophages in the subepithelial dome. At 7–21 dpf, increased PrP labelling was observed on the plasma membranes of FDCs in germinal centres of Peyer's patches from wild-type mice only, identifying FDCs as the first sites of PrP conversion and replication. Detection of PrP on extracellular vesicles displaying FAE enterocyte-derived A33 protein implied transport towards FDCs in association with FAE-derived vesicles. By 21 dpf, PrP was observed on the plasma membranes of neurons within neighbouring myenteric plexi. Together, these data identify a novel potential M cell-independent mechanism for prion transport, mediated by FAE enterocytes, which acts to initiate conversion and replication upon FDCs and subsequent infection of enteric nerves.

The Leiden Longevity Study consists of families that express extended survival across generations, decreased morbidity in middle-age, and beneficial metabolic profiles. To identify which pathways drive this complex phenotype of familial longevity and healthy aging, we performed a genome-wide gene expression study within this cohort to screen for mRNAs whose expression changes with age and associates with longevity. We first compared gene expression profiles from whole blood samples between 50 nonagenarians and 50 middle-aged controls, resulting in identification of 2,953 probes that associated with age. Next, we determined which of these probes associated with longevity by comparing the offspring of the nonagenarians (50 subjects) and the middle-aged controls. The expression of 360 probes was found to change differentially with age in members of the long-lived families. In a RT-qPCR replication experiment utilizing 312 controls, 332 offspring and 79 nonagenarians, we confirmed a nonagenarian specific expression profile for 21 genes out of 25 tested. Since only some of the offspring will have inherited the beneficial longevity profile from their long-lived parents, the contrast between offspring and controls is expected to be weak. Despite this dilution of the longevity effects, reduced expression levels of two genes, ASF1A and IL7R, involved in maintenance of chromatin structure and the immune system, associated with familial longevity already in middle-age. The size of this association increased when controls were compared to a subfraction of the offspring that had the highest probability to age healthily and become long-lived according to beneficial metabolic parameters. In conclusion, an ‘‘aging-signature’’ formed of 21 genes was identified, of which reduced expression of ASF1A and IL7R marked familial longevity already in middle-age. This indicates that expression changes of genes involved in metabolism, epigenetic control and immune function occur as a function of age, and some of these, like ASF1A and IL7R, represent early features of familial longevity and healthy ageing.

The β-amyloid precursor protein (APP), which is a key player in Alzheimer's disease, was recently reported to possess an Fe(II) binding site within its E2 domain which exhibits ferroxidase activity [Duce et al. 2010, Cell 142: 857]. The putative ligands of this site were compared to those in the ferroxidase site of ferritin. The activity was indirectly measured using transferrin, which scavenges the Fe(III) product of the reaction. A 22-residue synthetic peptide, named FD1, with the putative ferroxidase site of APP, and the E2 domain of APP were each reported to exhibit 40% of the ferroxidase activity of APP and of ceruloplasmin. It was also claimed that the ferroxidase activity of APP is inhibited by Zn(II) just as in ferritin. We measured the ferroxidase activity indirectly (i) by the incorporation of the Fe(III) product of the ferroxidase reaction into transferrin and directly (ii) by monitoring consumption of the substrate molecular oxygen. The results with the FD1 peptide were compared to the established ferroxidase activities of human H-chain ferritin and of ceruloplasmin. For FD1 we observed no activity above the background of non-enzymatic Fe(II) oxidation by molecular oxygen. Zn(II) binds to transferrin and diminishes its Fe(III) incorporation capacity and rate but it does not specifically bind to a putative ferroxidase site of FD1. Based on these results, and on comparison of the putative ligands of the ferroxidase site of APP with those of ferritin, we conclude that the previously reported results for ferroxidase activity of FD1 and – by implication – of APP should be re-evaluated.

Magnetic tweezers (MT) are a powerful tool for the study of DNA-enzyme interactions. Both the magnet-based manipulation and the camera-based detection used in MT are well suited for multiplexed measurements. Here, we systematically address challenges related to scaling of multiplexed magnetic tweezers (MMT) towards high levels of parallelization where large numbers of molecules (say 103) are addressed in the same amount of time required by a single-molecule measurement. We apply offline analysis of recorded images and show that this approach provides a scalable solution for parallel tracking of the xyz-positions of many beads simultaneously. We employ a large field-of-view imaging system to address many DNA-bead tethers in parallel. We model the 3D magnetic field generated by the magnets and derive the magnetic force experienced by DNA-bead tethers across the large field of view from first principles. We furthermore experimentally demonstrate that a DNA-bead tether subject to a rotating magnetic field describes a bicircular, Limaçon rotation pattern and that an analysis of this pattern simultaneously yields information about the force angle and the position of attachment of the DNA on the bead. Finally, we apply MMT in the high-throughput investigation of the distribution of the induced magnetic moment, the position of attachment of DNA on the beads, and DNA flexibility. The methods described herein pave the way to kilo-molecule level magnetic tweezers experiments.