Design, Synthesis, and Evaluation of New Hybrid Derivatives of 5,6-Dihydro-4H-pyrrolo[3,2,1-ij]quinolin-2(1H)-one as Potential Dual Inhibitors of Blood Coagulation Factors Xa and XIa.

Cardiovascular diseases caused by blood coagulation system disorders are one of the leading causes of morbidity and mortality in the world. Research shows that blood clotting factors are involved in these thrombotic processes. Among them, factor Xa occupies a key position in the blood coagulation cascade. Another coagulation factor, XIa, is also a promising target because its inhibition can suppress thrombosis with a limited contribution to normal hemostasis. In this regard, the development of dual inhibitors as new generation anticoagulants is an urgent problem. Here we report the synthesis and evaluation of novel potential dual inhibitors of coagulation factors Xa and XIa. Based on the principles of molecular design, we selected a series of compounds that combine in their structure fragments of pyrrolo[3,2,1-ij]quinolin-2-one and thiazole, connected through a hydrazine linker. The production of new hybrid molecules was carried out using a two-stage method. The reaction of 5,6-dihydropyrrolo[3,2,1-ij]quinoline-1,2-diones with thiosemicarbazide gave the corresponding hydrazinocarbothioamides. The reaction of the latter with DMAD led to the target methyl 2-(4-oxo-2-(2-(2-oxo-5,6-dihydro-4H-pyrrolo[3,2,1-ij]quinolin-1(2H)-ylidene)hydrazineyl)thiazol-5(4H)-ylidene)acetates in high yields. In vitro testing of the synthesized molecules revealed that ten of them showed high inhibition values for both the coagulation factors Xa and XIa, and the IC50 value for some compounds was also assessed. The resulting structures were also tested for their ability to inhibit thrombin.

New Hybrid Tetrahydropyrrolo[3,2,1-ij]quinolin-1-ylidene-2-thioxothiazolidin-4-ones as New Inhibitors of Factor Xa and Factor XIa: Design, Synthesis, and In Silico and Experimental Evaluation.

Despite extensive research in the field of thrombotic diseases, the prevention of blood clots remains an important area of study. Therefore, the development of new anticoagulant drugs with better therapeutic profiles and fewer side effects to combat thrombus formation is still needed. Herein, we report the synthesis and evaluation of novel pyrroloquinolinedione-based rhodanine derivatives, which were chosen from 24 developed derivatives by docking as potential molecules to inhibit the clotting factors Xa and XIa. For the synthesis of new hybrid derivatives of pyrrolo[3,2,1-ij]quinoline-2-one, we used a convenient structural modification of the tetrahydroquinoline fragment by varying the substituents in positions 2, 4, and 6. In addition, the design of target molecules was achieved by alkylating the amino group of the rhodanine fragment with propargyl bromide or by replacing the rhodanine fragment with 2-thioxoimidazolidin-4-one. The in vitro testing showed that eight derivatives are capable of inhibiting both coagulation factors, two compounds are selective inhibitors of factor Xa, and two compounds are selective inhibitors of factor XIa. Overall, these data indicate the potential anticoagulant activity of these molecules through the inhibition of the coagulation factors Xa and XIa.

Novel Inhibitors of 2'-O-Methyltransferase of the SARS-CoV-2 Coronavirus.

The COVID-19 pandemic is still affecting many people worldwide and causing a heavy burden to global health. To eliminate the disease, SARS-CoV-2, the virus responsible for the pandemic, can be targeted in several ways. One of them is to inhibit the 2'-O-methyltransferase (nsp16) enzyme that is crucial for effective translation of viral RNA and virus replication. For methylation of substrates, nsp16 utilizes S-adenosyl methionine (SAM). Binding of a small molecule in the protein site where SAM binds can disrupt the synthesis of viral proteins and, as a result, the replication of the virus. Here, we performed high-throughput docking into the SAM-binding site of nsp16 for almost 40 thousand structures, prepared for compounds from three libraries: Enamine Coronavirus Library, Enamine Nucleoside Mimetics Library, and Chemdiv Nucleoside Analogue Library. For the top scoring ligands, semi-empirical quantum-chemical calculations were performed, to better estimate protein-ligand binding enthalpy. Relying upon the calculated binding energies and predicted docking poses, we selected 21 compounds for experimental testing.

New Blood Coagulation Factor XIIa Inhibitors: Molecular Modeling, Synthesis, and Experimental Confirmation.

In the modern world, complications caused by disorders in the blood coagulation system are found in almost all areas of medicine. Thus, the development of new, more advanced drugs that can prevent pathological conditions without disrupting normal hemostasis is an urgent task. The blood coagulation factor XIIa is one of the most promising therapeutic targets for the development of anticoagulants based on its inhibitors. The initial stage of drug development is directly related to computational methods of searching for a lead compound. In this study, docking followed by quantum chemical calculations was used to search for noncovalent low-molecular-weight factor XIIa inhibitors in a focused library of druglike compounds. As a result of the study, four low-molecular-weight compounds were experimentally confirmed as factor XIIa inhibitors. Selectivity testing revealed that two of the identified factor XIIa inhibitors were selective over the coagulation factors Xa and XIa.

New Aspects of the Reaction of Thioacetamide and N-Substituted Maleimides.

N-Arylmaleimides are universal substrates for the synthesis of various heterocyclic compounds with a wide spectrum of biological activity. However, their reactions with thioacetamides have not been comprehensively studied. We studied the reactions of thioacetamide with N-arylmaleimides under various conditions. We established for the first time that three types of products: epithiopyrrolo[3,4-c]pyridines, pyrrolo[3,4-c]pyridines and 3,3'-thiobis(1-arylpyrrolidine-2,5-diones) can be obtained in different conditions. In all cases, two maleimide molecules are involved in the reaction. 3,3'-Thiobis(1-arylpyrrolidine-2,5-diones) are the major products when the reaction is conducted at boiling in acetic acid. When thioacetamide and N-arylmaleimide are kept in dioxane at 50 °C, epithiopyrrolo[3,4-c]pyridines can be isolated, which, when heated in dioxane, in acetic acid or in methanol in the presence of catalytic amounts of sodium methoxide, are converted into pyrrolo[3,4-c]pyridines by eliminating hydrogen sulfide. The reaction of thioacetamide and N-arylmaleimide in dioxane at boiling temperature with the portioned addition of N-arylmaleimide leads predominantly to the formation of pyrrolo[3,4-c]pyridines. The reaction of thioacetamide with N-alkylmaleimides under all the above conditions leads predominantly to the formation of the corresponding sulfides. The structure of the compounds obtained was characterized by a set of spectral analysis methods and X-ray diffraction (XRD) data.

New Chemicals Suppressing SARS-CoV-2 Replication in Cell Culture.

Candidates to being inhibitors of the main protease (Mpro) of SARS-CoV-2 were selected from the database of Voronezh State University using molecular modeling. The database contained approximately 19,000 compounds represented by more than 41,000 ligand conformers. These ligands were docked into Mpro using the SOL docking program. For one thousand ligands with best values of the SOL score, the protein-ligand binding enthalpy was calculated by the PM7 quantum-chemical method with the COSMO solvent model. Using the SOL score and the calculated protein-ligand binding enthalpies, eighteen compounds were selected for the experiments. Several of these inhibitors suppressed the replication of the coronavirus in cell culture, and we used the best three among them in the search for chemical analogs. Selection among analogs using the same procedure followed by experiments led to identification of seven inhibitors of the SARS-CoV-2 replication in cell culture with EC50 values at the micromolar level. The identified inhibitors belong to three chemical classes. The three inhibitors, 4,4-dimethyldithioquinoline derivatives, inhibit SARS-CoV-2 replication in Vero E6 cell culture just as effectively as the best published non-covalent inhibitors, and show low cytotoxicity. These results open up a possibility to develop antiviral drugs against the SARS-CoV-2 coronavirus.

Synthesis, Docking, and In Vitro Anticoagulant Activity Assay of Hybrid Derivatives of Pyrrolo[3,2,1-ij]Quinolin-2(1H)-one as New Inhibitors of Factor Xa and Factor XIa.

Coagulation factor Xa and factor XIa are proven to be convenient and crucial protein targets for treatment for thrombotic disorders and thereby their inhibitors can serve as effective anticoagulant drugs. In the present work, we focused on the structure-activity relationships of derivatives of pyrrolo[3,2,1-ij]quinolin-2(1H)-one and an evaluation of their activity against factor Xa and factor XIa. For this, docking-guided synthesis of nine compounds based on pyrrolo[3,2,1-ij]quinolin-2(1H)-one was carried out. For the synthesis of new hybrid hydropyrrolo[3,2,1-ij]quinolin-2(1H)-one derivatives, we used convenient structural modification of both the tetrahydro- and dihydroquinoline moiety by varying the substituents at the C6,8,9 positions. In vitro testing revealed that four derivatives were able to inhibit both coagulation factors and three compounds were selective factor XIa inhibitors. An IC50 value of 3.68 µM for was found for the best factor Xa inhibitor and 2 µM for the best factor XIa inhibitor.

Chemical Adjustment of Fibrinolysis.

Fibrinolysis is the process of the fibrin-platelet clot dissolution initiated after bleeding has been stopped. It is regulated by a cascade of proteolytic enzymes with plasmin at its core. In pathological cases, the balance of normal clot formation and dissolution is replaced by a too rapid lysis, leading to bleeding, or an insufficient one, leading to an increased thrombotic risk. The only approved therapy for emergency thrombus lysis in ischemic stroke is recombinant tissue plasminogen activator, though streptokinase or urokinase-type plasminogen activators could be used for other conditions. Low molecular weight compounds are of great interest for long-term correction of fibrinolysis dysfunctions. Their areas of application might go beyond the hematology field because the regulation of fibrinolysis could be important in many conditions, such as fibrosis. They enhance or weaken fibrinolysis without significant effects on other components of hemostasis. Here we will describe and discuss the main classes of these substances and their mechanisms of action. We will also explore avenues of research for the development of new drugs, with a focus on the use of computational models in this field.

Application of the docking program SOL for CSAR benchmark.

This paper is devoted to results obtained by the docking program SOL and the post-processing program DISCORE at the CSAR benchmark. SOL and DISCORE programs are described. SOL is the original docking program developed on the basis of the genetic algorithm, MMFF94 force field, rigid protein, precalculated energy grid including desolvation in the frame of simplified GB model, vdW, and electrostatic interactions and taking into account the ligand internal strain energy. An important SOL feature is the single- or multi-processor performance for up to hundreds of CPUs. DISCORE improves the binding energy scoring by the local energy optimization of the ligand docked pose and a simple linear regression on the base of available experimental data. The docking program SOL has demonstrated a good ability for correct ligand positioning in the active sites of the tested proteins in most cases of CSAR exercises. SOL and DISCORE have not demonstrated very exciting results on the protein-ligand binding free energy estimation. Nevertheless, for some target proteins, SOL and DISCORE were among the first in prediction of inhibition activity. Ways to improve SOL and DISCORE are discussed.

Atomistic Simulation of Physical Vapor Deposition of Optical Thin Films.

A review of the methods and results of atomistic modeling of the deposition of thin optical films and a calculation of their characteristics is presented. The simulation of various processes in a vacuum chamber, including target sputtering and the formation of film layers, is considered. Methods for calculating the structural, mechanical, optical, and electronic properties of thin optical films and film-forming materials are discussed. The application of these methods to studying the dependences of the characteristics of thin optical films on the main deposition parameters is considered. The simulation results are compared with experimental data.

Experimentally Validated Novel Factor XIIa Inhibitors Identified by Docking and Quantum Chemical Post-processing.

Antithrombotic agents based on factor XIIa inhibitors can become a new class of drugs to manage conditions associated with thrombosis. Herein, we report identification of two novel classes of factor XIIa inhibitors. The first one is triazolopyrimidine derivatives designed on the basis of the literature aminotriazole hit and identified using virtual screening of the focused library. The second class is a spirocyclic furo[3,4-c]pyrrole derivatives identified by virtual screening of a large chemical library of drug-like compounds performed in a previous study but confirmed in vitro here. In both cases, the prediction of inhibitory activity is based on the score of the SOL docking program, which uses the MMFF94 force field to calculate the binding energy. For the best ligands selected in virtual screening of the large chemical library, postprocessing with the PM7 semiempirical quantum-chemical method was used to calculate the enthalpy of protein-ligand binding to prioritize 16 compounds for testing in enzymatic assay, and one of them demonstrated micromolar activity. For triazolopyrimidine library, 21 compounds were prioritized for the testing based on docking scores, and visual inspection of docking poses. Of these, 4 compounds showed inhibition of factor XIIa at 30 µM.

On the origin of near-IR luminescence in Bi-doped materials (II). Subvalent monocation Bi⁺ and cluster Bi5³âº luminescence in AlCl3/ZnCl2/BiCl3 chloride glass.

Broadband NIR photoluminescence (from 1000 to 2500 nm) was observed from partially reduced AlCl3/ZnCl2/BiCl3 glass, containing subvalent bismuth species. The luminescence consists of three bands, assigned to Bi⁺ , Bi24⁺, and Bi5³âº ions. The physical and optical characteristics of these centers and possible contribution to NIR luminescence from bismuth-doped oxide glasses are discussed.

Quantum-Chemical Quasi-Docking for Molecular Dynamics Calculations.

The quantum quasi-docking procedure is used to compare the docking accuracies of two quantum-chemical semiempirical methods, namely, PM6-D3H4X and PM7. Quantum quasi-docking is an approximation to quantum docking. In quantum docking, it is necessary to search directly for the global minimum of the energy of the protein-ligand complex calculated by the quantum-chemical method. In quantum quasi-docking, firstly, we look for a wide spectrum of low-energy minima, calculated using the MMFF94 force field, and secondly, we recalculate the energies of all these minima using the quantum-chemical method, and among these recalculated energies we determine the lowest energy and the corresponding ligand position. Both PM6-D3H4X and PM7 are novel methods that describe well-dispersion interactions, hydrogen and halogen bonds. The PM6-D3H4X and PM7 methods are used with the COSMO implicit solvent model as it is implemented in the MOPAC program. The comparison is made for 25 high quality protein-ligand complexes. Firstly, the docking positioning accuracies have been compared, and we demonstrated that PM7+COSMO provides better positioning accuracy than PM6-D3H4X. Secondly, we found that PM7+COSMO demonstrates a much higher correlation between the calculated and measured protein-ligand binding enthalpies than PM6-D3H4X. For future quantum docking PM7+COSMO is preferable, but the COSMO model must be improved.

Molecular Dynamics Simulation of Laser Induced Heating of Silicon Dioxide Thin Films.

The full-atomistic classical molecular dynamics simulation of the laser heating of silicon dioxide thin films is performed. Both dense isotropic films and porous anisotropic films are investigated. It is assumed that heating occurs due to nodal structural defects, which are currently considered one of the possible causes of laser induced damage. It is revealed that heating to a temperature of 1000 K insignificantly affects the structure of the films and the concentration of point defects responsible for the radiation absorption. An increase in the heating temperature to 2000 K leads to the growth of the concentration of these defects. For "as deposited" films, this growth is greater in the case of a porous film deposited at a high deposition angle. Annealing of film reduces the difference in the concentration of laser induced defects in dense and porous films. The possible influence of optical active defects arising due to heating on the laser induced damage threshold is discussed.

Docking Paradigm in Drug Design.

Docking is in demand for the rational computer aided structure based drug design. A review of docking methods and programs is presented. Different types of docking programs are described. They include docking of non-covalent small ligands, protein-protein docking, supercomputer docking, quantum docking, the new generation of docking programs and the application of docking for covalent inhibitors discovery. Taking into account the threat of COVID-19, we present here a short review of docking applications to the discovery of inhibitors of SARS-CoV and SARS-CoV-2 target proteins, including our own result of the search for inhibitors of SARS-CoV-2 main protease using docking and quantum chemical post-processing. The conclusion is made that docking is extremely important in the fight against COVID-19 during the process of development of antivirus drugs having a direct action on SARS-CoV-2 target proteins.

Advances in Docking.

BACKGROUND: Design of small molecules which are able to bind to the protein responsible for a disease is the key step of the entire process of the new medicine discovery. Atomistic computer modeling can significantly improve effectiveness of such design. The accurate calculation of the free energy of binding a small molecule (a ligand) to the target protein is the most important problem of such modeling. Docking is one of the most popular molecular modeling methods for finding ligand binding poses in the target protein and calculating the protein-ligand binding energy. This energy is used for finding the most active compounds for the given target protein. This short review aims to give a concise description of distinctive features of docking programs focusing on computation methods and approximations influencing their accuracy. METHODS: This review is based on the peer-reviewed research literature including author's own publications. The main features of several representative docking programs are briefly described focusing on their characteristics influencing docking accuracy: force fields, energy calculations, solvent models, algorithms of the best ligand pose search, global and local optimizations, ligand and target protein flexibility, and the simplifications made for the docking accelerating. Apart from other recent reviews focused mainly on the performance of different docking programs, in this work, an attempt is made to extract the most important functional characteristics defining the docking accuracy. Also a roadmap for increasing the docking accuracy is proposed. This is based on the new generation of docking programs which have been realized recently. These programs and respective new global optimization algorithms are described shortly. RESULTS: Several popular conventional docking programs are considered. Their search of the best ligand pose is based explicitly or implicitly on the global optimization problem. Several algorithms are used to solve this problem, and among them, the heuristic genetic algorithm is distinguished by its popularity and an elaborate design. All conventional docking programs for their acceleration use the preliminary calculated grids of protein-ligand interaction potentials or preferable points of protein and ligand conjugation. These approaches and commonly used fitting parameters restrict strongly the docking accuracy. Solvent is considered in exceedingly simplified approaches in the course of the global optimization and the search for the best ligand poses. More accurate approaches on the base of implicit solvent models are used frequently for more careful binding energy calculations after docking. The new generation of docking programs are developed recently. They find the spectrum of low energy minima of a protein-ligand complex including the global minimum. These programs should be more accurate because they do not use a preliminary calculated grid of protein-ligand interaction potentials and other simplifications, the energy of any conformation of the molecular system is calculated in the frame of a given force field and there are no fitting parameters. A new docking algorithm is developed and fulfilled specially for the new docking programs. This algorithm allows docking a flexible ligand into a flexible protein with several dozen mobile atoms on the base of the global energy minimum search. Such docking results in improving the accuracy of ligand positioning in the docking process. The adequate choice of the method of molecular energy calculations also results in the better docking positioning accuracy. An advancement in the application of quantum chemistry methods to docking and scoring is revealed. CONCLUSION: The findings of this review confirm the great demand in docking programs for discovery of new medicine substances with the help of molecular modeling. New trends in docking programs design are revealed. These trends are focused on the increase of the docking accuracy at the expense of more accurate molecular energy calculations without any fitting parameters, including quantum-chemical methods and implicit solvent models, and by using new global optimization algorithms which make it possible to treat flexibility of ligands and mobility of protein atoms simultaneously. Finally, it is shown that all the necessary prerequisites for increasing the docking accuracy can be accomplished in practice.

New factor Xa inhibitors based on 1,2,3,4-tetrahydroquinoline developed by molecular modelling.

Factor Xa is a serine protease representing a crucial element in the coagulation process and an attractive target for anticoagulant therapy. At the present time there are several chemical classes of factor Xa inhibitors with proven activity. Furthermore, three factor Xa inhibitors have been approved for the medical use to date. However, therapy with these medications is accompanied by substantial adverse effects. In this background, the structure-based computational approach combining molecular docking and semiempirical quantum chemical calculations was applied for a search for new effective factor Xa inhibitors. We have undertaken a few virtual screening procedures to select potential candidates for synthesis and subsequent testing. The first screen of the focused library resulted in identifying 20 compounds among which 7 compounds showed the noticeable inhibition of factor Xa at maximal concentrations, allowed by solubility. The subsequent additional screens identified 20 additional candidates. Of these, 5 substances were shown to be capable of inhibiting factor Xa at 5 µM. The best two found 1,2,3,4-tetrahydroquinoline derivatives identified by means of modelling have demonstrated IC50 values in the micromolar range. One of them turned out to be selective factor Xa inhibitor over trypsin, factors IIa, IXa and XIa.

New generation of docking programs: Supercomputer validation of force fields and quantum-chemical methods for docking.

Discovery of new inhibitors of the protein associated with a given disease is the initial and most important stage of the whole process of the rational development of new pharmaceutical substances. New inhibitors block the active site of the target protein and the disease is cured. Computer-aided molecular modeling can considerably increase effectiveness of new inhibitors development. Reliable predictions of the target protein inhibition by a small molecule, ligand, is defined by the accuracy of docking programs. Such programs position a ligand in the target protein and estimate the protein-ligand binding energy. Positioning accuracy of modern docking programs is satisfactory. However, the accuracy of binding energy calculations is too low to predict good inhibitors. For effective application of docking programs to new inhibitors development the accuracy of binding energy calculations should be higher than 1kcal/mol. Reasons of limited accuracy of modern docking programs are discussed. One of the most important aspects limiting this accuracy is imperfection of protein-ligand energy calculations. Results of supercomputer validation of several force fields and quantum-chemical methods for docking are presented. The validation was performed by quasi-docking as follows. First, the low energy minima spectra of 16 protein-ligand complexes were found by exhaustive minima search in the MMFF94 force field. Second, energies of the lowest 8192 minima are recalculated with CHARMM force field and PM6-D3H4X and PM7 quantum-chemical methods for each complex. The analysis of minima energies reveals the docking positioning accuracies of the PM7 and PM6-D3H4X quantum-chemical methods and the CHARMM force field are close to one another and they are better than the positioning accuracy of the MMFF94 force field.

Evaluation of the novel algorithm of flexible ligand docking with moveable target-protein atoms.

We present the novel docking algorithm based on the Tensor Train decomposition and the TT-Cross global optimization. The algorithm is applied to the docking problem with flexible ligand and moveable protein atoms. The energy of the protein-ligand complex is calculated in the frame of the MMFF94 force field in vacuum. The grid of precalculated energy potentials of probe ligand atoms in the field of the target protein atoms is not used. The energy of the protein-ligand complex for any given configuration is computed directly with the MMFF94 force field without any fitting parameters. The conformation space of the system coordinates is formed by translations and rotations of the ligand as a whole, by the ligand torsions and also by Cartesian coordinates of the selected target protein atoms. Mobility of protein and ligand atoms is taken into account in the docking process simultaneously and equally. The algorithm is realized in the novel parallel docking SOL-P program and results of its performance for a set of 30 protein-ligand complexes are presented. Dependence of the docking positioning accuracy is investigated as a function of parameters of the docking algorithm and the number of protein moveable atoms. It is shown that mobility of the protein atoms improves docking positioning accuracy. The SOL-P program is able to perform docking of a flexible ligand into the active site of the target protein with several dozens of protein moveable atoms: the native crystallized ligand pose is correctly found as the global energy minimum in the search space with 157 dimensions using 4700 CPU ∗ h at the Lomonosov supercomputer.

Evaluation of Docking Target Functions by the Comprehensive Investigation of Protein-Ligand Energy Minima.

The adequate choice of the docking target function impacts the accuracy of the ligand positioning as well as the accuracy of the protein-ligand binding energy calculation. To evaluate a docking target function we compared positions of its minima with the experimentally known pose of the ligand in the protein active site. We evaluated five docking target functions based on either the MMFF94 force field or the PM7 quantum-chemical method with or without implicit solvent models: PCM, COSMO, and SGB. Each function was tested on the same set of 16 protein-ligand complexes. For exhaustive low-energy minima search the novel MPI parallelized docking program FLM and large supercomputer resources were used. Protein-ligand binding energies calculated using low-energy minima were compared with experimental values. It was demonstrated that the docking target function on the base of the MMFF94 force field in vacuo can be used for discovery of native or near native ligand positions by finding the low-energy local minima spectrum of the target function. The importance of solute-solvent interaction for the correct ligand positioning is demonstrated. It is shown that docking accuracy can be improved by replacement of the MMFF94 force field by the new semiempirical quantum-chemical PM7 method.