Twenty years ago, the landmark AM1 was introduced, and has since had an increasingly wide following among chemists due to its consistently good results and time-tested reliability--being presently available in countless computational quantum chemistry programs. However, semiempirical molecular orbital models still are of limited accuracy and need to be improved if the full potential of new linear scaling techniques, such as MOZYME and LocalSCF, is to be realized. Accordingly, in this article we present RM1 (Recife Model 1): a reparameterization of AM1. As before, the properties used in the parameterization procedure were: heats of formation, dipole moments, ionization potentials and geometric variables (bond lengths and angles). Considering that the vast majority of molecules of importance to life can be assembled by using only six elements: C, H, N, O, P, and S, and that by adding the halogens we can now build most molecules of importance to pharmaceutical research, our training set consisted of 1736 molecules, representative of organic and biochemistry, containing C, H, N, O, P, S, F, Cl, Br, and I atoms. Unlike AM1, and similar to PM3, all RM1 parameters have been optimized. For enthalpies of formation, dipole moments, ionization potentials, and interatomic distances, the average errors in RM1, for the 1736 molecules, are less than those for AM1, PM3, and PM5. Indeed, the average errors in kcal x mol(-1) of the enthalpies of formation for AM1, PM3, and PM5 are 11.15, 7.98, and 6.03, whereas for RM1 this value is 5.77. The errors, in Debye, of the dipole moments for AM1, PM3, PM5, and RM1 are, respectively, 0.37, 0.38, 0.50, and 0.34. Likewise, the respective errors for the ionization potentials, in eV, are 0.60, 0.55, 0.48, and 0.45, and the respective errors, in angstroms, for the interatomic distances are 0.036, 0.029, 0.037, and 0.027. The RM1 average error in bond angles of 6.82 degrees is only slightly higher than the AM1 figure of 5.88 degrees, and both are much smaller than the PM3 and PM5 figures of 6.98 degrees and 9.83 degrees, respectively. Moreover, a known error in PM3 nitrogen charges is corrected in RM1. Therefore, RM1 represents an improvement over AM1 and its similar successor PM3, and is probably very competitive with PM5, which is a somewhat different model, and not fully disclosed. RM1 possesses the same analytical construct and the same number of parameters for each atom as AM1, and, therefore, can be easily implemented in any software that already has AM1, not requiring any change in any line of code, with the sole exception of the values of the parameters themselves.
Our previously defined Sparkle model (Inorg. Chem. 2004, 43, 2346) has been reparameterized for Eu(III) as well as newly parameterized for Gd(III) and Tb(III). The parameterizations have been carried out in a much more extensive manner, aimed at producing a new, more accurate model called Sparkle/AM1, mainly for the vast majority of all Eu(III), Gd(III), and Tb(III) complexes, which possess oxygen or nitrogen as coordinating atoms. All such complexes, which comprise 80% of all geometries present in the Cambridge Structural Database for each of the three ions, were classified into seven groups. These were regarded as a "basis" of chemical ambiance around a lanthanide, which could span the various types of ligand environments the lanthanide ion could be subjected to in any arbitrary complex where the lanthanide ion is coordinated to nitrogen or oxygen atoms. From these seven groups, 15 complexes were selected, which were defined as the parameterization set and then were used with a numerical multidimensional nonlinear optimization to find the best parameter set for reproducing chemical properties. The new parameterizations yielded an unsigned mean error for all interatomic distances between the Eu(III) ion and the ligand atoms of the first sphere of coordination (for the 96 complexes considered in the present paper) of 0.09 A, an improvement over the value of 0.28 A for the previous model and the value of 0.68 A for the first model (Chem. Phys. Lett. 1994, 227, 349). Similar accuracies have been achieved for Gd(III) (0.07 A, 70 complexes) and Tb(III) (0.07 A, 42 complexes). Qualitative improvements have been obtained as well; nitrates now coordinate correctly as bidentate ligands. The results, therefore, indicate that Eu(III), Gd(III), and Tb(III) Sparkle/AM1 calculations possess geometry prediction accuracies for lanthanide complexes with oxygen or nitrogen atoms in the coordination polyhedron that are competitive with present day ab initio/effective core potential calculations, while being hundreds of times faster.
Applications of hybrid QM/MM methods range from reactions in active sites of small enzymes to multiple sites in large bioenergetic complexes. By combining the widely used molecular dynamics and visualization programs NAMD/VMD with the quantum chemistry packages ORCA, and MOPAC, we provide an integrated, comprehensive, customizable, and easy-to-use suite. Through the QwikMD interface, setup, execution, visualization, and analysis are streamlined for all levels of expertise.
In this study, we present some modifications in the semiempirical quantum chemistry MOPAC2009 code that accelerate single-point energy calculations (1SCF) of medium-size (up to 2500 atoms) molecular systems using GPU coprocessors and multithreaded shared-memory CPUs. Our modifications consisted of using a combination of highly optimized linear algebra libraries for both CPU (LAPACK and BLAS from Intel MKL) and GPU (MAGMA and CUBLAS) to hasten time-consuming parts of MOPAC such as the pseudodiagonalization, full diagonalization, and density matrix assembling. We have shown that it is possible to obtain large speedups just by using CPU serial linear algebra libraries in the MOPAC code. As a special case, we show a speedup of up to 14 times for a methanol simulation box containing 2400 atoms and 4800 basis functions, with even greater gains in performance when using multithreaded CPUs (2.1 times in relation to the single-threaded CPU code using linear algebra libraries) and GPUs (3.8 times). This degree of acceleration opens new perspectives for modeling larger structures which appear in inorganic chemistry (such as zeolites and MOFs), biochemistry (such as polysaccharides, small proteins, and DNA fragments), and materials science (such as nanotubes and fullerenes). In addition, we believe that this parallel (GPU-GPU) MOPAC code will make it feasible to use semiempirical methods in lengthy molecular simulations using both hybrid QM/MM and QM/QM potentials.
Herein, molecular modeling techniques were used with the main goal to obtain candidates from a drug database as potential targets to be used against SARS-CoV-2. This novel coronavirus, responsible by the COVID-19 outbreak since the end of 2019, became a challenge since there is not vaccine for this disease. The first step in this investigation was to solvate the isolated S-protein in water for molecular dynamics (MD) simulation, being observed a transition from "up" to "down" conformation of receptor-binding domain (RBD) of the S-protein with angle of 54.3 and 43.0 degrees, respectively. The RBD region was more exposed to the solvent and to the possible drugs due to its enhanced surface area. From the equilibrated MD structure, virtual screening by docking calculations were performed using a library contained 9091 FDA approved drugs. Among them, 24 best-scored ligands (14 traditional herbal isolate and 10 approved drugs) with the binding energy below -8.1 kcal/mol were selected as potential candidates to inhibit the SARS-CoV-2 S-protein, preventing the human cell infection and their replication. For instance, the ivermectin drug (present in our list of promise candidates) was recently used successful to control viral replication in vitro. MD simulations were performed for the three best ligands@S-protein complexes and the binding energies were calculated using the MM/ PBSA approach. Overall, it is highlighted an important strategy, some key residues, and chemical groups which may be considered on clinical trials for COVID-19 outbreak.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
customersupport@researchsolutions.com
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
This site is protected by reCAPTCHA and the Google Privacy Policy and Terms of Service apply.
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.