Search

The development of specific antiviral therapies targeting SARS-CoV-2 remains fundamental because of the continued high incidence of COVID-19 and limited accessibility to antivirals in some countries. In this context, dark chemical matter (DCM), a set of drug-like compounds with outstanding selectivity profiles that have never shown bioactivity despite being extensively assayed, appears to be an excellent starting point for drug development. Accordingly, in this study, we performed a high-throughput screening to identify inhibitors of the SARS-CoV-2 main protease (M pro ) using DCM compounds as ligands. Multiple receptors and two different docking scoring functions were employed to identify the best molecular docking poses. The selected structures were subjected to extensive conventional and Gaussian accelerated molecular dynamics. From the results, four compounds with the best molecular behavior and binding energy were selected for experimental testing, one of which presented inhibitory activity with a K i value of 48 ± 5 μM. Through virtual screening, we identified a significant starting point for drug development, shedding new light on DCM compounds.

SARS-CoV-2 must bind its principal receptor, ACE2, on the target cell to initiate infection. This interaction is largely driven by the receptor binding domain (RBD) of the viral Spike (S) protein. Accordingly, antiviral compounds that can block RBD/ACE2 interactions can constitute promising antiviral agents. To identify such molecules, we performed a virtual screening of the Selleck FDA approved drugs and the Selleck database of Natural Products using a multistep computational procedure. An initial set of candidates was identified from an ensemble docking process using representative structures determined from the analysis of four 3 μ s molecular dynamics trajectories of the RBD/ACE2 complex. Two procedures were used to construct an initial set of candidates including a standard and a pharmacophore guided docking procedure. The initial set was subsequently subjected to a multistep sieving process to reduce the number of candidates to be tested experimentally, using increasingly demanding computational procedures, including the calculation of the binding free energy computed using the MMPBSA and MMGBSA methods. After the sieving process, a final list of 10 candidates was proposed, compounds which were subsequently purchased and tested ex-vivo. The results identified estradiol cypionate and telmisartan as inhibitors of SARS-CoV-2 entry into cells. Our findings demonstrate that the methodology presented here enables the discovery of inhibitors targeting viruses for which high-resolution structures are available., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

In solid tumors, such as breast cancer, hypoxic microenvironment worsens patient prognoses. We have previously reported in MCF-7 breast cancer cells that, under hypoxic conditions, hydroxytyrosol (HT) downregulates the level of reactive oxygen species, reduces the expression of hypoxia inducible factor-1 α (HIF-1 α ), and, at high concentrations, can bind to the aryl hydrocarbon receptor (AhR). With this background, the present study investigated whether the most abundant extra virgin olive oil (EVOO) phenolic compound tyrosol (TYR), with a chemical structure similar to HT but with only one hydroxyl group, exerts comparable effects. Our results revealed that, although TYR did not show any antioxidant activity in hypoxic MCF-7 cells, it inhibited the PI3K/Akt/mTOR/S6 kinase (S6K) pathway and reduced the expression of HIF-1 α and some of its target genes. Besides, TYR showed a lower binding affinity with the cytosolic transcription factor AhR, and even reduced its transcriptional activity. Some of these results are positive to control tumor progression in a hypoxic environment; however, they are observed at doses unachievable with diet intake or nutraceutical presentations. Considering that EVOO phenols can have synergistic effects, a mixture of low doses of TYR and other phenols could be useful to achieve these beneficial effects.

Over 750 million cases of COVID-19, caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), have been reported since the onset of the global outbreak. The need for effective treatments has spurred intensive research for therapeutic agents based on pharmaceutical repositioning or natural products. In light of prior studies asserting the bioactivity of natural compounds of the autochthonous Peruvian flora, the present study focuses on the identification SARS-CoV-2 M pro main protease dimer inhibitors. To this end, a target-based virtual screening was performed over a representative set of Peruvian flora-derived natural compounds. The best poses obtained from the ensemble molecular docking process were selected. These structures were subjected to extensive molecular dynamics steps for the computation of binding free energies along the trajectory and evaluation of the stability of the complexes. The compounds exhibiting the best free energy behaviors were selected for in vitro testing, confirming the inhibitory activity of Hyperoside against M pro , with a K i value lower than 20 µM, presumably through allosteric modulation.

In the present work, the interactions of the novel kinase inhibitors BI-2536, Volasetib (BI-6727) and Ro-3280 with the pharmacological target PLK1 have been studied by fluorescence spectroscopy and molecular dynamics calculations. High Stern-Volmer constants were found in fluorescence experiments suggesting the formation of stable protein-ligand complexes. In addition, it was observed that the binding constant between BI-2536 and PLK1 increases about 100-fold in presence of the phosphopeptide Cdc25C-p that docks to the polo box domain of the protein and releases the kinase domain. All the determined binding constants are higher for the kinase inhibitors than for their competitor for the active center (ATP) being BI-2536 and Volasertib the inhibitors that showed more affinity for PLK1. Calculated binding free energies confirmed the higher affinity of PLK1 for BI-2536 and Volasertib than for ATP. The higher affinity of the inhibitors to PLK1 compared to ATP was mainly attributed to stronger van der Waals interactions. Results may help with the challenge of designing and developing new kinase inhibitors more effective in clinical cancer therapy., (Copyright © 2022. Published by Elsevier B.V.)

SARS-CoV-2 is a type of coronavirus responsible for the international outbreak of respiratory illness termed COVID-19 that forced the World Health Organization to declare a pandemic infectious disease situation of international concern at the beginning of 2020. The need for a swift response against COVID-19 prompted to consider different sources to identify bioactive compounds that can be used as therapeutic agents, including available drugs and natural products. Accordingly, this work reports the results of a virtual screening process aimed at identifying antiviral natural product inhibitors of the SARS-CoV-2 M pro viral protease. For this purpose, ca. 2000 compounds of the Selleck database of Natural Compounds were the subject of an ensemble docking process targeting the M pro protease. Molecules that showed binding to most of the protein conformations were retained for a further step that involved the computation of the binding free energy of the ligand-M pro complex along a molecular dynamics trajectory. The compounds that showed a smooth binding free energy behavior were selected for in vitro testing. From the resulting set of compounds, five compounds exhibited an antiviral profile, and they are disclosed in the present work.

Diverse computational methods to support fragment-based drug discovery (FBDD) are available in the literature. Despite their demonstrated efficacy in supporting FBDD campaigns, they exhibit some drawbacks such as protein denaturation or ligand aggregation that have not yet been clearly overcome in the framework of biomolecular simulations. In the present work, we discuss a systematic semi-automatic novel computational procedure, designed to surpass these difficulties. The method, named fragment dissolved Molecular Dynamics (fdMD), utilizes simulation boxes of solvated small fragments, adding a repulsive Lennard-Jones potential term to avoid aggregation, which can be easily used to solvate the targets of interest. This method has the advantage of solvating the target with a low number of ligands, thus preventing the denaturation of the target, while simultaneously generating a database of ligand-solvated boxes that can be used in further studies. A number of scripts are made available to analyze the results and obtain the descriptors proposed as a means to trustfully discard spurious binding sites. To test our method, four test cases of different complexity have been solvated with ligand boxes and four molecular dynamics runs of 200 ns length have been run for each system, which have been extended up to 1 μs when needed. The reported results point out that the selected number of replicas are enough to identify the correct binding sites irrespective of the initial structure, even in the case of proteins having several close binding sites for the same ligand. We also propose a set of descriptors to analyze the results, among which the average MMGBSA and the average KDEEP energies have emerged as the most robust ones.

Olive oil intake has been linked with a lower incidence of breast cancer. Hypoxic microenvironment in solid tumors, such as breast cancer, is known to play a crucial role in cancer progression and in the failure of anticancer treatments. HIF-1 is the foremost effector in hypoxic response, and given that hydroxytyrosol (HT) is one of the main bioactive compounds in olive oil, in this study we deepen into its modulatory role on HIF-1. Our results in MCF-7 breast cancer cells demonstrate that HT decreases HIF-1α protein, probably by downregulating oxidative stress and by inhibiting the PI3K/Akt/mTOR pathway. Strikingly, the expression of HIF-1 target genes does not show a parallel decrease. Particularly, adrenomedullin and vascular endothelial growth factor are up-regulated by high concentrations of HT even in HIF-1α silenced cells, pointing to HIF-1-independent mechanisms of regulation. In fact, we show, by in silico modelling and transcriptional analysis, that high doses of HT may act as an agonist of the aryl hydrocarbon receptor favoring the induction of these angiogenic genes. In conclusion, we suggest that the effect of HT in a hypoxic environment is largely affected by its concentration and involves both HIF-1 dependent and independent mechanisms.

Apoptosis is a key cell death pathway in mammalian cells. Understanding this process and its regulation has been a subject of study in the last three decades. Members of the Bcl-2 family of proteins are involved in the regulation of apoptosis through mitochondrial poration with the subsequent initiation of apoptosis. Deregulation of proapoptotic proteins contributes to the progression of many tumor processes. Understanding how these pore-forming Bcl-2 proteins Bak and Bax are activated is key to find new anticancer treatments. As no drug capable of activating Bak has been disclosed yet, the study of the structural features of BH3 peptides-known as Bak activators-relevant for binding along with its binding energy decomposition analysis, becomes essential for designing novel small-molecule mimics of BH3. Interestingly, a BH3 Bim analogue-inactivating Bak has recently been discovered, opening a question on the molecular features that determine the functions of BH3 peptides. Therefore, the present work is aimed at understanding the way BH3 peptides activate or inactivate Bak in order to identify differential structural features that can be used in drug design. For this purpose, complexes of Bak with an activator and an inhibitor have been subjected to a molecular dynamics study. Structural differences were assessed by means of the fluctuations of the corresponding principal components. Moreover, the MMPB/GBSA approach was used to compute the binding free energy of the diverse complexes to identify those residues of the BH3 peptide that exhibit the larger contributions to complex formation. The results obtained in this work show differences between activators and inhibitors, both in structural and energetic terms, which can be used in the design of new molecules that can activate or inactivate proapoptotic Bak.

Aim: Calmodulin interacts in many different ways with its ligands. We aim to shed light on its plasticity analyzing the changes followed by the linker region and the relative position of the lobes using conventional molecular dynamics, accelerated MD and scaled MD (sMD). Materials & methods: Three different structures of calmodulin are compared, obtaining a total of 2.5 μs of molecular dynamics, which have been analyzed using the principal component analysis and clustering methodologies. Results: sMD simulations reach conformations that conventional molecular dynamics is not able to, without compromising the stability of the protein. On the other hand, accelerated MD requires optimization of the setup parameters to be useful. Conclusion: sMD is useful to study flexible proteins, highlighting those factors that justify its promiscuity.

Hit-to-lead virtual screening frequently relies on a cascade of computational methods that starts with rapid calculations applied to a large number of compounds and ends with more expensive computations restricted to a subset of compounds that passed initial filters. This work focuses on set up protocols for alchemical free energy (AFE) scoring in the context of a Docking-MM/PBSA-AFE cascade. A dataset of 15 congeneric inhibitors of the ACK1 protein was used to evaluate the performance of AFE set up protocols that varied in the steps taken to prepare input files (using previously docked and best scored poses, manual selection of poses, manual placement of binding site water molecules). The main finding is that use of knowledge derived from X-ray structures to model binding modes, together with the manual placement of a bridging water molecule, improves the R2 from 0.45 ± 0.06 to 0.76 ± 0.02 and decreases the mean unsigned error from 2.11 ± 0.08 to 1.24 ± 0.04 kcal mol-1. By contrast a brute force automated protocol that increased the sampling time ten-fold lead to little improvements in accuracy. Besides, it is shown that for the present dataset hysteresis can be used to flag poses that need further attention even without prior knowledge of experimental binding affinities., Competing Interests: The authors have declared that no competing interests exist.

Computational and experimental studies unravel the structural and electronic properties of a novel supramolecular liquid crystal built through a hierarchical assembly process resulting in an H-bonded melamine rosette decorated with peripheral triphenylenes. The six-fold symmetry of the mesogen facilitates the formation of a highly organized hexagonal columnar mesophase stable at room temperature. X-ray diffraction and electron density maps confirm additional intra- and intercolumn segregation of functional subunits, and this paves the way for 1D charge transport. Indeed, hole mobility has been measured and found to be higher than for related mesogens. DFT calculations of HOMO and LUMO levels and parameters such as reorganization energy and transfer integral of the rosette structure have been achieved, and not only validate the columnar organization but also establish the way it translates into a favorable electronic architecture and molecular orbital interactions to promote charge carrier mobility., (© 2018 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Aim: Rescoring of docking-binding poses can significantly improve molecular docking results. Our aim was to evaluate postprocessing docking protocols in order to determine the most suitable methodology for the study of the binding of congeneric compounds to protein kinases., Materials & Methods: Diverse ligand-receptor poses generated after docking were submitted to different relaxation protocols. The Molecular Mechanics Poisson-Boltzmann (Generalized Born) Surface Area approach was applied for the evaluation of the binding affinity of complexes obtained. The performance of various Molecular Mechanics Poisson-Boltzmann (Generalized Born) Surface Area methodologies was compared., Results: The inclusion of a postprocessing protocol after docking enhances the quality of the results, although the best methodology is system dependent., Conclusion: An examination of the interactions established has allowed us to suggest useful modifications for the design of new type II inhibitors.

BI-2536 is a potent Polo-like kinase inhibitor which induces apoptosis in diverse human cancer cell lines. The binding affinity of BI-2536 for human serum albumin (HSA) protein may define its pharmacokinetic and pharmacodynamic profile. We have studied the binding of BI-2536 to HSA by means of different spectroscopic techniques and docking calculations. We have experimentally observed that the affinity of BI-2536 for HSA is higher than that of other common HSA binding drugs. Therefore, it can be postulated that the drug dose should be increased to achieve a certain concentration of free drug in plasma, although BI-2536 could also reach tumour tissues by uptaking HSA/BI-2536 complex. Only a single binding site on HSA has been observed for BI-2536 which seems to correspond to the subdomain IIA pocket. The formation of the HSA/BI-2536 complex is a spontaneous and entropy-driven process that does not cause a significant change of the secondary structure of the protein. Its endothermic character could be related to proton release. Thermodynamic analysis showed that the main protein-drug interactions are of the van der Waals type although the presence of amide and ether groups in BI-2536 could also allow H-bonding with some residues in the subdomain IIA pocket., (Copyright © 2017. Published by Elsevier B.V.)

Density functional theory calculations were carried out to investigate the evolvement of charge transport properties of a set of new discotic systems as a function of ring and heteroatom (B, Si, S, and Se) substitution on the basic structure of perylene. The replacement of six-membered rings by five-membered rings in the reference compound has shown a prominent effect on the electron reorganization energy that decreases ∼0.2 eV from perylene to the new carbon five-membered ring derivative. Heteroatom substitution with boron also revealed to lower the LUMO energy level and increase the electron affinity, therefore lowering the electron injection barrier compared to perylene. Since the rate of the charge transfer between two molecules in columnar discotic systems is strongly dependent on the orientation of the stacked cores, the total energy and transfer integral of a dimer as a disc is rotated with respect to the other along the stacking axis have been predicted. Aimed at obtaining a more realistic approach to the bulk structure, the molecular geometry of clusters made up of five discs was fully optimized, and charge transfer rate and mobilities were estimated for charge transport along a one dimensional pathway. Heteroatom substitution with selenium yields electron transfer integral values ∼0.3 eV with a relative disc orientation of 25°, which is the preferred angle according to the dimer energy profile. All the results indicate that the tetraselenium-substituted derivative, not synthetized so far, could be a promising candidate among those studied in this work for the fabrication of n-type semiconductors based on columnar discotic liquid crystals materials.

The effect of fluorine substitution on the molecular structure, crystal packing, and n-type semiconducting properties of a set of poly(arylene-ethynylene) polymers based on alternating thiadiazole and phenyl units linked through ethynylene groups has been studied by means of Density Functional Theory. As a result, an enlargement in the interplanar distance between cofacial polymer chains, as well as a decrease of the electronic coupling and electron mobility is predicted. On the other hand, fluorination could facilitate electron injection into the material. A polymer containing both alkoxy pendant chains and fluorine atoms is proposed as a compromise solution between efficiency of electron injection and charge transport within the material.

We present a density functional theory (DFT) study on charge-transport related properties in a series of discotic systems based on 1,3,5-triazine and tris[1,2,4]triazolo[1,3,5]triazine central cores as electron acceptor units, and phenyl-thiophene and N-carbazolyl-thiophene segments as electron donor units. The presence of both electron donor and acceptor moieties in the π-conjugated core could lead to new discotic liquid crystal (DLC) materials which are predicted to display ambipolar charge transport behavior in such a way that electrons could move through the central part of the next cores while holes mainly do through the peripheral groups. A significant increase in hole mobility when N-carbazolyl is present as an electron donor unit in the peripheral region is predicted. In addition, a detailed topological analysis of the electron charge density within the framework provided by Quantum Theory of Atoms in Molecules (QTAIM) has been performed in order to characterize intra- and intermolecular interactions in terms of hydrogen bonds and/or π···π stacking which contribute to the stabilization of the columnar stack and the helical self-assembly at the molecular scale.

B-Raf mutations are identified in 40-50% of patients with melanoma and among them, the substitution of valine for glutamic acid at position 600 ((V600E)B-Raf) is the most frequent. Treatment of these patients with B-Raf inhibitors has been associated with a clear clinical benefit. Unfortunately, multiple resistance mechanisms have been identified and new potent and selective inhibitors are currently needed. In this work, five different type II inhibitors, which bind (V600E)B-Raf in its DFG-out conformation, have been studied using molecular dynamics, free energy calculations and energy decomposition analysis. The ranking of calculated MM-PB/GBSA binding affinities is in good agreement with the experimentally measured ones. The per-residue decomposition of ΔGbinding, within the MM-GBSA approach, has been used to identify the key residues governing the allosteric binding of the studied compounds to the (V600E)B-Raf protein kinase. Results indicate that although van der Waals interactions are key determinants for binding, hydrogen bonds also play an important role. This work also provides a better structural understanding of the binding of DFG-out inhibitors to (V600E)B-Raf, which can be used in a further step for rational design of a new class of B-Raf potent inhibitors.

Targeting the BCL6 protein is a promising therapeutic strategy for the treatment of B cell lymphomas. One approach to treat these diseases consists of finding drug candidates able to disrupt the interactions established between BCL6 and its corepressors. Thus, this work presents a thorough comparative analysis of the interactions between the BCL6 BTB (bric-a-brac tramtrack broad complex) protein domain and its SMRT, NcoR and BCOR corepressor BBDs (BCL6 binding domain) through molecular dynamics. Moreover, a theoretical structure is presented and checked for the BCL6(BTB)-NcoR(BBD) complex. Considering the BBDs to be composed of 17 amino acids, our analyses show the region involving residues 4-15 of these 17 to play a main role in the protein-corepressor interactions. Particularly SER(11) seems to have a high relevance as it establishes specific bonds with BCL6(BTB) and is one of the only two residues sequence equivalent for the three studied corepressors. From this study, 14 pharmacophoric points have been proposed divided in two groups which coincide with residues 4-11 and 11-15, being SER(11) a hinge point. This finding suggests the possibility of searching for 2 small molecule inhibitors, mimicking 8 and 7 pharmacophoric points, respectively, which could incorporate a hydrogen donor pharmacophoric point mimicking SER(11) in any or both molecules. In short, the present work aims to contribute further knowledge in the modeling of drugs mimicking BCL6(BTB)-corepressor complexes., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Aimed to optimize the ratio accuracy/computational cost, in this work we study the performance of three different theoretical methodologies in the calculation of the optical bandgap for a test set made of a number of poly(aryl-ethynylene)s related polymers. Infinite, ideal polymer chains were first optimized by means of periodic calculations. Different length oligomers were afterward generated by direct replication of the corresponding periodic structure and their optical bandgaps were calculated by means of different time dependent-density functional theory (TD-DFT) methodologies. These results were fitted to an exponential function for each oligomer family in order to get a theoretical estimation of the optical bandgap for each polymer to be compared to the experimental reported values. The best result was obtained for TD-M06-2X yielding an average deviation of 3.4% with respect to the experimental values.

This work presents a Density Functional Theory (DFT) study on the charge transport related properties of two quinoidal diketopyrrolopyrrole (DPP) based systems. System A, recently synthesized, shows high efficiency as n-type organic semiconductor material while system B, not synthesized yet, has a linking benzothiadiazole (BT) unit between DPP moieties and would display an ambipolar character. The use of tuned, long-range corrected (LRC) functionals allows one to predict HOMO, LUMO, and charge transport properties for compound A in concordance with those experimentally observed. The use of BT building blocks allows for a conclusion that compound B is expected to display balanced and efficient charge injection along with high mobilities both for holes and electrons, which points to its potential to obtain high performances as an ambipolar semiconductor.

In the present study, a series of crystalline poly(arylene-ethynylene) copolymers containing phenylethynylene and 2,5-dialkoxy-phenylethynylene units together with 1,3,4-thiadiazole rings has been modeled by means of periodic calculations. Optimized three-dimensional polymeric structures show interchain distances that are consistent with the experimental values reported for a related polymer. It has also been observed that the presence of pendant alkoxy chains brings on both a further flattening and a separation of the coplanar chains. This fact is linked to a decrease of the interchain cofacial distance. The electron transport character of the polymer crystal structures was assessed through Marcus theory. Electronic coupling between neighboring polymer chains is most influenced by the presence of alkoxy chains giving rise to an expectable enhancement of the electron hopping mobility.

Insight into the electronic structure of disordered poly-2,5-bis(phenylethynyl)-1,3,4-thiadiazole in an amorphous region, in comparison to an ideal two-planar cofacial oligomer system, is pursued. The atomic structure of the amorphous polymer was obtained from classical molecular dynamics. It was subsequently used to calculate the electronic states and inter- and intrachain electronic coupling integrals using the density functional theory based charge patching method. The interchain electronic coupling integrals in the amorphous system were found to be an order of magnitude smaller than in the ordered system with similar distances between the chains. The results also suggest that the electronic structure of the whole system cannot be understood as a collection of the electronic structures of individual chains. The band gap of the whole system is significantly smaller than the band gaps of individual chains. This decrease originates from the disordered long range electrostatic potential created by the dipole moments of polymer repeat units, which should be minimized if one seeks good transport properties.

In the present work, we have studied from a theoretical perspective the geometry and electronic properties of the series of related compounds 2,5-bis(phenylethynyl)-1,3,4-thiadiazole, 2,5-bis(phenylethynyl)-1,3,4-oxadiazole, and 2,5-bis(phenylethynyl)-1,2,4-triazole as candidates for electron-conducting polymers and compounds with desirable (opto)electronic properties. The effect of the ethynyl group (-C[Triple Bond]C-) on the structure and electronic properties was also studied. The influence of planarity on electrical conductivity has been studied by a natural-bond-orbital analysis. The (opto)electronic properties and conducting capability were investigated through the highest occupied molecular orbital (HOMO)-lowest unoccupied molecular orbital (LUMO) gap, excitation energy, bond length alternation, LUMO energy, electron affinities, and intramolecular reorganization energy. Finally, the evolution of some properties such as optical bandgap and electron affinity with the increase of the number of repeat units in the oligomer chain has been checked.

2,5-Bis(phenylethynyl)-1,3,4-thiadiazole (PhEtTh) and 2,5-diphenyl-1,3,4-thiadiazole (PhTh) are expected to be building blocks for polymer materials that could be employed to conduct electricity due to their narrow highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) energy gaps. In this work, a theoretical, comparative study about the effect of the ethynyl group on the planarity and electrical conductivity of this kind of systems has been carried out. Thus, several ab initio (Hartree-Fock, Moller-Plesset) and DFT (B3LYP, B3PW91, M05, M05-2X) methods and basis sets (6-31G(*), 6-31G+G(**), 6-311G(**), cc-pVDZ, cc-pVTZ) have been tested. As a result, PhEtTh showed better properties for its use as electric conducting material relative to PhTh due to its smaller HOMO-LUMO gap, as well as its enhanced trend to retain the planarity provided the reduction in steric hindrances that the ethynyl group (-C[triple bond]C-) permits. Solvent effects were also modeled for ethanol and chloroform under the conductor-like polarizable continuum model approximation. Finally, electronic transitions in gas and solution phases were predicted by using TDDFT approximation in order to compare the theoretical lambda(max) with the experimental values reported in literature for both compounds.

We have explored a series of trisubstituted acridine-peptide conjugates for their ability to recognize and discriminate between DNA quadruplexes derived from the human telomere, and the c-kit and N-ras proto-oncogenes. Quadruplex affinity was measured as the peptide sequences were varied, together with their substitution position on the acridine, and the identity of the C-terminus (acid or amide). Surface plasmon resonance measurements revealed that all compounds bound to the human telomeric quadruplex with sub-micromolar affinity. Docking calculations from molecular modelling studies were used to model the effects of substituent orientation and peptide sequence. Modelling and experiment were in agreement that placement of the peptide over the face of the acridine is detrimental to binding affinity. The highest degrees of selectivity were observed towards the N-ras quadruplex by compounds capable of forming simultaneous contacts with their acridine and peptide moieties. The ligands that bound best displayed quadruplex affinities in the 1-5 nM range and at least 10-fold discrimination between the quadruplexes studied.

Apoptosis, also called programmed cell death, is a conserved mechanism inherent to all cells that sentences them to death when they receive the appropriate external stimuli. Inhibitor of apoptosis proteins (IAPs) are a family of regulatory proteins that suppress such cell death. XIAP is the most commonly studied member of the IAP family. It binds to and inhibits Caspases, an important family of apoptotic proteases. In addition, XIAP over-expression has been detected in numerous types of cancer. Smac/DIABLO, a mitochondrial protein that binds to IAPs and promotes Caspase activation, has the opposite action to XIAP and can be considered a key protein in the regulation of IAPs. Survivin, the smallest IAP protein, has received a lot of attention due to its specific expression in many cancer cell lines. It has been shown to interact with Smac/DIABLO, even though the structure of this complex has not yet been reported.We analysed the protein-protein interactions appearing in the Smac/DIABLO-XIAP and Smac/DIABLO-Survivin complexes fully, using molecular dynamics simulations. This information is a first step towards the design of Smac/DIABLO peptidomimetics that could be used as innovative therapeutic agents for the treatment of malignancy. Our results complement the experimental interactions described for the first complex and provide a detailed description for the second. We show that Smac/DIABLO interacts in a similar way with both targets through its amino terminal residues. In addition, we identify a pharmacophore formed by eight stable protein-protein interactions for the XIAP complex and seven stable protein-protein interactions for the Survivin complex, which describe the whole contact surface. This information is used to suggest the binding mode of embelin, the first non-peptidic inhibitor of XIAP, and two of its derivatives. Molecular docking and molecular dynamics simulations were also carried out to describe ligand and receptor flexibility. Finally, an MMGBSA protocol was used to obtain a more quantitative description of the binding in all the complexes studied., (2008 John Wiley & Sons, Ltd)