165 research outputs found
Functional and spatial rewiring jointly generate convergent-divergent units in self-organizing networks
Self-organization through adaptive rewiring of random neural networks
generates brain-like topologies comprising modular small-world structures with
rich club effects, merely as the product of optimizing the network topology. In
the nervous system, spatial organization is optimized no less by rewiring,
through minimizing wiring distance and maximizing spatially aligned wiring
layouts. We show that such spatial organization principles interact
constructively with adaptive rewiring, contributing to establish the networks'
connectedness and modular structures. We use an evolving neural network model
with weighted and directed connections, in which neural traffic flow is based
on consensus and advection dynamics, to show that wiring cost minimization
supports adaptive rewiring in creating convergent-divergent unit structures.
Convergent-divergent units consist of a convergent input-hub, connected to a
divergent output-hub via subnetworks of intermediate nodes, which may function
as the computational core of the unit. The prominence of minimizing wiring
distance in the dynamic evolution of the network determines the extent to which
the core is encapsulated from the rest of the network, i.e., the
context-sensitivity of its computations. This corresponds to the central role
convergent-divergent units play in establishing context-sensitivity in neuronal
information processing
Beyond sparse coding in V1
Growing evidence indicates that only a sparse subset from a pool of sensory
neurons is active for the encoding of visual stimuli at any instant in time.
Traditionally, to replicate such biological sparsity, generative models have
been using the norm as a penalty due to its convexity, which makes it
amenable to fast and simple algorithmic solvers. In this work, we use
biological vision as a test-bed and show that the soft thresholding operation
associated to the use of the norm is highly suboptimal compared to
other functions suited to approximating with
(including recently proposed Continuous Exact relaxations), both in terms of
performance and in the production of features that are akin to signatures of
the primary visual cortex. We show that sparsity produces a denser
code or employs a pool with more neurons, i.e. has a higher degree of
overcompleteness, in order to maintain the same reconstruction error as the
other methods considered. For all the penalty functions tested, a subset of the
neurons develop orientation selectivity similarly to V1 neurons. When their
code is sparse enough, the methods also develop receptive fields with varying
functionalities, another signature of V1. Compared to other methods, soft
thresholding achieves this level of sparsity at the expense of much degraded
reconstruction performance, that more likely than not is not acceptable in
biological vision. Our results indicate that V1 uses a sparsity inducing
regularization that is closer to the pseudo-norm rather than to the
norm
Tumor-immune metaphenotypes orchestrate an evolutionary bottleneck that promotes metabolic transformation
Introduction: Metabolism plays a complex role in the evolution of cancerous tumors, including inducing a multifaceted effect on the immune system to aid immune escape. Immune escape is, by definition, a collective phenomenon by requiring the presence of two cell types interacting in close proximity: tumor and immune. The microenvironmental context of these interactions is influenced by the dynamic process of blood vessel growth and remodelling, creating heterogeneous patches of well-vascularized tumor or acidic niches. Methods: Here, we present a multiscale mathematical model that captures the phenotypic, vascular, microenvironmental, and spatial heterogeneity which shapes acid-mediated invasion and immune escape over a biologically-realistic time scale. The model explores several immune escape mechanisms such as i) acid inactivation of immune cells, ii) competition for glucose, and iii) inhibitory immune checkpoint receptor expression (PD-L1). We also explore the efficacy of anti-PD-L1 and sodium bicarbonate buffer agents for treatment. To aid in understanding immune escape as a collective cellular phenomenon, we define immune escape in the context of six collective phenotypes (termed “meta-phenotypes”): Self-Acidify, Mooch Acid, PD-L1 Attack, Mooch PD-L1, Proliferate Fast, and Starve Glucose. Results: Fomenting a stronger immune response leads to initial benefits (additional cytotoxicity), but this advantage is offset by increased cell turnover that leads to accelerated evolution and the emergence of aggressive phenotypes. This creates a bimodal therapy landscape: either the immune system should be maximized for complete cure, or kept in check to avoid rapid evolution of invasive cells. These constraints are dependent on heterogeneity in vascular context, microenvironmental acidification, and the strength of immune response. Discussion: This model helps to untangle the key constraints on evolutionary costs and benefits of three key phenotypic axes on tumor invasion and treatment: acid-resistance, glycolysis, and PD-L1 expression. The benefits of concomitant anti-PD-L1 and buffer treatments is a promising treatment strategy to limit the adverse effects of immune escape
Genetic noise control via protein oligomerization
Gene expression in a cell entails random reaction events occurring over
disparate time scales. Thus, molecular noise that often results in phenotypic
and population-dynamic consequences sets a fundamental limit to biochemical
signaling. While there have been numerous studies correlating the architecture
of cellular reaction networks with noise tolerance, only a limited effort has
been made to understand the dynamic role of protein-protein interactions. Here
we have developed a fully stochastic model for the positive feedback control of
a single gene, as well as a pair of genes (toggle switch), integrating
quantitative results from previous in vivo and in vitro studies. We find that
the overall noise-level is reduced and the frequency content of the noise is
dramatically shifted to the physiologically irrelevant high-frequency regime in
the presence of protein dimerization. This is independent of the choice of
monomer or dimer as transcription factor and persists throughout the multiple
model topologies considered. For the toggle switch, we additionally find that
the presence of a protein dimer, either homodimer or heterodimer, may
significantly reduce its random switching rate. Hence, the dimer promotes the
robust function of bistable switches by preventing the uninduced (induced)
state from randomly being induced (uninduced). The specific binding between
regulatory proteins provides a buffer that may prevent the propagation of
fluctuations in genetic activity. The capacity of the buffer is a non-monotonic
function of association-dissociation rates. Since the protein oligomerization
per se does not require extra protein components to be expressed, it provides a
basis for the rapid control of intrinsic or extrinsic noise
The initial step of DNA hairpin folding: a kinetic analysis using fluorescence correlation spectroscopy
Conformational fluctuations of single-stranded DNA (ssDNA) oligonucleotides were studied in aqueous solution by monitoring contact-induced fluorescence quenching of the oxazine fluorophore MR121 by intrinsic guanosine residues (dG). We applied fluorescence correlation spectroscopy as well as steady-state and time-resolved fluorescence spectroscopy to analyze kinetics of DNA hairpin folding. We first characterized the reporter system by investigating bimolecular quenching interactions between MR121 and guanosine monophosphate in aqueous solution estimating rate constants, efficiency and stability for formation of quenched complexes. We then studied the kinetics of complex formation between MR121 and dG residues site-specifically incorporated in DNA hairpins. To uncover the initial steps of DNA hairpin folding we investigated complex formation in ssDNA carrying one or two complementary base pairs (dC–dG pairs) that could hybridize to form a short stem. Our data show that incorporation of a single dC–dG pair leads to non-exponential decays for opening and closing kinetics and reduces rate constants by one to two orders of magnitude. We found positive activation enthalpies independent of the number of dC–dG pairs. These results imply that the rate limiting step of DNA hairpin folding is not determined by loop dynamics, or by mismatches in the stem, but rather by interactions between stem and loop nucleotides
What drives the binding of minor groove-directed ligands to DNA hairpins?
Understanding the molecular basis of ligand–DNA-binding events, and its application to the rational design of novel drugs, requires knowledge of the structural features and forces that drive the corresponding recognition processes. Existing structural evidence on DNA complexation with classical minor groove-directed ligands and the corresponding studies of binding energetics have suggested that this type of binding can be described as a rigid-body association. In contrast, we show here that the binding-coupled conformational changes may be crucial for the interpretation of DNA (hairpin) association with a classical minor groove binder (netropsin). We found that, although the hairpin form is the only accessible state of ligand-free DNA, its association with the ligand may lead to its transition into a duplex conformation. It appears that formation of the fully ligated duplex from the ligand-free hairpin, occurring via two pathways, is enthalpically driven and accompanied by a significant contribution of the hydrophobic effect. Our thermodynamic and structure-based analysis, together with corresponding theoretical studies, shows that none of the predicted binding steps can be considered as a rigid-body association. In this light we anticipate our thermodynamic approach to be the basis of more sophisticated nucleic acid recognition mechanisms, which take into account the dynamic nature of both the nucleic acid and the ligand molecule
Complexity in the binding of minor groove agents: netropsin has two thermodynamically different DNA binding modes at a single site
Structural results with minor groove binding agents, such as netropsin, have provided detailed, atomic level views of DNA molecular recognition. Solution studies, however, indicate that there is complexity in the binding of minor groove agents to a single site. Netropsin, for example, has two DNA binding enthalpies in isothermal titration calorimetry (ITC) experiments that indicate the compound simultaneously forms two thermodynamically different complexes at a single AATT site. Two proposals for the origin of this unusual observation have been developed: (i) two different bound species of netropsin at single binding sites and (ii) a netropsin induced DNA hairpin to duplex transition. To develop a better understanding of DNA recognition complexity, the two proposals have been tested with several DNAs and the methods of mass spectrometry (MS), polyacrylamide gel electrophoresis (PAGE) and nuclear magnetic resonance spectroscopy in addition to ITC. All of the methods with all of the DNAs investigated clearly shows that netropsin forms two different complexes at AATT sites, and that the proposal for an induced hairpin to duplex transition in this system is incorrect
Morphology of the Homo naledi femora from Lesedi.
OBJECTIVES: The femoral remains recovered from the Lesedi Chamber are among the most complete South African fossil hominin femora discovered to date and offer new and valuable insights into the anatomy and variation of the bone in Homo naledi. While the femur is one of the best represented postcranial elements in the H. naledi assemblage from the Dinaledi Chamber, the fragmentary and commingled nature of the Dinaledi femoral remains has impeded the assessment of this element in its complete state. MATERIALS AND METHODS: Here we analyze and provide descriptions of three new relatively well-preserved femoral specimens of H. naledi from the Lesedi Chamber: U.W. 102a-001, U.W. 102a-003, and U.W. 102a-004. These femora are quantitatively and qualitatively compared to multiple extinct hominin femoral specimens, extant hominid taxa, and, where possible, each other. RESULTS: The Lesedi femora are morphologically similar to the Dinaledi femora for all overlapping regions, with differences limited to few traits of presently unknown significance. The Lesedi distal femur and mid-diaphysis preserve anatomy previously unidentified or unconfirmed in the species, including an anteroposteriorly expanded midshaft and anteriorly expanded patellar surface. The hypothesis that the Lesedi femoral sample may represent two individuals is supported. DISCUSSION: The Lesedi femora increase the range of variation of femoral morphology in H. naledi. Newly described features of the diaphysis and distal femur are either taxonomically uninformative or Homo-like. Overall, these three new femora are consistent with previous functional interpretations of the H. naledi lower limb as belonging to a species adapted for long distance walking and, possibly, running
Influence of substituent modifications on the binding of 2-amino-1,8-naphthyridines to cytosine opposite an AP site in DNA duplexes: thermodynamic characterization
Here, we report on a significant effect of substitutions on the binding affinity of a series of 2-amino-1,8-naphthyridines, i.e., 2-amino-1,8-naphthyridine (AND), 2-amino-7-methyl-1,8-naphthyridine (AMND), 2-amino-5,7-dimethyl-1,8-naphthyridine (ADMND) and 2-amino-5,6,7-trimethyl-1,8-naphthyridine (ATMND), all of which can bind to cytosine opposite an AP site in DNA duplexes. Fluorescence titration experiments show that the binding affinity for cytosine is effectively enhanced by the introduction of methyl groups to the naphthyridine ring, and the 1:1 binding constant (106 M−1) follows in the order of AND (0.30) < AMND (2.7) < ADMND (6.1) < ATMND (19) in solutions containing 110 mM Na+ (pH 7.0, at 20°C). The thermodynamic parameters obtained by isothermal titration calorimetry experiments indicate that the introduction of methyl groups effectively reduces the loss of binding entropy, which is indeed responsible for the increase in the binding affinity. The heat capacity change (ΔCp), as determined from temperature dependence of the binding enthalpy, is found to be significantly different between AND (−161 cal/mol K) and ATMND (−217 cal/mol K). The hydrophobic contribution appears to be a key force to explain the observed effect of substitutions on the binding affinity when the observed binding free energy (ΔGobs) is dissected into its component terms
Hexamer oligonucleotide topology and assembly under solution phase NMR and theoretical modeling scrutiny
The entire family of non-complementary hexamer oligodeoxyribonucleotides d(GCXYGC) (X and Y = A, G, C, or T) were assessed for topological indicators and equilibrium thermodynamics using a priori molecular modeling and solution phase NMR spectroscopy. Feasible modeled hairpin structures formed a basis from which solution structure and equilibria for each oligonucleotide were considered. 1H and 31P variable temperature (VT) and concentration dependent NMR data, NMR signal assignments and diffusion parameters led to d(GCGAGC) and d(GCGGGC) being understood as exceptions within the family in terms of self-association and topological character. A mean diffusion coefficient D298K = (2.0 ± 0.07) × 10-10 m2s-1 was evaluated across all hexamers except for d(GCGAGC) (D298K = 1.7 × 10-10 m2s-1) and d(GCGGGC) (D298K = 1.2 × 10-10 m2s-1). Melting under VT analysis (Tm = 323 K) combined with supporting NMR evidence confirmed d(GCGAGC) as the shortest tandem sheared GA mismatched duplex. Diffusion measurements were used to conclude that d(GCGGGC) preferentially exists as the shortest stable quadruplex structure. Thermodynamic analysis of all data led to the assertion that, with the exception of XY = GA and GG, the remaining non-complementary oligonucleotides adopt equilibria between monomer and duplex, contributed to largely by monomer random-coil forms. Contrastingly d(GCGAGC) showed preference for tandem sheared GA mismatch duplex formation with an association constant K = 3.9 × 105 M-1. No direct evidence was acquired for hairpin formation in any instance although its potential existence is considered possible for d(GCGAGC) on the basis of molecular modeling studies
- …