How 'Protein-Docking' Translates into the New Emerging Field of Docking Small Molecules to Nucleic Acids?

In this review, we retraced the '40-year evolution' of molecular docking algorithms. Over the course of the years, their development allowed to progress from the so-called 'rigid-docking' searching methods to the more sophisticated 'semi-flexible' and 'flexible docking' algorithms. Together with the advancement of computing architecture and power, molecular docking's applications also exponentially increased, from a single-ligand binding calculation to large screening and polypharmacology profiles. Recently targeting nucleic acids with small molecules has emerged as a valuable therapeutic strategy especially for cancer treatment, along with bacterial and viral infections. For example, therapeutic intervention at the mRNA level allows to overcome the problematic of undruggable proteins without modifying the genome. Despite the promising therapeutic potential of nucleic acids, molecular docking programs have been optimized mostly for proteins. Here, we have analyzed literature data on nucleic acid to benchmark some of the widely used docking programs. Finally, the comparison between proteins and nucleic acid targets docking highlighted similarity and differences, which are intrinsically related to their chemical and structural nature.

Cooperative binding of the tandem WW domains of PLEKHA7 to PDZD11 promotes conformation-dependent interaction with tetraspanin 33.

Pleckstrin homology domain-containing A7 (PLEKHA7) is a cytoplasmic protein at adherens junctions that has been implicated in hypertension, glaucoma, and responses to Staphylococcus aureus α-toxin. Complex formation between PLEKHA7, PDZ domain-containing 11 (PDZD11), tetraspanin 33, and the α-toxin receptor ADAM metallopeptidase domain 10 (ADAM10) promotes junctional clustering of ADAM10 and α-toxin-mediated pore formation. However, how the N-terminal region of PDZD11 interacts with the N-terminal tandem WW domains of PLEKHA7 and how this interaction promotes tetraspanin 33 binding to the WW1 domain is unclear. Here, we used site-directed mutagenesis, glutathione S-transferase pulldown experiments, immunofluorescence, molecular modeling, and docking experiments to characterize the mechanisms driving these interactions. We found that Asp-30 of WW1 and His-75 of WW2 interact through a hydrogen bond and, together with Thr-35 of WW1, form a binding pocket that accommodates a polyproline stretch within the N-terminal PDZD11 region. By strengthening the interactions of the ternary complex, the WW2 domain stabilized the WW1 domain and cooperatively promoted the interaction with PDZD11. Modeling results indicated that, in turn, PDZD11 binding induces a conformational rearrangement, which strengthens the ternary complex, and contributes to enlarging a "hydrophobic hot spot" region on the WW1 domain. The last two lipophilic residues of tetraspanin 33, Trp-283 and Tyr-282, were required for its interaction with PLEKHA7. Docking of the tetraspanin 33 C terminus revealed that it fits into the hydrophobic hot spot region of the accessible surface of WW1. We conclude that communication between the two tandem WW domains of PLEKHA7 and the PLEKHA7-PDZD11 interaction modulate the ligand-binding properties of PLEKHA7.

Structure-based drug design and in vitro testing reveal new inhibitors of enoyl-acyl carrier protein reductases.

The need for new antibacterial agents is increasingly becoming of great importance as bacterial resistance to current drugs is quickly spreading. Enoyl-acyl carrier protein reductases (FabI) are important enzymes for fatty acid biosynthesis in bacteria and other micro-organisms. In this project, we conducted structure-based virtual screening against the FabI enzyme, and accordingly, 37 compounds were selected for experimental testing. Interestingly, five compounds were able to demonstrate antimicrobial effect with variable inhibition activity against various strains of bacteria and fungi. Minimum inhibitory concentrations of the active compounds were determined and showed to be in low to medium micromolar range. Subsequently, enzyme inhibition assay was carried out for our five antimicrobial hits to confirm their biological target and determine their IC50 values. Three of these tested compounds exhibited inhibition activity for the FabI enzyme where our best hit MN02 had an IC50 value of 7.8 µM. Furthermore, MN02 is a small bisphenolic compound that is predicted to have all required features to firmly bind with the target enzyme. To sum up, hits discovered in this work can act as a good starting point for the future development of new and potent antimicrobial agents.

Targeting RNA structure in SMN2 reverses spinal muscular atrophy molecular phenotypes.

Modification of SMN2 exon 7 (E7) splicing is a validated therapeutic strategy against spinal muscular atrophy (SMA). However, a target-based approach to identify small-molecule E7 splicing modifiers has not been attempted, which could reveal novel therapies with improved mechanistic insight. Here, we chose as a target the stem-loop RNA structure TSL2, which overlaps with the 5' splicing site of E7. A small-molecule TSL2-binding compound, homocarbonyltopsentin (PK4C9), was identified that increases E7 splicing to therapeutic levels and rescues downstream molecular alterations in SMA cells. High-resolution NMR combined with molecular modelling revealed that PK4C9 binds to pentaloop conformations of TSL2 and promotes a shift to triloop conformations that display enhanced E7 splicing. Collectively, our study validates TSL2 as a target for small-molecule drug discovery in SMA, identifies a novel mechanism of action for an E7 splicing modifier, and sets a precedent for other splicing-mediated diseases where RNA structure could be similarly targeted.

Modeling and resistant alleles explain the selectivity of antimalarial compound 49c towards apicomplexan aspartyl proteases.

Toxoplasma gondii aspartyl protease 3 (TgASP3) phylogenetically clusters with Plasmodium falciparum Plasmepsins IX and X (PfPMIX, PfPMX). These proteases are essential for parasite survival, acting as key maturases for secreted proteins implicated in invasion and egress. A potent antimalarial peptidomimetic inhibitor (49c) originally developed against Plasmepsin II selectively targets TgASP3, PfPMIX, and PfPMX To unravel the molecular basis for the selectivity of 49c, we constructed homology models of PfPMIX, PfPMX, and TgASP3 that were first validated by identifying the determinants of microneme and rhoptry substrate recognition. The flap and flap-like structures of several reported Plasmepsins are highly flexible and critically modulate the access to the binding cavity. Molecular docking of 49c to TgASP3, PfPMIX, and PfPMX models predicted that the conserved phenylalanine residues in the flap, F344, F291, and F305, respectively, account for the sensitivity toward 49c. Concordantly, phenylalanine mutations in the flap of the three proteases increase twofold to 15-fold the IC50 values of 49c. Compellingly the selection of mutagenized T. gondii resistant strains to 49c reproducibly converted F344 to a cysteine residue.