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

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.

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.

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.

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.

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.

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.

Structure, dynamics and thermodynamics of the human centrin 2/hSfi1 complex.

Centrin, an EF-hand calcium-binding protein, has been shown to be involved in the duplication of centrosomes, and Sfi1 (Suppressor of fermentation-induced loss of stress resistance protein 1) is one of its centrosomal targets. There are three isoforms of human centrin, but here we only considered centrin 2 (HsCen2). This protein has the ability to bind to any of the approximately 25 repeats of human Sfi1 (hSfi1) with more or less affinity. In this study, we mainly focused on the 17th repeat (R17-hSfi1-20), which presents the highest level of similarity with a well-studied 17-residue peptide (P17-XPC) from human xeroderma pigmentosum complementation group C protein, another centrin target for DNA repair. The only known structure of HsCen2 was resolved in complex with P17-XPC. The 20-residue peptide R17-hSfi1-20 exhibits the motif L8L4W1, which is the reverse of the XPC motif, W1L4L8. Consequently, the dipole of the helix formed by this motif has a reverse orientation. We wished to ascertain the impact of this reversal on the structure, dynamics and affinity of centrin. To address this question, we determined the structure of C-HsCen2 [the C-terminal domain of HsCen2 (T94-Y172)] in complex with R17-hSfi1-20 and monitored its dynamics by NMR, after having verified that the N-terminal domain of HsCen2 does not interact with the peptide. The structure shows that the binding mode is similar to that of P17-XPC. However, we observed a 2 -A translation of the R17-hSfi1-20 helix along its axis, inducing less anchorage in the protein and the disruption of a hydrogen bond between a tryptophan residue in the peptide and a well-conserved nearby glutamate in C-HsCen2. NMR dynamic studies of the complex strongly suggested the existence of an unusual calcium secondary binding mode in calcium-binding loop III, made possible by the uncommon residue composition of this loop. The secondary metal site is only populated at high calcium concentration and depends on the type of bound ligand.

What determines the degree of compactness of a calcium-binding protein?

The EF-hand calcium-binding proteins may exist either in an extended or a compact conformation. This conformation is sometimes correlated with the function of the calcium-binding protein. For those proteins whose structure and function are known, calcium sensors are usually extended and calcium buffers compact; hence, there is interest in predicting the form of the protein starting from its sequence. In the present study, we used two different procedures: one that already exists in the literature, the sosuidumbbell algorithm, mainly based on the charges of the two EF-hand domains, and the other comprising a novel procedure that is based on linker average hydrophilicity. The linker consists of the residues that connect the domains. The two procedures were tested on 17 known-structure calcium-binding proteins and then applied to 59 unknown-structure centrins. The sosuidumbbell algorithm yielded the correct conformations for only 15 of the known-structure proteins and predicted that all centrins should be in a closed form. The linker average hydrophilicity procedure discriminated well between all the extended and non-extended forms of the known-structure calcium-binding proteins, and its prediction concerning centrins reflected well their phylogenetic classification. The linker average hydrophilicity criterion is a simple and powerful means to discriminate between extended and non-extended forms of calcium-binding proteins. What is remarkable is that only a few residues that constitute the linker (between 2 and 20 in our tested sample of proteins) are responsible for the form of the calcium-binding protein, showing that this form is mainly governed by short-range interactions.

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.

Expanding the chemical diversity of CK2 inhibitors.

None of the already described CK2 inhibitors did fulfill the requirements for successful clinical settings. In order to find innovative CK2 inhibitors based on new scaffolds, we have performed a high-throughput screening of diverse chemical libraries. We report here the identification and characterization of several classes of new inhibitors. Whereas some share characteristics of previously known CK2 inhibitors, others are chemically unrelated and may represent new opportunities for the development of better CK2 inhibitors. By combining structure-activity relationships with a docking procedure, we were able to determine the binding mode of these inhibitors. Interestingly, beside the identification of several nanomolar ATP-competitive inhibitors, one class of chemical inhibitors displays a non-ATP competitive mode of inhibition, a feature that suggests that CK2 possess distinct druggable binding sites. For the most promising inhibitors, selectivity profiling was performed. We also provide evidence that some chemical compounds are inhibiting CK2 in living cells. Finally, the collected data allowed us to draw the rules about the chemical requirements for CK2 inhibition both in vitro and in a cellular context.

Nucleosome chiral transition under positive torsional stress in single chromatin fibers.

Using magnetic tweezers to investigate the mechanical response of single chromatin fibers, we show that fibers submitted to large positive torsion transiently trap positive turns at a rate of one turn per nucleosome. A comparison with the response of fibers of tetrasomes (the [H3-H4](2) tetramer bound with approximately 50 bp of DNA) obtained by depletion of H2A-H2B dimers suggests that the trapping reflects a nucleosome chiral transition to a metastable form built on the previously documented right-handed tetrasome. In view of its low energy, <8 kT, we propose that this transition is physiologically relevant and serves to break the docking of the dimers on the tetramer that in the absence of other factors exerts a strong block against elongation of transcription by the main RNA polymerase.

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.