A method for validating Rent's rule for technological and biological networks

Rent’s rule is empirical power law introduced in an effort to describe and optimize the wiring complexity of computer logic graphs. It is known that brain and neuronal networks also obey Rent’s rule, which is consistent with the idea that wiring costs play a fundamental role in brain evolution and development. Here we propose a method to validate this power law for a certain range of network partitions. This method is based on the bifurcation phenomenon that appears when the network is subjected to random alterations preserving its degree distribution. It has been tested on a set of VLSI circuits and real networks, including biological and technological ones. We also analyzed the effect of different types of random alterations on the Rentian scaling in order to test the influence of the degree distribution. There are network architectures quite sensitive to these randomization procedures with significant increases in the values of the Rent exponents

Uncertainty-aware predictions of molecular X-ray absorption spectra using neural network ensembles

As machine learning (ML) methods continue to be applied to a broad scope of problems in the physical sciences, uncertainty quantification is becoming correspondingly more important for their robust application. Uncertainty aware machine learning methods have been used in select applications, but largely for scalar properties. In this work, we showcase an exemplary study in which neural network ensembles are used to predict the X-ray absorption spectra of small molecules, as well as their point-wise uncertainty, from local atomic environments. The performance of the resulting surrogate clearly demonstrates quantitative correlation between errors relative to ground truth and the predicted uncertainty estimates. Significantly, the model provides an upper bound on the expected error. Specifically, an important quality of this uncertainty-aware model is that it can indicate when the model is predicting on out-of-sample data. This allows for its integration with large scale sampling of structures together with active learning or other techniques for structure refinement. Additionally, our models can be generalized to larger molecules than those used for training, and also successfully track uncertainty due to random distortions in test molecules. While we demonstrate this workflow on a specific example, ensemble learning is completely general. We believe it could have significant impact on ML-enabled forward modeling of a broad array of molecular and materials properties.Comment: 24 pages, 16 figure

Efficient Physical Embedding of Topologically Complex Information Processing Networks in Brains and Computer Circuits

Nervous systems are information processing networks that evolved by natural selection, whereas very large scale integrated (VLSI) computer circuits have evolved by commercially driven technology development. Here we follow historic intuition that all physical information processing systems will share key organizational properties, such as modularity, that generally confer adaptivity of function. It has long been observed that modular VLSI circuits demonstrate an isometric scaling relationship between the number of processing elements and the number of connections, known as Rent's rule, which is related to the dimensionality of the circuit's interconnect topology and its logical capacity. We show that human brain structural networks, and the nervous system of the nematode C. elegans, also obey Rent's rule, and exhibit some degree of hierarchical modularity. We further show that the estimated Rent exponent of human brain networks, derived from MRI data, can explain the allometric scaling relations between gray and white matter volumes across a wide range of mammalian species, again suggesting that these principles of nervous system design are highly conserved. For each of these fractal modular networks, the dimensionality of the interconnect topology was greater than the 2 or 3 Euclidean dimensions of the space in which it was embedded. This relatively high complexity entailed extra cost in physical wiring: although all networks were economically or cost-efficiently wired they did not strictly minimize wiring costs. Artificial and biological information processing systems both may evolve to optimize a trade-off between physical cost and topological complexity, resulting in the emergence of homologous principles of economical, fractal and modular design across many different kinds of nervous and computational networks

Validação de heterogeneidade estrutural em dados de Crio-ME por comitês de agrupadores

Orientadores: Fernando José Von Zuben, Rodrigo Villares PortugalDissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Elétrica e de ComputaçãoResumo: Análise de Partículas Isoladas é uma técnica que permite o estudo da estrutura tridimensional de proteínas e outros complexos macromoleculares de interesse biológico. Seus dados primários consistem em imagens de microscopia eletrônica de transmissão de múltiplas cópias da molécula em orientações aleatórias. Tais imagens são bastante ruidosas devido à baixa dose de elétrons utilizada. Reconstruções 3D podem ser obtidas combinando-se muitas imagens de partículas em orientações similares e estimando seus ângulos relativos. Entretanto, estados conformacionais heterogêneos frequentemente coexistem na amostra, porque os complexos moleculares podem ser flexíveis e também interagir com outras partículas. Heterogeneidade representa um desafio na reconstrução de modelos 3D confiáveis e degrada a resolução dos mesmos. Entre os algoritmos mais populares usados para classificação estrutural estão o agrupamento por k-médias, agrupamento hierárquico, mapas autoorganizáveis e estimadores de máxima verossimilhança. Tais abordagens estão geralmente entrelaçadas à reconstrução dos modelos 3D. No entanto, trabalhos recentes indicam ser possível inferir informações a respeito da estrutura das moléculas diretamente do conjunto de projeções 2D. Dentre estas descobertas, está a relação entre a variabilidade estrutural e manifolds em um espaço de atributos multidimensional. Esta dissertação investiga se um comitê de algoritmos de não-supervisionados é capaz de separar tais "manifolds conformacionais". Métodos de "consenso" tendem a fornecer classificação mais precisa e podem alcançar performance satisfatória em uma ampla gama de conjuntos de dados, se comparados a algoritmos individuais. Nós investigamos o comportamento de seis algoritmos de agrupamento, tanto individualmente quanto combinados em comitês, para a tarefa de classificação de heterogeneidade conformacional. A abordagem proposta foi testada em conjuntos sintéticos e reais contendo misturas de imagens de projeção da proteína Mm-cpn nos estados "aberto" e "fechado". Demonstra-se que comitês de agrupadores podem fornecer informações úteis na validação de particionamentos estruturais independetemente de algoritmos de reconstrução 3DAbstract: Single Particle Analysis is a technique that allows the study of the three-dimensional structure of proteins and other macromolecular assemblies of biological interest. Its primary data consists of transmission electron microscopy images from multiple copies of the molecule in random orientations. Such images are very noisy due to the low electron dose employed. Reconstruction of the macromolecule can be obtained by averaging many images of particles in similar orientations and estimating their relative angles. However, heterogeneous conformational states often co-exist in the sample, because the molecular complexes can be flexible and may also interact with other particles. Heterogeneity poses a challenge to the reconstruction of reliable 3D models and degrades their resolution. Among the most popular algorithms used for structural classification are k-means clustering, hierarchical clustering, self-organizing maps and maximum-likelihood estimators. Such approaches are usually interlaced with the reconstructions of the 3D models. Nevertheless, recent works indicate that it is possible to infer information about the structure of the molecules directly from the dataset of 2D projections. Among these findings is the relationship between structural variability and manifolds in a multidimensional feature space. This dissertation investigates whether an ensemble of unsupervised classification algorithms is able to separate these "conformational manifolds". Ensemble or "consensus" methods tend to provide more accurate classification and may achieve satisfactory performance across a wide range of datasets, when compared with individual algorithms. We investigate the behavior of six clustering algorithms both individually and combined in ensembles for the task of structural heterogeneity classification. The approach was tested on synthetic and real datasets containing a mixture of images from the Mm-cpn chaperonin in the "open" and "closed" states. It is shown that cluster ensembles can provide useful information in validating the structural partitionings independently of 3D reconstruction methodsMestradoEngenharia de ComputaçãoMestre em Engenharia Elétric

Comparison of methodologies to estimate state-of-health of commercial Li-ion cells from electrochemical frequency response data

Various impedance-based and nonlinear frequency response-based methods for determining the state-of-health (SOH) of commercial lithium-ion cells are evaluated. Frequency response-based measurements provide a spectral representation of dynamics of underlying physicochemical processes in the cell, giving evidence about its internal physical state. The investigated methods can be carried out more rapidly than controlled full discharge and thus constitute prospectively more efficient measurement procedures to determine the SOH of aged lithium-ion cells. We systematically investigate direct use of electrochemical impedance spectroscopy (EIS) data, equivalent circuit fits to EIS, distribution of relaxation times analysis on EIS, and nonlinear frequency response analysis. SOH prediction models are developed by correlating key parameters of each method with conventional capacity measurement (i.e., current integration). The practical feasibility, reliability and uncertainty of each of the established SOH models are considered: all models show average RMS error in the range 0.75%–1.5% SOH units, attributable principally to cell-to-cell variation. Methods based on processed data (equivalent circuit, distribution of relaxation times) are more experimentally and numerically demanding but show lower average uncertainties and may offer more flexibility for future application

Rapid UHPLC-MS metabolite profiling and phenotypic assays reveal genotypic impacts of nitrogen supplementation in oats

IntroductionOats (Avena sativa L.) are a whole grain cereal recognised for their health benefits and which are cultivated largely in temperate regions providing both a source of food for humans and animals, as well as being used in cosmetics and as a potential treatment for a number of diseases. Oats are known as being a cereal source high in dietary fibre (e.g. β-glucans), as well as being high in antioxidants, minerals and vitamins. Recently, oats have been gaining increased global attention due to their large number of beneficial health effects. Consumption of oats has been proven to lower blood LDL cholesterol levels and blood pressure, thus reducing the risk of heart disease, as well as reducing blood-sugar and insulin levels.ObjectivesOats are seen as a low input cereal. Current agricultural guidelines on nitrogen application are believed to be suboptimal and only consider the effect of nitrogen on grain yield. It is important to understand the role of both variety and of crop management in determining nutritional quality of oats. In this study the response of yield, grain quality and grain metabolites to increasing nitrogen application to levels greater than current guidelines were investigated.MethodsFour winter oat varieties (Mascani, Tardis, Balado and Gerald) were grown in a replicated nitrogen response trial consisting of a no added nitrogen control and four added nitrogen treatments between 50 and 200 kg N ha-1 in a randomised split-plot design. Grain yield, milling quality traits, β-glucan, total protein and oil content were assessed. The de-hulled oats (groats) were also subjected to a rapid Ultra High Performance Liquid Chromatography-Mass Spectrometry (UHPLC-MS) metabolomic screening approach.ResultsApplication of nitrogen had a significant effect on grain yield but there was no significant difference between the response of the four varieties. Grain quality traits however displayed significant differences both between varieties and nitrogen application level. β-glucan content significantly increased with nitrogen application. The UHPLC-MS approach has provided a rapid, sub 15 min per sample, metabolite profiling method that is repeatable and appropriate for the screening of large numbers of cereal samples. The method captured a wide range of compounds, inclusive of primary metabolites such as the amino acids, organic acids, vitamins and lipids, as well as a number of key secondary metabolites, including the avenanthramides, caffeic acid, and sinapic acid and its derivatives and was able to identify distinct metabolic phenotypes for the varieties studied. Amino acid metabolism was massively upregulated by nitrogen supplementation as were total protein levels, whilst the levels of organic acids were decreased, likely due to them acting as a carbon skeleton source. Several TCA cycle intermediates were also impacted, potentially indicating increased TCA cycle turn over, thus providing the plant with a source of energy and reductant power to aid elevated nitrogen assimilation. Elevated nitrogen availability was also directed towards the increased production of nitrogen containing phospholipids. A number of both positive and negative impacts on the metabolism of phenolic compounds that have influence upon the health beneficial value of oats and their products were also observed.ConclusionsAlthough the developed method has broad applicability as a rapid screening method or a rapid metabolite profiling method and in this study has provided valuable metabolic insights, it still must be considered that much greater confidence in metabolite identification, as well as quantitative precision, will be gained by the application of higher resolution chromatography methods, although at a large expense to sample throughput. Follow up studies will apply higher resolution GC (gas chromatography) and LC (reversed phase and HILIC) approaches, oats will be also analysed from across multiple growth locations and growth seasons, effectively providing a cross validation for the results obtained within this preliminary study. It will also be fascinating to perform more controlled experiments with sampling of green tissues, as well as oat grains, throughout the plants and grains development, to reveal greater insight of carbon and nitrogen metabolism balance, as well as resource partitioning into lipid and secondary metabolism

Continuum solvation models: Dissecting the free energy of solvation

The most usual self-consistent reaction field (SCRF) continuum models for the description of solvation within the quantum mechanical (QM) framework are reviewed, trying to emphasize their common roots as well as the inherent approximations assumed in the calculation of the free energy of solvation. Particular attention is also paid to the specific features involved in the development of current state-of-the-art QM SCRF continuum models. This is used to discuss the need to maintain a close correspondence between each SCRF formalism and the specific details entailing its parametrization, as well as the need to be cautious in analyzing the balance between electrostatic and non-electrostatic contributions to the solvation free energy between different SCRF models. Finally, special emphasis is given to the post-processing of the free energy of solvation to derive parameters providing a compact picture of the ability of a molecule to interact with different solvents, which can be of particular interest in biopharmaceutical studies