Insights from incorporating quantum computing into drug design workflows.

MOTIVATION: While many quantum computing (QC) methods promise theoretical advantages over classical counterparts, quantum hardware remains limited. Exploiting near-term QC in computer-aided drug design (CADD) thus requires judicious partitioning between classical and quantum calculations. RESULTS: We present HypaCADD, a hybrid classical-quantum workflow for finding ligands binding to proteins, while accounting for genetic mutations. We explicitly identify modules of our drug-design workflow currently amenable to replacement by QC: non-intuitively, we identify the mutation-impact predictor as the best candidate. HypaCADD thus combines classical docking and molecular dynamics with quantum machine learning (QML) to infer the impact of mutations. We present a case study with the coronavirus (SARS-CoV-2) protease and associated mutants. We map a classical machine-learning module onto QC, using a neural network constructed from qubit-rotation gates. We have implemented this in simulation and on two commercial quantum computers. We find that the QML models can perform on par with, if not better than, classical baselines. In summary, HypaCADD offers a successful strategy for leveraging QC for CADD. AVAILABILITY AND IMPLEMENTATION: Jupyter Notebooks with Python code are freely available for academic use on GitHub: https://www.github.com/hypahub/hypacadd_notebook. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

LongISLND: in silico sequencing of lengthy and noisy datatypes.

LongISLND is a software package designed to simulate sequencing data according to the characteristics of third generation, single-molecule sequencing technologies. The general software architecture is easily extendable, as demonstrated by the emulation of Pacific Biosciences (PacBio) multi-pass sequencing with P5 and P6 chemistries, producing data in FASTQ, H5, and the latest PacBio BAM format. We demonstrate its utility by downstream processing with consensus building and variant calling. AVAILABILITY AND IMPLEMENTATION: LongISLND is implemented in Java and available at http://bioinform.github.io/longislnd CONTACT: hugo.lam@roche.comSupplementary information: Supplementary data are available at Bioinformatics online.

Theory of the magnetic and metal-insulator transitions in RNiO3 bulk and layered structures.

A slave rotor--Hartree-Fock formalism is presented for studying the properties of the p-d model describing perovskite transition metal oxides, and a flexible and efficient numerical formalism is developed for its solution. The methodology is shown to yield, within a unified formulation, the significant aspects of the rare-earth nickelate phase diagram, including the paramagnetic metal state observed for the LaNiO3 and the correct ground-state magnetic order of insulating compounds. It is then used to elucidate ground state changes occurring as morphology is varied from bulk to strained and unstrained thin-film form. For ultrathin films, epitaxial strain and charge transfer to the apical out-of-plane oxygen sites are shown to have significant impact on the phase diagram.

Gene Expression Scoring of Immune Activity Levels for Precision Use of Hydrocortisone in Vasodilatory Shock.

PURPOSE: Among patients with vasodilatory shock, gene expression scores may identify different immune states. We aimed to test whether such scores are robust in identifying patients' immune state and predicting response to hydrocortisone treatment in vasodilatory shock. MATERIALS AND METHODS: We selected genes to generate continuous scores to define previously established subclasses of sepsis. We used these scores to identify a patient's immune state. We evaluated the potential for these states to assess the differential effect of hydrocortisone in two randomized clinical trials of hydrocortisone versus placebo in vasodilatory shock. RESULTS: We initially identified genes associated with immune-adaptive, immune-innate, immune-coagulant functions. From these genes, 15 were most relevant to generate expression scores related to each of the functions. These scores were used to identify patients as immune-adaptive prevalent (IA-P) and immune-innate prevalent (IN-P). In IA-P patients, hydrocortisone therapy increased 28-day mortality in both trials (43.3% vs 14.7%, P = 0.028) and (57.1% vs 0.0%, P = 0.99). In IN-P patients, this effect was numerically reversed. CONCLUSIONS: Gene expression scores identified the immune state of vasodilatory shock patients, one of which (IA-P) identified those who may be harmed by hydrocortisone. Gene expression scores may help advance the field of personalized medicine.

High-spin polaron in lightly doped CuO2 planes.

We derive and investigate numerically a minimal yet detailed spin-polaron model that describes lightly doped CuO2 layers. The low-energy physics of a hole is studied by total-spin-resolved exact diagonalization on clusters of up to 32 CuO2 unit cells, revealing features missed by previous studies. In particular, spin-polaron states with total spin 3/2 are the lowest eigenstates in some regions of the Brillouin zone. In these regions, and also at other points, the quasiparticle weight is identically zero indicating orthogonal states to those represented in the one electron Green's function. This highlights the importance of the proper treatment of spin fluctuations in the many-body background.

Deep convolutional neural networks for accurate somatic mutation detection.

Accurate detection of somatic mutations is still a challenge in cancer analysis. Here we present NeuSomatic, the first convolutional neural network approach for somatic mutation detection, which significantly outperforms previous methods on different sequencing platforms, sequencing strategies, and tumor purities. NeuSomatic summarizes sequence alignments into small matrices and incorporates more than a hundred features to capture mutation signals effectively. It can be used universally as a stand-alone somatic mutation detection method or with an ensemble of existing methods to achieve the highest accuracy.

Comprehensive variation discovery in single human genomes.

Complete knowledge of the genetic variation in individual human genomes is a crucial foundation for understanding the etiology of disease. Genetic variation is typically characterized by sequencing individual genomes and comparing reads to a reference. Existing methods do an excellent job of detecting variants in approximately 90% of the human genome; however, calling variants in the remaining 10% of the genome (largely low-complexity sequence and segmental duplications) is challenging. To improve variant calling, we developed a new algorithm, DISCOVAR, and examined its performance on improved, low-cost sequence data. Using a newly created reference set of variants from the finished sequence of 103 randomly chosen fosmids, we find that some standard variant call sets miss up to 25% of variants. We show that the combination of new methods and improved data increases sensitivity by several fold, with the greatest impact in challenging regions of the human genome.