Machine learning (ML) methods offer a promising route to the construction of universal molecular potentials with high accuracy and low computational cost. It is becoming evident that integrating physical principles into these models, or utilizing them in a Δ-ML scheme, significantly enhances their robustness and transferability. This paper introduces PM6-ML, a Δ-ML method that synergizes the semiempirical quantum-mechanical (SQM) method PM6 with a state-of-the-art ML potential applied as a universal correction. The method demonstrates superior performance over standalone SQM and ML approaches and covers a broader chemical space than its predecessors. It is scalable to systems with thousands of atoms, which makes it applicable to large biomolecular systems. Extensive benchmarking confirms PM6-ML's accuracy and robustness. Its practical application is facilitated by a direct interface to MOPAC. The code and parameters are available at https://github.com/Honza-R/mopac-ml.

• Publikační typ
• časopisecké články MeSH

Schistosomiasis, caused by a parasitic blood fluke of the genus Schistosoma, is a global health problem for which new chemotherapeutic options are needed. We explored the scaffold of gallinamide A, a natural peptidic metabolite of marine cyanobacteria that has previously been shown to inhibit cathepsin L-type proteases. We screened a library of 19 synthetic gallinamide A analogs and identified nanomolar inhibitors of the cathepsin B-type protease SmCB1, which is a drug target for the treatment of schistosomiasis mansoni. Against cultured S. mansoni schistosomula and adult worms, many of the gallinamides generated a range of deleterious phenotypic responses. Imaging with a fluorescent-activity-based probe derived from gallinamide A demonstrated that SmCB1 is the primary target for gallinamides in the parasite. Furthermore, we solved the high-resolution crystal structures of SmCB1 in complex with gallinamide A and its two analogs and describe the acrylamide covalent warhead and binding mode in the active site. Quantum chemical calculations evaluated the contribution of individual positions in the peptidomimetic scaffold to the inhibition of the target and demonstrated the importance of the P1' and P2 positions. Our study introduces gallinamides as a powerful chemotype that can be exploited for the development of novel antischistosomal chemotherapeutics.

• Klíčová slova
• Schistosoma mansoni, acrylamide inhibitor, cathepsin B, cysteine protease, drug target, parasite,
• MeSH
• kathepsin B * antagonisté a inhibitory   metabolismus   MeSH
• krystalografie rentgenová MeSH
• molekulární modely MeSH
• Schistosoma mansoni * enzymologie   účinky léků   MeSH
• schistosomicidy farmakologie   chemie   MeSH
• vazba proteinů MeSH
• zvířata MeSH
• Check Tag
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• kathepsin B * MeSH
• schistosomicidy MeSH

Accurate estimation of protein-ligand binding affinity is the cornerstone of computer-aided drug design. We present a universal physics-based scoring function, named SQM2.20, addressing key terms of binding free energy using semiempirical quantum-mechanical computational methods. SQM2.20 incorporates the latest methodological advances while remaining computationally efficient even for systems with thousands of atoms. To validate it rigorously, we have compiled and made available the PL-REX benchmark dataset consisting of high-resolution crystal structures and reliable experimental affinities for ten diverse protein targets. Comparative assessments demonstrate that SQM2.20 outperforms other scoring methods and reaches a level of accuracy similar to much more expensive DFT calculations. In the PL-REX dataset, it achieves excellent correlation with experimental data (average R2 = 0.69) and exhibits consistent performance across all targets. In contrast to DFT, SQM2.20 provides affinity predictions in minutes, making it suitable for practical applications in hit identification or lead optimization.

• MeSH
• ligandy MeSH
• proteiny * metabolismus   MeSH
• racionální návrh léčiv * MeSH
• termodynamika MeSH
• vazba proteinů MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• ligandy MeSH
• proteiny * MeSH

In this theoretical study, we set out to demonstrate the substitution effect of PEDOT analogues on planarity as an intrinsic indicator for electronic performance. We perform a quantum mechanical (DFT) study of PEDOT and analogous model systems and demonstrate the usefulness of the ωB97X-V functional to simulate chalcogen bonds and other noncovalent interactions. We confirm that the chalcogen bond stabilizes the planar conformation and further visualize its presence via the electrostatic potential surface. In comparison to the prevalent B3LYP, we gain 4-fold savings in computational time and simulate model systems of up to a dodecamer. Implications for design of conductive polymers can be drawn from the results, and an example for self-doped polymers is presented where modulation of the strength of the chalcogen bond plays a significant role.

• Publikační typ
• časopisecké články MeSH

Spectroscopic characteristics of Me3Si-H···Y complexes (Y = ICF3, BrCN, and HCN) containing a hydridic hydrogen were determined experimentally by low-temperature IR experiments based on the direct spectral measurement of supersonically expanded intermediates on a cold substrate or by the technique of argon-matrix isolation as well as computationally at harmonic and one-dimensional anharmonic levels. The computations were based on DFT-D, MP2, MP2-F12, and CCSD(T)-F12 levels using various extended AO basis sets. The formation of all complexes related to the redshift of the Si-H stretching frequency upon complex formation was accompanied by an increase in its intensity. Similar results were obtained for another 10 electron acceptors of different types, positive σ-, π-, and p-holes and cations. The formation of HBe-H···Y complexes, studied only computationally and again containing a hydridic hydrogen, was characterized by the blueshift of the Be-H stretching frequency upon complexation accompanied by an increase in its intensity. The spectral shifts and stabilization energies obtained for all presently studied hydridic H-bonded complexes were comparable to those in protonic H-bonded complexes, which has prompted us to propose a modification of the existing IUPAC definition of H-bonding that covers, besides the classical protonic form, the non-classical hydridic and dihydrogen forms.

• Publikační typ
• časopisecké články MeSH

Cathepsin K (CatK) is a target for the treatment of osteoporosis, arthritis, and bone metastasis. Peptidomimetics with a cyanohydrazide warhead represent a new class of highly potent CatK inhibitors; however, their binding mechanism is unknown. We investigated two model cyanohydrazide inhibitors with differently positioned warheads: an azadipeptide nitrile Gü1303 and a 3-cyano-3-aza-β-amino acid Gü2602. Crystal structures of their covalent complexes were determined with mature CatK as well as a zymogen-like activation intermediate of CatK. Binding mode analysis, together with quantum chemical calculations, revealed that the extraordinary picomolar potency of Gü2602 is entropically favoured by its conformational flexibility at the nonprimed-primed subsites boundary. Furthermore, we demonstrated by live cell imaging that cyanohydrazides effectively target mature CatK in osteosarcoma cells. Cyanohydrazides also suppressed the maturation of CatK by inhibiting the autoactivation of the CatK zymogen. Our results provide structural insights for the rational design of cyanohydrazide inhibitors of CatK as potential drugs.

• Klíčová slova
• Cathepsin K, azadipeptide nitrile, cyanohydrazide warhead, protease inhibitor, structure,
• MeSH
• hydraziny chemie   farmakologie   MeSH
• inhibitory proteas chemie   farmakologie   MeSH
• kathepsin K antagonisté a inhibitory   metabolismus   MeSH
• lidé MeSH
• molekulární modely MeSH
• molekulární struktura MeSH
• nádorové buňky kultivované MeSH
• nitrily chemie   farmakologie   MeSH
• vztah mezi dávkou a účinkem léčiva MeSH
• vztahy mezi strukturou a aktivitou MeSH
• Check Tag
• lidé MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• CTSK protein, human MeSH Prohlížeč
• hydraziny MeSH
• inhibitory proteas MeSH
• kathepsin K MeSH
• nitrily MeSH

Community efforts in the computational molecular sciences (CMS) are evolving toward modular, open, and interoperable interfaces that work with existing community codes to provide more functionality and composability than could be achieved with a single program. The Quantum Chemistry Common Driver and Databases (QCDB) project provides such capability through an application programming interface (API) that facilitates interoperability across multiple quantum chemistry software packages. In tandem with the Molecular Sciences Software Institute and their Quantum Chemistry Archive ecosystem, the unique functionalities of several CMS programs are integrated, including CFOUR, GAMESS, NWChem, OpenMM, Psi4, Qcore, TeraChem, and Turbomole, to provide common computational functions, i.e., energy, gradient, and Hessian computations as well as molecular properties such as atomic charges and vibrational frequency analysis. Both standard users and power users benefit from adopting these APIs as they lower the language barrier of input styles and enable a standard layout of variables and data. These designs allow end-to-end interoperable programming of complex computations and provide best practices options by default.

• Publikační typ
• časopisecké články MeSH

Schistosomiasis, a parasitic disease caused by blood flukes of the genus Schistosoma, is a global health problem with over 200 million people infected. Treatment relies on just one drug, and new chemotherapies are needed. Schistosoma mansoni cathepsin B1 (SmCB1) is a critical peptidase for the digestion of host blood proteins and a validated drug target. We screened a library of peptidomimetic vinyl sulfones against SmCB1 and identified the most potent SmCB1 inhibitors reported to date that are active in the subnanomolar range with second order rate constants (k2nd) of ∼2 × 105 M-1 s-1. High resolution crystal structures of the two best inhibitors in complex with SmCB1 were determined. Quantum chemical calculations of their respective binding modes identified critical hot spot interactions in the S1' and S2 subsites. The most potent inhibitor targets the S1' subsite with an N-hydroxysulfonic amide moiety and displays favorable functional properties, including bioactivity against the pathogen, selectivity for SmCB1 over human cathepsin B, and reasonable metabolic stability. Our results provide structural insights for the rational design of next-generation SmCB1 inhibitors as potential drugs to treat schistosomiasis.

• Klíčová slova
• Schistosoma mansoni, cathepsin B, cysteine peptidase, drug target, parasite, vinyl sulfone inhibitor,
• MeSH
• kathepsin B * MeSH
• lidé MeSH
• Schistosoma mansoni MeSH
• schistosomóza * farmakoterapie   MeSH
• sulfony farmakologie   MeSH
• zvířata MeSH
• Check Tag
• lidé MeSH
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Research Support, N.I.H., Extramural MeSH
• Názvy látek
• divinyl sulfone MeSH Prohlížeč
• kathepsin B * MeSH
• sulfony MeSH

We report the generation of gas-phase riboguanosine radicals that were tagged at ribose with a fixed-charge 6-(trimethylammonium)hexane-1-aminocarbonyl group. The radical generation relied on electron transfer from fluoranthene anion to noncovalent dibenzocrown-ether dication complexes which formed nucleoside cation radicals upon one-electron reduction and crown-ether ligand loss. The cation radicals were characterized by collision-induced dissociation (CID), photodissociation (UVPD), and UV-vis action spectroscopy. Identification of charge-tagged guanosine radicals was challenging because of spontaneous dissociations by loss of a hydrogen atom and guanine that occurred upon storing the ions in the ion trap without further excitation. The loss of H proceeded from an exchangeable position on N-7 in guanine that was established by deuterium labeling and was the lowest energy dissociation of the guanosine radicals according to transition-state energy calculations. Rate constant measurements revealed an inverse isotope effect on the loss of either hydrogen or deuterium with rate constants kH = 0.25-0.26 s-1 and kD = 0.39-0.54 s-1. We used time-dependent density functional theory calculations, including thermal vibronic effects, to predict the absorption spectra of several protomeric radical isomers. The calculated spectra of low-energy N-7-H guanine-radical tautomers closely matched the action spectra. Transition-state-theory calculations of the rate constants for the loss of H-7 and guanine agreed with the experimental rate constants for a narrow range of ion effective temperatures. Our calculations suggest that the observed inverse isotope effect does not arise from the isotope-dependent differences in the transition-state energies. Instead, it may be caused by the dynamics of post-transition-state complexes preceding the product separation.

• Publikační typ
• časopisecké články MeSH

Influenza A virus (IAV) encodes a polymerase composed of three subunits: PA, with endonuclease activity, PB1 with polymerase activity and PB2 with host RNA five-prime cap binding site. Their cooperation and stepwise activation include a process called cap-snatching, which is a crucial step in the IAV life cycle. Reproduction of IAV can be blocked by disrupting the interaction between the PB2 domain and the five-prime cap. An inhibitor of this interaction called pimodivir (VX-787) recently entered the third phase of clinical trial; however, several mutations in PB2 that cause resistance to pimodivir were observed. First major mutation, F404Y, causing resistance was identified during preclinical testing, next the mutation M431I was identified in patients during the second phase of clinical trials. The mutation H357N was identified during testing of IAV strains at Centers for Disease Control and Prevention. We set out to provide a structural and thermodynamic analysis of the interactions between cap-binding domain of PB2 wild-type and PB2 variants bearing these mutations and pimodivir. Here we present four crystal structures of PB2-WT, PB2-F404Y, PB2-M431I and PB2-H357N in complex with pimodivir. We have thermodynamically analysed all PB2 variants and proposed the effect of these mutations on thermodynamic parameters of these interactions and pimodivir resistance development. These data will contribute to understanding the effect of these missense mutations to the resistance development and help to design next generation inhibitors.

• Klíčová slova
• VX-787, antivirals, influenza A polymerase, pimodivir, resistance,
• MeSH
• krystalografie rentgenová MeSH
• kvantová teorie MeSH
• molekulární modely MeSH
• mutace genetika   MeSH
• mutantní proteiny metabolismus   MeSH
• podjednotky proteinů antagonisté a inhibitory   chemie   metabolismus   MeSH
• proteinové domény MeSH
• pyridiny chemie   farmakologie   MeSH
• pyrimidiny chemie   farmakologie   MeSH
• pyrroly chemie   farmakologie   MeSH
• RNA-dependentní RNA-polymerasa antagonisté a inhibitory   chemie   metabolismus   MeSH
• termodynamika MeSH
• virová léková rezistence účinky léků   MeSH
• virové proteiny antagonisté a inhibitory   chemie   metabolismus   MeSH
• virus chřipky A účinky léků   enzymologie   MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• mutantní proteiny MeSH
• PB2 protein, Influenzavirus A MeSH Prohlížeč
• pimodivir MeSH Prohlížeč
• podjednotky proteinů MeSH
• pyridiny MeSH
• pyrimidiny MeSH
• pyrroly MeSH
• RNA-dependentní RNA-polymerasa MeSH
• virové proteiny MeSH