Machine learning predictions of MHC-II specificities reveal alternative binding mode of class II epitopes.

CD4+ T cells orchestrate the adaptive immune response against pathogens and cancer by recognizing epitopes presented on class II major histocompatibility complex (MHC-II) molecules. The high polymorphism of MHC-II genes represents an important hurdle toward accurate prediction and identification of CD4+ T cell epitopes. Here we collected and curated a dataset of 627,013 unique MHC-II ligands identified by mass spectrometry. This enabled us to precisely determine the binding motifs of 88 MHC-II alleles across humans, mice, cattle, and chickens. Analysis of these binding specificities combined with X-ray crystallography refined our understanding of the molecular determinants of MHC-II motifs and revealed a widespread reverse-binding mode in HLA-DP ligands. We then developed a machine-learning framework to accurately predict binding specificities and ligands of any MHC-II allele. This tool improves and expands predictions of CD4+ T cell epitopes and enables us to discover viral and bacterial epitopes following the aforementioned reverse-binding mode.

Structures of calmodulin-melittin complexes show multiple binding modes lacking classical anchoring interactions.

Calmodulin (CaM) is a Ca2+ sensor protein found in all eukaryotic cells that regulates a large number of target proteins in a Ca2+ concentration-dependent manner. As a transient-type hub protein, it recognizes linear motifs of its targets, though for the Ca2+-dependent binding, no consensus sequence was identified. Its complex with melittin, a major component of bee venom, is often used as a model system of protein-protein complexes. Yet, the structural aspects of the binding are not well understood, as only diverse, low-resolution data are available concerning the association. We present the crystal structure of melittin in complex with Ca2+-saturated CaMs from two, evolutionarily distant species, Homo sapiens and Plasmodium falciparum, representing three binding modes of the peptide. Results-augmented by molecular dynamics simulations-indicate that multiple binding modes can exist for CaM-melittin complexes, as an intrinsic characteristic of the binding. While the helical structure of melittin remains, swapping of its salt bridges and partial unfolding of its C-terminal segment can occur. In contrast to the classical way of target recognition by CaM, we found that different sets of residues can anchor at the hydrophobic pockets of CaM, which were considered as main recognition sites. Finally, the nanomolar binding affinity of the CaM-melittin complex is created by an ensemble of arrangements of similar stability-tight binding is achieved not by optimized specific interactions but by simultaneously satisfying less optimal interaction patterns in co-existing different conformers.

Favipiravir Analogues as Inhibitors of SARS-CoV-2 RNA-Dependent RNA Polymerase, Combined Quantum Chemical Modeling, Quantitative Structure-Property Relationship, and Molecular Docking Study.

Our study was motivated by the urgent need to develop or improve antivirals for effective therapy targeting RNA viruses. We hypothesized that analogues of favipiravir (FVP), an inhibitor of RNA-dependent RNA polymerase (RdRp), could provide more effective nucleic acid recognition and binding processes while reducing side effects such as cardiotoxicity, hepatotoxicity, teratogenicity, and embryotoxicity. We proposed a set of FVP analogues together with their forms of triphosphate as new SARS-CoV-2 RdRp inhibitors. The main aim of our study was to investigate changes in the mechanism and binding capacity resulting from these modifications. Using three different approaches, QTAIM, QSPR, and MD, the differences in the reactivity, toxicity, binding efficiency, and ability to be incorporated by RdRp were assessed. Two new quantum chemical reactivity descriptors, the relative electro-donating and electro-accepting power, were defined and successfully applied. Moreover, a new quantitative method for comparing binding modes was developed based on mathematical metrics and an atypical radar plot. These methods provide deep insight into the set of desirable properties responsible for inhibiting RdRp, allowing ligands to be conveniently screened. The proposed modification of the FVP structure seems to improve its binding ability and enhance the productive mode of binding. In particular, two of the FVP analogues (the trifluoro- and cyano-) bind very strongly to the RNA template, RNA primer, cofactors, and RdRp, and thus may constitute a very good alternative to FVP.

Elucidating the Role of Noncovalent Interactions in Favipiravir, a Drug Active against Various Human RNA Viruses; a 1H-14N NQDR/Periodic DFT/QTAIM/RDS/3D Hirshfeld Surfaces Combined Study.

Favipiravir (6-fluoro-3-hydroxypyrazine-2-carboxamide, FPV), an active pharmaceutical component of the drug discovered and registered in March 2014 in Japan under the name Avigan, with an indication for pandemic influenza, has been studied. The study of this compound was prompted by the idea that effective processes of recognition and binding of FPV to the nucleic acid are affected predominantly by the propensity to form intra- and intermolecular interactions. Three nuclear quadrupole resonance experimental techniques, namely 1H-14N cross-relaxation, multiple frequency sweeps, and two-frequency irradiation, followed by solid-state computational modelling (density functional theory supplemented by the quantum theory of atoms in molecules, 3D Hirshfeld Surfaces, and reduced density gradient) approaches were applied. The complete NQR spectrum consisting of nine lines indicating the presence of three chemically inequivalent nitrogen sites in the FPV molecule was detected, and the assignment of lines to particular sites was performed. The description of the nearest vicinity of all three nitrogen atoms was used to characterize the nature of the intermolecular interactions from the perspective of the local single atoms and to draw some conclusions on the nature of the interactions required for effective recognition and binding. The propensity to form the electrostatic N-H···O, N-H···N, and C-H···O intermolecular hydrogen bonds competitive with two intramolecular hydrogen bonds, strong O-H···O and very weak N-H···N, closing the 5-member ring and stiffening the structure, as well as π···π and F···F dispersive interactions, were analysed in detail. The hypothesis regarding the similarity of the interaction pattern in the solid and the RNA template was verified. It was discovered that the -NH2 group in the crystal participates in intermolecular hydrogen bonds N-H···N and N-H···O, in the precatalytic state only in N-H···O, while in the active state in N-H···N and N-H···O hydrogen bonds, which is of importance to link FVP to the RNA template. Our study elucidates the binding modes of FVP (in crystal, precatalytic, and active forms) in detail and should guide the design of more potent analogues targeting SARS-CoV-2. Strong direct binding of FVP-RTP to both the active site and cofactor discovered by us suggests a possible alternative, allosteric mechanism of FVP action, which may explain the scattering of the results of clinical trials or the synergistic effect observed in combined treatment against SARS-CoV-2.

Dinuclear platinum(II) complexes as the pattern for phosphate backbone binding: a new perspective for recognition of binding modes to DNA.

The mechanism of action of most approved drugs in use today is based on their binding to specific proteins or DNA. One of the achievements of this research is a new perspective for recognition of binding modes to DNA by monitoring of changes in measured and stoichiometric values of absorbance at 260 nm. UV-Vis and IR spectroscopy, gel electrophoresis and docking study were used for investigation of binding properties of three dinuclear platinum(II) complexes containing different pyridine-based bridging ligands, [{Pt(en)Cl}2(µ-4,4'-bipy)]Cl2·2H2O (Pt1), [{Pt(en)Cl}2(µ-bpa)]Cl2·4H2O (Pt2) and [{Pt(en)Cl}2(µ-bpe)]Cl2·4H2O (Pt3) to DNA (4,4'-bipy, bpa and bpe are 4,4'-bipyridine, 1,2-bis(4-pyridyl)ethane and 1,2-bis(4-pyridyl)ethene, respectively). In contrast to the system with well-known intercalated ligand (EtBr), covalently bound ligand (cis-Pt) and with minor groove binder (Hoechst 33258), which do not have significant differences in measured and stoichiometric values, the most pronounced deviations are recorded for two dinuclear platinum(II) complexes (Pt1 and Pt2), as a consequence of complex binding to the phosphate backbone and bending of DNA helix. The hydrolysis of complexes and changes in DNA conformation were also analysed as phenomena that may have an impact on the changes in absorbance.

Outliers in SAR and QSAR: 4. effects of allosteric protein-ligand interactions on the classical quantitative structure-activity relationships.

Effects of allosteric interactions on the classical structure-activity relationship (SAR) and quantitative SAR (QSAR) have been investigated. Apprehending the outliers in SAR and QSAR studies can improve the quality, predictability, and use of QSAR in designing unknown compounds in drug discovery research. We explored allosteric protein-ligand interactions as a possible source of outliers in SAR/QSAR. We used glycogen phosphorylase as an example of a protein that has an allosteric site. Examination of the ligand-bound x-ray crystal structures of glycogen phosphorylase revealed that many inhibitors bound at more than one binding site. The results of QSAR analyses of the inhibitors included a QSAR that recognized an outlier bound at a distinctive allosteric binding site. The case provided an example of constructive use of QSAR identifying outliers with alternative binding modes. Other allosteric QSARs that captured our attention were the inverted parabola/bilinear QSARs. The x-ray crystal structures and the QSAR analyses indicated that the inverted parabola QSARs could be associated with the conformational changes in the allosteric interactions. Our results showed that the normal parabola, as well as the inverted parabola QSARs, can describe the allosteric interactions. Examination of the ligand-bound X-ray crystal structures of glycogen phosphorylase revealed that many inhibitors bound at more than one binding site. The results of QSAR analyses of the inhibitors included a QSAR that recognized an outlier bound at a distinctive allosteric binding site.

Indole-Based Tubulin Inhibitors: Binding Modes and SARs Investigations.

Tubulin inhibitors can interfere with normal cell mitosis and inhibit cell proliferation through interfering with the normal structure and function of microtubules, forming spindle filaments. Indole, as a privileged pharmacological skeleton, has been widely used in anti-cancer inhibitors. A variety of alkaloids containing an indole core obtained from natural sources have been proven to inhibit tubulin polymerization, and an ever-increasing number of synthetic indole-based tubulin inhibitors have been reported. Among these, several kinds of indole-based derivatives, such as TMP analogues, aroylindoles, arylthioindoles, fused indole, carbazoles, azacarbolines, alkaloid nortopsentin analogues and bis-indole derivatives, have shown good inhibition activities towards tubulin polymerization. The binding modes and SARs investigations of synthetic indole derivatives, along with a brief mechanism on their anti-tubulin activity, are presented in this review.

Machine learning predicts nucleosome binding modes of transcription factors.

BACKGROUND: Most transcription factors (TFs) compete with nucleosomes to gain access to their cognate binding sites. Recent studies have identified several TF-nucleosome interaction modes including end binding (EB), oriented binding, periodic binding, dyad binding, groove binding, and gyre spanning. However, there are substantial experimental challenges in measuring nucleosome binding modes for thousands of TFs in different species. RESULTS: We present a computational prediction of the binding modes based on TF protein sequences. With a nested cross-validation procedure, our model outperforms several fine-tuned off-the-shelf machine learning (ML) methods in the multi-label classification task. Our binary classifier for the EB mode performs better than these ML methods with the area under precision-recall curve achieving 75%. The end preference of most TFs is consistent with low nucleosome occupancy around their binding site in GM12878 cells. The nucleosome occupancy data is used as an alternative dataset to confirm the superiority of our EB classifier. CONCLUSIONS: We develop the first ML-based approach for efficient and comprehensive analysis of nucleosome binding modes of TFs.

What coronavirus 3C-like protease tells us: From structure, substrate selectivity, to inhibitor design.

The emergence of a variety of coronaviruses (CoVs) in the last decades has posed huge threats to human health. Especially, the ongoing pandemic of coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has led to more than 70 million infections and over 1.6 million of deaths worldwide in the past few months. None of the efficacious antiviral agents against human CoVs have been approved yet. 3C-like protease (3CLpro ) is an attractive target for antiviral intervention due to its essential role in processing polyproteins translated from viral RNA, and its conserved structural feature and substrate specificity among CoVs in spite of the sequence variation. This review focuses on all available crystal structures of 12 CoV 3CLpro s and their inhibitors, and intends to provide a comprehensive understanding of this protease from multiple aspects including its structural features, substrate specificity, inhibitor binding modes, and more importantly, to recapitulate the similarity and diversity among different CoV 3CLpro s and the structure-activity relationship of various types of inhibitors. Such an attempt could gain a deep insight into the inhibition mechanisms and drive future structure-based drug discovery targeting 3CLpro s.

Superprotonic Conductivity of MOF-808 Achieved by Controlling the Binding Mode of Grafted Sulfamate.

A metal-organic framework (MOF) having superprotonic conductivity, MOF-808, is prepared by modulating the binding mode of the sulfamate (SA) moieties grafted onto the metal clusters. The activation of the SA-grafted MOF-808 at 150 °C changes the binding mode of the grafted SA from monodentate to bridging bidentate, thus converting the neutral amido (-S-NH2 ) moiety of the grafted SA to the more acidic cationic sulfiliminium (-S=NH2+ ) moiety. Further, the acidic sulfiliminium moiety of MOF-808-4SA-150 results in more efficient proton conduction than the amido moiety of MOF-808-4SA-60. At 60 °C and 95 % relative humidity, MOF-808-4SA-150 is found to have a proton conductivity of 7.89×10-2  S cm-1 , which is more than 30-times higher than that of MOF-808-4SA-60. Moreover, this superprotonic conductivity is well maintained over 1000 cycles of conductivity measurements and for similar cyclic measurements each day for seven days.

Binding Modes of Salicylic Acids to Titanium Oxide Molecular Surfaces.

A set of titanium oxide clusters (TOCs) comprised of 4 to 16 Ti atoms are synthesized with substituted salicylates (SSAs). The interfacial coordination environment of these SSA/Ti oxide hybrids are surveyed and found to be limited to four binding modes, with the bridging chelate mode being the most common one. The SSA-functionalized TOCs show strong visible light absorption properties. The contribution of the SSAs in the frontier orbitals of the TOCs are analyzed by using TD-DFT calculations based on the molecular geometries determined by X-ray diffraction. For TOCs of relatively high O/Ti ratio, the SSAs narrow the band gap of the TOCs by contributing solely to the HOMOs. Both binding modes and locations of the SSAs are important for the roles of SSAs in changing the HOMOs and thereby the absorption onsets.

Discovery of isoquinolinone and naphthyridinone-based inhibitors of poly(ADP-ribose) polymerase-1 (PARP1) as anticancer agents: Structure activity relationship and preclinical characterization.

The exploitation of GLU988 and LYS903 residues in PARP1 as targets to design isoquinolinone (I & II) and naphthyridinone (III) analogues is described. Compounds of structure I have good biochemical and cellular potency but suffered from inferior PK. Constraining the linear propylene linker of structure I into a cyclopentene ring (II) offered improved PK parameters, while maintaining potency for PARP1. Finally, to avoid potential issues that may arise from the presence of an anilinic moiety, the nitrogen substituent on the isoquinolinone ring was incorporated as part of the bicyclic ring. This afforded a naphthyridinone scaffold, as shown in structure III. Further optimization of naphthyridinone series led to identification of a novel and highly potent PARP1 inhibitor 34, which was further characterized as preclinical candidate molecule. Compound 34 is orally bioavailable and displayed favorable pharmacokinetic (PK) properties. Compound 34 demonstrated remarkable antitumor efficacy both as a single-agent as well as in combination with chemotherapeutic agents in the BRCA1 mutant MDA-MB-436 breast cancer xenograft model. Additionally, compound 34 also potentiated the effect of agents such as temozolomide in breast cancer, pancreatic cancer and Ewing's sarcoma models.

Recent advances in the development of ligands specifically targeting telomeric multimeric G-quadruplexes.

Long human telomeric DNA sequence could form higher-order G-quadruplex structures, namely telomeric multimeric G-quadruplexes. The formation of telomeric multimeric G-quadruplexes has been demonstrated. Several efforts have been devoted to the development of ligands targeting telomeric multimeric G-quadruplexes in recent years. The reported ligands specifically targeting telomeric multimeric G-quadruplexes exhibited either high anticancer activity with effective stabilization ability or distinct fluorescence responses to telomeric multimeric G-quadruplexes. In this review, the ligands including three types of small molecules are discussed which focus on their structural features and binding modes. Furthermore, we put forward the promising prospects and current challenges. This presented review might provide new strategies to exploit more sophisticated ligands targeting telomeric multimeric G-quadruplexes.

Interactions of mono spermine porphyrin derivative with DNAs.

In this work, we have characterized the interactions of monospermine porphyrin derivative with calf thymus DNA (ct-DNA) and poly (dG-dC)2 in both B and Z conformation. By several spectroscopic techniques (UV-vis, electronic circular dichroism and resonance light scattering), the binding modes of monospermine porphyrin derivative with different DNA sequences have been elucidated. In the presence of ct-DNA, the porphyrin binds along the external double helix as well as in the presence of B conformation of poly (dG-dC)2 . Whilst when the Z form of the poly (dG-dC)2 is induced, a slight intercalation of the porphyrin between the basis has been detected.

Binding modes identification through molecular dynamic simulations: A case study with carnosine enantiomers and the Teicoplanin A2-2-based chiral stationary phase.

In the present study, an in silico methodology able to define the binding modes adopted by carnosine enantiomers in the setting of the chiral recognition process is described. The inter- and intramolecular forces involved in the enantioseparation process with the Teicoplanin A2-2 chiral selector and carnosine as model compound are successfully identified. This approach fully rationalizes, at a molecular level, the (S) < (R) enantiomeric elution order obtained under reversed-phase conditions. Consistent explanations were achieved by managing molecular dynamics results with advanced techniques of data analysis. As a result, the time-dependent identification of all the interactions simultaneously occurring in the chiral selector-enantiomeric analyte binding process was obtained. Accordingly, it was found that only (R)-carnosine is able to engage a stabilizing charge-charge interaction through its ionized imidazole ring with the carboxylate counter-part on the chiral selector. Instead, (S)-carnosine establishes intramolecular contacts between its ionized functional groups, that limit its conformational freedom and impair the association with the chiral selector unit.

X-ray Structure-Based Chemoinformatic Analysis Identifies Promiscuous Ligands Binding to Proteins from Different Classes with Varying Shapes.

(1) Background: Compounds with multitarget activity are of interest in basic research to explore molecular foundations of promiscuous binding and in drug discovery as agents eliciting polypharmacological effects. Our study has aimed to systematically identify compounds that form complexes with proteins from distinct classes and compare their bioactive conformations and molecular properties. (2) Methods: A large-scale computational investigation was carried out that combined the analysis of complex X-ray structures, ligand binding modes, compound activity data, and various molecular properties. (3) Results: A total of 515 ligands with multitarget activity were identified that included 70 organic compounds binding to proteins from different classes. These multiclass ligands (MCLs) were often flexible and surprisingly hydrophilic. Moreover, they displayed a wide spectrum of binding modes. In different target structure environments, binding shapes of MCLs were often similar, but also distinct. (4) Conclusions: Combined structural and activity data analysis identified compounds with activity against proteins with distinct structures and functions. MCLs were found to have greatly varying shape similarity when binding to different protein classes. Hence, there were no apparent canonical binding shapes indicating multitarget activity. Rather, conformational versatility characterized MCL binding.

Functionalization of Titanium Oxide Cluster Ti17 O24 (Oi C3 H7 )20 with Catechols: Structures and Ligand-Exchange Reactivities.

The functionalization of Ti17 O24 (Oi C3 H7 )20 (Ti17 ) with substituted catechols is studied by using crystallography, Raman spectroscopy, and stopped-flow kinetics. The knowledges on the number of accessible functionalities, their exact location correlated with their Raman assignment, and the kinetic parameters are acquired. A catecholate ligand binds to a five-coordinated surface Ti of Ti17 (denoted as Tia site) adopting the mono-protonated, chelate-bidentate binding mode, whereas it binds to a six-coordinated surface-Ti (denoted as Tib site) adopting the mono-dentate mode. With low numbers of equivalents of added catechols the Tia sites show higher reactivity than the Tib sites toward functionalization. Two binding modes may co-exist and equilibrate in solution. Our results also imply that at most eight of the twenty Oi C3 H7 ligands of Ti17 are exchangeable without damage of the core structure. The kinetic studies point out that the ligand-exchange reaction is second order and occurs very fast. The current findings are helpful for the controlled functionalization of Ti17 and other Ti oxide clusters, and the further application of them as building blocks in supramolecular chemistry for the assembly of well-defined organic-inorganic hybrid materials.

A Metadynamics-Based Protocol for the Determination of GPCR-Ligand Binding Modes.

G protein-coupled receptors (GPCRs) are a main drug target and therefore a hot topic in pharmaceutical research. One important prerequisite to understand how a certain ligand affects a GPCR is precise knowledge about its binding mode and the specific underlying interactions. If no crystal structure of the respective complex is available, computational methods can be used to deduce the binding site. One of them are metadynamics simulations which have the advantage of an enhanced sampling compared to conventional molecular dynamics simulations. However, the enhanced sampling of higher-energy states hampers identification of the preferred binding mode. Here, we present a novel protocol based on clustering of multiple walker metadynamics simulations which allows identifying the preferential binding mode from such conformational ensembles. We tested this strategy for three different model systems namely the histamine H1 receptor in combination with its physiological ligand histamine, as well as the ß 2 adrenoceptor with its agonist adrenaline and its antagonist alprenolol. For all three systems, the proposed protocol was able to reproduce the correct binding mode known from the literature suggesting that the approach can more generally be applied to the prediction of GPCR ligand binding in future.

Naphthalene Diimides as Multimodal G-Quadruplex-Selective Ligands.

G-quadruplexes are four-stranded nucleic acids structures that can form in guanine-rich sequences. Following the observation that G-quadruplexes are particularly abundant in genomic regions related to cancer, such as telomeres and oncogenes promoters, several G-quadruplex-binding molecules have been developed for therapeutic purposes. Among them, naphthalene diimide derivatives have reported versatility, consistent selectivity and high affinity toward the G-quadruplex structures. In this review, we present the chemical features, synthesis and peculiar optoelectronic properties (absorption, emission, redox) that make naphtalene diimides so versatile for biomedical applications. We present the latest developments on naphthalene diimides as G-quadruplex ligands, focusing on their ability to bind G-quadruplexes at telomeres and oncogene promoters with consequent anticancer activity. Their different binding modes (reversible versus irreversible/covalent) towards G-quadruplexes and their additional use as antimicrobial agents are also presented and discussed.

Reduced state transition barrier of CDK6 from open to closed state induced by Thr177 phosphorylation and its implication in binding modes of inhibitors.

BACKGROUND: CDK6 is considered as a highly validated anticancer drug target due to its essential role in regulating cell cycle progression at G1 restriction point. Activation of CDK6 requires the phosphorylation of Thr177 on A-loop, but the structural insights of the activation mechanism remain unclear. METHODS: Herein, all-atoms molecular dynamics (MD) simulations were used to study the effects of Thr177 phosphorylation on the dynamic structure of CDK6-Vcyclin complex. RESULTS: MD results indicated that the free energy barrier of the transition from open to closed state decreased ~47.2% after Thr177 phosphorylation. Key steps along the state transition process were obtained from a cluster analysis. Binding preference of ten different inhibitors to open or closed state were also investigated through molecular docking along with MD simulations methods. CONCLUSIONS: Our results indicated that Thr177 phosphorylation increased the flexibility around the ATP-binding pocket. The transition of the ATP-binding pocket between open and closed states should be considered for understanding the binding of CDK6 inhibitors. GENERAL SIGNIFICANCE: This work could deepen the understanding of CDKs activation mechanism, and provide useful information for the discovery of new CDKs inhibitors with high affinity and specificity.