MM-ISMSA: An Ultrafast and Accurate Scoring Function for Protein–Protein Docking

Acquired Immune Deficiency Syndrome correlation with SARS-CoV-2 N genotypes

BACKGROUND

Epigenetics and clinical observations referring to Betacoronavirus lead to the conjecture that Sarbecovirus may have the ability to infect lymphocytes using a different way than the spike protein. In addition to inducing the death of lymphocytes, thus drastically reducing their population and causing a serious immune deficiency, allows it to remain hidden for long periods of latency using them as a viral reservoir in what is named Long-Covid Disease. Exploring possibilities, the hypothesis is focused on that N protein may be the key of infecting lymphocytes.

METHOD

The present article exhibits a computational assay for the latest complete sequences reported to GISAID, correlating N genotypes with an enhancement in the affinity of the complex that causes immune deficiency in order to determine a good docking with the N protein and some receptors in lymphocytes.

RESULTS

A novel high-interaction coupling of N-RBD and CD147 is presented as the main way of infecting lymphocytes, allowing to define those genotypes involved in their affinity enhancement.

CONCLUSION

The hypothesis is consistent with the mutagenic deriving observed on the in-silico assay, which reveals that genotypes N/120 and N/152 are determinant to reduce the Immune Response of the host infecting lymphocytes, allowing the virus persists indefinitely and causing an Acquire Immune Deficiency Syndrome.

Flexible structural arrangement and DNA-binding properties of protein p6 from Bacillus subtillis phage φ29

The genome-organizing protein p6 of Bacillus subtilis bacteriophage φ29 plays an essential role in viral development by activating the initiation of DNA replication and participating in the early-to-late transcriptional switch. These activities require the formation of a nucleoprotein complex in which the DNA adopts a right-handed superhelix wrapping around a multimeric p6 scaffold, restraining positive supercoiling and compacting the viral genome. Due to the absence of homologous structures, prior attempts to unveil p6's structural architecture failed. Here, we employed AlphaFold2 to engineer rational p6 constructs yielding crystals for three-dimensional structure determination. Our findings reveal a novel fold adopted by p6 that sheds light on its self-association mechanism and its interaction with DNA. By means of protein-DNA docking and molecular dynamic simulations, we have generated a comprehensive structural model for the nucleoprotein complex that consistently aligns with its established biochemical and thermodynamic parameters. Besides, through analytical ultracentrifugation, we have confirmed the hydrodynamic properties of the nucleocomplex, further validating in solution our proposed model. Importantly, the disclosed structure not only provides a highly accurate explanation for previously experimental data accumulated over decades, but also enhances our holistic understanding of the structural and functional attributes of protein p6 during φ29 infection.

Good Vibrations Report on the DNA Quadruplex Binding of an Excited State Amplified Ruthenium Polypyridyl IR Probe

The nitrile containing Ru(II)polypyridyl complex [Ru(phen)2(11,12-dCN-dppz)]2+ (1) is shown to act as a sensitive infrared probe of G-quadruplex (G4) structures. UV-visible absorption spectroscopy reveals enantiomer sensitive binding for the hybrid htel(K) and antiparallel htel(Na) G4s formed by the human telomer sequence d[AG3(TTAG3)3]. Time-resolved infrared (TRIR) of 1 upon 400 nm excitation indicates dominant interactions with the guanine bases in the case of Λ-1/htel(K), Δ-1/htel(K), and Λ-1/htel(Na) binding, whereas Δ-1/htel(Na) binding is associated with interactions with thymine and adenine bases in the loop. The intense nitrile transient at 2232 cm-1 undergoes a linear shift to lower frequency as the solution hydrogen bonding environment decreases in DMSO/water mixtures. This shift is used as a sensitive reporter of the nitrile environment within the binding pocket. The lifetime of 1 in D2O (ca. 100 ps) is found to increase upon DNA binding, and monitoring of the nitrile and ligand transients as well as the diagnostic DNA bleach bands shows that this increase is related to greater protection from the solvent environment. Molecular dynamics simulations together with binding energy calculations identify the most favorable binding site for each system, which are in excellent agreement with the observed TRIR solution study. This study shows the power of combining the environmental sensitivity of an infrared (IR) probe in its excited state with the TRIR DNA "site effect" to gain important information about the binding site of photoactive agents and points to the potential of such amplified IR probes as sensitive reporters of biological environments.

Molecular basis of the final step of cell division in Streptococcus pneumoniae

Bacterial cell-wall hydrolases must be tightly regulated during bacterial cell division to prevent aberrant cell lysis and to allow final separation of viable daughter cells. In a multidisciplinary work, we disclose the molecular dialogue between the cell-wall hydrolase LytB, wall teichoic acids, and the eukaryotic-like protein kinase StkP in Streptococcus pneumoniae. After characterizing the peptidoglycan recognition mode by the catalytic domain of LytB, we further demonstrate that LytB possesses a modular organization allowing the specific binding to wall teichoic acids and to the protein kinase StkP. Structural and cellular studies notably reveal that the temporal and spatial localization of LytB is governed by the interaction between specific modules of LytB and the final PASTA domain of StkP. Our data collectively provide a comprehensive understanding of how LytB performs final separation of daughter cells and highlights the regulatory role of eukaryotic-like kinases on lytic machineries in the last step of cell division in streptococci.

Molecular dynamics simulations reveal the impact of NUDT15 R139C and R139H variants in structural conformation and dynamics

NUDT15, also known as MTH2, is a member of the NUDIX protein family that catalyzes the hydrolysis of nucleotides and deoxynucleotides, as well as thioguanine analogues. NUDT15 has been reported as a DNA sanitizer in humans, and more recent studies have shown that some genetic variants are related to a poor prognosis in neoplastic and immunologic diseases treated with thioguanine drugs. Despite this, the role of NUDT15 in physiology and molecular biology is quite unclear, as is the mechanism of action of this enzyme. The existence of clinically relevant variants has prompted the study of these enzymes, whose capacity to bind and hydrolyze thioguanine nucleotides is still poorly understood. By using a combination of biomolecular modeling techniques and molecular dynamics, we have studied the monomeric wild type NUDT15 as well as two important variants, R139C and R139H. Our findings reveal not only how nucleotide binding stabilizes the enzyme but also how two loops are responsible for keeping the enzyme in a packed, close conformation. Mutations in α2 helix affect a network of hydrophobic and π-interactions that enclose the active site. This knowledge contributes to the understanding of NUDT15 structural dynamics and will be valuable for the design of new chemical probes and drugs targeting this protein.Communicated by Ramaswamy H. Sarma.

Chemical modulation of microtubule structure through the laulimalide/peloruside site

Taxanes are microtubule-stabilizing agents used in the treatment of many solid tumors, but they often involve side effects affecting the peripheral nervous system. It has been proposed that this could be related to structural modifications on the filament upon drug binding. Alternatively, laulimalide and peloruside bind to a different site also inducing stabilization, but they have not been exploited in clinics. Here, we use a combination of the parental natural compounds and derived analogs to unravel the stabilization mechanism through this site. These drugs settle lateral interactions without engaging the M loop, which is part of the key and lock involved in the inter-protofilament contacts. Importantly, these drugs can modulate the angle between protofilaments, producing microtubules of different diameters. Among the compounds studied, we have found some showing low cytotoxicity and able to induce stabilization without compromising microtubule native structure. This opens the window of new applications for microtubule-stabilizing agents beyond cancer treatment.

Allosteric Inhibition of Neutral Sphingomyelinase 2 (nSMase2) by DPTIP: From Antiflaviviral Activity to Deciphering Its Binding Site through In Silico Studies and Experimental Validation

Flavivirus comprises globally emerging and re-emerging pathogens such as Zika virus (ZIKV), Dengue virus (DENV), and West Nile virus (WNV), among others. Although some vaccines are available, there is an unmet medical need as no effective antiviral treatment has been approved for flaviviral infections. The development of host-directed antivirals (HDAs) targeting host factors that are essential for viral replication cycle offers the opportunity for the development of broad-spectrum antivirals. In the case of flaviviruses, recent studies have revealed that neutral sphingomyelinase 2, (nSMase2), involved in lipid metabolism, plays a key role in WNV and ZIKV infection. As a proof of concept, we have determined the antiviral activity of the non-competitive nSMase2 inhibitor DPTIP against WNV and ZIKV virus. DPTIP showed potent antiviral activity with EC₅₀ values of 0.26 µM and 1.56 µM for WNV and ZIKV, respectively. In order to unravel the allosteric binding site of DPTIP in nSMase2 and the details of the interaction, computational studies have been carried out. These studies have revealed that DPTIP could block the DK switch in nSMase2. Moreover, the analysis of the residues contributing to the binding identified His463 as a crucial residue. Interestingly, the inhibitory activity of DPTIP on the H463A mutant protein supported our hypothesis. Thus, an allosteric cavity in nSMase2 has been identified that can be exploited for the development of new inhibitors with anti-flaviviral activity.

Glycan-Adhering Lectins and Experimental Evaluation of a Lectin FimH Inhibitor in Enterohemorrhagic Escherichia coli (EHEC) O157:H7 Strain EDL933

In this study, we tried to develop a FimH inhibitor that inhibits adhesion of enterohemorrhagic Escherichia coli (EHEC) on the epithelium of human intestine during the initial stage of infections. Using a T7 phage display method with a reference strain, EHEC EDL933, FimH was selected as an adherent lectin to GM1a and Gb3 glycans. In order to detect the ligand binding domain (LBD) of FimH, we used a docking simulation and found three binding site sequences of FimH, i.e., P1, P2, and P3. Among Gb3 mimic peptides, P2 was found to have the strongest binding strength. Moreover, in vitro treatment with peptide P2 inhibited binding activity in a concentration-dependent manner. Furthermore, we conducted confirmation experiments through several strains isolated from patients in Korea, EHEC NCCP15736, NCCP15737, and NCCP15739. In addition, we analyzed the evolutionary characteristics of the predicted FimH lectin-like adhesins to construct a lectin-glycan interaction (LGI). We selected 70 recently differentiated strains from the phylogenetic tree of 2240 strains with Shiga toxin in their genome. We can infer EHEC strains dynamically evolved but FimH was conserved during the evolution time according to the phylogenetic tree. Furthermore, FimH could be a reliable candidate of drug target in terms of evolution. We examined how pathogen lectins interact with host glycans early in infection in EDL933 as well as several field strains and confirmed that glycan-like peptides worked as an initial infection inhibitor.

Unravelling the Turn-On Fluorescence Mechanism of a Fluorescein-Based Probe in GABAA Receptors

GABAA (γ-aminobutyric acid type A) receptors are ligand-gated ion channels mediating fast inhibitory transmission in the mammalian brain. Here we report the molecular and electronic mechanism governing the turn-on emission of a fluorescein-based imaging probe able to target the human GABAA receptor. Multiscale calculations evidence a drastic conformational change of the probe from folded in solution to extended upon binding to the receptor. Intramolecular ππ-stacking interactions present in the folded probe are responsible for quenching fluorescence in solution. In contrast, unfolding within the GABAA receptor changes the nature of the bright excited state triggering emission. Remarkably, this turn-on effect only manifests for the dianionic prototropic form of the imaging probe, which is found to be the strongest binder to the GABAA receptor. This study is expected to assist the design of new photoactivatable screening tools for allosteric modulators of the GABAA receptor.

Unravelling the Turn‐On Fluorescence Mechanism of a Fluorescein‐Based Probe in GABA A Receptors

GABA_A (γ‐aminobutyric acid type A) receptors are ligand‐gated ion channels mediating fast inhibitory transmission in the mammalian brain. Here we report the molecular and electronic mechanism governing the turn‐on emission of a fluorescein‐based imaging probe able to target the human GABA_A receptor. Multiscale calculations evidence a drastic conformational change of the probe from folded in solution to extended upon binding to the receptor. Intramolecular ππ‐stacking interactions present in the folded probe are responsible for quenching fluorescence in solution. In contrast, unfolding within the GABA_A receptor changes the nature of the bright excited state triggering emission. Remarkably, this turn‐on effect only manifests for the dianionic prototropic form of the imaging probe, which is found to be the strongest binder to the GABA_A receptor. This study is expected to assist the design of new photoactivatable screening tools for allosteric modulators of the GABA_A receptor.

Insight into the mechanism of molecular recognition between human Integrin-Linked Kinase and Cpd22 and its implication at atomic level

Pseudokinases have received increasing attention over the past decade because of their role in different physiological phenomena. Although pseudokinases lack several active-site residues, thereby hindering their catalytic activity, recent discoveries have shown that these proteins can play a role in intracellular signaling thanks to their non-catalytic functions. Integrin-linked kinase (ILK) was discovered more than two decades ago and was subsequently validated as a promising target for neoplastic diseases. Since then, only a few small-molecule inhibitors have been described, with the V-shaped pyrazole Cpd22 being the most interesting and characterized. However, little is known about its detailed mechanism of action at atomic level. In this study, using a combination of computational chemistry methods including PELE calculations, docking, molecular dynamics and experimental surface plasmon resonance, we were able to prove the direct binding of this molecule to ILK, thus providing the basis of its molecular recognition by the protein and the effect over its architecture. Our breakthroughs show that Cpd22 binding stabilizes the ILK domain by binding to the pseudo-active site in a similar way to the ATP, possibly modulating its scaffolding properties as pseudokinase. Moreover, our results explain the experimental observations obtained during Cpd22 development, thus paving the way to the development of new chemical probes and potential drugs.

Graphical abstract

A Computer-Driven Scaffold-Hopping Approach Generating New PTP1B Inhibitors from the Pyrrolo[1,2-a]quinoxaline Core

Protein tyrosine phosphatase 1B (PTP1B) is a very promising target for the treatment of metabolic disorders such as type II diabetes mellitus. Although it was validated as a promising target for this disease more than 30 years ago, as yet there is no drug in advanced clinical trials, and its biochemical mechanism and functions are still being studied. In the present study, based on our experience generating PTP1B inhibitors, we have developed and implemented a scaffold-hopping approach to vary the pyrrole ring of the pyrrolo[1,2-a]quinoxaline core, supported by extensive computational techniques aimed to explain the molecular interaction with PTP1B. Using a combination of docking, molecular dynamics and end-point free-energy calculations, we have rationally designed a hypothesis for new PTP1B inhibitors, supporting their recognition mechanism at a molecular level. After the design phase, we were able to easily synthesize proposed candidates and their evaluation against PTP1B was found to be in good concordance with our predictions. Moreover, the best candidates exhibited glucose uptake increments in cellulo model, thus confirming their utility for PTP1B inhibition and validating this approach for inhibitors design and molecules thus obtained.

Current trends in computer aided drug design and a highlight of drugs discovered via computational techniques: A review

Computer-aided drug design (CADD) is one of the pivotal approaches to contemporary pre-clinical drug discovery, and various computational techniques and software programs are typically used in combination, in a bid to achieve the desired outcome. Several approved drugs have been developed with the aid of CADD. On SciFinder®, we evaluated more than 600 publications through systematic searching and refining, using the terms, virtual screening; software methods; computational studies and publication year, in order to obtain data concerning particular aspects of CADD. The primary focus of this review was on the databases screened, virtual screening and/or molecular docking software program used. Furthermore, we evaluated the studies that subsequently performed molecular dynamics (MD) simulations and we reviewed the software programs applied, the application of density functional theory (DFT) calculations and experimental assays. To represent the latest trends, the most recent data obtained was between 2015 and 2020, consequently the most frequently employed techniques and software programs were recorded. Among these, the ZINC database was the most widely preferred with an average use of 31.2%. Structure-based virtual screening (SBVS) was the most prominently used type of virtual screening and it accounted for an average of 57.6%, with AutoDock being the preferred virtual screening/molecular docking program with 41.8% usage. Following the screening process, 38.5% of the studies performed MD simulations to complement the virtual screening and GROMACS with 39.3% usage, was the popular MD software program. Among the computational techniques, DFT was the least applied whereby it only accounts for 0.02% average use. An average of 36.5% of the studies included reports on experimental evaluations following virtual screening. Ultimately, since the inception and application of CADD in pre-clinical drug discovery, more than 70 approved drugs have been discovered, and this number is steadily increasing over time.

Efficient Dimerization Disruption of Leishmania infantum Trypanothione Reductase by Triazole-phenyl-thiazoles

Inhibition of Leishmania infantum trypanothione disulfide reductase (LiTryR) by disruption of its homodimeric interface has proved to be an alternative and unexploited strategy in the search for novel antileishmanial agents. Proof of concept was first obtained by peptides and peptidomimetics. Building on previously reported dimerization disruptors containing an imidazole-phenyl-thiazole scaffold, we now report a new 1,2,3-triazole-based chemotype that yields noncompetitive, slow-binding inhibitors of LiTryR. Several compounds bearing (poly)aromatic substituents dramatically improve the ability to disrupt LiTryR dimerization relative to reference imidazoles. Molecular modeling studies identified an almost unexplored hydrophobic region at the interfacial domain as the putative binding site for these compounds. A subsequent structure-based design led to a symmetrical triazole analogue that displayed even more potent inhibitory activity over LiTryR and enhanced leishmanicidal activity. Remarkably, several of these novel triazole-bearing compounds were able to kill both extracellular and intracellular parasites in cell cultures.

Binding of non-canonical peptidoglycan controls Vibrio cholerae broad spectrum racemase activity

•

Cell wall muropeptides bind Vibrio cholerae's broad spectrum racemase BsrV and inhibit its activity.

•

BsrV activity is regulated by polymeric peptidoglycan.

•

Muropeptides modified with non-canonical D-amino acids show the strongest inhibitory activity of BsrV.

Broad-spectrum amino acid racemases (Bsrs) enable bacteria to generate non-canonical D-amino acids (NCDAAs), whose roles and impact on microbial physiology, including modulation of cell wall structure and dissolution of biofilms, are just beginning to be appreciated. Here we used a diverse array of structural, biochemical and molecular simulation studies to define and characterize how BsrV is post-translationally regulated. We discovered that contrary to Vibrio cholerae alanine racemase AlrV highly compacted active site, BsrV’s is broader and can be occupied by cell wall stem peptides. We found that peptidoglycan peptides modified with NCDAAs are better stabilized by BsrV’s catalytic cavity and show better inhibitory capacity than canonical muropeptides. Notably, BsrV binding and inhibition can be recapitulated by undigested peptidoglycan sacculi as it exists in the cell. Docking simulations of BsrV binding the peptidoglycan polymer generate a model where the peptide stems are perfectly accommodated and stabilized within each of the dimeŕs active sites. Taking these biochemical and structural data together, we propose that inhibition of BsrV by peptidoglycan peptides underlies a negative regulatory mechanism to avoid excessive NCDAA production. Our results collectively open the door to use “à la carte” synthetic peptides as a tool to modulate DAAs production of Bsr enzymes.

Structural Cues for Understanding eEF1A2 Moonlighting

Spontaneous mutations in the EEF1A2 gene cause epilepsy and severe neurological disabilities in children. The crystal structure of eEF1A2 protein purified from rabbit skeletal muscle reveals a post-translationally modified dimer that provides information about the sites of interaction with numerous binding partners, including itself, and maps these mutations onto the dimer and tetramer interfaces. The spatial locations of the side chain carboxylates of Glu301 and Glu374, to which phosphatidylethanolamine is uniquely attached via an amide bond, define the anchoring points of eEF1A2 to cellular membranes and interorganellar membrane contact sites. Additional bioinformatic and molecular modeling results provide novel structural insight into the demonstrated binding of eEF1A2 to SH3 domains, the common MAPK docking groove, filamentous actin, and phosphatidylinositol-4 kinase IIIβ. In this new light, the role of eEF1A2 as an ancient, multifaceted, and articulated G protein at the crossroads of autophagy, oncogenesis and viral replication appears very distant from the "canonical" one of delivering aminoacyl-tRNAs to the ribosome that has dominated the scene and much of the thinking for many decades.

Spiropyran Meets Guanine Quadruplexes: Isomerization Mechanism and DNA Binding Modes of Quinolizidine-Substituted Spiropyran Probes

The recent delivery of a fluorescent quinolizidine-substituted spiropyran, which is able to switch in vivo and bind to guanine quadruplexes (G4) at physiological pH values, urged us to elucidate its molecular switching and binding mechanism. Combining multiscale dynamical studies and accurate quantum chemical calculations, we show that, both in water and in the G4 environment, the switching of the spiropyran ring is not promoted by an initial protonation event-as expected by the effect of low pH solutions-but that the deprotonated merocyanine form is an intermediate of the reaction leading to the protonated open species. Additionally, we investigate the binding of both deprotonated and protonated open forms of merocyanine to c-MYC G4s. Both species bind to G4s albeit with different hydrogen-bond patterns and provide distinct rotamers around the exocyclic double bond of the merocyanine forms. Altogether, our study sheds light on the pharmacophoric points for the binding of these probes to DNA, and thereby, contributes to future developments of new G4 binders of the remarkable family of quinolizidine-substituted spiropyrans.

Identification of the naphthoquinone derivative inhibitors binding site in heat shock protein 90: an induced-fit docking, molecular dynamics and 3D-QSAR study

The combination of molecular modeling methods to identify the putative binding site of inhibitors constitutes an important tool in drug discovery. In this work, we used these analyses to understand the potent inhibitory effect of naphthoquinone derivatives on heat shock protein 90 (Hsp90), one of the proteins involved in many types of cancer. Molecular docking results indicated that some favorable interactions of key amino acid residues at the binding site of Hsp90 with these quinones would be responsible for the inhibition of Hsp90 activity. Molecular docking and molecular dynamics simulation were carried out to further understand the binding modes and the interactions between the protein and these inhibitors. The main residues of the internal cavity were Val136, Phe138, Tyr139, Val150, Trp162 and Val186. The high concordance between the docking results and 3D-QSAR contour maps gives us helpful information about the environment of the binding site. Our results provide the bases for a rational modification of new molecules based in quinone scaffold, in order to design more potent Hsp90 inhibitors, which would exhibit highly potent antitumor activity.Communicated by Ramaswamy H. Sarma.

Distinct binding of cetirizine enantiomers to human serum albumin and the human histamine receptor H1

Cetirizine, a major metabolite of hydroxyzine, became a marketed second-generation H₁ antihistamine that is orally active and has a rapid onset of action, long duration of effects and a very good safety record at recommended doses. The approved drug is a racemic mixture of (S)-cetirizine and (R)-cetirizine, the latter being the levorotary enantiomer that also exists in the market as a third-generation, non-sedating and highly selective antihistamine. Both enantiomers bind tightly to the human histamine H₁ receptor (hH₁R) and behave as inverse agonists but the affinity and residence time of (R)-cetirizine are greater than those of (S)-cetirizine. In blood plasma, cetirizine exists in the zwitterionic form and more than 90% of the circulating drug is bound to human serum albumin (HSA), which acts as an inactive reservoir. Independent X-ray crystallographic work has solved the structure of the hH₁R:doxepin complex and has identified two drug-binding sites for cetirizine on equine serum albumin (ESA). Given this background, we decided to model a membrane-embedded hH₁R in complex with either (R)- or (S)-cetirizine and also the complexes of both ESA and HSA with these two enantiomeric drugs to analyze possible differences in binding modes between enantiomers and also among targets. The ensuing molecular dynamics simulations in explicit solvent and additional computational chemistry calculations provided structural and energetic information about all of these complexes that is normally beyond current experimental possibilities. Overall, we found very good agreement between our binding energy estimates and extant biochemical and pharmacological evidence. A much higher degree of solvent exposure in the cetirizine-binding site(s) of HSA and ESA relative to the more occluded orthosteric binding site in hH₁R is translated into larger positional fluctuations and considerably lower affinities for these two nonspecific targets. Whereas it is demonstrated that the two known pockets in ESA provide enough stability for cetirizine binding, only one such site does so in HSA due to a number of amino acid replacements. At the histamine-binding site in hH₁R, the distinct interactions established between the phenyl and chlorophenyl moieties of the two enantiomers with the amino acids lining up the pocket and between their free carboxylates and Lys179 in the second extracellular loop account for the improved pharmacological profile of (R)-cetirizine.

Reaction mechanism of nucleoside 2'-deoxyribosyltransferases: free-energy landscape supports an oxocarbenium ion as the reaction intermediate

Insight into the catalytic mechanism of Lactobacillus leichmannii nucleoside 2'-deoxyribosyltransferase (LlNDT) has been gained by calculating a quantum mechanics-molecular mechanics (QM/MM) free-energy landscape of the reaction within the enzyme active site. Our results support an oxocarbenium species as the reaction intermediate and thus an SN1 reaction mechanism in this family of bacterial enzymes. Our mechanistic proposal is validated by comparing experimental kinetic data on the impact of the single amino acid replacements Tyr7, Glu98 and Met125 with Ala, Asp and Ala/norLeu, respectively, and accounts for the specificity shown by this enzyme on a non-natural substrate. This work broadens our understanding of enzymatic C-N bond cleavage and C-N bond formation.

Binding site plasticity in viral PPxY Late domain recognition by the third WW domain of human NEDD4

The recognition of PPxY viral Late domains by the third WW domain of the HECT-E3 ubiquitin ligase NEDD4 (hNEDD4-WW3) is essential for the completion of the budding process of numerous enveloped viruses, including Ebola, Marburg, HTLV1 or Rabies. hNEDD4-WW3 has been validated as a promising target for the development of novel host-oriented broad spectrum antivirals. Nonetheless, finding inhibitors with good properties as therapeutic agents remains a challenge since the key determinants of binding affinity and specificity are still poorly understood. We present here a detailed structural and thermodynamic study of the interactions of hNEDD4-WW3 with viral Late domains combining isothermal titration calorimetry, NMR structural determination and molecular dynamics simulations. Structural and energetic differences in Late domain recognition reveal a highly plastic hNEDD4-WW3 binding site that can accommodate PPxY-containing ligands with varying orientations. These orientations are mostly determined by specific conformations adopted by residues I859 and T866. Our results suggest a conformational selection mechanism, extensive to other WW domains, and highlight the functional relevance of hNEDD4-WW3 domain conformational flexibility at the binding interface, which emerges as a key element to consider in the search for potent and selective inhibitors of therapeutic interest.

Pyrrolopyrimidine vs Imidazole-Phenyl-Thiazole Scaffolds in Nonpeptidic Dimerization Inhibitors of Leishmania infantum Trypanothione Reductase

Disruption of protein-protein interactions of essential oligomeric enzymes by small molecules represents a significant challenge. We recently reported some linear and cyclic peptides derived from an α-helical region present in the homodimeric interface of Leishmania infantum trypanothione reductase ( Li-TryR) that showed potent effects on both dimerization and redox activity of this essential enzyme. Here, we describe our first steps toward the design of nonpeptidic small-molecule Li-TryR dimerization disruptors using a proteomimetic approach. The pyrrolopyrimidine and the 5-6-5 imidazole-phenyl-thiazole α-helix-mimetic scaffolds were suitably decorated with substituents that could mimic three key residues (K, Q, and I) of the linear peptide prototype (PKIIQSVGIS-Nle-K-Nle). Extensive optimization of previously described synthetic methodologies was required. A library of 15 compounds bearing different hydrophobic alkyl and aromatic substituents was synthesized. The imidazole-phenyl-thiazole-based analogues outperformed the pyrrolopyrimidine-based derivatives in both inhibiting the enzyme and killing extracellular and intracellular parasites in cell culture. The most active imidazole-phenyl-thiazole compounds 3e and 3f inhibit Li-TryR and prevent growth of the parasites at low micromolar concentrations similar to those required by the peptide prototype. The intrinsic fluorescence of these compounds inside the parasites visually demonstrates their good permeability in comparison with previous peptide-based Li-TryR dimerization disruptors.

Viral engagement with host receptors blocked by a novel class of tryptophan dendrimers that targets the 5-fold-axis of the enterovirus-A71 capsid

Enterovirus A71 (EV-A71) is a non-polio neurotropic enterovirus with pandemic potential. There are no antiviral agents approved to prevent or treat EV-A71 infections. We here report on the molecular mechanism by which a novel class of tryptophan dendrimers inhibits (at low nanomolar to high picomolar concentration) EV-A71 replication in vitro. A lead compound in the series (MADAL385) prevents binding and internalization of the virus but does not, unlike classical capsid binders, stabilize the particle. By means of resistance selection, reverse genetics and cryo-EM, we map the binding region of MADAL385 to the 5-fold vertex of the viral capsid and demonstrate that a single molecule binds to each vertex. By interacting with this region, MADAL385 prevents the interaction of the virus with its cellular receptors PSGL1 and heparan sulfate, thereby blocking the attachment of EV-A71 to the host cells.

Enterovirus A71 (EV-A71) is the virus responsible for most of the severe forms of hand, foot and mouth disease (HFMD) associated with neurological involvement and mortality in young children under the age of 5. Seasonal outbreaks of HFMD -with a 2–3 years epidemic cycle- are recurring around the world, especially in the Asia-Pacific region. To date, no antiviral agent has been approved for the treatment of EV-A71 infections. Here, we report on a recently uncovered class of tryptophan dendrimers with an extraordinary antiviral activity in vitro against circulating EV-A71 clinical isolates. Mode of action studies revealed that this class of compounds targets the 5-fold vertex of EV-A71, in turn blocking receptor binding. Our finding may open an entirely novel line of research and largely aid in anti-enterovirus drug development.

Active learning strategies with COMBINE analysis: new tricks for an old dog

The COMBINE method was designed to study congeneric series of compounds including structural information of ligand-protein complexes. Although very successful, the method has not received the same level of attention than other alternatives to study Quantitative Structure Active Relationships (QSAR) mainly because lack of ways to measure the uncertainty of the predictions and the need for large datasets. Active learning, a semi-supervised learning approach that makes use of uncertainty to enhance models' performance while reducing the size of the training sets, has been used in this work to address both problems. We propose two estimators of uncertainty: the pool of regressors and the distance to the training set. The performance of the methods has been evaluated by testing the resulting active learning workflows in 3 diverse datasets: HIV-1 protease inhibitors, Taxol-derivatives and BRD4 inhibitors. The proposed strategies were successful in 80% of the cases for the taxol-derivatives and BRD4 inhibitors, while outperformed random selection in the case of the HIV-1 protease inhibitors time-split. Our results suggest that AL-COMBINE might be an effective way of producing consistently superior QSAR models with a limited number of samples.

Affinity Calculations of Cyclodextrin Host-Guest Complexes: Assessment of Strengths and Weaknesses of End-Point Free Energy Methods

The end-point methods like MM/PBSA or MM/GBSA estimate the free energy of a biomolecule by combining its molecular mechanics energy with solvation free energy and entropy terms. On the one hand, their performance largely depends on the particular system of interest, and despite numerous attempts to improve their reliability that have resulted in many variants, there is still no clear alternative to improve their accuracy. On the other hand, the relatively small cyclodextrin host-guest complexes, for which high-quality binding calorimetric data are usually available, are becoming reference models for testing the accuracy of free energy methods. In this work, we further assess the performance of various MM/PBSA-like approaches as applied to cyclodextrin complexes. To this end, we select a set of complexes between β-cyclodextrin and 57 small organic molecules that has been previously studied with the binding energy distribution analysis method in combination with an implicit solvent model ( Wickstrom, L.; He, P.; Gallicchio, E.; Levy, R. M. J. Chem. Theory Comput. 2013 , 9 , 3136 - 3150 ). For each complex, a conventional 1.0 μs molecular dynamics simulation in explicit solvent is performed. Then we employ semiempirical quantum chemical calculations, several variants of the MM-PB(GB)SA methods, entropy estimations, etc., to assess the reliability of the end-point affinity calculations. The best end-point protocol in this study, which combines DFTB3 energies with entropy corrections, yields estimations of the binding free energies that still have substantial errors (RMSE = 2.2 kcal/mol), but it exhibits a good prediction capacity in terms of ligand ranking ( R² = 0.66) that is close to or even better than that of rigorous free energy methodologies. Our results can be helpful to discriminate between the intrinsic limitations of the end-point methods and other sources of error, such as the underlying energy and continuum solvation methods.

In silico-designed mutations increase variable new-antigen receptor single-domain antibodies for VEGF165 neutralization

The stability, binding, and tissue penetration of variable new-antigen receptor (VNAR) single-domain antibodies have been tested as part of an investigation into their ability to serve as novel therapeutics. V13 is a VNAR that recognizes vascular endothelial growth factor 165 (VEGF₁₆₅). In the present study V13 was used as a parental molecule into which we introduced mutations designed in silico. Two of the designed VNAR mutants were expressed, and their ability to recognize VEGF₁₆₅ was assessed in vitro and in vivo. One mutation (Pro98Tyr) was designed to increase VEGF₁₆₅ recognition, while the other (Arg97Ala) was designed to inhibit VEGF₁₆₅ binding. Compared to parental V13, the Pro98Tyr mutant showed enhanced VEGF₁₆₅ recognition and neutralization, as indicated by inhibition of angiogenesis and tumor growth. This molecule thus appears to have therapeutic potential for neutralizing VEGF₁₆₅ in cancer treatment.

Investigation of the complexation between cyclodextrins and medetomidine enantiomers by capillary electrophoresis, NMR spectroscopy and molecular modeling

The migration order of the enantiomers of medetomidine in the presence of cyclodextrins studied by capillary electrophoresis in phosphate buffer, pH 2.5, depended on the cavity size and the substitution pattern of the cyclodextrins. Opposite migration order was observed in the presence of β-cyclodextrin (β-CD) and γ-cyclodextrin (γ-CD) as well as randomly sulfated β-CD (S-β-CD) and heptakis(6-O-sulfo)-β-CD (HS-β-CD). This could be rationalized by the fact that dexmedetomidine formed more stable complexes with β-CD and S-β-CD, while levomedetomidine interacted stronger with γ-CD and HS-β-CD. The structure of the complexes was derived from rotating frame nuclear Overhauser (ROESY) experiments for β-CD, γ-CD and HS-β-CD. In the case of the native CDs, the phenyl ring of medetomidine entered the cavity through the wider secondary rim of the CDs, whereas the protonated imidazole ring was positioned inside the CD cavity interacting with the sulfate groups of HS-β-CD. Furthermore, molecular dynamics calculations also suggested opposite affinities of the medetomidine enantiomers toward β-CD and γ-CD.

A Series of Enthalpically Optimized Docetaxel Analogues Exhibiting Enhanced Antitumor Activity and Water Solubility

A dual-purpose strategy aimed at enhancing the binding affinity for microtubules and improving the water solubility of docetaxel led to the design and synthesis of a series of C-2- and C-3'-modified analogues. Both aims were realized when the C-3' phenyl group present in docetaxel was replaced with a propargyl alcohol. The resulting compound, 3f, was able to overcome drug resistance in cultured P-gp-overexpressing tumor cells and showed greater activity than docetaxel against drug-resistant A2780/AD ovarian cancer xenografts in mice. In addition, the considerably lower hydrophobicity of 3f relative to both docetaxel and paclitaxel led to better aqueous solubility. A molecular model of tubulin-bound 3f revealed novel hydrogen-bonding interactions between the propargyl alcohol and the polar environment provided by the side chains of Ser236, Glu27, and Arg320.

High-affinity ligands of the colchicine domain in tubulin based on a structure-guided design

Microtubule-targeting agents that bind at the colchicine-site of tubulin are of particular interest in antitumoral therapy due to their dual mechanism of action as antimitotics and vascular disrupting agents. Cyclohexanediones derivatives have been described as a new family of colchicine-domain binders with an association constant to tubulin similar to that of colchicine. Here, the high-resolution structures of tubulin in complex with cyclohexanediones TUB015 and TUB075 were solved by X-ray crystallography. A detailed analysis of the tubulin-TUB075 interaction by means of computational affinity maps allowed the identification of two additional regions at the binding site that were addressed with the design and synthesis of a new series of cyclohexanediones with a distal 2-substituted benzofurane. These new compounds showed potent antiproliferative activity with IC50 values in the nM range, arrested cell cycle progression at the G2/M phase and induced apoptosis at sub μM concentrations. Moreover, they caused the destruction of a preformed vascular network in vitro and inhibited the migration of endothelial cells at non-toxic concentrations. Finally, these compounds displayed high affinity for tubulin as substantiated by a K b value of 2.87 × 108 M-1 which, to the best of our knowledge, represents the highest binding constant measured to date for a colchicine-domain ligand.

Conformational mimetics of the α-methyl chalcone TUB091 binding tubulin: Design, synthesis and antiproliferative activity

Based on the conformation of the α-methyl chalcone TUB091 in its complex with tubulin, a series of conformational mimetics have been designed and synthesized where the methyl group of the chalcone has been fused to phenyl ring B resulting in 1,2,3,4-tetrahydronaphthalen-2-yl aryl ketones. Among the synthesized compounds, the 5-amino-6-methoxy derivative, with a similar substitution pattern to that of TUB091, showed antiproliferative activity around 20 nM against tumor and endothelial cells. Tubulin binding experiments confirmed its binding to tubulin at the colchicine site with a Kb of 2.4 × 10⁶ M^-1 resulting in the inhibition of the in vitro assembly of purified tubulin. Moreover, based on the recently reported complex of combretastatin A4 (CA4) with tubulin, a comparative analysis of the binding mode of CA4 and the α-methyl chalcone to tubulin has been performed.

Structural rationale for the chiral separation and migration order reversal of clenpenterol enantiomers in capillary electrophoresis using two different β-cyclodextrins

NMR spectroscopy experiments, molecular dynamics simulations, and theoretical chemistry calculations provide insight into the structural and energetic determinants of the distinct binding of clenpenterol enantiomers to two cyclodextrins and the migration order reversal of their respective inclusion complexes in capillary electrophoresis.

Making sense of the past: hyperstability of ancestral thioredoxins explained by free energy simulations

Thioredoxin (Trx), a small and globular protein, is present in all kinds of organisms, from Archea to higher mammals. Throughout evolution, the Trx sequence has undergone subtle modifications to adapt to varying environmental conditions. The high degree of sequence conservation makes Trx very amenable to ancestral protein reconstruction techniques. In this work, we address the study of the structural and energetic determinants of thermostability in E. coli Trx using a dataset of mutations inspired by ancestral reconstruction. We compute, from first principles, the expected contribution of 19 different amino acid substitutions to the stability (ΔΔG) and the melting temperature (ΔT_m) of the protein. We also describe the specific changes in structure and protein dynamics responsible for the stabilizing or destabilizing effects of these mutations. Our results point to local and independent changes for most of the variants. Our predictions are accurate enough to substantiate the proposal of new hypotheses regarding evolutionary relationships between mutations, as in the case of T89R, P68A and G74S or K90L and F102A, and reach beyond the initial set to suggest improved variants, such as K90I or K90Y.

Development of inhibitors of receptor protein tyrosine phosphatase β/ζ (PTPRZ1) as candidates for CNS disorders

A new series of blood-brain barrier permeable molecules designed to mimic the activity of Pleiotrophin in the CNS has been designed and synthesized. These compounds exert their action by interacting with the intracellular domain PD1 of the Protein Tyrosine-Phosphatase Receptor Z1 (PTPRZ1), and inhibiting its tyrosine phosphatase activity. The most potent compounds 10a and 12b (IC₅₀ = 0,1 μM) significantly increase the phosphorylation of key tyrosine residues of PTPRZ1 substrates involved in neuronal survival and differentiation, and display protective effects against amphetamine-induced toxicity. Docking and molecular dynamics experiments have been used to analyze the binding mode and to explain the observed selectivity against PTP1B. An In vivo experiment has demonstrated that 10a can cross the BBB, thus promoting the possibility of moving forward these candidates for the development of drugs for the treatment of CNS disorders, such as drug addiction and neurodegenerative diseases.

Combined molecular modelling and 3D-QSAR study for understanding the inhibition of NQO1 by heterocyclic quinone derivatives

A combination of three-dimensional quantitative structure-activity relationship (3D-QSAR), and molecular modelling methods were used to understand the potent inhibitory NAD(P)H:quinone oxidoreductase 1 (NQO1) activity of a set of 52 heterocyclic quinones. Molecular docking results indicated that some favourable interactions of key amino acid residues at the binding site of NQO1 with these quinones would be responsible for an improvement of the NQO1 activity of these compounds. The main interactions involved are hydrogen bond of the amino group of residue Tyr128, π-stacking interactions with Phe106 and Phe178, and electrostatic interactions with flavin adenine dinucleotide (FADH) cofactor. Three models were prepared by 3D-QSAR analysis. The models derived from Model I and Model III, shown leave-one-out cross-validation correlation coefficients (q²_LOO ) of .75 and .73 as well as conventional correlation coefficients (R² ) of .93 and .95, respectively. In addition, the external predictive abilities of these models were evaluated using a test set, producing the predicted correlation coefficients (r²_pred ) of .76 and .74, respectively. The good concordance between the docking results and 3D-QSAR contour maps provides helpful information about a rational modification of new molecules based in quinone scaffold, in order to design more potent NQO1 inhibitors, which would exhibit highly potent antitumor activity.

2'-Deoxyribosyltransferase from Leishmania mexicana, an efficient biocatalyst for one-pot, one-step synthesis of nucleosides from poorly soluble purine bases

Processes catalyzed by enzymes offer numerous advantages over chemical methods although in many occasions the stability of the biocatalysts becomes a serious concern. Traditionally, synthesis of nucleosides using poorly water-soluble purine bases, such as guanine, xanthine, or hypoxanthine, requires alkaline pH and/or high temperatures in order to solubilize the substrate. In this work, we demonstrate that the 2'-deoxyribosyltransferase from Leishmania mexicana (LmPDT) exhibits an unusually high activity and stability under alkaline conditions (pH 8-10) across a broad range of temperatures (30-70 °C) and ionic strengths (0-500 mM NaCl). Conversely, analysis of the crystal structure of LmPDT together with comparisons with hexameric, bacterial homologues revealed the importance of the relationships between the oligomeric state and the active site architecture within this family of enzymes. Moreover, molecular dynamics and docking approaches provided structural insights into the substrate-binding mode. Biochemical characterization of LmPDT identifies the enzyme as a type I NDT (PDT), exhibiting excellent activity, with specific activity values 100- and 4000-fold higher than the ones reported for other PDTs. Interestingly, LmPDT remained stable during 36 h at different pH values at 40 °C. In order to explore the potential of LmPDT as an industrial biocatalyst, enzymatic production of several natural and non-natural therapeutic nucleosides, such as vidarabine (ara A), didanosine (ddI), ddG, or 2'-fluoro-2'-deoxyguanosine, was carried out using poorly water-soluble purines. Noteworthy, this is the first time that the enzymatic synthesis of 2'-fluoro-2'-deoxyguanosine, ara G, and ara H by a 2'-deoxyribosyltransferase is reported.

First example of peptides targeting the dimer interface of Leishmania infantum trypanothione reductase with potent in vitro antileishmanial activity. Eur J Med Chem 2017;135:49-59. [PMID: 28431354 DOI: 10.1016/j.ejmech.2017.04.020]

A series of 9-mer and 13-mer amide-bridged cyclic peptides derived from the linear prototype Ac-PKIIQSVGIS-Nle-K-Nle-NH₂ (Toro et al. ChemBioChem2013) has been designed and synthesized by introduction of the lactam between amino acid side chains that are separated by one helical turn (i, i+4). All of these compounds were tested in vitro as both dimerization and enzyme inhibitors of Leishmania infantum trypanothione reductase (Li-TryR). Three of the 13-mer cyclic peptide derivatives (3, 4 and 6) inhibited the oxidoreductase activity of Li-TryR in the low micromolar range and they also disrupted enzyme dimerization. Cyclic analogues 3 and 4 were more resistant to proteases than was the linear prototype. Furthermore, the most potent TryR inhibitors in the linear and cyclic series displayed potent in vitro activity against Leishmania infantum upon conjugation with cationic cell-penetrating peptides. To date, these conjugated peptides can be considered the first example of TryR dimerization inhibitors that are active in cell culture.

Structural Determinants of the Dictyostatin Chemotype for Tubulin Binding Affinity and Antitumor Activity Against Taxane- and Epothilone-Resistant Cancer Cells

A combined biochemical, structural, and cell biology characterization of dictyostatin is described, which enables an improved understanding of the structural determinants responsible for the high-affinity binding of this anticancer agent to the taxane site in microtubules (MTs). The study reveals that this macrolide is highly optimized for MT binding and that only a few of the structural modifications featured in a library of synthetic analogues resulted in small gains in binding affinity. The high efficiency of the dictyostatin chemotype in overcoming various kinds of clinically relevant resistance mechanisms highlights its potential for therapeutic development for the treatment of drug-resistant tumors. A structural explanation is advanced to account for the synergy observed between dictyostatin and taxanes on the basis of their differential effects on the MT lattice. The X-ray crystal structure of a tubulin-dictyostatin complex and additional molecular modeling have allowed the rationalization of the structure-activity relationships for a set of synthetic dictyostatin analogues, including the highly active hybrid 12 with discodermolide. Altogether, the work reported here is anticipated to facilitate the improved design and synthesis of more efficacious dictyostatin analogues and hybrids with other MT-stabilizing agents.

Computational Approaches to Toll-Like Receptor 4 Modulation

Toll-like receptor 4 (TLR4), along with its accessory protein myeloid differentiation factor 2 (MD-2), builds a heterodimeric complex that specifically recognizes lipopolysaccharides (LPS), which are present on the cell wall of Gram-negative bacteria, activating the innate immune response. Some TLR4 modulators are undergoing preclinical and clinical evaluation for the treatment of sepsis, inflammatory diseases, cancer and rheumatoid arthritis. Since the relatively recent elucidation of the X-ray crystallographic structure of the extracellular domain of TLR4, research around this fascinating receptor has risen to a new level, and thus, new perspectives have been opened. In particular, diverse computational techniques have been applied to decipher some of the basis at the atomic level regarding the mechanism of functioning and the ligand recognition processes involving the TLR4/MD-2 system at the atomic level. This review summarizes the reported molecular modeling and computational studies that have recently provided insights into the mechanism regulating the activation/inactivation of the TLR4/MD-2 system receptor and the key interactions modulating the molecular recognition process by agonist and antagonist ligands. These studies have contributed to the design and the discovery of novel small molecules with promising activity as TLR4 modulators.

Activation of Human Toll-like Receptor 4 (TLR4)·Myeloid Differentiation Factor 2 (MD-2) by Hypoacylated Lipopolysaccharide from a Clinical Isolate of Burkholderia cenocepacia

Lung infection by Burkholderia species, in particular Burkholderia cenocepacia, accelerates tissue damage and increases post-lung transplant mortality in cystic fibrosis patients. Host-microbe interplay largely depends on interactions between pathogen-specific molecules and innate immune receptors such as Toll-like receptor 4 (TLR4), which recognizes the lipid A moiety of the bacterial lipopolysaccharide (LPS). The human TLR4·myeloid differentiation factor 2 (MD-2) LPS receptor complex is strongly activated by hexa-acylated lipid A and poorly activated by underacylated lipid A. Here, we report that B. cenocepacia LPS strongly activates human TLR4·MD-2 despite its lipid A having only five acyl chains. Furthermore, we show that aminoarabinose residues in lipid A contribute to TLR4-lipid A interactions, and experiments in a mouse model of LPS-induced endotoxic shock confirmed the proinflammatory potential of B. cenocepacia penta-acylated lipid A. Molecular modeling combined with mutagenesis of TLR4-MD-2 interactive surfaces suggests that longer acyl chains and the aminoarabinose residues in the B. cenocepacia lipid A allow exposure of the fifth acyl chain on the surface of MD-2 enabling interactions with TLR4 and its dimerization. Our results provide a molecular model for activation of the human TLR4·MD-2 complex by penta-acylated lipid A explaining the ability of hypoacylated B. cenocepacia LPS to promote proinflammatory responses associated with the severe pathogenicity of this opportunistic bacterium.

Development of urea and thiourea kynurenamine derivatives: synthesis, molecular modeling, and biological evaluation as nitric oxide synthase inhibitors

Herein we describe the synthesis of a new family of kynurenamine derivatives with a urea or thiourea moiety, together with their in vitro biological evaluation as inhibitors of both neuronal and inducible nitric oxide synthases (nNOS and iNOS, respectively), enzymes responsible for the biogenesis of NO. These compounds were synthesized from a 5-substituted-2-nitrophenyl vinyl ketone scaffold in a five-step procedure with moderate to high chemical yields. In general, the assayed compounds show greater inhibition of iNOS than of nNOS, with 1-[3-(2-amino-5-chlorophenyl)-3-oxopropyl]-3-ethylurea (compound 5 n) being the most potent iNOS inhibitor in the series and the most iNOS/nNOS-selective compound. In this regard, we performed molecular modeling studies to propose a binding mode for this family of compounds to both enzymes and, thereby, to elucidate the differential molecular features that could explain the observed selectivity between iNOS and nNOS.