Extensive targeting of chemical space at the prime side of ketoamide inhibitors of rhomboid proteases by branched substituents empowers their selectivity and potency.

Rhomboid intramembrane serine proteases have been implicated in several pathologies, and emerge as attractive pharmacological target candidates. The most potent and selective rhomboid inhibitors available to date are peptidyl α-ketoamides, but their selectivity for diverse rhomboid proteases and strategies to modulate it in relevant contexts are poorly understood. This gap, together with the lack of suitable in vitro models, hinders ketoamide development for relevant eukaryotic rhomboid enzymes. Here we explore the structure-activity relationship principles of rhomboid inhibiting ketoamides by medicinal chemistry and enzymatic in vitro and in-cell assays with recombinant rhomboid proteases GlpG, human mitochondrial rhomboid PARL and human RHBDL2. We use X-ray crystallography in lipidic cubic phase to understand the binding mode of one of the best ketoamide inhibitors synthesized here containing a branched terminal substituent bound to GlpG. In addition, to extend the interpretation of the co-crystal structure, we use quantum mechanical calculations and quantify the relative importance of interactions along the inhibitor molecule. These combined experimental analyses implicates that more extensive exploration of chemical space at the prime side is unexpectedly powerful for the selectivity of rhomboid inhibiting ketoamides. Together with variations in the peptide sequence at the non-prime side, or its non-peptidic alternatives, this strategy enables targeted tailoring of potent and selective ketoamides towards diverse rhomboid proteases including disease-relevant ones such as PARL and RHBDL2.

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.

Structural insights into the interaction between adenovirus C5 hexon and human lactoferrin.

Adenovirus (AdV) infection of the respiratory epithelium is common but poorly understood. Human AdV species C types, such as HAdV-C5, utilize the Coxsackie-adenovirus receptor (CAR) for attachment and subsequently integrins for entry. CAR and integrins are however located deep within the tight junctions in the mucosa where they would not be easily accessible. Recently, a model for CAR-independent AdV entry was proposed. In this model, human lactoferrin (hLF), an innate immune protein, aids the viral uptake into epithelial cells by mediating interactions between the major capsid protein, hexon, and yet unknown host cellular receptor(s). However, a detailed understanding of the molecular interactions driving this mechanism is lacking. Here, we present a new cryo-EM structure of HAdV-5C hexon at high resolution alongside a hybrid structure of HAdV-5C hexon complexed with human lactoferrin (hLF). These structures reveal the molecular determinants of the interaction between hLF and HAdV-C5 hexon. hLF engages hexon primarily via its N-terminal lactoferricin (Lfcin) region, interacting with hexon's hypervariable region 1 (HVR-1). Mutational analyses pinpoint critical Lfcin contacts and also identify additional regions within hLF that critically contribute to hexon binding. Our study sheds more light on the intricate mechanism by which HAdV-C5 utilizes soluble hLF/Lfcin for cellular entry. These findings hold promise for advancing gene therapy applications and inform vaccine development. IMPORTANCE: Our study delves into the structural aspects of adenovirus (AdV) infections, specifically HAdV-C5 in the respiratory epithelium. It uncovers the molecular details of a novel pathway where human lactoferrin (hLF) interacts with the major capsid protein, hexon, facilitating viral entry, and bypassing traditional receptors such as CAR and integrins. The study's cryo-EM structures reveal how hLF engages hexon, primarily through its N-terminal lactoferricin (Lfcin) region and hexon's hypervariable region 1 (HVR-1). Mutational analyses identify critical Lfcin contacts and other regions within hLF vital for hexon binding. This structural insight sheds light on HAdV-C5's mechanism of utilizing soluble hLF/Lfcin for cellular entry, holding promise for gene therapy and vaccine development advancements in adenovirus research.

Real-time Imaging of Nascent DNA in Live Cells by Monitoring the Fluorescence Lifetime of DNA-Incorporated Thiazole Orange-Modified Nucleotides.

A modified 2'-deoxycytidine triphosphate derivative (dCTO TP) bearing a thiazole orange moiety tethered via an oligoethylene glycol linker was designed and synthesized. The nucleotide was incorporated into DNA by DNA polymerases in vitro as well as in live cells. Upon incorporation of dCTO TP into DNA, the thiazole orange moiety exhibited a fluorescence lifetime that differed significantly from the non-incorporated (i.e. free and non-covalently intercalated) forms of dCTO TP. When dCTO TP was delivered into live U-2 OS cells using a synthetic nucleoside triphosphate transporter, it allowed us to distinguish and monitor cells that were actively synthesizing DNA in real time, from the very first moments after the treatment. We anticipate that this probe could be used to study chromatin organization and dynamics.

Powerful Avidity with a Limited Valency for Virus-Attachment Blockers on DC-SIGN: Combining Chelation and Statistical Rebinding with Structural Plasticity of the Receptor.

The C-type lectin receptor DC-SIGN has been highlighted as the coreceptor for the spike protein of the SARS-CoV-2 virus. A multivalent glycomimetic ligand, Polyman26, has been found to inhibit DC-SIGN-dependent trans-infection of SARS-CoV-2. The molecular details underlying avidity generation in such systems remain poorly characterized. In an effort to dissect the contribution of the known multivalent effects - chelation, clustering, and statistical rebinding - we studied a series of dendrimer constructs related to Polyman26 with a rod core rationally designed to engage simultaneously two binding sites of the tetrameric DC-SIGN. Binding properties of these compounds have been studied with a range of biophysical techniques, including recently developed surface plasmon resonance oriented-surface methodology. Using molecular modeling we addressed, for the first time, the impact of the carbohydrate recognition domains' flexibility of the DC-SIGN tetramer on the compounds' avidity. We were able to gain deeper insight into the role of different binding modes, which in combination produce a construct with a nanomolar affinity despite a limited valency. This multifaceted experimental-theoretical approach provides detailed understanding of multivalent ligand/multimeric protein interactions which can lead to future predictions. This work opens the way to the development of new virus attachment blockers adapted to different C-type lectin receptors of viruses.

Helquat dyes targeting G-quadruplexes as a new class of anti-HIV-1 inhibitors.

The secondary structure of nucleic acids containing quartets of guanines, termed G-quadruplexes, is known to regulate the transcription of many genes. Several G-quadruplexes can be formed in the HIV-1 long terminal repeat promoter region and their stabilization results in the inhibition of HIV-1 replication. Here, we identified helquat-based compounds as a new class of anti-HIV-1 inhibitors that inhibit HIV-1 replication at the stage of reverse transcription and provirus expression. Using Taq polymerase stop and FRET melting assays, we have demonstrated their ability to stabilize G-quadruplexes in the HIV-1 long-terminal repeat sequence. Moreover, these compounds were not binding to the general G-rich region, but rather to G-quadruplex-forming regions. Finally, docking and molecular dynamics calculations indicate that the structure of the helquat core greatly affects the binding mode to the individual G-quadruplexes. Our findings can provide useful information for the further rational design of inhibitors targeting G-quadruplexes in HIV-1.

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.

Targeting the insulin receptor with hormone and peptide dimers.

Insulin is a key hormone involved in the regulation of overall energetic homeostasis of the organism. The dimeric character of the receptor for insulin evokes ideas about its activation or inhibition with peptide dimers that could either trigger or block the structural transition of the insulin receptor, leading to its activation. Herewith, we present the chemical engineering and biological characterization of several series of insulin dimers or dimers of specific peptides that should be able to bind receptors for insulin or insulin growth factor 1. The hormones or peptides in the dimers were interconnected with different linkers, consisting of triazole moieties and 3, 6, 8, 11, or 23 polyethylene glycol units. The prepared dimers were weaker in binding to insulin receptors than human insulin. However, some of the insulin dimers showed preferential binding specificity toward the isoform A of the insulin receptor, and the insulin dimers also stimulated the insulin receptor more strongly than would be consistent with their binding affinities. Our results suggest that designing insulin dimers may be a promising strategy for modulating the ability of the hormone to activate the receptor or to alter its specificity toward insulin receptor isoforms.

In Vitro Evolution Reveals Noncationic Protein-RNA Interaction Mediated by Metal Ions.

RNA-peptide/protein interactions have been of utmost importance to life since its earliest forms, reaching even before the last universal common ancestor (LUCA). However, the ancient molecular mechanisms behind this key biological interaction remain enigmatic because extant RNA-protein interactions rely heavily on positively charged and aromatic amino acids that were absent (or heavily under-represented) in the early pre-LUCA evolutionary period. Here, an RNA-binding variant of the ribosomal uL11 C-terminal domain was selected from an approximately 1010 library of partially randomized sequences, all composed of ten prebiotically plausible canonical amino acids. The selected variant binds to the cognate RNA with a similar overall affinity although it is less structured in the unbound form than the wild-type protein domain. The variant complex association and dissociation are both slower than for the wild-type, implying different mechanistic processes involved. The profile of the wild-type and mutant complex stabilities along with molecular dynamics simulations uncovers qualitative differences in the interaction modes. In the absence of positively charged and aromatic residues, the mutant uL11 domain uses ion bridging (K+/Mg2+) interactions between the RNA sugar-phosphate backbone and glutamic acid residues as an alternative source of stabilization. This study presents experimental support to provide a new perspective on how early protein-RNA interactions evolved, where the lack of aromatic/basic residues may have been compensated by acidic residues plus metal ions.

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.

Engineering the Ligand Specificity of the Human Galectin-1 by Incorporation of Tryptophan Analogues.

Galectin-1 is a ß-galactoside-binding lectin with manifold biological functions. A single tryptophan residue (W68) in its carbohydrate binding site plays a major role in ligand binding and is highly conserved among galectins. To fine tune galectin-1 specificity, we introduced several non-canonical tryptophan analogues at this position of human galectin-1 and analyzed the resulting variants using glycan microarrays. Two variants containing 7-azatryptophan and 7-fluorotryptophan showed a reduced affinity for 3'-sulfated oligosaccharides. Their interaction with different ligands was further analyzed by fluorescence polarization competition assay. Using molecular modeling we provide structural clues that the change in affinities comes from modulated interactions and solvation patterns. Thus, we show that the introduction of subtle atomic mutations in the ligand binding site of galectin-1 is an attractive approach for fine-tuning its interactions with different ligands.

3H-Pyrazolo[4,3-f]quinoline-Based Kinase Inhibitors Inhibit the Proliferation of Acute Myeloid Leukemia Cells In Vivo.

The 3H-pyrazolo[4,3-f]quinoline moiety has been recently shown to be a privileged kinase inhibitor core with potent activities against acute myeloid leukemia (AML) cell lines in vitro. Herein, various 3H-pyrazolo[4,3-f]quinoline-containing compounds were rapidly assembled via the Doebner-Povarov multicomponent reaction from the readily available 5-aminoindazole, ketones, and heteroaromatic aldehydes in good yields. The most active compounds potently inhibit the recombinant FLT3 kinase and its mutant forms with nanomolar IC50 values. Docking studies with the FLT3 kinase showed a type I binding mode, where the 3H-pyrazolo group interacts with Cys694 in the hinge region. The compounds blocked the proliferation of AML cell lines harboring oncogenic FLT3-ITD mutations with remarkable IC50 values, which were comparable to the approved FLT3 inhibitor quizartinib. The compounds also inhibited the growth of leukemia in a mouse-disseminated AML model, and hence, the novel 3H-pyrazolo[4,3-f]quinoline-containing kinase inhibitors are potential lead compounds to develop into anticancer agents, especially for kinase-driven cancers.

Structural determinants for subnanomolar inhibition of the secreted aspartic protease Sapp1p from Candida parapsilosis.

Pathogenic Candida albicans yeasts frequently cause infections in hospitals. Antifungal drugs lose effectiveness due to other Candida species and resistance. New medications are thus required. Secreted aspartic protease of C. parapsilosis (Sapp1p) is a promising target. We have thus solved the crystal structures of Sapp1p complexed to four peptidomimetic inhibitors. Three potent inhibitors (Ki: 0.1, 0.4, 6.6 nM) resembled pepstatin A (Ki: 0.3 nM), a general aspartic protease inhibitor, in terms of their interactions with Sapp1p. However, the weaker inhibitor (Ki: 14.6 nM) formed fewer nonpolar contacts with Sapp1p, similarly to the smaller HIV protease inhibitor ritonavir (Ki: 1.9 µM), which, moreover, formed fewer H-bonds. The analyses have revealed the structural determinants of the subnanomolar inhibition of C. parapsilosis aspartic protease. Because of the high similarity between Saps from different Candida species, these results can further be used for the design of potent and specific Sap inhibitor-based antimycotic drugs.

Multipodal insulin mimetics built on adamantane or proline scaffolds.

Multi-orthogonal molecular scaffolds can be applied as core structures of bioactive compounds. Here, we prepared four tri-orthogonal scaffolds based on adamantane or proline skeletons. The scaffolds were used for the solid-phase synthesis of model insulin mimetics bearing two different peptides on the scaffolds. We found that adamantane-derived compounds bind to the insulin receptor more effectively (Kd value of 0.5 µM) than proline-derived compounds (Kd values of 15-38 µM) bearing the same peptides. Molecular dynamics simulations suggest that spacers between peptides and central scaffolds can provide greater flexibility that can contribute to increased binding affinity. Molecular modeling showed possible binding modes of mimetics to the insulin receptor. Our data show that the structure of the central scaffold and flexibility of attached peptides in this type of compound are important and that different scaffolds should be considered when designing peptide hormone mimetics.

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.

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 (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.

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.

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.

A rapid synthesis of low-nanomolar divalent LecA inhibitors in four linear steps from d-galactose pentaacetate.

Chronic infections with Pseudomonas aeruginosa are associated with the formation of bacterial biofilms. The tetrameric P. aeruginosa lectin LecA is a virulence factor and an anti-biofilm drug target. Increasing the overall binding affinity by multivalent presentation of binding epitopes can enhance the weak carbohydrate-ligand interactions. Low-nanomolar divalent LecA ligands/inhibitors with up to 260-fold valency-normalized potency boost and excellent selectivity over human galectin-1 were synthesized from d-galactose pentaacetate and benzaldehyde-based linkers in four linear steps.

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.