Comprehensive analysis of intramolecular G-quadruplex structures: furthering the understanding of their formalism.

G-quadruplexes (G4) are helical structures found in guanine-rich DNA or RNA sequences. Generally, their formalism is based on a few dozen structures, which can produce some inconsistencies or incompleteness. Using the website ASC-G4, we analyzed the structures of 333 intramolecular G4s, of all types, which allowed us to clarify some key concepts and present new information. To each of the eight distinguishable topologies corresponds a groove-width signature and a predominant glycosidic configuration (gc) pattern governed by the directions of the strands. The relative orientations of the stacking guanines within the strands, which we quantified and related to their vertical gc successions, determine the twist and tilt of the helices. The latter impact the minimum groove widths, which represent the space available for lateral ligand binding. The G4 four helices have similar twists, even when these twists are irregular, meaning that they have various angles along the strands. Despite its importance, the vertical gc succession has no strict one-to-one relationship with the topology, which explains the discrepancy between some topologies and their corresponding circular dichroism spectra. This study allowed us to introduce the new concept of platypus G4s, which are structures with properties corresponding to several topologies.

ASC-G4, an algorithm to calculate advanced structural characteristics of G-quadruplexes.

ASC-G4 is an algorithm for the calculation of the advanced structural characteristics of G-quadruplexes (G4). It allows the unambiguous determination of the intramolecular G4 topology, based on the oriented strand numbering. It also resolves the ambiguity in the determination of the guanine glycosidic configuration. With this algorithm, we showed that the use of the C3' or C5' atoms to calculate the groove width in G4 is more appropriate than the P atoms and that the groove width does not always reflect the space available within the groove. For the latter, the minimum groove width is more appropriate. The application of ASC-G4 to 207 G4 structures guided the choices made for the calculations. A website based on ASC-G4 (http://tiny.cc/ASC-G4) was created, where the user uploads his G4 structure and gets its topology, the types of its loops and their lengths, the presence of snapbacks and bulges, the distribution of guanines in the tetrads and strands, the glycosidic configuration of these guanines, their rise, the groove widths, the minimum groove widths, the tilt and twist angles, the backbone dihedral angles, etc. It also provides a large number of atom-atom and atom-plane distances that are relevant to evaluating the quality of the structure.

Optimization of G-Quadruplex Ligands through a SAR Study Combining Parallel Synthesis and Screening of Cationic Bis(acylhydrazones).

G-quadruplexes (G4s), secondary structures adopted by guanine-rich DNA and RNA sequences, are implicated in numerous biological processes and have been suggested as potential drug targets. Accordingly, there is an increasing interest in developing high-throughput methods that allow the generation of congeneric series of G4-targeting molecules ("ligands") and investigating their interactions with the targets. We have developed an operationally simple method of parallel synthesis to generate "ready-to-screen" libraries of cationic acylhydrazones, a motif that we have previously identified as a promising scaffold for potent, biologically active G4 ligands. Combined with well-established screening techniques, such as fluorescence melting, this method enables the rapid synthesis and screening of combinatorial libraries of potential G4 ligands. Following this protocol, we synthesized a combinatorial library of 90 bis(acylhydrazones) and screened it against five different nucleic acid structures. This way, we were able to analyze the structure-activity relationships within this series of G4 ligands, and identified three novel promising ligands whose interactions with G4-DNAs of different topologies were studied in detail by a combination of several biophysical techniques, including native mass spectrometry, and molecular modeling.

Corrigendum: Deciphering the mechanism of inhibition of SERCA1a by sarcolipin using molecular simulations.

[This corrects the article DOI: 10.3389/fmolb.2020.606254.].

New potent dual inhibitors of CK2 and Pim kinases: discovery and structural insights.

Protein kinase casein kinase 2 (CK2) is a serine/threonine kinase with evidence of implication in growth dysregulation and apoptosis resistance, making it a relevant target for cancer therapy. Several CK2 inhibitors have been developed showing variable efficiency, emphasizing the need to expand the chemical diversity of those inhibitors. We report the identification and characterization of 2,8-difurandicarboxylic acid derivatives as a new class of nanomolar ATP-competitive inhibitors. Selectivity profiling pointed out proviral insertion Moloney virus kinases (Pim kinases) as the only other kinases that are significantly inhibited. By combining structure-activity relationship analysis with structural determination, we were able to determine the binding mode of these inhibitors for both kinases and to explain their strong inhibitory potency. Essential chemical features necessary for activity on both kinases were then identified. The described compounds are not cell permeable: however, they could provide a lead for developing novel inhibitors usable also in vivo. Given the similar but not redundant pathophysiological functions of CK2 and Pim family members, such inhibitors would provide new attractive leads for targeted cancer therapy. This work highlights that 2 functionally related kinases from different kinome branches display exquisite sensitivity to a common inhibitor.

How to choose relevant multiple receptor conformations for virtual screening: a test case of Cdk2 and normal mode analysis.

Better treatment of protein flexibility is essential in structure-based drug design projects such as virtual screening and protein-ligand docking. Diversity in ligand-binding mechanisms and receptor conformational changes makes it difficult to treat dynamic features of the receptor during the docking simulation. Thus, the use of pregenerated multiple receptor conformations is applied today in virtual screening studies. However, generation of a small relevant set of receptor conformations remains challenging. To address this problem, we propose a new protocol for the generation of multiple receptor conformations via normal mode analysis and for the selection of several receptor conformations suitable for docking/virtual screening. We validated this protocol on cyclin-dependent kinase 2, which possesses a binding site located at the interface between two subdomains and is known to undergo significant conformational changes in the active site region upon ligand binding. We believe that the suggested rules for the choice of suitable receptor conformations can be applied to other targets when dealing with in silico screening on flexible receptors.

Deciphering the Mechanism of Inhibition of SERCA1a by Sarcolipin Using Molecular Simulations.

SERCA1a is an ATPase calcium pump that transports Ca2+ from the cytoplasm to the sarco/endoplasmic reticulum lumen. Sarcolipin (SLN), a transmembrane peptide, regulates the activity of SERCA1a by decreasing its Ca2+ transport rate, but its mechanism of action is still not well-understood. To decipher this mechanism, we have performed normal mode analysis in the all-atom model, with the SERCA1a-SLN complex, or the isolated SERCA1a, embedded in an explicit membrane. The comparison of the results allowed us to provide an explanation at the atomic level for the action of SLN that is in good agreement with experimental observations. In our analyses, the presence of SLN locally perturbs the TM6 transmembrane helix and as a consequence modifies the position of D800, one of the key metal-chelating residues. Additionally, it reduces the flexibility of the gating residues, V304, and E309 in TM4, at the entrance of the Ca2+ binding sites, which would decrease the affinity for Ca2+. Unexpectedly, SLN has also an effect on the ATP binding site more than 35 Å away, due to the straightening of TM5, a long helix considered as the spine of the protein. The straightening of TM5 modifies the structure of the P-N linker that sits above it, and which comprises the 351DKTG354 conserved motif, resulting in an increase of the distance between ATP and the phosphorylation site. As a consequence, the turn-over rate could be affected. All this gives SERCA1a the propensity to go toward a Ca2+ low-affinity E2-like state in the presence of SLN and toward a Ca2+ high-affinity E1-like state in the absence of SLN. In addition to a general mechanism of inhibition of SERCA1a regulatory peptides, this study also provides an insight into the conformational transition between the E2 and E1 states.

Salicylaldehyde derivatives as new protein kinase CK2 inhibitors.

Protein kinase CK2 is a Ser/Thr kinase, with a constitutive activity, that is considered as a promising target for cancer therapy. The currently available CK2 inhibitors lack the potency and the pharmacological properties necessary to be suitable and successful in clinical settings. We report the development of new potent CK2 inhibitors from salicylaldehyde derivatives identified by automated screening of a proprietary small-molecule library. Docking simulations and analysis of the structure-activity relationship for the hits allowed to determine their binding modes on CK2, and to carry out the optimization of their structures. This strategy led to the discovery of potent CK2 inhibitors with novel structures, one of which was able to inhibit CK2 activity in living cells and promote tumor cell death. The essential features required for potent CK2 inhibitory activity of this class of compounds are discussed.

Identification of New Potent Acetylcholinesterase Inhibitors Using Virtual Screening and in vitro Approaches.

Acetylcholinesterase (AChE) is currently the most favorable target for the symptomatic treatment and reduction of Alzheimer's disease (AD). In order to identify new potent inhibitors of this enzyme, we describe herein a new structure-based virtual screening (SBVS) using the Institut Curie-CNRS chemical library (ICCL), which contained at the screening date 14307 compounds. The strategy undertaken in this work consisted of the use of several docking programs in SBVS calculations followed by the application of a consensus method (vSDC) and a scrupulous visual analysis. It allowed us to obtain a high degree of success, with a yield of almost 86 %, since 12 hits were identified among only 14 molecules tested in vitro. Still more remarkably, 6 of these hits were more active than galantamine, the reference inhibitor. These hits were predicted to have good ADMET properties. The two most promising compounds can serve as leads for AD treatment.

CO migration pathways in cytochrome P450cam studied by molecular dynamics simulations.

Previous laser flash photolysis investigations between 100 and 300 K have shown that the kinetics of CO rebinding with cytochrome P450(cam)(camphor) consist of up to four different processes revealing a complex internal dynamics after ligand dissociation. In the present work, molecular dynamics simulations were undertaken on the ternary complex P450(cam)(cam)(CO) to explore the CO migration pathways, monitor the internal cavities of the protein, and localize the CO docking sites. One trajectory of 1 nsec with the protein in a water box and 36 trajectories of 1 nsec in the vacuum were calculated. In each trajectory, the protein contained only one CO ligand on which no constraints were applied. The simulations were performed at 200, 300, and 320 K. The results indicate the presence of seven CO docking sites, mainly hydrophobic, located in the same moiety of the protein. Two of them coincide with xenon binding sites identified by crystallography. The protein matrix exhibits eight persistent internal cavities, four of which corresponding to the ligand docking sites. In addition, it was observed that water molecules entering the protein were mainly attracted into the polar pockets, far away from the CO docking sites. Finally, the identified CO migration pathways provide a consistent interpretation of the experimental rebinding kinetics.

Efficient conformational sampling and weak scoring in docking programs? Strategy of the wisdom of crowds.

BACKGROUND: In drug design, an efficient structure-based optimization of a ligand needs the precise knowledge of the protein-ligand interactions. In the absence of experimental information, docking programs are necessary for ligand positioning, and the choice of a reliable program is essential for the success of such an optimization. The performances of four popular docking programs, Gold, Glide, Surflex and FlexX, were investigated using 100 crystal structures of complexes taken from the Directory of Useful Decoys-Enhanced database. RESULTS: The ligand conformational sampling was rather efficient, with a correct pose found for a maximum of 84 complexes, obtained by Surflex. However, the ranking of the correct poses was not as efficient, with a maximum of 68 top-rank or 75 top-4 rank correct poses given by Glidescore. No relationship was found between either the sampling or the scoring performance of the four programs and the properties of either the targets or the small molecules, except for the number of ligand rotatable bonds. As well, no exploitable relationship was found between each program performance in docking and in virtual screening; a wrong top-rank pose may obtain a good score that allows it to be ranked among the most active compounds and vice versa. Also, to improve the results of docking, the strengths of the programs were combined either by using a rescoring procedure or the United Subset Consensus (USC). Oddly, positioning with Surflex and rescoring with Glidescore did not improve the results. However, USC based on docking allowed us to obtain a correct pose in the top-4 rank for 87 complexes. Finally, nine complexes were scrutinized, because a correct pose was found by at least one program but poorly ranked by all four programs. Contrarily to what was expected, except for one case, this was not due to weaknesses of the scoring functions. CONCLUSIONS: We conclude that the scoring functions should be improved to detect the correct poses, but sometimes their failure may be due to other varied considerations. To increase the chances of success, we recommend to use several programs and combine their results. Graphical abstract Summary of the results obtained by semi-rigid docking of crystallographic ligands. The docking was done on 100 protein-ligand X-ray structures, taken from the DUD-E database, and using four programs, Glide, Gold, Surflex and FlexX. Based on the docking results, we applied our United Subset Consensus method (USC), for which only the top4-rank poses are relevant. The number of complexes for which the best pose is correct, is represented by the gray boxes, the blue and red boxes correspond to the number of complexes with a correct pose ranked as the top 1 or within the top 4. A pose is considered correct when its root-mean-square deviation from the crystal structure is less than 2 Å.

Structure-based mutational analysis of ICAT residues mediating negative regulation of ß-catenin co-transcriptional activity.

ICAT (Inhibitor of ß-CAtenin and TCF) is a small acidic protein that negatively regulates ß-catenin co-transcriptional activity by competing with TCF/LEF factors in their binding to ß-catenin superhelical core. In melanoma cells, ICAT competes with LEF1 to negatively regulate the M-MITF and NEDD9 target genes. The structure of ICAT consists of two domains: the 3-helix bundle N-terminal domain binds to ß-catenin Armadillo (Arm) repeats 10-12 and the C-terminal tail binds to Arm repeats 5-9. To elucidate the structural mechanisms governing ICAT/ß-catenin interactions in melanoma cells, three ICAT residues Y15, K19 and V22 in the N-terminal domain, contacting hydrophobic ß-catenin residue F660, were mutated and interaction was assessed by immunoprecipitation. Despite the moderate hydrophobicity of the contact, its removal completely abolished the interaction. In the ICAT C-terminal tail consensus sequence, neutralization of the electrostatic interactions between residues D66, E75 and ß-catenin residues K435, K312, coupled to deletion of the hydrophobic contact between F71 and ß-catenin R386, markedly reduced, but failed to abolish the ICAT-mediated negative regulation of M-MITF and NEDD9 promoters. We conclude that in melanoma cells, anchoring of ICAT N-terminal domain to ß-catenin through the hook made by residue F660, trapped in the pincers formed by ICAT residues Y15 and V22, is crucial for stabilizing the ICAT/ß-catenin complex. This is a prerequisite for binding of the consensus peptide to Arm repeats 5-9 and competition with LEF1. Differences between ICAT and LEF1 in their affinity for ß-catenin may rely on the absence in ICAT of hydrophilic residues between D66 and F71.

Internal cavities and ligand passageways in human hemoglobin characterized by molecular dynamics simulations.

Molecular dynamics simulations of the unliganded T state of human hemoglobin showed the existence of a spontaneous, very wide cavity on the distal side of the alpha subunit. This cavity consists of three tunnels spreading from the vicinity of the iron atom (the ligand binding site) to the surface of the subunit, constituting possible passageways for the entrance of the ligand. A fourth passageway was characterized due to the trajectory of water molecules entering or leaving the heme pocket. Analogous passages were observed in the beta subunits. They all appear and disappear dynamically, although some parts of them are more persistent along the trajectories. The most persistent regions within these tunnels correspond to all the xenon docking sites of human cytoglobin and to some of those of sperm whale and horse heart myoglobins and group I truncated hemoglobins.

Benchmark of four popular virtual screening programs: construction of the active/decoy dataset remains a major determinant of measured performance.

BACKGROUND: In a structure-based virtual screening, the choice of the docking program is essential for the success of a hit identification. Benchmarks are meant to help in guiding this choice, especially when undertaken on a large variety of protein targets. Here, the performance of four popular virtual screening programs, Gold, Glide, Surflex and FlexX, is compared using the Directory of Useful Decoys-Enhanced database (DUD-E), which includes 102 targets with an average of 224 ligands per target and 50 decoys per ligand, generated to avoid biases in the benchmarking. Then, a relationship between these program performances and the properties of the targets or the small molecules was investigated. RESULTS: The comparison was based on two metrics, with three different parameters each. The BEDROC scores with α = 80.5, indicated that, on the overall database, Glide succeeded (score > 0.5) for 30 targets, Gold for 27, FlexX for 14 and Surflex for 11. The performance did not depend on the hydrophobicity nor the openness of the protein cavities, neither on the families to which the proteins belong. However, despite the care in the construction of the DUD-E database, the small differences that remain between the actives and the decoys likely explain the successes of Gold, Surflex and FlexX. Moreover, the similarity between the actives of a target and its crystal structure ligand seems to be at the basis of the good performance of Glide. When all targets with significant biases are removed from the benchmarking, a subset of 47 targets remains, for which Glide succeeded for only 5 targets, Gold for 4 and FlexX and Surflex for 2. CONCLUSION: The performance dramatic drop of all four programs when the biases are removed shows that we should beware of virtual screening benchmarks, because good performances may be due to wrong reasons. Therefore, benchmarking would hardly provide guidelines for virtual screening experiments, despite the tendency that is maintained, i.e., Glide and Gold display better performance than FlexX and Surflex. We recommend to always use several programs and combine their results. Graphical AbstractSummary of the results obtained by virtual screening with the four programs, Glide, Gold, Surflex and FlexX, on the 102 targets of the DUD-E database. The percentage of targets with successful results, i.e., with BDEROC(α = 80.5) > 0.5, when the entire database is considered are in Blue, and when targets with biased chemical libraries are removed are in Red.

vSDC: a method to improve early recognition in virtual screening when limited experimental resources are available.

BACKGROUND: In drug design, one may be confronted to the problem of finding hits for targets for which no small inhibiting molecules are known and only low-throughput experiments are available (like ITC or NMR studies), two common difficulties encountered in a typical academic setting. Using a virtual screening strategy like docking can alleviate some of the problems and save a considerable amount of time by selecting only top-ranking molecules, but only if the method is very efficient, i.e. when a good proportion of actives are found in the 1-10 % best ranked molecules. RESULTS: The use of several programs (in our study, Gold, Surflex, FlexX and Glide were considered) shows a divergence of the results, which presents a difficulty in guiding the experiments. To overcome this divergence and increase the yield of the virtual screening, we created the standard deviation consensus (SDC) and variable SDC (vSDC) methods, consisting of the intersection of molecule sets from several virtual screening programs, based on the standard deviations of their ranking distributions. CONCLUSIONS: SDC allowed us to find hits for two new protein targets by testing only 9 and 11 small molecules from a chemical library of circa 15,000 compounds. Furthermore, vSDC, when applied to the 102 proteins of the DUD-E benchmarking database, succeeded in finding more hits than any of the four isolated programs for 13-60 % of the targets. In addition, when only 10 molecules of each of the 102 chemical libraries were considered, vSDC performed better in the number of hits found, with an improvement of 6-24 % over the 10 best-ranked molecules given by the individual docking programs.Graphical abstractIn drug design, for a given target and a given chemical library, the results obtained with different virtual screening programs are divergent. So how to rationally guide the experimental tests, especially when only a few number of experiments can be made? The variable Standard Deviation Consensus (vSDC) method was developed to answer this issue. Left panel the vSDC principle consists of intersecting molecule sets, chosen on the basis of the standard deviations of their ranking distributions, obtained from various virtual screening programs. In this study Glide, Gold, FlexX and Surflex were used and tested on the 102 targets of the DUD-E database. Right panel Comparison of the average percentage of hits found with vSDC and each of the four programs, when only 10 molecules from each of the 102 chemical libraries of the DUD-E database were considered. On average, vSDC was capable of finding 38 % of the findable hits, against 34 % for Glide, 32 % for Gold, 16 % for FlexX and 14 % for Surflex, showing that with vSDC, it was possible to overcome the unpredictability of the virtual screening results and to improve them.

Integrating three views of Arf1 activation dynamics.

The proteins Arno and Gea2 of the Sec7 family can promote GDP-GTP exchange on Arf1, a small GTP-binding protein, which coordinates coated vesicle formation for protein transport within the cell. Crystal structures of the essential Sec7 domain (Sec7d) of Gea2 in the free and Arf1-bound forms suggest that conformational dynamics of the Sec7d as well as those of the G-protein play a role in nucleotide exchange. Starting from a set of complementary crystal structures, we compared the collective movements of unbound Gea2 and Arno Sec7 domains, Arf1-GDP, and the Arf1-Gea2(Sec7d) nucleotide-free complex using normal modes analyses. In all unbound Sec7d analyses, significant low-energy movements were found to lead to closure of the hydrophobic groove towards the form seen in the Arf1-Gea2(Sec7d) complex, suggesting that groove closure is a general feature of the Sec7 family. Low-energy movements in Arf1-GDP implicate critical switch 1 and 2 residues which are coupled to modifications in the myristoylated N-terminal-helix binding site at the other end of the "interswitch" beta hairpin. It is suggested that Sec7d groove closure upon docking of the two molecules may permit extraction of switch 1 from Arf1-GDP and prepare the complex for movement of the interswitch, which is central to the membrane-linked exchange activity. Large-scale collective movements in the Arf1-Sec7d complex appear to participate in the insertion of the Sec7d Glu finger into the GDP binding site to promote actual nucleotide release.

Molecular dynamics assignment of NMR correlation times to specific motions in a "basket-handle porphyrin" heme.

Iron (II) basket-handle porphyrins (BHP) are a series of encumbered heme models designed several years ago to mimic the ligand binding site of hemoproteins. Contrary to expectations, kinetic investigations have revealed that the k(on) rates for CO and/or O2 binding were only marginally affected by the assumed central steric hindrance of the iron atom. Thus, it was hypothesized that the internal dynamics of the molecule might be at the origin of the poor steric protection. To address this issue, measurements of nuclear magnetic resonance relaxation rates, fluorescence anisotropy experiments, and molecular dynamics simulations were undertaken. The size of BHP is small enough to allow the simulation in explicit chloroform with an almost complete sampling of the conformational space. The order parameters calculated from the MD trajectory compare well with the NMR experimental data and the predicted rotational correlation time corresponding to the Brownian motion of the molecule is in good agreement with the fluorescence measurements. Moreover, combining the results obtained using the three techniques allows the attribution of each internal NMR correlation time to a particular internal motion, revealing that even such medium-sized molecules are able to display quite complex internal dynamics. In particular, the handle phenyls that were assumed to sandwich the porphyrin have in fact a vanishing probability to be found in the proximity of the iron atom. They are therefore unable to reduce ligand accessibility significantly, which may explain the behavior of the k(on) rates.

Molecular dynamics of the FixJ receiver domain: movement of the beta4-alpha4 loop correlates with the in and out flip of Phe101.

FixJ is a two-domain response regulator involved in nitrogen fixation in Sinorhizobium meliloti. Recent X-ray characterization of both the native (unphosphorylated) and the active (phosphorylated) states of the protein identify conformational changes of the beta4-alpha4 loop and the conserved residue Phe101 as the key switches in activation. These structures also allowed investigation of the transition between conformations of this two-component regulatory receiver domain by molecular dynamics simulations. The path for the conformational change was studied with a distance constraint directing the system from one state to the other. The simulations provide evidence for a correlation between the conformation of the beta4-alpha4 loop and the orientation of the residue Phe101. A model presenting the sequence of events during the activation/deactivation process is discussed.

An atomistic view of microtubule stabilization by GTP.

A microtubule is a dynamic system formed of αß-tubulins. The presence of nonhydrolyzable guanosine-5'-triphosphate (GTP)/guanosine diphosphate (GDP) on the ß-tubulins provokes microtubule polymerization/depolymerization. Despite the large number of experimental studies of this dynamical process, its mechanism is still unclear. To provide insights into this mechanism we studied the first depolymerization steps of GDP/GTP-bound microtubules by normal-mode analysis with the all-atom model. We also constructed a depolymerizing microtubule and compared it to cryo-electron microscopy tomograms (cyro-ET). The results show that during depolymerization, the protofilaments not only curve but twist to weaken their lateral interactions. These interactions are stabilized by GTP, but not evenly. Not all of the interface residues are of equal importance: five of them, belonging to the H2-S3 loop, play a special role; acting as a lock whose key is the γ-phosphate of GTP. Sequence alignments of several tubulins confirm the importance of these residues.

Novel 8-arylated purines as inhibitors of glycogen synthase kinase.

A series of 8-arylated purine derivatives bearing either an aniline or an alkyl amide at position 6 were found to inhibit glycogen synthase kinase-3, with good selectivity over ten kinases. Molecular modeling studies indicated that the most active compounds (8a and 8e), adopt a planar conformation, close to the shape of AMPPNP in the crystal structure of GSK-3. These compounds are stabilized by hydrophobic contacts between the 8-aromatic group and the protein adenine pocket and by electrostatic contacts.