Perspective on the Current State-of-the-Art of Quantum Computing for Drug Discovery Applications

Future Potential of Quantum Computing and Simulations in Biological Science

The review article presents the recent progress in quantum computing and simulation within the field of biological sciences. The article is designed mainly into two portions: quantum computing and quantum simulation. In the first part, significant aspects of quantum computing was illustrated, such as quantum hardware, quantum RAM and big data, modern quantum processors, qubit, superposition effect in quantum computation, quantum interference, quantum entanglement, and quantum logic gates. Simultaneously, in the second part, vital features of the quantum simulation was illustrated, such as the quantum simulator, algorithms used in quantum simulations, and the use of quantum simulation in biological science. Finally, the review provides exceptional views to future researchers about different aspects of quantum simulation in biological science.

Ab initio extended Hubbard model of short polyenes for efficient quantum computing

We propose introducing an extended Hubbard Hamiltonian derived via the ab initio downfolding method, which was originally formulated for periodic materials, toward efficient quantum computing of molecular electronic structure calculations. By utilizing this method, the first-principles Hamiltonian of chemical systems can be coarse-grained by eliminating the electronic degrees of freedom in higher energy space and reducing the number of terms of electron repulsion integral from O(N4) to O(N2). Our approach is validated numerically on the vertical excitation energies and excitation characters of ethylene, butadiene, and hexatriene. The dynamical electron correlation is incorporated within the framework of the constrained random phase approximation in advance of quantum computations, and the constructed models capture the trend of experimental and high-level quantum chemical calculation results. As expected, the L1-norm of the fermion-to-qubit mapped model Hamiltonians is significantly lower than that of conventional ab initio Hamiltonians, suggesting improved scalability of quantum computing. Those numerical outcomes and the results of the simulation of excited-state sampling demonstrate that the ab initio extended Hubbard Hamiltonian may hold significant potential for quantum chemical calculations using quantum computers.

Revolutionizing Healthcare: The Emerging Role of Quantum Computing in Enhancing Medical Technology and Treatment

The healthcare sector faces complex challenges that call for innovative solutions to improve diagnostic accuracy, treatment efficacy, and data management. Quantum computing, with its unique capabilities, holds the potential to revolutionize various aspects of healthcare. This narrative review critically examines the existing literature on the application of quantum computing in healthcare, focusing on its utility in enhancing diagnostics, data processing, and treatment planning. Quantum computing's ability to handle large, complex datasets more efficiently than classical computers can significantly impact domains such as genomics, medical imaging, and personalized medicine. Quantum algorithms can accelerate the identification of genetic markers associated with diseases, facilitate the analysis of medical images, and optimize treatment plans based on individual genetic profiles. Moreover, quantum cryptography offers a robust security solution for safeguarding sensitive patient data, a critical need as healthcare increasingly relies on digital platforms. Despite the promising outlook, the integration of quantum computing into healthcare faces technical, ethical, and regulatory challenges. The delicate nature of quantum hardware, the need for error correction, and the scalability of quantum systems pose barriers to widespread adoption. Additionally, concerns around patient privacy and data security, as well as the need for updated regulatory frameworks, must be addressed. Ongoing research and collaborative efforts involving researchers, healthcare providers, and technology developers are crucial to overcoming these hurdles and realizing the full potential of quantum computing in transforming healthcare. As quantum computing continues to evolve, its impact on the future of healthcare could be profound, leading to earlier disease detection, more personalized treatments, and improved patient outcomes. For instance, quantum computing has already been applied to enhance drug discovery processes, with companies like D-Wave Systems (Burnaby, Canada) demonstrating faster molecular simulations for pharmaceutical research and IBM's (Armonk, USA) quantum systems being used to model chemical reactions for new drug development.

Quantum Molecular Docking with a Quantum-Inspired Algorithm

Molecular docking (MD) is a crucial task in drug design, which predicts the position, orientation, and conformation of the ligand when it is bound to a target protein. It can be interpreted as a combinatorial optimization problem, where quantum annealing (QA) has shown a promising advantage for solving combinatorial optimization. In this work, we propose a novel quantum molecular docking (QMD) approach based on a QA-inspired algorithm. We construct two binary encoding methods to efficiently discretize the degrees of freedom with an exponentially reduced number of bits and propose a smoothing filter to rescale the rugged objective function. We propose a new quantum-inspired algorithm, hopscotch simulated bifurcation (hSB), showing great advantages in optimizing over extremely rugged energy landscapes. This hSB can be applied to any formulation of an objective function under binary variables. An adaptive local continuous search is also introduced for further optimization of the discretized solution from hSB. Concerning the stability of docking, we propose a perturbation detection method to help rank the candidate poses. We demonstrate our approach on a typical data set. QMD has shown advantages over the search-based Autodock Vina and the deep-learning DIFFDOCK in both redocking and self-docking scenarios. These results indicate that quantum-inspired algorithms can be applied to solve practical problems in drug discovery even before quantum hardware become mature.

Quantum computing in bioinformatics: a systematic review mapping

The field of quantum computing (QC) is expanding, with efforts being made to apply it to areas previously covered by classical algorithms and methods. Bioinformatics is one such domain that is developing in terms of QC. This article offers a broad mapping review of methods and algorithms of QC in bioinformatics, marking the first of its kind. It presents an overview of the domain and aids researchers in identifying further research directions in the early stages of this field of knowledge. The work presented here shows the current state-of-the-art solutions, focuses on general future directions, and highlights the limitations of current methods. The gathered data includes a comprehensive list of identified methods along with descriptions, classifications, and elaborations of their advantages and disadvantages. Results are presented not just in a descriptive table but also in an aggregated and visual format.

A hybrid quantum computing pipeline for real world drug discovery

Quantum computing, with its superior computational capabilities compared to classical approaches, holds the potential to revolutionize numerous scientific domains, including pharmaceuticals. However, the application of quantum computing for drug discovery has primarily been limited to proof-of-concept studies, which often fail to capture the intricacies of real-world drug development challenges. In this study, we diverge from conventional investigations by developing a hybrid quantum computing pipeline tailored to address genuine drug design problems. Our approach underscores the application of quantum computation in drug discovery and propels it towards more scalable system. We specifically construct our versatile quantum computing pipeline to address two critical tasks in drug discovery: the precise determination of Gibbs free energy profiles for prodrug activation involving covalent bond cleavage, and the accurate simulation of covalent bond interactions. This work serves as a pioneering effort in benchmarking quantum computing against veritable scenarios encountered in drug design, especially the covalent bonding issue present in both of the case studies, thereby transitioning from theoretical models to tangible applications. Our results demonstrate the potential of a quantum computing pipeline for integration into real world drug design workflows.

From GPUs to AI and quantum: three waves of acceleration in bioinformatics

The enormous growth in the amount of data generated by the life sciences is continuously shifting the field from model-driven science towards data-driven science. The need for efficient processing has led to the adoption of massively parallel accelerators such as graphics processing units (GPUs). Consequently, the development of bioinformatics methods nowadays often heavily depends on the effective use of these powerful technologies. Furthermore, progress in computational techniques and architectures continues to be highly dynamic, involving novel deep neural network models and artificial intelligence (AI) accelerators, and potentially quantum processing units in the future. These are expected to be disruptive for the life sciences as a whole and for drug discovery in particular. Here, we identify three waves of acceleration and their applications in a bioinformatics context: (i) GPU computing, (ii) AI and (iii) next-generation quantum computers.

A Vision for the Future of Multiscale Modeling

The rise of modern computer science enabled physical chemistry to make enormous progresses in understanding and harnessing natural and artificial phenomena. Nevertheless, despite the advances achieved over past decades, computational resources are still insufficient to thoroughly simulate extended systems from first principles. Indeed, countless biological, catalytic and photophysical processes require ab initio treatments to be properly described, but the breadth of length and time scales involved makes it practically unfeasible. A way to address these issues is to couple theories and algorithms working at different scales by dividing the system into domains treated at different levels of approximation, ranging from quantum mechanics to classical molecular dynamics, even including continuum electrodynamics. This approach is known as multiscale modeling and its use over the past 60 years has led to remarkable results. Considering the rapid advances in theory, algorithm design, and computing power, we believe multiscale modeling will massively grow into a dominant research methodology in the forthcoming years. Hereby we describe the main approaches developed within its realm, highlighting their achievements and current drawbacks, eventually proposing a plausible direction for future developments considering also the emergence of new computational techniques such as machine learning and quantum computing. We then discuss how advanced multiscale modeling methods could be exploited to address critical scientific challenges, focusing on the simulation of complex light-harvesting processes, such as natural photosynthesis. While doing so, we suggest a cutting-edge computational paradigm consisting in performing simultaneous multiscale calculations on a system allowing the various domains, treated with appropriate accuracy, to move and extend while they properly interact with each other. Although this vision is very ambitious, we believe the quick development of computer science will lead to both massive improvements and widespread use of these techniques, resulting in enormous progresses in physical chemistry and, eventually, in our society.

Hückel molecular orbital theory on a quantum computer: A scalable system-agnostic variational implementation with compact encoding

Hückel molecular orbital (HMO) theory provides a semi-empirical treatment of the electronic structure in conjugated π-electronic systems. A scalable system-agnostic execution of HMO theory on a quantum computer is reported here based on a variational quantum deflation (VQD) algorithm for excited state quantum simulation. A compact encoding scheme is proposed here that provides an exponential advantage over the direct mapping and allows for quantum simulation of the HMO model for systems with up to 2n conjugated centers with n qubits. The transformation of the Hückel Hamiltonian to qubit space is achieved by two different strategies: an iterative refinement transformation and the Frobenius-inner-product-based transformation. These methods are tested on a series of linear, cyclic, and hetero-nuclear conjugated π-electronic systems. The molecular orbital energy levels and wavefunctions from the quantum simulation are in excellent agreement with the exact classical results. However, the higher excited states of large systems are found to suffer from error accumulation in the VQD simulation. This is mitigated by formulating a variant of VQD that exploits the symmetry of the Hamiltonian. This strategy has been successfully demonstrated for the quantum simulation of C60 fullerene containing 680 Pauli strings encoded on six qubits. The methods developed in this work are easily adaptable to similar problems of different complexity in other fields of research.

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to prevent unphysical rearrangements during geometry relaxation. A simple fixed-atom or "coordinate-lock" approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

Scalable narrow linewidth high power laser for barium ion optical qubits

The linewidth of a laser plays a pivotal role in ensuring the high fidelity of ion trap quantum processors and optical clocks. As quantum computing endeavors scale up in qubit number, the demand for higher laser power with ultra-narrow linewidth becomes imperative, and leveraging fiber amplifiers emerges as a promising approach to meet these requirements. This study explores the effectiveness of thulium-doped fiber amplifiers (TDFAs) as a viable solution for addressing optical qubit transitions in trapped barium ion qubits. We demonstrate that by performing high-fidelity gates on the qubit while introducing minimal intensity noise, TDFAs do not significantly broaden the linewidth of the seed lasers. We employed a Voigt fitting scheme in conjunction with a delayed self-heterodyne method to accurately measure the linewidth independently, corroborating our findings through quadrupole spectroscopy with trapped barium ions. Our results show linewidth values of 160 ± 15 Hz and 156 ± 16 Hz, respectively, using these two methods, underscoring the reliability of our measurement techniques. The slight variation within the error-bars of the two methods can be attributed to factors such as amplified spontaneous emission in the TDFA or the influence of 1/f noise within the heterodyne setup delay line. These contribute to advancing our understanding of laser linewidth control in the context of ion trap quantum computing as well as stretching the availability of narrow linewidth, high-power tunable lasers beyond the C-band.

Solvent Distribution Effects on Quantum Chemical Calculations with Quantum Computers

We present a combination of three-dimensional reference interaction site model self-consistent field (3D-RISM-SCF) theory and the variational quantum eigensolver (VQE) to consider the solvent distribution effects within the framework of quantum-classical hybrid computing. The present method, 3D-RISM-VQE, does not include any statistical errors from the solvent configuration sampling owing to the analytical treatment of the statistical solvent distribution. We apply 3D-RISM-VQE to compute the spatial distribution functions of solvent water around a water molecule, the potential and Helmholtz energy curves of NaCl, and to analyze the Helmholtz energy component and related properties of H2O and NH4+. Moreover, we utilize 3D-RISM-VQE to analyze the extent to which solvent effects alter the efficiency of quantum calculations compared with calculations in the gas phase using the L1-norms of molecular electronic Hamiltonians. Our results demonstrate that the efficiency of quantum chemical calculations on a quantum computer in solution is virtually the same as that in the gas phase.

Distributed Implementation of Full Configuration Interaction for One Trillion Determinants

Full configuration interaction (FCI) can provide an exact molecular ground-state energy within a given basis set and serve as a benchmark for approximate methods in quantum chemical calculations, including the emerging variational quantum eigensolver. However, its exponential computational and memory requirements easily exceed the capability of a single server and limit its applicability to large molecules. In this paper, we present a distributed FCI implementation employing a hybrid parallelization scheme with multithreading and multiprocessing to expand FCI's applicability. We optimize this scheme to minimize the bottlenecks arising from interprocess communications and interthread data management. Our implementation achieves higher scalability than the naive combination of prior works and successfully calculates the exact energy of C3H8/STO-3G with 1.31 trillion determinants, which is the largest FCI calculation to the best of our knowledge. Furthermore, we provide a comprehensive list of FCI results with 136 combinations of molecules and basis sets for future evaluation and development of approximate methods.

Quantum Computing in the Next-Generation Computational Biology Landscape: From Protein Folding to Molecular Dynamics

Modern biological science is trying to solve the fundamental complex problems of molecular biology, which include protein folding, drug discovery, simulation of macromolecular structure, genome assembly, and many more. Currently, quantum computing (QC), a rapidly emerging technology exploiting quantum mechanical phenomena, has developed to address current significant physical, chemical, biological issues, and complex questions. The present review discusses quantum computing technology and its status in solving molecular biology problems, especially in the next-generation computational biology scenario. First, the article explained the basic concept of quantum computing, the functioning of quantum systems where information is stored as qubits, and data storage capacity using quantum gates. Second, the review discussed quantum computing components, such as quantum hardware, quantum processors, and quantum annealing. At the same time, article also discussed quantum algorithms, such as the grover search algorithm and discrete and factorization algorithms. Furthermore, the article discussed the different applications of quantum computing to understand the next-generation biological problems, such as simulation and modeling of biological macromolecules, computational biology problems, data analysis in bioinformatics, protein folding, molecular biology problems, modeling of gene regulatory networks, drug discovery and development, mechano-biology, and RNA folding. Finally, the article represented different probable prospects of quantum computing in molecular biology.

Applications of noisy quantum computing and quantum error mitigation to "adamantaneland": a benchmarking study for quantum chemistry

The field of quantum computing has the potential to transform quantum chemistry. The variational quantum eigensolver (VQE) algorithm has allowed quantum computing to be applied to chemical problems in the noisy intermediate-scale quantum (NISQ) era. Applications of VQE have generally focused on predicting absolute energies instead of chemical properties that are relative energy differences and that are most interesting to chemists studying a chemical problem. We address this shortcoming by constructing a molecular benchmark data set in this work containing isomers of C10H16 and carbocationic rearrangements of C10H15+, calculated at a high-level of theory. Using the data set, we compared noiseless VQE simulations to conventionally performed density functional and wavefunction theory-based methods to understand the quality of results. We also investigated the effectiveness of a quantum state tomography-based error mitigation technique in applications of VQE under noise (simulated and real). Our findings reveal that the use of quantum error mitigation is crucial in the NISQ era and advantageous to yield almost noiseless quality results.

Visioning synthetic futures for yeast research within the context of current global techno-political trends

Yeast research is entering into a new period of scholarship, with new scientific tools, new questions to ask and new issues to consider. The politics of emerging and critical technology can no longer be separated from the pursuit of basic science in fields, such as synthetic biology and engineering biology. Given the intensifying race for technological leadership, yeast research is likely to attract significant investment from government, and that it offers huge opportunities to the curious minded from a basic research standpoint. This article provides an overview of new directions in yeast research with a focus on Saccharomyces cerevisiae, and places these trends in their geopolitical context. At the highest level, yeast research is situated within the ongoing convergence of the life sciences with the information sciences. This convergent effect is most strongly pronounced in areas of AI-enabled tools for the life sciences, and the creation of synthetic genomes, minimal genomes, pan-genomes, neochromosomes and metagenomes using computer-assisted design tools and methodologies. Synthetic yeast futures encompass basic and applied science questions that will be of intense interest to government and nongovernment funding sources. It is essential for the yeast research community to map and understand the context of their research to ensure their collaborations turn global challenges into research opportunities.

Selectivity Mechanism of Pyrrolopyridone Analogues Targeting Bromodomain 2 of Bromodomain-Containing Protein 4 from Molecular Dynamics Simulations

Bromodomain and extra-terminal domain (BET) proteins play an important role in epigenetic regulation and are linked to several diseases; therefore, they are interesting targets. BET has two bromodomains: bromodomain 1 (BD1) and BD2. Selective targeting of BD1 or BD2 may produce different activities and greater effects than pan-BD inhibitors. However, the selective mechanism of the specific core must be studied at the atomic level. This study determined the effectiveness of pyrrolopyridone analogues to selectively inhibit BD2 using a pan-BD inhibitor (ABBV-075) and a selective-BD2 inhibitor (ABBV-744). Molecular dynamics simulations and calculations of binding free energies were used to systematically study the selectivity of BD2 inhibition by the pyrrolopyridone analogues. Overall, the pyrrolopyridone analogue inhibitors targeting BD2 interacted mainly with the following amino acid pairs between bromodomain-containing protein 4 (BRD4)-BD1 and BRD4-BD2 complexes: I146/V439, N140/N433, D144/H437, P82/P375, V87/V380, D88/D381, and Y139/Y432. The pyrrolopyridone analogues targeting BRD4-BD2 were divided into five regions based on selectivity mechanism. These results suggest that the R3 and R5 regions of pyrrolopyridone analogues can be modified to improve the selectivity between BRD4-BD1 and BRD4-BD2. The selectivity of BD2 inhibition by pyrrolopyridone analogues can be used to design novel BD2 inhibitors based on a pyrrolopyridone core.

Quantum Computer Simulation of Protein Protonation

In attempts to simulate the protonation of proteins, a major challenge is that the number of protonation states grows rapidly as a function (2N) of the number of protonation sites (N). Expression on the free energy of the protonation state as an N-site Ising model ─ using an empirical Generalized-Born model ─ allows a quantum computer to efficiently determine the important states at a given pH value and subsequently reconstruct the pH titration process at all sites. Compared with the exact results painstakingly obtained with classical computers, the results obtained using quantum computers show good agreement for staphylococcal nuclease and excellent agreement for α-lactalbumin. This work illustrates the effectiveness of quantum computers in sampling important physical states, which may be useful in attacking challenging biomolecular problems.