Nature-Inspired Gallinamides Are Potent Antischistosomal Agents: Inhibition of the Cathepsin B1 Protease Target and Binding Mode Analysis

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.

SQM2.20: Semiempirical quantum-mechanical scoring function yields DFT-quality protein-ligand binding affinity predictions in minutes

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.

Chlorinated polyhedral selenaboranes revisited by joint experimental/computational efforts: the formation of closo-1-SeB9Cl9 and the crystal structure of closo-SeB11Cl11

The recent success in the formation of chlorinated telluraboranes and the reactivities of pnictogenaboranes prompted us to re-examine the vacuum co-pyrolysis of B2Cl4 with Se2Cl2 at various molar ratios and temperatures in order to search for the generation of other polyhedral selenaboranes than closo-SeB5Cl5 (1a) and closo-SeB11Cl11 (1b), the latter being observed earlier. Interestingly, a new compound with the elemental composition SeB9Cl9 (2) was detected, this time by high- and low-resolution mass spectrometry. Further characterization by 1- and 2-D 11B-NMR spectroscopy suggests that 2 should adopt a closed bicapped square-antiprismatic geometry with selenium at the apical position. Moreover, vacuum sublimation gave suitable crystals of 1b, which were subjected to single-crystal X-ray structure determination. Crystallographic data analysis confirmed that 1b, consistent with its 26 skeletal electron count, adopts a distorted icosahedral structure close to the symmetry of C5v. Computations at the DFT-D3 level have revealed that 33% of the total computed binding motifs in the grown 1b crystals are due to the very strong chalcogen bonding. Moreover, SAPT decomposition has shown that the bonding motifs in the crystals are stabilized mainly by dispersion and electrostatic terms. Homodecoupling and high resolution 11B NMR and 77Se NMR experiments have resolved both coupling constants 1J(11B11B) and 1J(77Se11B) as well as the 77Se chemical shift of 1a and 1b, which are in reasonable agreement with the corresponding computed values. The computed 11B chemical shifts of 2 were determined by the well-established DFT/GIAO/NMR structural tool based on its B3LYP/6-311+G** internal coordinates. They agree well with the experimental values and provide a good representation of the molecular structure of 2 in solution. The extraordinary downfield 11B NMR chemical shift of B(10) in 2 has been ascribed to the intensive paramagnetic contribution to the shielding tensor in this bicapped square-antiprismatic motif. Calculations of the synproportionation free energies of smaller (n - 1) closo-selenaboranes with larger-sized (n + 1) ones support the extraordinary stability of octahedral, bicapped square-antiprismatic and icosahedral closo motifs in the SeBnCln family (n = 4-12).

Boron-based octahedral dication experimentally detected: DFT surface confirms its availability

Borane and heteroborane clusters have been known as neutral or anionic species. In contrast to them, several ten-vertex monocationic nido and closo dicarbaborane-based systems have recently emerged from the reaction of the parent bicapped-square antiprismatic dicarbaboranes with N-heterocyclic carbenes followed by the protonization of the corresponding nido intermediates. The expansion of these efforts has afforded the very first closo-dicationic octahedral phosphahexaborane along with new closo-monocationic pnictogenahexaboranes of the same shapes. All are the products of the one-pot procedure that consists in the reaction of the same carbenes with the parent closo-1,2-Pn₂B₄Br₄ (Pn = As, P). Whereas in the case of phosphorus such a monocation appears to be a mixture of stable intermediates, and arsenahexaboranyl monocation has occurred as the final product, all of them without using any subsequent reaction. The well-established DFT/ZORA/NMR approach has unambiguously confirmed the existence of these species in solution, and computed electrostatic potentials have revealed the delocalization of the positive charge in these monocations and in the very first dication, namely within the octahedral shapes in both cases.

DFT Surface Infers Ten-Vertex Cationic Carboranes from the Corresponding Neutral closo Ten-Vertex Family: The Computed Background Confirming Their Experimental Availability

Modern computational protocols based on the density functional theory (DFT) infer that polyhedral closo ten-vertex carboranes are key starting stationary states in obtaining ten-vertex cationic carboranes. The rearrangement of the bicapped square polyhedra into decaborane-like shapes with open hexagons in boat conformations is caused by attacks of N-heterocyclic carbenes (NHCs) on the closo motifs. Single-point computations on the stationary points found during computational examinations of the reaction pathways have clearly shown that taking the "experimental" NHCs into account requires the use of dispersion correction. Further examination has revealed that for the purposes of the description of reaction pathways in their entirety, i.e., together with all transition states and intermediates, a simplified model of NHCs is sufficient. Many of such transition states resemble in their shapes those that dictate Z-rearrangement among various isomers of closo ten-vertex carboranes. Computational results are in very good agreement with the experimental findings obtained earlier.

Strategies for the Design of PEDOT Analogues Unraveled: the Use of Chalcogen Bonds and σ-Holes

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.

The Telluraboranes closo-TeB5Cl5 and closo-TeB11Cl11 with Exceptionally Long Body Diagonals: Synthetic and Bonding Motifs for Innovative Octahedral and Icosahedral Geometries

Six-vertex closo-TeB5Cl5 (1) and twelve-vertex closo-TeB11Cl11 (2) telluraboranes have been prepared via co-pyrolysis of B2Cl4 with TeCl4in vacuo at temperatures between 360ºC and 400°C. Both compounds are sublimable, off-white solids, and they have been characterized by one- and two-dimensional 11B-NMR and high-resolution mass spectroscopy. Both ab initio/GIAO/NMR and DFT/ZORA/NMR computations support octahedral and icosahedral geometries for 1 and 2, respectively, as expected due to their closo-electron counts. The octahedral structure of 1 has been confirmed by single-crystal X-ray diffraction on an incommensurately modulated crystal. The corresponding bonding properties have been analyzed in terms of the intrinsic bond orbital (IBO) approach. 1 is the first example of a polyhedral telluraborane with a cluster size smaller than 11 vertices.

B-H⋯π and C-H⋯π interactions in protein-ligand complexes: carbonic anhydrase II inhibition by carborane sulfonamides

Among non-covalent interactions, B-H⋯π and C-H⋯π hydrogen bonding is rather weak and less studied. Nevertheless, since both can affect the energetics of protein-ligand binding, their understanding is an important prerequisite for reliable predictions of affinities. Through a combination of high-resolution X-ray crystallography and quantum-chemical calculations on carbonic anhydrase II/carborane-based inhibitor systems, this paper provides the first example of B-H⋯π hydrogen bonding in a protein-ligand complex. It shows that the B-H⋯π interaction is stabilized by dispersion, followed by electrostatics. Furthermore, it demonstrates that the similar C-H⋯π interaction is twice as strong, with a slightly smaller contribution of dispersion and a slightly higher contribution of electrostatics. Such a detailed insight will facilitate the rational design of future protein ligands, controlling these types of non-covalent interactions.

Reactivity of Perhalogenated Octahedral Phospha- and Arsaboranes toward THF: A Joint Experimental/Computational Study

Reactions of the perhalogenated polyhedral pnictogenaboranes closo-1,2-Pn₂B₄Hal₄ (Pn = P, As; Hal = Cl, Br) with Lewis bases are presently being studied with a focus on rationalizing the sites of nucleophilic attacks on clusters bearing σ-holes. These σ-holes are localized both on pnictogens and, for Hal = Br, on bromine atoms, as revealed by electrostatic potential (ESP) and intrinsic bond orbital (IBO) analyses. Surprisingly, the attack of the cyclic ether THF on closo-1,2-Pn₂B₄Br₄ does not occur on the site with the largest positive partial charge, centered in the middle of the pnictogen-pnictogen vector. Instead, presumably promoted by the positivated bromine substituents, THF inserts into the boron-bromine bonds of the negatively charged boron atoms opposite to the pnictogen atoms to form 4-(4-bromobut-1-oxy)-closo-1,2-Pn₂B₄Br₃ (1-PB and 1-AsB) and 4,6-(4-bromobut-1-oxy)₂-closo-1,2-Pn₂B₄Br₂ (2-PB and 2-AsB). ¹¹B and ³¹P chemical shift computations at various levels support the assignments of the signals, which reflect the correctness of the molecular geometries in solutions. The Lewis-acidic perchlorinated analogues closo-1,2-P₂B₄Cl₄, closo-1,2-As₂B₄Cl₄, and the mixed closo-1,2-AsPB₄Cl₄ bear negative charges. These negative charges are revealed by the V_s,max values when computing the electrostatic potentials both on the boron and the chlorine atoms. Due to this negative charge, the analogues do not react with THF unless they are heated above 66 °C, where they slowly decompose to borate esters B(OR)₃ without the formation of concrete intermediates. The evaluation of ³¹P NMR data of 1-PB has allowed the experimental determination of the coupling constant ¹J(³¹P(1), ³¹P(2)) = |143| Hz in a closo-diphosphaborane for the first time, which agrees well with the computed value of -178 Hz. The pioneering joint experimental vs computational interpretation of ³¹P NMR spectra in the area of boron cluster chemistry was decisive for the structural characterization of 1-PB and 2-PB.

Highly potent inhibitors of cathepsin K with a differently positioned cyanohydrazide warhead: structural analysis of binding mode to mature and zymogen-like enzymes

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.

Transformation of various multicenter bondings within bicapped-square antiprismatic motifs: Z-rearrangement

Reported herein are mutual rearrangements in the whole series of seven bicapped-square antiprismatic closo-C₂B₈H₁₀ by means of high-quality computations that disprove the earlier postulated dsd (diamond-square-diamond) scheme for these isomerizations. The experimentally existing closo-1,2-C₂B₈H₁₀ was able to be converted to 1,6-, and 1,10-isomers by pyrolysis, and the dsd (diamond-square-diamond) mechanism was offered as an explanation of these processes. However, these computations disprove the postulated dsd scheme for these isomerizations that take place in the ten-vertex closo series. Experimentally observed thermal rearrangements, both in the parent and substituted closo-1,2-C₂B₈H₁₀, closo-1-CB₉H₁₀^-, and closo-B₁₀H₁₀^2-, indirectly support these refined computations. All these processes are based on the new concept of the so-called Z-mechanism, being consistent with a transition state of a boat shape with an open hexagonal belt that results from the initial breakage of three bonds. Such bond breakings and the consequent bond formations bring to mind the shape of the letter Z. In effect, the pattern of multicenter bonding shifts from reactant through a transition state to product. The molecular rearrangements that are available experimentally favour either the axial or equatorial isomers, and this ratio depends on temperature and the type of cluster and its substitution.

Access to cationic polyhedral carboranes via dynamic cage surgery with N-heterocyclic carbenes

Polyhedral boranes and heteroboranes appear almost exclusively as neutral or anionic species, while the cationic ones are protonated at exoskeletal heteroatoms or they are instable. Here we report the reactivity of 10-vertex closo-dicarbadecaboranes with one or two equivalents of N-heterocyclic carbene to 10-vertex nido mono- and/or bis-carbene adducts, respectively. These complexes easily undergo a reaction with HCl to give cages of stable and water soluble 10-vertex nido-type cations with protonation in the form of a BHB bridge or 10-vertex closo-type cations containing one carbene ligand when originating from closo-1,10-dicarbadecaborane. The reaction of a 10-vertex nido mono-carbene adduct with phosphorus trichloride gives nido-11-vertex 2-phospha-7,8-dicarbaundecaborane, which undergoes an oxidation of the phosphorus atom to P = O, while the product of a bis-carbene adduct reaction is best described as a distorted C₂B₆H₈ fragment bridged by the (BH)₂PCl₂⁺ moiety.

In comparison to their neutral or anionic counterparts, examples of cationic boron clusters remain scarce. Here, the authors prepare a variety of cationic polyhedral boranes by reacting closo-10-vertex carboranes with N-heterocyclic carbenes; the resulting open-cage cationic nido- arachno- or closo- derivatives are water soluble, which may enable unprecedented applications.

The Structure-Based Design of SARS-CoV-2 nsp14 Methyltransferase Ligands Yields Nanomolar Inhibitors

In this study, we have focused on the structure-based design of the inhibitors of one of the two SARS-CoV-2 methyltransferases (MTases), nsp14. This MTase catalyzes the transfer of the methyl group from S-adenosyl-l-methionine (SAM) to cap the guanosine triphosphate moiety of the newly synthesized viral RNA, yielding the methylated capped RNA and S-adenosyl-l-homocysteine (SAH). As the crystal structure of SARS-CoV-2 nsp14 is unknown, we have taken advantage of its high homology to SARS-CoV nsp14 and prepared its homology model, which has allowed us to identify novel SAH derivatives modified at the adenine nucleobase as inhibitors of this important viral target. We have synthesized and tested the designed compounds in vitro and shown that these derivatives exert unprecedented inhibitory activity against this crucial enzyme. The docking studies nicely explain the contribution of an aromatic part attached by a linker to the position 7 of the 7-deaza analogues of SAH.

Thiaborane Icosahedral Barrier Increased by the Functionalization of all Terminal Hydrogens in closo-1-SB11H11

The electrophilic substitution of icosahedral closo-1-SB₁₁H₁₁ with methyl iodide has resulted in two B-functionalized thiaboranes, 7,12-I₂-2,3,4,5,6,8,9,10,11-(CH₃)₉-1-closo-SB₁₁ and 7,8,12-I₃-2,3,4,5,6,9,10,11-(CH₃)₈-closo-1-SB₁₁, with the former being significantly predominant. These two icosahedral thiaboranes are the first cases of polysubstituted polyhedral boron clusters with another vertex that differs from B and C. Such polyfunctionalizations have increased the earlier observed thiaborane icosahedral barrier, not exhibiting any reactivity toward bases, unlike the parent thiaborane. The search for methylation pathways has revealed that the complete B₁₁-methylation is impossible, like in the case of decaborane(14), where this seems to be a result of the positively charged upper parts of these two molecules.

Druggable Hot Spots in the Schistosomiasis Cathepsin B1 Target Identified by Functional and Binding Mode Analysis of Potent Vinyl Sulfone Inhibitors

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 (k_2nd) of ∼2 × 10⁵ 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.

Structural and Thermodynamic Analysis of the Resistance Development to Pimodivir (VX-787), the Clinical Inhibitor of Cap Binding to PB2 Subunit of Influenza A Polymerase

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.

Azanitrile Inhibitors of the SmCB1 Protease Target Are Lethal to Schistosoma mansoni: Structural and Mechanistic Insights into Chemotype Reactivity

Azapeptide nitriles are postulated to reversibly covalently react with the active-site cysteine residue of cysteine proteases and form isothiosemicarbazide adducts. We investigated the interaction of azadipeptide nitriles with the cathepsin B1 drug target (SmCB1) from Schistosoma mansoni, a pathogen that causes the global neglected disease schistosomiasis. Azadipeptide nitriles were superior inhibitors of SmCB1 over their parent carba analogs. We determined the crystal structure of SmCB1 in complex with an azadipeptide nitrile and analyzed the reaction mechanism using quantum chemical calculations. The data demonstrate that azadipeptide nitriles, in contrast to their carba counterparts, undergo a change from E- to Z-configuration upon binding, which gives rise to a highly favorable energy profile of noncovalent and covalent complex formation. Finally, azadipeptide nitriles were considerably more lethal than their carba analogs against the schistosome pathogen in culture, supporting the further development of this chemotype as a treatment for schistosomiasis.

Benchmark Data Sets of Boron Cluster Dihydrogen Bonding for the Validation of Approximate Computational Methods

The success of approximate computational methods, such as molecular mechanics, or dispersion-corrected density functional theory, in the description of non-covalent interactions relies on accurate parameterizations. Benchmark data sets are thus required. This area is well developed for organic molecules and biomolecules but practically non-existent for boron clusters, which have been gaining in importance in modern drug as well as material design. To fill this gap, we have introduced two data sets featuring the most common non-covalent interaction of boron clusters, the dihydrogen bond, and calculated reference interaction energies at the "golden standard" CCSD(T)/CBS level. The boron clusters studied interact with formamide, methanol, water and methane at various distances and in two geometrical arrangements. The performance of the tested approximate methods is variable and recommendations for further use are given.

Bromination Mechanism of closo-1,2-C2 B10 H12 and the Structure of the Resulting 9-Br-closo-1,2-C2 B10 H11 Determined by Gas Electron Diffraction

9-Br-closo-1,2-C₂ B₁₀ H₁₁ has been prepared and its gas-phase structure has been examined by means of gas electron diffraction. The structure of the carbaborane core is similar to the structure of the parent compound, which is of C_2v symmetry. A DFT-based search for the corresponding reaction pathway of the bromination of closo-1,2-C₂ B₁₀ H₁₂ revealed that the catalytic amount of aluminum reduces the barrier of the initial attack of the bromination agent toward the negatively charged part of the icosahedral carbaborane, i. e., the first transition state, from about 40 to about 27 kcalmol^-1 . The Br-Br bond is weakened by an intermediate binding to the large π-hole on the aluminum atom of AlBr₃ , which is the driving force for the AlBr₃ -catalyzed bromination.

SQM/COSMO Scoring Function: Reliable Quantum-Mechanical Tool for Sampling and Ranking in Structure-Based Drug Design

Invited for this month's cover is the group of Prof. Pavel Hobza, Czech Academy of Sciences, Prague. The cover picture shows a powerful automated quantum mechanics based SQM/COSMO approach to protein-ligand scoring. It comprises thorough preparation of ligand structures, extensive generation of binding complexes, fast geometry relaxation and reliable affinity prediction. Read the full text of the Minireview at 10.1002/cplu.202000120.

Unraveling the anti-influenza effect of flavonoids: Experimental validation of luteolin and its congeners as potent influenza endonuclease inhibitors

The biological effects of flavonoids on mammal cells are diverse, ranging from scavenging free radicals and anti-cancer activity to anti-influenza activity. Despite appreciable effort to understand the anti-influenza activity of flavonoids, there is no clear consensus about their precise mode-of-action at a cellular level. Here, we report the development and validation of a screening assay based on AlphaScreen technology and illustrate its application for determination of the inhibitory potency of a large set of polyols against PA N-terminal domain (PA-Nter) of influenza RNA-dependent RNA polymerase featuring endonuclease activity. The most potent inhibitors we identified were luteolin with an IC₅₀ of 72 ± 2 nM and its 8-C-glucoside orientin with an IC₅₀ of 43 ± 2 nM. Submicromolar inhibitors were also evaluated by an in vitro endonuclease activity assay using single-stranded DNA, and the results were in full agreement with data from the competitive AlphaScreen assay. Using X-ray crystallography, we analyzed structures of the PA-Nter in complex with luteolin at 2.0 Å resolution and quambalarine B at 2.5 Å resolution, which clearly revealed the binding pose of these polyols coordinated to two manganese ions in the endonuclease active site. Using two distinct assays along with the structural work, we have presumably identified and characterized the molecular mode-of-action of flavonoids in influenza-infected cells.

SQM/COSMO Scoring Function: Reliable Quantum-Mechanical Tool for Sampling and Ranking in Structure-Based Drug Design

Quantum mechanical (QM) methods have been gaining importance in structure-based drug design where a reliable description of protein-ligand interactions is of utmost significance. However, strategies i. e. QM/MM, fragmentation or semiempirical (SQM) methods had to be pursued to overcome the unfavorable scaling of QM methods. Various SQM-based approaches have significantly contributed to the accuracy of docking and improvement of lead compounds. Parametrizations of SQM and implicit solvent methods in our laboratory have been instrumental to obtain a reliable SQM-based scoring function. The experience gained in its application for activity ranking of ligands binding to tens of protein targets resulted in setting up a faster SQM/COSMO scoring approach, which outperforms standard scoring methods in native pose identification for two dozen protein targets with ten thousand poses. Recently, SQM/COSMO was effectively applied in a proof-of-concept study of enrichment in virtual screening. Due to its superior performance, feasibility and chemical generality, we propose the SQM/COSMO approach as an efficient tool in structure-based drug design.

Face-Fusion of Icosahedral Boron Hydride Increases Affinity to γ-Cyclodextrin: closo,closo-[B21 H18 ]- as an Anion with Very Low Free Energy of Dehydration

The supramolecular recognition of closo,closo-[B21 H18 ]- by cyclodextrins (CDs) has been studied in aqueous solution by isothermal titration calorimetry and nuclear magnetic resonance spectroscopy. These solution studies follow up on previous mass-spectrometric measurements and computations, which indicated the formation and stability of CD ⋅ B21 H18- complexes in the gas phase. The thermodynamic signature of solution-phase binding is exceptional, the association constant for the γ-CD complex with B21 H18- reaches 1.8×106  M-1 , which is on the same order of magnitude as the so far highest observed value for the complex between γ-CD and a metallacarborane. The nature of the intermolecular interaction is also examined by quantum-mechanical computational protocols. These suggest that the desolvation penalty, which is particularly low for the B21 H18- anion, is the decisive factor for its high binding strength. The results further suggest that the elliptical macropolyhedral boron hydride is another example of a CD binder, whose extraordinary binding affinity is driven by the chaotropic effect, which describes the intrinsic affinity of large polarizable and weakly solvated chaotropic anions to hydrophobic cavities and surfaces in aqueous solution.

Complexation and stability of the fungicide penconazole in the presence of zinc and copper ions

RATIONALE

For the risk assessment of penconazole in the environment and the evaluation of the possible consequences of its use, it is important to determine how its reactivity and degradation are influenced by metals commonly found in nature, such as copper and zinc.

METHODS

Changes in the reactivity of penconazole in the presence of zinc/copper ions were studied using electrospray ionisation mass spectrometry and density functional theory calculations.

RESULTS

Many penconazole complexes with copper and zinc ions were created; a comparison of the elements showed that a few complexes were formed analogously (doubly charged complexes with four penconazole molecules, singly charged complexes with chlorine as a counterion and singly charged complexes with deprotonated penconazole as a counterion). The metal complexes with different structures indicated different reactivity of penconazole with copper and zinc.

CONCLUSIONS

The experimental and computational approaches have revealed different changes in the structure of penconazole. In the Zn(II) complex, penconazole deprotonated to stabilise the bond to Zn(II). In the Cu(II) complex, it loses one chlorine atom, creates an additional ring between the triazole ring and the phenyl ring, and/or creates a double bond in the short aliphatic chain.

Biomimetic Macrocyclic Inhibitors of Human Cathepsin D: Structure-Activity Relationship and Binding Mode Analysis

Human cathepsin D (CatD), a pepsin-family aspartic protease, plays an important role in tumor progression and metastasis. Here, we report the development of biomimetic inhibitors of CatD as novel tools for regulation of this therapeutic target. We designed a macrocyclic scaffold to mimic the spatial conformation of the minimal pseudo-dipeptide binding motif of pepstatin A, a microbial oligopeptide inhibitor, in the CatD active site. A library of more than 30 macrocyclic peptidomimetic inhibitors was employed for scaffold optimization, mapping of subsite interactions, and profiling of inhibitor selectivity. Furthermore, we solved high-resolution crystal structures of three macrocyclic inhibitors with low nanomolar or subnanomolar potency in complex with CatD and determined their binding mode using quantum chemical calculations. The study provides a new structural template and functional profile that can be exploited for design of potential chemotherapeutics that specifically inhibit CatD and related aspartic proteases.

Thiaboranes on Both Sides of the Icosahedral Barrier: Retaining and Breaking the Barrier with Carbon Functionalities

The concept of icosahedral barrier has been expanded from the chemistry of carbaboranes to the area of thiaboranes. Both representatives of this barrier, i. e., closo-1,2-C₂ B₁₀ H₁₂ and closo-1-SB₁₁ H₁₁ , are similar in their electron distribution, which is dominated by positive charge in the midpoint of the C-C vector and on the sulfur atom with experimentally determined dipole moments of 4.50 D and 3.64 D, respectively. This is a driving force for their reactivity as exemplified by their reactions with different carbon functionalities. Icosahedral closo-1-SB₁₁ H₁₁ reacts both with an electron sextet containing carbon (in the form of N-heterocyclic carbenes), reported earlier, and with methyl iodide with an electron octet on the carbon. The latter reaction provides hexamethylated thiaborane on the basis of methylation so far unknown in this area of heteroborane chemistry. The computations of the heat of formation (ΔH_f ²⁹⁸ ) make it possible to estimate the height of the barrier as well as to propose closo-thiaboranes beyond the barrier. Eleven and twelve vertex thiaboranes with nido electron count are known experimentally for breaking the barrier. These computations also suggest that the larger nido-thiaboranes are promising candidates for the corresponding experimental availability, i. e., the ΔH_f ²⁹⁸ of a 13-vertex nido-thiaborane cluster has been computed to be more negative than that of the well-known nido-SB₁₀ H₁₁ ^- cluster (-6.7 and -5.6 kcal mol^-1 per vertex, respectively).

Surface termination of MgB2 unveiled by a combination of adsorption experiments and theoretical calculations

Superconductivity in polycrystalline and thin-film MgB2 is strongly affected by the termination of its surface, but a reliable determination of the surface termination is still a challenging task of surface chemistry. Here, the surface properties of superconducting MgB2 were investigated using a combination of inverse gas chromatography and van der Waals corrected density functional theory calculations. The dispersive surface energy was measured as a function of the surface coverage and its value (58 mJ m-2 to 48 mJ m-2) was verified by high-level non-local EXX + RPA calculations, which predicted that the dispersive contribution to the cleavage energy was 56 mJ m-2. The isosteric adsorption enthalpies of cyclohexane, dioxane, acetone and acetonitrile molecules were measured on an MgB2 sample and compared to the DFT calculated enthalpies for the Mg-terminated MgB2, B-terminated MgB2 and MgO(001) surfaces. The close agreement between theory and experiment for the Mg-terminated surface suggested that the magnesium termination is the dominant surface phase of MgB2. Thus, combining inverse gas chromatography experiments with theoretical calculations may provide information about the surface termination.

Thiaborane clusters with an exoskeletal B-H group

The thiaboranes closo-1-SB11H11 (1a) and 12-I-closo-1-SB11H10 (1b) react with 4-(dimethylamino)pyridine under inert conditions upon the formation of the nido-type thiaboranes 9-B{(4-Me2N)C5NH4}2(H)-7-SB10H11 and 9-B{(4-Me2N)C5NH4}2(H)-5-I-7-SB10H10 containing an exoskeletal B-H group. The same type of B-H moiety is also stabilised by one bipyridine molecule in a chelating fashion. These complexes are unstable in solution, and in air and hydrolyse to monodeboronated ionic compounds having [nido-7-SB10H11]- or [5-I-nido-7-SB10H10]- anions which are also products of the reactions of 1a and 1b with other N-bases such as pyridine, ammonia and DABCO. The extrusion of one boron and one sulphur atom takes place when 1a reacts with 2,6-di-tert-butylpyridine to yield decaborane.

Investigation of Thiaborane closo- nido Conversion Pathways Promoted by N-Heterocyclic Carbenes

The 12-X- closo-SB₁₁H₁₀ (X = H or I) thiaboranes react with one or two molar equivalents of various N-heterocyclic carbenes (NHCs) to give the deprotonated 12-vertex species of [12-X-SB₁₁H₉·NHC]^-[NHC-H]⁺composition as kinetic products. The use of one molar equivalent of a sterically more hindered NHC reactant leads to the formation of 12-X-SB₁₁H₁₀·NHC adducts with a heavily distorted cage and the nido electron count. Further reaction of 12-I-SB₁₁H₁₀·NHC to deboronated 12-X-SB₁₀H₉·NHC proceeds in acetone to complete the closo- nido reaction pathway under the thermodynamic control. The structures of all compounds have been investigated by NMR spectroscopy and diffraction techniques. The results are supported by theoretical methods.

Synthesis of closo-1,2-H2C2B8Me8 and 1,2-H2C2B8Me7X (X = I and OTf) Dicarbaboranes and Their Rearrangement Reactions

Methyl-camouflaged dicarbaboranes closo-1,2- and 1,10-H₂C₂B₈Me₈ have been prepared in high yields either from nido-5,6-H₂C₂B₈H₁₀ or closo-1,2-H₂C₂B₈H₈ via electrophilic methylation reactions and cluster-rearrangement methods. Prepared were also monosubstituted derivatives of general formulation closo-H₂C₂B₈Me₇-X (X = I or OTf). The permethylated compounds exhibit extreme air stability in comparison to unprotected counterparts as a consequence of rigid, egg-shaped hydrocarbon structures incorporating inner C₂B₈ carborane scaffolding. The structures of all compounds isolated were confirmed unambiguously by multinuclear (¹¹B, ¹H, ¹³C, and ¹⁹F) NMR measurements, supported by X-ray diffraction analyses and geometry optimization methods on several compounds.

A theoretical analysis of the structure and properties of B26H30 isomers. Consequences to the laser and semiconductor doping capabilities of large borane clusters

The most stable isomer B₂₆H₃₀ with its predicted applications.

Quantitative syntheses of permethylated closo-1,10-R2C2B8Me8 (R = H, Me) carboranes. Egg-shaped hydrocarbons on the Frontier between inorganic and organic chemistry

Electrophilic methylation of the closo-1,10-R2C2B8H8 (1) (R = H or Me) dicarbaboranes at higher temperatures or thermal rearrangement of the 1,6-R2C2B8Me8 (3) compounds at 400-500 °C generated the B-permethylated derivatives closo-1,10-R2C2B8Me8 (2) in quantitative (>95%) yields. The compounds exhibit extreme air stability as a consequence of a rigid, egg shaped hydrocarbon structures incorporating inner 1,10-C2B8 carborane core.

Outerly functionalized and non-functionalized boron clusters intercalated into layered hydroxides with different modes of binding: materials for superacid storage

Two binary boron hydrides (NH4)2B10H10 and Na2B12H12 and mono- and dicarboxy p- and m-carboranes (namely, 1-(COOH)-closo-1,7-C2B10H11, 1,12-(COOH)2-closo-1,12-C2B10H10 and 1,7-(COOH)2-closo-1,7-C2B10H10) were intercalated into ZnAl-layered double hydroxides (ZnAl-LDH) and into Zn5(OH)8(NO3)2·2H2O. The formed compounds were characterized using elemental analysis, thermogravimetry analysis, X-ray powder diffraction, infrared spectroscopy and solid state NMR. All the intercalated boron compounds are present in the interlayer space of the layered hosts as anions. It is presumed that in the case of B10H102-, B12H122- and 1,12-(COO)2-closo-1,12-C2B10H102-, the guest molecules form a monolayer, whereas in the case of 1-(COO)-closo-1,7-C2B10H111- and 1,7-(COO)2-closo-1,7-C2B10H102- a bilayer arrangement is more probable. In the case of 1,7-(COO)2-closo-1,7-C2B10H102-, the guest molecules are strongly interdigitated resulting in lowering of the interlayer distance. Two different modes of binding were found. Whereas the carboxylate derivatives of p- and m-carboranes are bonded through classical hydrogen bonds, the corresponding parent borane anions interact with the host structures by mainly dihydrogen bonding. In effect, both kinds of hydrogen bonding are mainly of an electrostatic nature. The dihydrogen bond is detected, e.g. in crystal engineering, and represents a driving force for interactions of boranes with biomolecules. Since the latter dicarboxylic acids were found to be superacids, their interactions with the host structures should be stronger than in the case of the benzoic and terephthalic acid intercalates.

Methyl camouflage in the ten-vertex closo-dicarbaborane(10) series. Isolation of closo-1,6-R2C2B8Me8 (R = H and Me) and their monosubstituted analogues

Reported are procedures leading to the first types of methyl camouflaged dicarbadecaboranes with fewer than eleven vertices. The compounds contain the closo-1,6-C2B8 scaffolding inside the egg-shaped hepta - decamethyl sheath, which imparts unusually high air and solvolytic stability to all of these compounds.

Chalcogen Bonding in Protein-Ligand Complexes: PDB Survey and Quantum Mechanical Calculations

A chalcogen bond is a nonclassical noncovalent interaction which can stabilise small-molecule crystals as well as protein structures. Here, we systematically explore the stabilising potential of chalcogen bonding in protein-ligand complexes in the Protein Data Bank (PDB). We have found that a large fraction (23 %) of complexes with a S/Se-containing ligand feature close S/Se⋅⋅⋅O/N/S contacts. Eleven non-redundant representative potential S/Se⋅⋅⋅O chalcogen-bond motifs were selected and truncated to model systems and seven more model systems were prepared by S-to-Se substitution. These systems were then subjected to analysis by quantum chemical (QM) methods-electrostatic potential, geometry optimisation or interaction energy calculations, including solvent effects. The QM calculations indicate that chalcogen bonding does indeed play a dominant role in stabilising some of the interaction motifs studied. We thus advocate further exploration of chalcogen bonding with the aim of potential future use in structure-based drug design.

Ranking Power of the SQM/COSMO Scoring Function on Carbonic Anhydrase II-Inhibitor Complexes

Accurate prediction of protein-ligand binding affinities is essential for hit-to-lead optimization and virtual screening. The reliability of scoring functions can be improved by including quantum effects. Here, we demonstrate the ranking power of the semiempirical quantum mechanics (SQM)/implicit solvent (COSMO) scoring function by using a challenging set of 10 inhibitors binding to carbonic anhydrase II through Zn²⁺ in the active site. This new dataset consists of the high-resolution (1.1-1.4 Å) crystal structures and experimentally determined inhibitory constant (K_i ) values. It allows for evaluation of the common approximations, such as representing the solvent implicitly or by using a single target conformation combined with a set of ligand docking poses. SQM/COSMO attained a good correlation of R² of 0.56-0.77 with the experimental inhibitory activities, benefiting from careful handling of both noncovalent interactions (e.g. charge transfer) and solvation. This proof-of-concept study of SQM/COSMO ranking for metalloprotein-ligand systems demonstrates its potential for hit-to-lead applications.

A systematic examination of classical and multi-center bonding in heteroborane clusters

A systematic revision of bonding types on a broad series of heteroboranes covering closo, nido, arachno and hypho architectures with incorporated tetrel, pnictogen or chalcogen heterovertices.

Various types of non-covalent interactions contributing towards crystal packing of halogenated diphospha-dicarbaborane with an open pentagonal belt

We have synthesized and crystalized 3-Cl-10-I-nido-7,8,9,11-P₂C₂B₇H₇. Quantum chemical calculations have demonstrated that the obtained crystal structure is stabilized by hydrogen, dihydrogen and pnictogen bonds.

Mimicking of cyproconazole behavior in the presence of Cu and Zn

RATIONALE

The frequently used pesticide cyproconazole (Cyp) interacts with the essential elements commonly present in soil which play important roles in various enzymatic processes. These interactions predetermine the Cyp degradation pathways. We suggest a simple experimental and theoretical approach for the prediction of pesticide behavior.

METHODS

Cu/Cyp complexes are explored because of the typical Cu(II) reduction in complexes. Its level and the stability of the Cu-ligand bond depend on the type and the number of the surrounding ligands. Zn/Cyp complexes were compared as it is not expected that Zn(II) will reduce. The complexations were studied by means of electrospray ionization ion trap mass spectrometry and MS/MS collision-induced dissociations with comparative and explicative density functional theory calculations.

RESULTS

The Cyp ligand allows both Cu(II) reduction as well as, in specific cases, it protects the higher Cu oxidation state. The reduction is observed in the complexes with solely neutral Cyp where the number of ligands is below 3; a higher number protects the Cu(II) state. The metal atom binds to Cyp via N2 of the triazole ring as well as via π-electrons of the benzene ring; additional stabilization brings an interaction with the deprotonated OH group.

CONCLUSIONS

The character of Cyp interactions with doubly charged metals (Cu(II), Zn(II)) clarified the creation of Cyp metabolites. The phenyl and triazole rings are bound to the metal cation and enable access for the isopropyl ring to be cleaved leaving the common metabolite (CAS Number: 58905-19-4).

A novel stibacarbaborane cluster with adjacent antimony atoms exhibiting unique pnictogen bond formation that dominates its crystal packing

We have prepared nido-7,8,9,11-Sb₂C₂B₇H₉, the first cluster with simultaneous Sb-B, Sb-C and Sb-Sb atom pairs with interatomic separations with magnitudes that approach the respective sums of covalent radii. However, the length of the Sb-Sb separation in this cluster is slightly less than the sum of the covalent radii. Quantum chemical analysis has revealed that the crystal packing of nido-7,8,9,11-Sb₂C₂B₇H₉ is predominantly dictated by pnictogen (Pn) bonding, an unconventional σ-hole interaction. Indeed, the interaction energy of a very strong Sb₂H-B Pn-bond in the nido-7,8,9,11-Sb₂C₂B₇H₉ dimer exceeds -6.0 kcal mol^-1. This is a very large value and is comparable to the strengths of known Pn-bonds in Cl₃Pnπ complexes (Pn = As, Sb).

Superior Performance of the SQM/COSMO Scoring Functions in Native Pose Recognition of Diverse Protein-Ligand Complexes in Cognate Docking

General and reliable description of structures and energetics in protein-ligand (PL) binding using the docking/scoring methodology has until now been elusive. We address this urgent deficiency of scoring functions (SFs) by the systematic development of corrected semiempirical quantum mechanical (SQM) methods, which correctly describe all types of noncovalent interactions and are fast enough to treat systems of thousands of atoms. Two most accurate SQM methods, PM6-D3H4X and SCC-DFTB3-D3H4X, are coupled with the conductor-like screening model (COSMO) implicit solvation model in so-called "SQM/COSMO" SFs and have shown unique recognition of native ligand poses in cognate docking in four challenging PL systems, including metalloprotein. Here, we apply the two SQM/COSMO SFs to 17 diverse PL complexes and compare their performance with four widely used classical SFs (Glide XP, AutoDock4, AutoDock Vina, and UCSF Dock). We observe superior performance of the SQM/COSMO SFs and identify challenging systems. This method, due to its generality, comparability across the chemical space, and lack of need for any system-specific parameters, gives promise of becoming, after comprehensive large-scale testing in the near future, a useful computational tool in structure-based drug design and serving as a reference method for the development of other SFs.

B–H⋯π: a nonclassical hydrogen bond or dispersion contact?

Quantum mechanical calculations disprove the attractive electrostatic nature of B–H⋯π motif and define it as dispersion-driven contact.

Binary twinned-icosahedral [B21H18]− interacts with cyclodextrins as a precedent for its complexation with other organic motifs

A macropolyhedral boron hydride anion with two counterions can form stable complexes with β- and γ-cyclodextrin in the gas phase.