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 ℓ1\ell_1 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 $\ell_1$ 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 $\ell_1$ norm is highly suboptimal compared to other functions suited to approximating $\ell_q$ with $0 \leq q < 1$ (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 $\ell_1$ 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 $\ell_0$ pseudo-norm rather than to the $\ell_1$ norm

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

RNAstructure: software for RNA secondary structure prediction and analysis

<p>Abstract</p> <p>Background</p> <p>To understand an RNA sequence's mechanism of action, the structure must be known. Furthermore, target RNA structure is an important consideration in the design of small interfering RNAs and antisense DNA oligonucleotides. RNA secondary structure prediction, using thermodynamics, can be used to develop hypotheses about the structure of an RNA sequence.</p> <p>Results</p> <p>RNAstructure is a software package for RNA secondary structure prediction and analysis. It uses thermodynamics and utilizes the most recent set of nearest neighbor parameters from the Turner group. It includes methods for secondary structure prediction (using several algorithms), prediction of base pair probabilities, bimolecular structure prediction, and prediction of a structure common to two sequences. This contribution describes new extensions to the package, including a library of C++ classes for incorporation into other programs, a user-friendly graphical user interface written in JAVA, and new Unix-style text interfaces. The original graphical user interface for Microsoft Windows is still maintained.</p> <p>Conclusion</p> <p>The extensions to RNAstructure serve to make RNA secondary structure prediction user-friendly. The package is available for download from the Mathews lab homepage at <url>http://rna.urmc.rochester.edu/RNAstructure.html</url>.</p