wfl Python toolkit for creating machine learning interatomic potentials and related atomistic simulation workflows.

Predictive atomistic simulations are increasingly employed for data intensive high throughput studies that take advantage of constantly growing computational resources. To handle the sheer number of individual calculations that are needed in such studies, workflow management packages for atomistic simulations have been developed for a rapidly growing user base. These packages are predominantly designed to handle computationally heavy ab initio calculations, usually with a focus on data provenance and reproducibility. However, in related simulation communities, e.g., the developers of machine learning interatomic potentials (MLIPs), the computational requirements are somewhat different: the types, sizes, and numbers of computational tasks are more diverse and, therefore, require additional ways of parallelization and local or remote execution for optimal efficiency. In this work, we present the atomistic simulation and MLIP fitting workflow management package wfl and Python remote execution package ExPyRe to meet these requirements. With wfl and ExPyRe, versatile atomic simulation environment based workflows that perform diverse procedures can be written. This capability is based on a low-level developer-oriented framework, which can be utilized to construct high level functionality for user-friendly programs. Such high level capabilities to automate machine learning interatomic potential fitting procedures are already incorporated in wfl, which we use to showcase its capabilities in this work. We believe that wfl fills an important niche in several growing simulation communities and will aid the development of efficient custom computational tasks.

Cancer lineage-specific regulation of YAP responsive elements revealed through large-scale functional epigenomic screens.

YAP is a key transcriptional co-activator of TEADs, it regulates cell growth and is frequently activated in cancer. In Malignant Pleural Mesothelioma (MPM), YAP is activated by loss-of-function mutations in upstream components of the Hippo pathway, while, in Uveal Melanoma (UM), YAP is activated in a Hippo-independent manner. To date, it is unclear if and how the different oncogenic lesions activating YAP impact its oncogenic program, which is particularly relevant for designing selective anti-cancer therapies. Here we show that, despite YAP being essential in both MPM and UM, its interaction with TEAD is unexpectedly dispensable in UM, limiting the applicability of TEAD inhibitors in this cancer type. Systematic functional interrogation of YAP regulatory elements in both cancer types reveals convergent regulation of broad oncogenic drivers in both MPM and UM, but also strikingly selective programs. Our work reveals unanticipated lineage-specific features of the YAP regulatory network that provide important insights to guide the design of tailored therapeutic strategies to inhibit YAP signaling across different cancer types.

Identification of the Transcription Factor ATF3 as a Direct and Indirect Regulator of the LDLR.

Coronary artery disease (CAD) is a complex, multifactorial disease caused, in particular, by inflammation and cholesterol metabolism. At the molecular level, the role of tissue-specific signaling pathways leading to CAD is still largely unexplored. This study relied on two main resources: (1) genes with impact on atherosclerosis/CAD, and (2) liver-specific transcriptome analyses from human and mouse studies. The transcription factor activating transcription factor 3 (ATF3) was identified as a key regulator of a liver network relevant to atherosclerosis and linked to inflammation and cholesterol metabolism. ATF3 was predicted to be a direct and indirect (via MAF BZIP Transcription Factor F (MAFF)) regulator of low-density lipoprotein receptor (LDLR). Chromatin immunoprecipitation DNA sequencing (ChIP-seq) data from human liver cells revealed an ATF3 binding motif in the promoter regions of MAFF and LDLR. siRNA knockdown of ATF3 in human Hep3B liver cells significantly upregulated LDLR expression (p < 0.01). Inflammation induced by lipopolysaccharide (LPS) stimulation resulted in significant upregulation of ATF3 (p < 0.01) and subsequent downregulation of LDLR (p < 0.001). Liver-specific expression data from human CAD patients undergoing coronary artery bypass grafting (CABG) surgery (STARNET) and mouse models (HMDP) confirmed the regulatory role of ATF3 in the homeostasis of cholesterol metabolism. This study suggests that ATF3 might be a promising treatment candidate for lowering LDL cholesterol and reducing cardiovascular risk.

A Hybrid Machine Learning Approach for Structure Stability Prediction in Molecular Co-crystal Screenings.

Co-crystals are a highly interesting material class as varying their components and stoichiometry in principle allows tuning supramolecular assemblies toward desired physical properties. The in silico prediction of co-crystal structures represents a daunting task, however, as they span a vast search space and usually feature large unit cells. This requires theoretical models that are accurate and fast to evaluate, a combination that can in principle be accomplished by modern machine-learned (ML) potentials trained on first-principles data. Crucially, these ML potentials need to account for the description of long-range interactions, which are essential for the stability and structure of molecular crystals. In this contribution, we present a strategy for developing Δ-ML potentials for co-crystals, which use a physical baseline model to describe long-range interactions. The applicability of this approach is demonstrated for co-crystals of variable composition consisting of an active pharmaceutical ingredient and various co-formers. We find that the Δ-ML approach offers a strong and consistent improvement over the density functional tight binding baseline. Importantly, this even holds true when extrapolating beyond the scope of the training set, for instance in molecular dynamics simulations under ambient conditions.

Vascular Tissue Specific miRNA Profiles Reveal Novel Correlations with Risk Factors in Coronary Artery Disease.

Cardiovascular disease (CVD) is the leading cause of morbidity and mortality worldwide. Non-coding RNAs have already been linked to CVD development and progression. While microRNAs (miRs) have been well studied in blood samples, there is little data on tissue-specific miRs in cardiovascular relevant tissues and their relation to cardiovascular risk factors. Tissue-specific miRs derived from Arteria mammaria interna (IMA) from 192 coronary artery disease (CAD) patients undergoing coronary artery bypass grafting (CABG) were analyzed. The aims of the study were 1) to establish a reference atlas which can be utilized for identification of novel diagnostic biomarkers and potential therapeutic targets, and 2) to relate these miRs to cardiovascular risk factors. Overall, 393 individual miRs showed sufficient expression levels and passed quality control for further analysis. We identified 17 miRs-miR-10b-3p, miR-10-5p, miR-17-3p, miR-21-5p, miR-151a-5p, miR-181a-5p, miR-185-5p, miR-194-5p, miR-199a-3p, miR-199b-3p, miR-212-3p, miR-363-3p, miR-548d-5p, miR-744-5p, miR-3117-3p, miR-5683 and miR-5701-significantly correlated with cardiovascular risk factors (correlation coefficient >0.2 in both directions, p-value (p < 0.006, false discovery rate (FDR) <0.05). Of particular interest, miR-5701 was positively correlated with hypertension, hypercholesterolemia, and diabetes. In addition, we found that miR-629-5p and miR-98-5p were significantly correlated with acute myocardial infarction. We provide a first atlas of miR profiles in IMA samples from CAD patients. In perspective, these miRs might play an important role in improved risk assessment, mechanistic disease understanding and local therapy of CAD.

Systematic dissection of transcriptional regulatory networks by genome-scale and single-cell CRISPR screens.

Millions of putative transcriptional regulatory elements (TREs) have been cataloged in the human genome, yet their functional relevance in specific pathophysiological settings remains to be determined. This is critical to understand how oncogenic transcription factors (TFs) engage specific TREs to impose transcriptional programs underlying malignant phenotypes. Here, we combine cutting edge CRISPR screens and epigenomic profiling to functionally survey ≈15,000 TREs engaged by estrogen receptor (ER). We show that ER exerts its oncogenic role in breast cancer by engaging TREs enriched in GATA3, TFAP2C, and H3K27Ac signal. These TREs control critical downstream TFs, among which TFAP2C plays an essential role in ER-driven cell proliferation. Together, our work reveals novel insights into a critical oncogenic transcription program and provides a framework to map regulatory networks, enabling to dissect the function of the noncoding genome of cancer cells.

Data-efficient machine learning for molecular crystal structure prediction.

The combination of modern machine learning (ML) approaches with high-quality data from quantum mechanical (QM) calculations can yield models with an unrivalled accuracy/cost ratio. However, such methods are ultimately limited by the computational effort required to produce the reference data. In particular, reference calculations for periodic systems with many atoms can become prohibitively expensive for higher levels of theory. This trade-off is critical in the context of organic crystal structure prediction (CSP). Here, a data-efficient ML approach would be highly desirable, since screening a huge space of possible polymorphs in a narrow energy range requires the assessment of a large number of trial structures with high accuracy. In this contribution, we present tailored Δ-ML models that allow screening a wide range of crystal candidates while adequately describing the subtle interplay between intermolecular interactions such as H-bonding and many-body dispersion effects. This is achieved by enhancing a physics-based description of long-range interactions at the density functional tight binding (DFTB) level-for which an efficient implementation is available-with a short-range ML model trained on high-quality first-principles reference data. The presented workflow is broadly applicable to different molecular materials, without the need for a single periodic calculation at the reference level of theory. We show that this even allows the use of wavefunction methods in CSP.

Mapping Materials and Molecules.

The visualization of data is indispensable in scientific research, from the early stages when human insight forms to the final step of communicating results. In computational physics, chemistry and materials science, it can be as simple as making a scatter plot or as straightforward as looking through the snapshots of atomic positions manually. However, as a result of the "big data" revolution, these conventional approaches are often inadequate. The widespread adoption of high-throughput computation for materials discovery and the associated community-wide repositories have given rise to data sets that contain an enormous number of compounds and atomic configurations. A typical data set contains thousands to millions of atomic structures, along with a diverse range of properties such as formation energies, band gaps, or bioactivities.It would thus be desirable to have a data-driven and automated framework for visualizing and analyzing such structural data sets. The key idea is to construct a low-dimensional representation of the data, which facilitates navigation, reveals underlying patterns, and helps to identify data points with unusual attributes. Such data-intensive maps, often employing machine learning methods, are appearing more and more frequently in the literature. However, to the wider community, it is not always transparent how these maps are made and how they should be interpreted. Furthermore, while these maps undoubtedly serve a decorative purpose in academic publications, it is not always apparent what extra information can be garnered from reading or making them.This Account attempts to answer such questions. We start with a concise summary of the theory of representing chemical environments, followed by the introduction of a simple yet practical conceptual approach for generating structure maps in a generic and automated manner. Such analysis and mapping is made nearly effortless by employing the newly developed software tool ASAP. To showcase the applicability to a wide variety of systems in chemistry and materials science, we provide several illustrative examples, including crystalline and amorphous materials, interfaces, and organic molecules. In these examples, the maps not only help to sift through large data sets but also reveal hidden patterns that could be easily missed using conventional analyses.The explosion in the amount of computed information in chemistry and materials science has made visualization into a science in itself. Not only have we benefited from exploiting these visualization methods in previous works, we also believe that the automated mapping of data sets will in turn stimulate further creativity and exploration, as well as ultimately feed back into future advances in the respective fields.

A Trickster in Disguise: Hyaluronan's Ambivalent Roles in the Matrix.

Hyaluronan (HA) is a simple but diverse glycosaminoglycan. It plays a major role in aging, cellular senescence, cancer, and tissue homeostasis. In which way HA affects the surrounding tissues greatly depends on the molecular weight of HA. Whereas high molecular weight HA is associated with homeostasis and protective effects, HA fragments tend to be linked to the pathologic state. Furthermore, the interaction of HA with its binding partners, the hyaladherins, such as CD44, is essential for sustaining tissue integrity and is likewise related to cancer. The naked mole rat, a rodent species, possesses a special form of very high molecular weight (vHMW) HA, which is associated with the extraordinary cancer resistance and longevity of those animals. This review addresses HA and its diverse facets: from HA synthesis to degradation, from oligomeric HA to vHMW-HA and from its beneficial properties to the involvement in pathologies. We further discuss the functions of HA in the naked mole rat and compare them to human conditions. Though intensively researched, this simple polymer bears some secrets that may hold the key for a better understanding of cellular processes and the development of diseases, such as cancer.