First-principles calculation of DNA looping in tethered particle experiments

We calculate the probability of DNA loop formation mediated by regulatory proteins such as Lac repressor (LacI), using a mathematical model of DNA elasticity. Our model is adapted to calculating quantities directly observable in Tethered Particle Motion (TPM) experiments, and it accounts for all the entropic forces present in such experiments. Our model has no free parameters; it characterizes DNA elasticity using information obtained in other kinds of experiments. [...] We show how to compute both the "looping J factor" (or equivalently, the looping free energy) for various DNA construct geometries and LacI concentrations, as well as the detailed probability density function of bead excursions. We also show how to extract the same quantities from recent experimental data on tethered particle motion, and then compare to our model's predictions. [...] Our model successfully reproduces the detailed distributions of bead excursion, including their surprising three-peak structure, without any fit parameters and without invoking any alternative conformation of the LacI tetramer. Indeed, the model qualitatively reproduces the observed dependence of these distributions on tether length (e.g., phasing) and on LacI concentration (titration). However, for short DNA loops (around 95 basepairs) the experiments show more looping than is predicted by the harmonic-elasticity model, echoing other recent experimental results. Because the experiments we study are done in vitro, this anomalously high looping cannot be rationalized as resulting from the presence of DNA-bending proteins or other cellular machinery. We also show that it is unlikely to be the result of a hypothetical "open" conformation of the LacI tetramer.Comment: See the supplement at http://www.physics.upenn.edu/~pcn/Ms/TowlesEtalSuppl.pdf . This revised version accepted for publication at Physical Biolog

Concentration and Length Dependence of DNA Looping in Transcriptional Regulation

In many cases, transcriptional regulation involves the binding of transcription factors at sites on the DNA that are not immediately adjacent to the promoter of interest. This action at a distance is often mediated by the formation of DNA loops: Binding at two or more sites on the DNA results in the formation of a loop, which can bring the transcription factor into the immediate neighborhood of the relevant promoter. These processes are important in settings ranging from the historic bacterial examples (bacterial metabolism and the lytic-lysogeny decision in bacteriophage), to the modern concept of gene regulation to regulatory processes central to pattern formation during development of multicellular organisms. Though there have been a variety of insights into the combinatorial aspects of transcriptional control, the mechanism of DNA looping as an agent of combinatorial control in both prokaryotes and eukaryotes remains unclear. We use single-molecule techniques to dissect DNA looping in the lac operon. In particular, we measure the propensity for DNA looping by the Lac repressor as a function of the concentration of repressor protein and as a function of the distance between repressor binding sites. As with earlier single-molecule studies, we find (at least) two distinct looped states and demonstrate that the presence of these two states depends both upon the concentration of repressor protein and the distance between the two repressor binding sites. We find that loops form even at interoperator spacings considerably shorter than the DNA persistence length, without the intervention of any other proteins to prebend the DNA. The concentration measurements also permit us to use a simple statistical mechanical model of DNA loop formation to determine the free energy of DNA looping, or equivalently, the J-factor for looping

Structure-based stabilization of insulin as a therapeutic protein assembly via enhanced aromatic-aromatic interactions

Key contributions to protein structure and stability are provided by weakly polar interactions, which arise from asymmetric electronic distributions within amino acids and peptide bonds. Of particular interest are aromatic side chains whose directional π-systems commonly stabilize protein interiors and interfaces. Here, we consider aromatic-aromatic interactions within a model protein assembly: the dimer interface of insulin. Semi-classical simulations of aromatic-aromatic interactions at this interface suggested that substitution of residue TyrB26 by Trp would preserve native structure while enhancing dimerization (and hence hexamer stability). The crystal structure of a [TrpB26]insulin analog (determined as a T3Rf3 zinc hexamer at a resolution of 2.25 Å) was observed to be essentially identical to that of WT insulin. Remarkably and yet in general accordance with theoretical expectations, spectroscopic studies demonstrated a 150-fold increase in the in vitro lifetime of the variant hexamer, a critical pharmacokinetic parameter influencing design of long-acting formulations. Functional studies in diabetic rats indeed revealed prolonged action following subcutaneous injection. The potency of the TrpB26-modified analog was equal to or greater than an unmodified control. Thus, exploiting a general quantum-chemical feature of protein structure and stability, our results exemplify a mechanism-based approach to the optimization of a therapeutic protein assembly

Missense-depleted regions in population exomes implicate ras superfamily nucleotide-binding protein alteration in patients with brain malformation.

Genomic sequence interpretation can miss clinically relevant missense variants for several reasons. Rare missense variants are numerous in the exome and difficult to prioritise. Affected genes may also not have existing disease association. To improve variant prioritisation, we leverage population exome data to identify intragenic missense-depleted regions (MDRs) genome-wide that may be important in disease. We then use missense depletion analyses to help prioritise undiagnosed disease exome variants. We demonstrate application of this strategy to identify a novel gene association for human brain malformation. We identified de novo missense variants that affect the GDP/GTP-binding site of ARF1 in three unrelated patients. Corresponding functional analysis suggests ARF1 GDP/GTP-activation is affected by the specific missense mutations associated with heterotopia. These findings expand the genetic pathway underpinning neurologic disease that classically includes FLNA. ARF1 along with ARFGEF2 add further evidence implicating ARF/GEFs in the brain. Using functional ontology, top MDR-containing genes were highly enriched for nucleotide-binding function, suggesting these may be candidates for human disease. Routine consideration of MDR in the interpretation of exome data for rare diseases may help identify strong genetic factors for many severe conditions, infertility/reduction in reproductive capability, and embryonic conditions contributing to preterm loss

Analysis of Nematic Liquid Crystals with Disclination Lines

We investigate the structure of nematic liquid crystal thin films described by the Landau--de Gennes tensor-valued order parameter with Dirichlet boundary conditions of nonzero degree. We prove that as the elasticity constant goes to zero a limiting uniaxial texture forms with disclination lines corresponding to a finite number of defects, all of degree 1/2 or all of degree -1/2. We also state a result on the limiting behavior of minimizers of the Chern-Simons-Higgs model without magnetic field that follows from a similar proof.Comment: 40 pages, 1 figur

NEXUS/Physics: An interdisciplinary repurposing of physics for biologists

In response to increasing calls for the reform of the undergraduate science curriculum for life science majors and pre-medical students (Bio2010, Scientific Foundations for Future Physicians, Vision & Change), an interdisciplinary team has created NEXUS/Physics: a repurposing of an introductory physics curriculum for the life sciences. The curriculum interacts strongly and supportively with introductory biology and chemistry courses taken by life sciences students, with the goal of helping students build general, multi-discipline scientific competencies. In order to do this, our two-semester NEXUS/Physics course sequence is positioned as a second year course so students will have had some exposure to basic concepts in biology and chemistry. NEXUS/Physics stresses interdisciplinary examples and the content differs markedly from traditional introductory physics to facilitate this. It extends the discussion of energy to include interatomic potentials and chemical reactions, the discussion of thermodynamics to include enthalpy and Gibbs free energy, and includes a serious discussion of random vs. coherent motion including diffusion. The development of instructional materials is coordinated with careful education research. Both the new content and the results of the research are described in a series of papers for which this paper serves as an overview and context.Comment: 12 page

Extensive HST Ultraviolet Spectra and Multi-wavelength Observations of SN 2014J in M82 Indicate Reddening and Circumstellar Scattering by Typical Dust

SN 2014J in M82 is the closest detected Type Ia supernova (SN Ia) in at least 28 years and perhaps in 410 years. Despite its small distance of 3.3 Mpc, SN 2014J is surprisingly faint, peaking at V = 10.6 mag, and assuming a typical SN Ia luminosity, we infer an observed visual extinction of A_V = 2.0 +/- 0.1 mag. But this picture, with R_V = 1.6 +/- 0.2, is too simple to account for all observations. We combine 10 epochs (spanning a month) of HST/STIS ultraviolet through near-infrared spectroscopy with HST/WFC3, KAIT, and FanCam photometry from the optical to the infrared and 9 epochs of high-resolution TRES spectroscopy to investigate the sources of extinction and reddening for SN 2014J. We argue that the wide range of observed properties for SN 2014J is caused by a combination of dust reddening, likely originating in the interstellar medium of M82, and scattering off circumstellar material. For this model, roughly half of the extinction is caused by reddening from typical dust (E(B-V ) = 0.45 mag and R_V = 2.6) and roughly half by scattering off LMC-like dust in the circumstellar environment of SN 2014J.Comment: 17 pages (excluding references and tables), 15 figures, accepted to MNRAS. A high-resolution HST image of SN 2014J in M82 is available upon reques

Multiscale, multimodal analysis of tumor heterogeneity in IDH1 mutant vs wild-type diffuse gliomas.

Glioma is recognized to be a highly heterogeneous CNS malignancy, whose diverse cellular composition and cellular interactions have not been well characterized. To gain new clinical- and biological-insights into the genetically-bifurcated IDH1 mutant (mt) vs wildtype (wt) forms of glioma, we integrated data from protein, genomic and MR imaging from 20 treatment-naïve glioma cases and 16 recurrent GBM cases. Multiplexed immunofluorescence (MxIF) was used to generate single cell data for 43 protein markers representing all cancer hallmarks, Genomic sequencing (exome and RNA (normal and tumor) and magnetic resonance imaging (MRI) quantitative features (protocols were T1-post, FLAIR and ADC) from whole tumor, peritumoral edema and enhancing core vs equivalent normal region were also collected from patients. Based on MxIF analysis, 85,767 cells (glioma cases) and 56,304 cells (GBM cases) were used to generate cell-level data for 24 biomarkers. K-means clustering was used to generate 7 distinct groups of cells with divergent biomarker profiles and deconvolution was used to assign RNA data into three classes. Spatial and molecular heterogeneity metrics were generated for the cell data. All features were compared between IDH mt and IDHwt patients and were finally combined to provide a holistic/integrated comparison. Protein expression by hallmark was generally lower in the IDHmt vs wt patients. Molecular and spatial heterogeneity scores for angiogenesis and cell invasion also differed between IDHmt and wt gliomas irrespective of prior treatment and tumor grade; these differences also persisted in the MR imaging features of peritumoral edema and contrast enhancement volumes. A coherent picture of enhanced angiogenesis in IDHwt tumors was derived from multiple platforms (genomic, proteomic and imaging) and scales from individual proteins to cell clusters and heterogeneity, as well as bulk tumor RNA and imaging features. Longer overall survival for IDH1mt glioma patients may reflect mutation-driven alterations in cellular, molecular, and spatial heterogeneity which manifest in discernable radiological manifestations

Identification and inference in a simultaneous equation under alternative information sets and sampling schemes

In simple static linear simultaneous equation models, the empirical distributions of IV and OLS are examined under alternative sampling schemes and compared with their first-order asymptotic approximations. We demonstrate that the limiting distribution of consistent IV is not affected by conditioning on exogenous regressors, whereas that of inconsistent OLS is. The OLS asymptotic and simulated actual variances are shown to diminish by extending the set of exogenous variables kept fixed in sampling, whereas such an extension disrupts the distribution of IV and deteriorates the accuracy of its standard asymptotic approximation, not only when instruments are weak. Against this background, the consequences for the identification of parameters of interest are examined for a setting in which (in practice often incredible) assumptions regarding the zero correlation between instruments and disturbances are replaced by (generally more credible) interval assumptions on the correlation between endogenous regressor and disturbance. This yields OLS-based modified confidence intervals, which are usually conservative, as is established by simulation. Often they compare favourably with IV-based intervals and accentuate their frailty. The latter is demonstrated in an empirical illustration

Prospective memory: age related change is influenced by APOE genotype

Non-focal prospective memory (PM) is sensitive to age-related decline; an additional impairment in focal PM is characteristic of mild stage Alzheimer’s disease. This research explored whether, by mid-adulthood, the distinct demands of focal and non-focal PM expose differences in carriers of an APOE ε4 allele, a genetic risk factor for Alzheimer’s disease. Thirty-three young and 55 mid-age adults, differentiated by APOE genotype, completed a category-decision task with a concurrent focal or non-focal PM demand. Only mid-age ε4 carriers showed a cost of carrying a focal PM intention. In addition, mid-age ε4 carriers showed a significantly greater cost of carrying a non-focal PM intention than young ε4 carriers, supporting a profile of accelerated aging. Consistency in the profile of cost differences observed in mid-age ε4 carriers and pathological aging may indicate premature vulnerability. Future research correlating a shift in PM performance with early genotype differences in brain-based markers of decline is important