Assessing Early Cardiac Outflow Tract Adaptive Responses Through Combined Experimental-Computational Manipulations

Mechanical forces are essential for proper growth and remodeling of the primitive pharyngeal arch arteries (PAAs) into the great vessels of the heart. Despite general acknowledgement of a hemodynamic-malformation link, the direct correlation between hemodynamics and PAA morphogenesis remains poorly understood. The elusiveness is largely due to difficulty in performing isolated hemodynamic perturbations and quantifying changes in-vivo. Previous in-vivo arch artery occlusion/ablation experiments either did not isolate the effects of hemodynamics, did not analyze the results in a 3D context or did not consider the effects of varying degrees of occlusion. Here, we overcome these limitations by combining minimally invasive occlusion experiments in the avian embryo with 3D anatomical models of development and in-silico testing of experimental phenomenon. We detail morphological and hemodynamic changes 24 hours post vessel occlusion. 3D anatomical models showed that occlusion geometries had more circular cross-sectional areas and more elongated arches than their control counterparts. Computational fluid dynamics revealed a marked change in wall shear stress-morphology trends. Instantaneous (in-silico) occlusion models provided mechanistic insights into the dynamic vessel adaptation process, predicting pressure-area trends for a number of experimental occlusion arches. We follow the propagation of small defects in a single embryo Hamburger Hamilton (HH) Stage 18 embryo to a more serious defect in an HH29 embryo. Results demonstrate that hemodynamic perturbation of the presumptive aortic arch, through varying degrees of vessel occlusion, overrides natural growth mechanisms and prevents it from becoming the dominant arch of the aorta.

Insights into the dynamic trajectories of protein filament division revealed by numerical investigation into the mathematical model of pure fragmentation

The dynamics by which polymeric protein filaments divide in the presence of negligible growth, for example due to the depletion of free monomeric precursors, can be described by the universal mathematical equations of ‘pure fragmentation’. The rates of fragmentation reactions reflect the stability of the protein filaments towards breakage, which is of importance in biology and biomedicine for instance in governing the creation of amyloid seeds and the propagation of prions. Here, we devised from mathematical theory inversion formulae to recover the division rates and division kernel information from time-dependent experimental measurements of filament size distribution. The numerical approach to systematically analyze the behaviour of pure fragmentation trajectories was also developed. We illustrate how these formulae can be used, provide some insights on their robustness, and show how they inform the design of experiments to measure fibril fragmentation dynamics. These advances are made possible by our central theoretical result on how the length distribution profile of the solution to the pure fragmentation equation aligns with a steady distribution profile for large times.

Amyloid fibrils are fibrillar protein structures involved in many neurodegenerative illnesses, such as Parkinson’s disease or Alzheimer’s disease. To propagate in disease, these misfolded protein aggregates must grow and divide to proliferate. Therefore, the intrinsic characteristics of their division, including the division rate and the pattern of division in terms of whether the fibrils are likely to break in the middle or at the edges, impact the disease aetiology. Here, we discovered mathematical formulae that can be used to directly extract the fibril division characteristics from recent experiments data obtained from time-dependent fibril length distribution measurements. We explain how these formulae can be used, and prove the robustness of the division rate formula where small errors in the measurement leads to small errors in the division rate. We also demonstrate that the mathematical formula is not robust enough to precisely decipher the pattern of division in the data, and suggest instead new future experimental design with short time measurements in experiments starting with fibril suspensions where all fibrils have similar size, which would be suitable to provide improved estimates.

PPIDomainMiner: Inferring domain-domain interactions from multiple sources of protein-protein interactions

Many biological processes are mediated by protein-protein interactions (PPIs). Because protein domains are the building blocks of proteins, PPIs likely rely on domain-domain interactions (DDIs). Several attempts exist to infer DDIs from PPI networks but the produced datasets are heterogeneous and sometimes not accessible, while the PPI interactome data keeps growing.

We describe a new computational approach called “PPIDM” (Protein-Protein Interactions Domain Miner) for inferring DDIs using multiple sources of PPIs. The approach is an extension of our previously described “CODAC” (Computational Discovery of Direct Associations using Common neighbors) method for inferring new edges in a tripartite graph. The PPIDM method has been applied to seven widely used PPI resources, using as “Gold-Standard” a set of DDIs extracted from 3D structural databases. Overall, PPIDM has produced a dataset of 84,552 non-redundant DDIs. Statistical significance (p-value) is calculated for each source of PPI and used to classify the PPIDM DDIs in Gold (9,175 DDIs), Silver (24,934 DDIs) and Bronze (50,443 DDIs) categories. Dataset comparison reveals that PPIDM has inferred from the 2017 releases of PPI sources about 46% of the DDIs present in the 2020 release of the 3did database, not counting the DDIs present in the Gold-Standard. The PPIDM dataset contains 10,229 DDIs that are consistent with more than 13,300 PPIs extracted from the IMEx database, and nearly 23,300 DDIs (27.5%) that are consistent with more than 214,000 human PPIs extracted from the STRING database. Examples of newly inferred DDIs covering more than 10 PPIs in the IMEx database are provided.

Further exploitation of the PPIDM DDI reservoir includes the inventory of possible partners of a protein of interest and characterization of protein interactions at the domain level in combination with other methods. The result is publicly available at http://ppidm.loria.fr/.

We revisit at a large scale the question of inferring DDIs from PPIs. Compared to previous studies, we take a unified approach accross multiple sources of PPIs. This approach is a method for inferring new edges in a tripartite graph setting and can be compared to link prediction approaches in knowledge graphs. Aggregation of several sources is performed using an optimized weighted average of the individual scores calculated in each source. A huge dataset of over 84K DDIs is produced which far exceeds the previous datasets. We show that a significant portion of the PPIDM dataset covers a large number of PPIs from curated (IMEx) or non curated (STRING) databases. Such a reservoir of DDIs deserves further exploration and can be combined with high-throughput methods such as cross-linking mass spectrometry to identify plausible protein partners of proteins of interest.

Stepwise target controllability identifies dysregulations of macrophage networks in multiple sclerosis

Identifying the nodes able to drive the state of a network is crucial to understand, and eventually control, biological systems. Despite recent advances, such identification remains difficult because of the huge number of equivalent controllable configurations, even in relatively simple networks. Based on the evidence that in many applications it is essential to test the ability of individual nodes to control a specific target subset, we develop a fast and principled method to identify controllable driver-target configurations in sparse and directed networks. We demonstrate our approach on simulated networks and experimental gene networks to characterize macrophage dysregulation in human subjects with multiple sclerosis.

Combining magnetoencephalography with magnetic resonance imaging enhances learning of surrogate-biomarkers

Electrophysiological methods, that is M/EEG, provide unique views into brain health. Yet, when building predictive models from brain data, it is often unclear how electrophysiology should be combined with other neuroimaging methods. Information can be redundant, useful common representations of multimodal data may not be obvious and multimodal data collection can be medically contraindicated, which reduces applicability. Here, we propose a multimodal model to robustly combine MEG, MRI and fMRI for prediction. We focus on age prediction as a surrogate biomarker in 674 subjects from the Cam-CAN dataset. Strikingly, MEG, fMRI and MRI showed additive effects supporting distinct brain-behavior associations. Moreover, the contribution of MEG was best explained by cortical power spectra between 8 and 30 Hz. Finally, we demonstrate that the model preserves benefits of stacking when some data is missing. The proposed framework, hence, enables multimodal learning for a wide range of biomarkers from diverse types of brain signals.

Monitoring of cardiovascular physiology augmented by a patient-specific biomechanical model during general anesthesia. A proof of concept study

During general anesthesia (GA), direct analysis of arterial pressure or aortic flow waveforms may be inconclusive in complex situations. Patient-specific biomechanical models, based on data obtained during GA and capable to perform fast simulations of cardiac cycles, have the potential to augment hemodynamic monitoring. Such models allow to simulate Pressure-Volume (PV) loops and estimate functional indicators of cardiovascular (CV) system, e.g. ventricular-arterial coupling (Vva), cardiac efficiency (CE) or myocardial contractility, evolving throughout GA. In this prospective observational study, we created patient-specific biomechanical models of heart and vasculature of a reduced geometric complexity for n = 45 patients undergoing GA, while using transthoracic echocardiography and aortic pressure and flow signals acquired in the beginning of GA (baseline condition). If intraoperative hypotension (IOH) appeared, diluted norepinephrine (NOR) was administered and the model readjusted according to the measured aortic pressure and flow signals. Such patients were a posteriori assigned into a so-called hypotensive group. The accuracy of simulated mean aortic pressure (MAP) and stroke volume (SV) at baseline were in accordance with the guidelines for the validation of new devices or reference measurement methods in all patients. After NOR administration in the hypotensive group, the percentage of concordance with 10% exclusion zone between measurement and simulation was >95% for both MAP and SV. The modeling results showed a decreased Vva (0.64±0.37 vs 0.88±0.43; p = 0.039) and an increased CE (0.8±0.1 vs 0.73±0.11; p = 0.042) in hypotensive vs normotensive patients. Furthermore, Vva increased by 92±101%, CE decreased by 13±11% (p < 0.001 for both) and contractility increased by 14±11% (p = 0.002) in the hypotensive group post-NOR administration. In this work we demonstrated the application of fast-running patient-specific biophysical models to estimate PV loops and functional indicators of CV system using clinical data available during GA. The work paves the way for model-augmented hemodynamic monitoring at operating theatres or intensive care units to enhance the information on patient-specific physiology.

Enhanced production of heterologous proteins by a synthetic microbial community: Conditions and trade-offs

Synthetic microbial consortia have been increasingly utilized in biotechnology and experimental evidence shows that suitably engineered consortia can outperform individual species in the synthesis of valuable products. Despite significant achievements, though, a quantitative understanding of the conditions that make this possible, and of the trade-offs due to the concurrent growth of multiple species, is still limited. In this work, we contribute to filling this gap by the investigation of a known prototypical synthetic consortium. A first E. coli strain, producing a heterologous protein, is sided by a second E. coli strain engineered to scavenge toxic byproducts, thus favoring the growth of the producer at the expense of diverting part of the resources to the growth of the cleaner. The simplicity of the consortium is ideal to perform an in depth-analysis and draw conclusions of more general interest. We develop a coarse-grained mathematical model that quantitatively accounts for literature data from different key growth phenotypes. Based on this, assuming growth in chemostat, we first investigate the conditions enabling stable coexistence of both strains and the effect of the metabolic load due to heterologous protein production. In these conditions, we establish when and to what extent the consortium outperforms the producer alone in terms of productivity. Finally, we show in chemostat as well as in a fed-batch scenario that gain in productivity comes at the price of a reduced yield, reflecting at the level of the consortium resource allocation trade-offs that are well-known for individual species.

In nature, microorganisms occur in communities comprising a variety of mutually interacting species. Established through evolution, these interactions allow for the survival and growth of microorganisms in their natural environment, and give rise to complex dynamics that could not be exhibited by any of the species in isolation. The richness of microbial community dynamics has been leveraged to outperform individual species in biotechnological production processes and other processes of high societal value. Yet, in view of their complexity, natural communities are difficult to study and control. In order to overcome these issues, a rapidly growing research field concerns the rational design and engineering of synthetic microbial consortia. Despite the great potential of synthetic microbial consortia, and significant efforts devoted to their mathematical modelling and analysis, a detailed understanding of how enhanced production can be achieved, and at what cost, is still unavailable. In this work, based on a quantitative model of a prototypical synthetic microbial consortium, we determine precise conditions under which a consortium outperforms individual species in the production of a recombinant protein. Moreover, we identify the inherent trade-offs between productivity and efficiency of substrate utilization.

Atlases of cognition with large-scale human brain mapping

To map the neural substrate of mental function, cognitive neuroimaging relies on controlled psychological manipulations that engage brain systems associated with specific cognitive processes. In order to build comprehensive atlases of cognitive function in the brain, it must assemble maps for many different cognitive processes, which often evoke overlapping patterns of activation. Such data aggregation faces contrasting goals: on the one hand finding correspondences across vastly different cognitive experiments, while on the other hand precisely describing the function of any given brain region. Here we introduce a new analysis framework that tackles these difficulties and thereby enables the generation of brain atlases for cognitive function. The approach leverages ontologies of cognitive concepts and multi-label brain decoding to map the neural substrate of these concepts. We demonstrate the approach by building an atlas of functional brain organization based on 30 diverse functional neuroimaging studies, totaling 196 different experimental conditions. Unlike conventional brain mapping, this functional atlas supports robust reverse inference: predicting the mental processes from brain activity in the regions delineated by the atlas. To establish that this reverse inference is indeed governed by the corresponding concepts, and not idiosyncrasies of experimental designs, we show that it can accurately decode the cognitive concepts recruited in new tasks. These results demonstrate that aggregating independent task-fMRI studies can provide a more precise global atlas of selective associations between brain and cognition.

Cognitive neuroscience uses neuroimaging to identify brain systems engaged in specific cognitive tasks. However, linking unequivocally brain systems with cognitive functions is difficult: each task probes only a small number of facets of cognition, while brain systems are often engaged in many tasks. We develop a new approach to generate a functional atlas of cognition, demonstrating brain systems selectively associated with specific cognitive functions. This approach relies upon an ontology that defines specific cognitive functions and the relations between them, along with an analysis scheme tailored to this ontology. Using a database of thirty neuroimaging studies, we show that this approach provides a highly-specific atlas of mental functions, and that it can decode the mental processes engaged in new tasks.

Traceability, reproducibility and wiki-exploration for "à-la-carte" reconstructions of genome-scale metabolic models

Genome-scale metabolic models have become the tool of choice for the global analysis of microorganism metabolism, and their reconstruction has attained high standards of quality and reliability. Improvements in this area have been accompanied by the development of some major platforms and databases, and an explosion of individual bioinformatics methods. Consequently, many recent models result from "à la carte" pipelines, combining the use of platforms, individual tools and biological expertise to enhance the quality of the reconstruction. Although very useful, introducing heterogeneous tools, that hardly interact with each other, causes loss of traceability and reproducibility in the reconstruction process. This represents a real obstacle, especially when considering less studied species whose metabolic reconstruction can greatly benefit from the comparison to good quality models of related organisms. This work proposes an adaptable workspace, AuReMe, for sustainable reconstructions or improvements of genome-scale metabolic models involving personalized pipelines. At each step, relevant information related to the modifications brought to the model by a method is stored. This ensures that the process is reproducible and documented regardless of the combination of tools used. Additionally, the workspace establishes a way to browse metabolic models and their metadata through the automatic generation of ad-hoc local wikis dedicated to monitoring and facilitating the process of reconstruction. AuReMe supports exploration and semantic query based on RDF databases. We illustrate how this workspace allowed handling, in an integrated way, the metabolic reconstructions of non-model organisms such as an extremophile bacterium or eukaryote algae. Among relevant applications, the latter reconstruction led to putative evolutionary insights of a metabolic pathway.