Search

The occurrence of a G-triplex folding intermediate of thrombin binding aptamer (TBA) has been recently predicted by metadynamics calculations, and experimentally supported by Nuclear Magnetic Resonance (NMR), Circular Dichroism (CD) and Differential Scanning Calorimetry (DSC) data collected on a 3′ end TBA-truncated 11-mer oligonucleotide (11-mer-3′-t-TBA). Here we present the solution structure of 11-mer-3′-t-TBA in the presence of potassium ions. This structure is the first experimental example of a G-triplex folding, where a network of Hoogsteen-like hydrogen bonds stabilizes six guanines to form two G:G:G triad planes. The G-triplex folding of 11-mer-3′-t-TBA is stabilized by the potassium ion and destabilized by increasing the temperature. The superimposition of the experimental structure with that predicted by metadynamics shows a great similarity, with only significant differences involving two loops. These new structural data show that 11-mer-3′-t-TBA assumes a G-triplex DNA conformation as its stable form, reinforcing the idea that G-triplex folding intermediates may occur in vivo in human guanine-rich sequences. NMR and CD screening of eight different constructs obtained by removing from one to four bases at either the 3′ and the 5′ ends show that only the 11-mer-3′-t-TBA yields a relatively stable G-triplex., Nucleic Acids Research, 42 (21), ISSN:1362-4962, ISSN:0301-5610

Bile acids are extensively investigated for their potential in the treatment of human disorders. The liver X receptors (LXRs), activated by oxysterols and by a secondary bile acid named hyodeoxycholic acid (HDCA), have been found essential in the regulation of lipid homeostasis in mammals. Unfortunately, LXRα activates lipogenic enzymes causing accumulation of lipid in the liver. In addition to LXRs, HDCA has been also shown to function as ligand for GPBAR1, a G protein coupled receptor for secondary bile acids whose activation represents a promising approach to liver steatosis. In the present study, we report a library of HDCA derivatives endowed with modulatory activity on the two receptors. The lead optimization of HDCA moiety was rationally driven by the structural information on the binding site of the two targets and results from pharmacological characterization allowed the identification of hyodeoxycholane derivatives with selective agonistic activity toward LXRα and GPBAR1 and notably to the identification of the first example of potent dual LXRα/GPBAR1 agonists. The new chemical entities might hold utility in the treatment of dyslipidemic disorders.

Bile acids are signaling molecules interacting with the nuclear receptor FXR and the G-protein coupled receptor 1 (GP-BAR1/TGR5). GP-BAR1 is a promising pharmacological target for the treatment of steatohepatitis, type 2 diabetes, and obesity. Endogenous bile acids and currently available semisynthetic bile acids are poorly selective toward GP-BAR1 and FXR. Thus, in the present study we have investigated around the structure of UDCA, a clinically used bile acid devoid of FXR agonist activity, to develop a large family of side chain modified 3α,7β-dihydroxyl cholanoids that selectively activate GP-BAR1. In vivo and in vitro pharmacological evaluation demonstrated that administration of compound 16 selectively increases the expression of pro-glucagon 1, a GP-BAR1 target, in the small intestine, while it had no effect on FXR target genes in the liver. Further, compound 16 results in a significant reshaping of bile acid pool in a rodent model of cholestasis. These data demonstrate that UDCA is a useful scaffold to generate novel and selective steroidal ligands for GP-BAR1.

A detailed description of the events ruling ligand/protein interaction and an accurate estimation of the drug affinity to its target is of great help in speeding drug discovery strategies. We have developed a metadynamics-based approach, named funnel metadynamics, that allows the ligand to enhance the sampling of the target binding sites and its solvated states. This method leads to an efficient characterization of the binding free-energy surface and an accurate calculation of the absolute protein–ligand binding free energy. We illustrate our protocol in two systems, benzamidine/trypsin and SC-558/cyclooxygenase 2. In both cases, the X-ray conformation has been found as the lowest free-energy pose, and the computed protein–ligand binding free energy in good agreement with experiments. Furthermore, funnel metadynamics unveils important information about the binding process, such as the presence of alternative binding modes and the role of waters. The results achieved at an affordable computational cost make funnel metadynamics a valuable method for drug discovery and for dealing with a variety of problems in chemistry, physics, and material science.

An exhaustive description of the molecular recognition mechanism between a ligand and its biological target is of great value because it provides the opportunity for an exogenous control of the related process. Very often this aim can be pursued using high resolution structures of the complex in combination with inexpensive computational protocols such as docking algorithms. Unfortunately, in many other cases a number of factors, like protein flexibility or solvent effects, increase the degree of complexity of ligand/protein interaction and these standard techniques are no longer sufficient to describe the binding event. We have experienced and tested these limits in the present study in which we have developed and revealed the mechanism of binding of a new series of potent inhibitors of Adenosine Deaminase. We have first performed a large number of docking calculations, which unfortunately failed to yield reliable results due to the dynamical character of the enzyme and the complex role of the solvent. Thus, we have stepped up the computational strategy using a protocol based on metadynamics. Our approach has allowed dealing with protein motion and solvation during ligand binding and finally identifying the lowest energy binding modes of the most potent compound of the series, 4-decyl-pyrazolo[1,5- a ]pyrimidin-7-one.

The widely used nonsteroidal anti-inflammatory drugs block the cyclooxygenase enzymes (COXs) and are clinically used for the treatment of inflammation, pain, and cancers. A selective inhibition of the different isoforms, particularly COX-2, is desirable, and consequently a deeper understanding of the molecular basis of selective inhibition is of great demand. Using an advanced computational technique we have simulated the full dissociation process of a highly potent and selective inhibitor, SC-558, in both COX-1 and COX-2. We have found a previously unreported alternative binding mode in COX-2 explaining the time-dependent inhibition exhibited by this class of inhibitors and consequently their long residence time inside this isoform. Our metadynamics-based approach allows us to illuminate the highly dynamical character of the ligand/protein recognition process, thus explaining a wealth of experimental data and paving the way to an innovative strategy for designing new COX inhibitors with tuned selectivity.

NMDA receptors are glutamate-gated ion channels (iGluRs) that are involved in several important physiological functions such as neuronal development, synaptic plasticity, learning, and memory. Among iGluRs, NMDA receptors have been perhaps the most actively investigated for their role in chronic neurodegeneration such as Alzheimer's, Parkinson's, and Huntington's diseases. Recent studies have shown that the NTD of subunit NR2B modulates ion channel gating through the binding of allosteric modulators such as the prototypical compound ifenprodil. In the present paper, the construction of a three-dimensional model for the NR2B modulatory domain is described and docking calculations allow, for the first time, definition of the ifenprodil binding pose at an atomic level and fully explain all the available structure-activity relationships. Moreover, in an attempt to add further insight into the ifenprodil mechanism of action, as it is not completely clear if it binds and stabilizes an open or a closed conformation of the NR2B modulatory domain, a matter, which is fundamental for the rational design of NMDA antagonists, MD simulations followed by an MM-PBSA analysis were performed. These calculations reveal that the closed conformation of the R1-R2 domain, rather than the open, constitutes the high affinity binding site for ifenprodil and that a profound stabilization of the closed conformation upon ifenprodil binding occurs. Thus, for a rational design and/or for virtual screening experiments, the closed conformation of the R1-R2 domain should be taken into account and our 3D model can provide valuable hints for the design of NR2B-selective antagonists. © 2007 Wiley-VCH Verlag GmbH & Co. KGaA.

Adenylate Kinase (AK) is a signal transducing protein that regulates cellular energy homeostasis balancing between different conformations. An alteration of its activity can lead to severe pathologies such as heart failure, cancer and neurodegenerative diseases. A comprehensive elucidation of the large-scale conformational motions that rule the functional mechanism of this enzyme is of great value to guide rationally the development of new medications. Here using a metadynamics-based computational protocol we elucidate the thermodynamics and structural properties underlying the AK functional transitions. The free energy estimation of the conformational motions of the enzyme allows characterizing the sequence of events that regulate its action. We reveal the atomistic details of the most relevant enzyme states, identifying residues such as Arg119 and Lys13, which play a key role during the conformational transitions and represent druggable spots to design enzyme inhibitors. Our study offers tools that open new areas of investigation on large-scale motion in proteins., Scientific Reports, 5, ISSN:2045-2322

We are grateful to Prof. Bernard Knoops and Andre Clippe from University of Louvain-la-Neuve, Belgium, for the generous gifts of biological materials to express the human peroxiredoxin in Escherichia coli. We thank Acellera Company for providing of a complementary 6 month ACEMD multi-GPU software licence and for the generous gift of a Nvidia Tesla K20 supercomputing unit; International audience; One of the intrinsic properties of proteins is their capacity to interact selectively with other molecules in their environment, inducing many chemical equilibria each differentiated by the mutual affinities of the components. A comprehensive understanding of these molecular binding processes at atomistic resolution requires formally the complete description of the system dynamics and statistics at the relevant time scales. While solution NMR observables are averaged over different time scales, from picosecond to second, recent new molecular dynamics protocols accelerated considerably the simulation time of realistic model systems. Based on known ligands recently discovered either by crystallography or NMR for the human peroxiredoxin 5, their affinities were for the first time accurately evaluated at atomistic resolution comparing absolute binding free-energy estimated by funnel-metadynamics simulations and solution NMR experiments. In particular, free-energy calculations are demonstrated to discriminate two closely related ligands as pyrocatechol and 4-methylpyrocathecol separated just by 1 kcal/mol in aqueous solution. The results provide a new experimental and theoretical basis for the estimation of ligand-protein affinities.

Glycoproteins are often recognized as not-self molecules by antibodies triggering the onset of severe autoimmune diseases such as Multiple Sclerosis (MS). Thus, the development of antigen-mimicking biomarkers represents an attractive strategy for an early diagnosis of the disease. An example is the synthetic glycopeptide CSF114(Glc), which was designed and tested as MS biomarker and whose clinical application was limited by its reduced ability to detect autoantibodies in MS patients. In the attempt to improve the efficacy of CSF114(Glc), we have characterized all the events leading to the final binding of the biomarker to the autoantibody using atomistic simulations, ESR and NMR experiments. The glycosydic moiety plays a primary role in the whole process. In particular, in an environment mimicking that used in the clinical tests the glycopeptide assumes a α-helix structure that is functional for the interaction with the antibody. In this conformation CSF114(Glc) binds the monoclonal antibody mAb8-18C5 similarly to the myelin oligodendrocyte glycoprotein MOG, which is a known MS auto-antigen, thus explaining its diagnostic activity. Our study offers new molecular bases to design more effective biomarkers and provides a most valid protocol to investigate other systems where the environment effect is determinant for the biological activity.

During the chromatographic purification of organic extracts obtained from Plakortis cfr. lita we obtained three highly degraded steroid derivatives, the known incisterol A2 (1) and the new incisterols A5 (2) and A6 (3). The new compounds were characterized basing on NMR and MS evidences along with comparison with model compounds. Incisterol A5 proved to bear a 17S-ethyl-15E,18-diene (incisterol numbering system) side chain, found for the first time in a marine organism. The new incisterols A5 and A6 proved to be potent inducers of transactivation of the pregnane X receptor (PXR) and they also stimulate the expression of CYP7A4 and MDR1 with a potency comparable to that of Rifaximin. These observations prompt to consider incisterols A5 and A6 as new potent agonists of PXR. On the other hand, the 17R-methyl analogue incisterol A2 shows only a poor PXR agonist activity. Molecular docking simulations elucidated the binding mechanism of the active incisterols in the ligand binding domain of PXR.

Intervention in integrin-mediated cell adhesion and integrin signaling pathways is an ongoing area of research in medicinal chemistry and drug development. One key element in integrin-ligand interaction is the coordination of the bivalent cation at the metal ion-dependent adhesion site (MIDAS) by a carboxylic acid function, a consistent feature of all integrin ligands. With the exception of the recently discovered hydroxamic acids, all bioisosteric attempts to replace the carboxylic acid of integrin ligands failed. We report that phosphinates as well as monomethyl phosphonates represent excellent isosters, when introduced into integrin antagonists for the platelet integrin αIIbβ3. The novel inhibitors exhibit in vitro and ex vivo activities in the low nanomolar range. Steric and charge requirements of the MIDAS region were unraveled, thus paving the way for an in silico prediction of ligand activity and in turn the rational design of the next generation of integrin antagonists.

3-Substituted 6-phenyl-4H-imidazo[1,5-a][1,4]benzodiazepines and related compounds were synthesized as central benzodiazepine receptor (CBR) ligands. Most of the compounds showed high affinity for bovine and human CBR, their K(i) values spanning from the low nanomolar to the submicromolar range. In particular, imidazoester 5f was able to promote a massive flow of (36)Cl(-) in rat cerebrocortical synaptoneurosomes overlapping its efficacy profile with that of a typical full agonist. Compound 5f was then examined in mice for its pharmacological effects where it proved to be a safe anxiolytic agent devoid of the unpleasant myorelaxant and amnesic effects of the classical 1,4-benzodiazepines. Moreover, the selectivity of some selected compounds has been assessed in recombinant α(1)β(2)γ(2)L, α(2)β(1)γ(2)L, and α(5)β(2)γ(2)L human GABA(A) receptors. Finally, some compounds were submitted to molecular docking calculations along with molecular dynamics simulations in the Cromer's GABA(A) homology model.

Activated leukocyte cell adhesion molecule (ALCAM) plays a relevant role in tumor biology and progression. Our previous studies showed that ALCAM is expressed at the surface of epithelial ovarian cancer (EOC) cells and is released in a soluble form by ADAM-17-mediated shedding. This process is relevant to EOC cell motility and invasiveness, which is reduced by nonspecific inhibitors of ADAM-17. For this reason, ADAM-17 may represent a new useful target in anticancer therapy. Herein, we report the synthesis and biological evaluation of new ADAM-17 inhibitors containing an arylsulfonamidic scaffold. Among the new potential inhibitors, two very promising compounds 17 and 18 were discovered, with a nanomolar activity for ADAM-17 isolated enzyme. These compounds proved to be also the most potent in inhibiting soluble ALCAM release in cancer cells, showing a nanomolar activity on A2774 and SKOV3 cell lines.

Here we report the synthesis of a number of compounds structurally related to arginine methyltransferase inhibitor 1 (AMI-1). The structural alterations that we made included: 1) the substitution of the sulfonic groups with the bioisosteric carboxylic groups; 2) the replacement of the ureidic function with a bis-amidic moiety; 3) the introduction of a N-containing basic moiety; and 4) the positional isomerization of the aminohydroxynaphthoic moiety. We have assessed the biological activity of these compounds against a panel of arginine methyltransferases (fungal RmtA, hPRMT1, hCARM1, hPRMT3, hPRMT6) and a lysine methyltransferase (SET7/9) using histone and nonhistone proteins as substrates. Molecular modeling studies for a deep binding-mode analysis of test compounds were also performed. The bis-carboxylic acid derivatives 1 b and 7 b emerged as the most effective PRMT inhibitors, both in vitro and in vivo, being comparable or even better than the reference compound (AMI-1) and practically inactive against the lysine methyltransferase SET7/9.

Adenosine induces glioma cell proliferation by means of an antiapoptotic effect, which is blocked by cotreatment with selective A(3) AR antagonists. In this study, a novel series of N(2)-substituted pyrazolo[3,4-d]pyrimidines 2a-u was developed as highly potent and selective A(3) AR antagonists. The most performing compounds were derivatives 2a (R(1) = CH(3) and R(2) = COC(6)H(5); K(i) 334, 728, and 0.60 nM at the human A(1), A(2A), and A(3) ARs, respectively) and 2b (R(1) = CH(3) and R(2) = COC(6)H(4)-4-OCH(3); K(i) 1037, 3179, and 0.18 nM at the human A(1), A(2A), and A(3) ARs, respectively), which counteracted the effect of the A(3) AR agonists Cl-IB-MECA and IB-MECA on human glioma U87MG cell proliferation. This effect was concentration-dependent, with IC(50) values comparable to A(3) AR binding affinity values of 2a and 2b, thereby suggesting that their effects were receptor-mediated. Furthermore, the antiproliferative activity of the new compounds was demonstrated to be mediated by the block of A(3) AR agonist activation of intracellular kinases ERK 1/2.

The study of DNA G-quadruplex stabilizers has enjoyed a great momentum in the late years due to their application as anticancer agents. The recognition of the grooves of these structural motifs is expected to result in a higher degree of selectivity over other DNA structures. Therefore, to achieve an enhanced knowledge on the structural and conformational requisites for quadruplex groove recognition, distamycin A, the only compound for which a pure groove binding has been proven, has been chemically modified. Surprisingly, structural and thermodynamic studies revealed that the absence of Coulombic interactions results in an unprecedented binding position in which both the groove and the 3' end of the DNA are occupied. This further contribution adds another piece to the so far elusive puzzle of the recognition between ligands and DNA quadruplexes and will serve as a platform for a rational design of new groove binders.

Ethyl 8-fluoro-6-(4-nitrophenyl)- and ethyl 8-fluoro-6-(3-nitrophenyl)-4 H-imidazo[1,5-a][1,4]benzodiazepine 3-carboxylate 6 and 7 were synthesized as central benzodiazepine receptor (CBR) ligands and tested for their ability to displace [(3)H]flumazenil from bovine and human cortical brain membranes. Both compounds showed high affinity for bovine and human CBR. In particular, compound 7 emerged as the most interesting compound, having a partial agonist profile in vitro while possessing useful activity in various animal models of anxiety. In accordance with its partial agonist profile, compound 7 was devoid of typical benzodiazepine side effects. The homology model of the GABA A receptor developed by Cromer et al. was used to assess the binding modes of ligands 6 and 7. From our docking results, the partial agonist activity elicited by compound 7 is likely to be due to the 3'-nitro substituent, which is in the appropriate position to interact with Thr193 of the gamma 2-subunit by means of a hydrogen bond.

A number of 1,2,4-oxadiazol-5-yl-acetic acids and oxazol-4-yl-acetic acids were synthesized and tested for their ability to inhibit aldose reductase (ALR2). The oxadiazole derivatives, 7c, 7f, 7i, and 8h, 8i, proved to be the most active compounds, exhibiting inhibitory levels in the submicromolar range. In this series, the phenyl group turned out to be the preferred substitution pattern, as its lengthening to a benzyl moiety determined a general reduction of the inhibitory potency. The lead compound, 2-[3-(4-methoxyphenyl)-1,2,4-oxadiazol-5-yl]acetic acid, 7c, showed an excellent in vivo activity, proving to prevent cataract development in severely galactosemic rats when administered as an eye-drop solution in the precorneal region of the animals. Computational studies on the ALR2 inhibitors were performed to rationalize the structure-activity relationships observed and to provide the basis for further structure-guided design of novel ALR2 inhibitors.

The Cdc25 dual specificity phosphatases have central roles in coordinating cellular signalling processes and cell proliferation. It has been reported that an improper amplification or activation of these enzymes is a distinctive feature of a number of human cancers, including breast cancers. Thus, the inhibition of Cdc25 phosphatases might provide a novel approach for the discovery of new and selective antitumor agents. By using the crystal structure of the catalytic domain of Cdc25B, structural models for the interaction of various Cdc25B inhibitors (1-13) with the enzyme were generated by computational docking. The parallel use of two efficient and predictive docking programs, AutoDock and GOLD, allowed mutual validation of the predicted binding poses. To evaluate their quality, the models were validated with known structure-activity relationships and site-directed mutagenesis data. The results provide an improved basis for structure-based ligand design and suggest a possible explanation for the inhibition mechanism of the examined Cdc25B ligands. We suggest that the recurring motif of a tight interaction between the inhibitor and the two arginine residues, 482 and 544, is of prime importance for reversible enzyme inhibition. In contrast, the irreversible inhibition mechanism of 1-4 seems to be associated with the close vicinity of the quinone ring and the Cys473 catalytic thiolate. We believe that this extensive study might provide useful hints to guide the development of new potent Cdc25B inhibitors as novel anticancer drugs. © 2006 Wiley-VCH Verlag GmbH & Co. KGaA.

The PLUMED consortium unifies developers and contributors to PLUMED, an open-source library for enhanced-sampling, free-energy calculations and the analysis of molecular dynamics simulations. Here, we outline our efforts to promote transparency and reproducibility by disseminating protocols for enhanced-sampling molecular simulations.