ترغب بنشر مسار تعليمي؟ اضغط هنا

Novel algorithms and high-performance cloud computing enable efficient fully quantum mechanical protein-ligand scoring

48   0   0.0 ( 0 )
 نشر من قبل Narbe Mardirossian
 تاريخ النشر 2020
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

Ranking the binding of small molecules to protein receptors through physics-based computation remains challenging. Though inroads have been made using free energy methods, these fail when the underlying classical mechanical force fields are insufficient. In principle, a more accurate approach is provided by quantum mechanical density functional theory (DFT) scoring, but even with approximations, this has yet to become practical on drug discovery-relevant timescales and resources. Here, we describe how to overcome this barrier using algorithms for DFT calculations that scale on widely available cloud architectures, enabling full density functional theory, without approximations, to be applied to protein-ligand complexes with approximately 2500 atoms in tens of minutes. Applying this to a realistic example of 22 ligands binding to MCL1 reveals that density functional scoring outperforms classical free energy perturbation theory for this system. This raises the possibility of broadly applying fully quantum mechanical scoring to real-world drug discovery pipelines.



قيم البحث

اقرأ أيضاً

We predict the enhanced light harvesting of a protein-pigment complex when assembled to a quantum dot (QD) antenna. Our prototypical nanoassembly setup is composed of a Fenna-Mattews-Olson system hosting 8 Bacteriochlorophyll (BChl) a dyes, and a nea r-infrared emitting CdSe$_x$Te$_{(1-x)}$/ZnS alloy-core/shell nanocrystal. BChl a has two wide windows of poor absorption in the green and orange-red bands, precisely where most of the sunlight energy lies. The selected QD is able to collect sunlight efficiently in a broader band and funnel its energy by a (non-radiative) Forster resonance energy transfer mechanism to the dyes embedded in the protein. By virtue of the coupling between the QD and the dyes, the nanoassembly absorption is dramatically improved in the poor absorption window of the BChl a.
We describe the development of a scientific cloud computing (SCC) platform that offers high performance computation capability. The platform consists of a scientific virtual machine prototype containing a UNIX operating system and several materials s cience codes, together with essential interface tools (an SCC toolset) that offers functionality comparable to local compute clusters. In particular, our SCC toolset provides automatic creation of virtual clusters for parallel computing, including tools for execution and monitoring performance, as well as efficient I/O utilities that enable seamless connections to and from the cloud. Our SCC platform is optimized for the Amazon Elastic Compute Cloud (EC2). We present benchmarks for prototypical scientific applications and demonstrate performance comparable to local compute clusters. To facilitate code execution and provide user-friendly access, we have also integrated cloud computing capability in a JAVA-based GUI. Our SCC platform may be an alternative to traditional HPC resources for materials science or quantum chemistry applications.
Computational drug discovery provides an efficient tool helping large scale lead molecules screening. One of the major tasks of lead discovery is identifying molecules with promising binding affinities towards a target, a protein in general. The accu racies of current scoring functions which are used to predict the binding affinity are not satisfactory enough. Thus, machine learning (ML) or deep learning (DL) based methods have been developed recently to improve the scoring functions. In this study, a deep convolutional neural network (CNN) model (called OnionNet) is introduced and the features are based on rotation-free element-pair specific contacts between ligands and protein atoms, and the contacts were further grouped in different distance ranges to cover both the local and non-local interaction information between the ligand and the protein. The prediction power of the model is evaluated and compared with other scoring functions using the comparative assessment of scoring functions (CASF-2013) benchmark and the v2016 core set of PDBbind database. When compared to a previous CNN-based scoring function, our model shows improvements of 0.08 and 0.16 in the correlations (R) and standard deviations (SD) of regression, respectively, between the predicted binding affinities and the experimental measured binding affinities. The robustness of the model is further explored by predicting the binding affinities of the complexes generated from docking simulations instead of experimentally determined PDB structures.
The conventional technique for solving the equations of quantum chemistry (of solid state) is extended unconventionally to the structures possessing certain symmetries. This proposal concerns changing the way for selection of occupied orbitals, allow ing, in turn, to release the unoccupied electronic states located lower than the ground state Fermi level of a specific system. Such states can be treated as spectral holes. Application of this technique, in particular, when calculating the electronic structure of the HTSC-compound $YBa_2Cu_3O_{7-x}$ (0<x<1) results in the following. The spectral holes of high spatial localization are found. These spatial spectral holes are located, mainly at the $P_y$-orbitals of the apex oxygens. These orbitals overlap and form linear chains which are parallel to the known Cu(1)-O chains, disappearing when x is closed to 1. One can suppose that the linear chains of the overlapping hole states form a superconducting channel. Some other parameters closely related to the critical characteristics of HTSC-materials are also calculated. The calculations show that the superconducting channel is broken when the oxygen chain atoms O(1) are removed (x>0). One could easily connect the obtained results to the high-temperature superconductivity of Littles linear chains, as well as Ginzburgs two-dimensional layers and even to BCS-model.
We present the implementation of an electronic-structure approach dedicated to ionization dynamics of molecules interacting with x-ray free-electron laser (XFEL) pulses. In our scheme, molecular orbitals for molecular core-hole states are represented by linear combination of numerical atomic orbitals that are solutions of corresponding atomic core-hole states. We demonstrate that our scheme efficiently calculates all possible multiple-hole configurations of molecules formed during XFEL pulses. The present method is suitable to investigate x-ray multiphoton multiple ionization dynamics and accompanying nuclear dynamics, providing essential information on the chemical dynamics relevant for high-intensity x-ray imaging.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا