Optimizing potentiometric ionophore and electrode design for environmental on-site control of antibiotic drugs: Application to sulfamethoxazole

Potentiometric sensors are typically unable to carry out on-site monitoring of environmental drug contaminants because of their high limits of detection (LODs). Designing a novel ligand material for the target analyte and managing the composition of the internal reference solution have been the strategies employed here to produce for the first time a potentiometric-based direct reading method for an environmental drug contaminant. This concept has been applied to sulfamethoxazole (SMX), one of the many antibiotics used in aquaculture practices that may occur in environmental waters. The novel ligand has been produced by imprinting SMX on the surface of graphitic carbon nanostructures (CN) < 500 nm. The imprinted carbon nanostructures (ICN) were dispersed in plasticizer and entrapped in a PVC matrix that included (or not) a small amount of a lipophilic additive. The membrane composition was optimized on solid-contact electrodes, allowing near-Nernstian responses down to 5.2 μg/mL and detecting 1.6 μg/mL. The membranes offered good selectivity against most of the ionic compounds in environmental water. The best membrane cocktail was applied on the smaller end of a 1000 μL micropipette tip made of polypropylene. The tip was then filled with inner reference solution containing SMX and chlorate (as interfering compound). The corresponding concentrations were studied for 1 × 10−5 to 1 × 10−10 and 1 × 10−3 to 1 × 10−8 mol/L. The best condition allowed the detection of 5.92 ng/L (or 2.3 × 10−8 mol/L) SMX for a sub-Nernstian slope of −40.3 mV/decade from 5.0 × 10−8 to 2.4 × 10−5 mol/L.

The design, synthesis and characterization of new building blocks for the preparation of molecule-based magnetic materials /

Two new families of building blocks have been prepared and fully characterized and their coordination chemistry exploited for the preparation of molecule-based magnetic materials. The first class of compounds were prepared by exploiting the chemistry of 3,3'-diamino-2,2'-bipyridine together with 2-pyridine carbonyl chloride or 2-pyridine aldehyde. Two new ligands, 2,2'-bipyridine-3,3'-[2-pyridinecarboxamide] (Li, 2.3) and N'-6/s(2-pyridylmethyl) [2,2'bipyridine]-3,3'-diimine (L2, 2.7), were prepared and characterized. For ligand L4, two copper(II) coordination compounds were isolated with stoichiometrics [Cu2(Li)(hfac)2] (2.4) and [Cu(Li)Cl2] (2.5). The molecular structures of both complexes were determined by X-ray crystallography. In both complexes the ligand is in the dianionic form and coordinates the divalent Cu(II) ions via one amido and two pyridine nitrogen donor atoms. In (2.4), the coordination geometry around both Cu11 ions is best described as distorted trigonal bipyramidal where the remaining two coordination sites are satisfied by hfac counterions. In (2.5), both Cu(II) ions adopt a (4+1) distorted square pyramidal geometry. One copper forms a longer apical bond to an adjacent carbonyl oxygen atom, whereas the second copper is chelated to a neighboring Cu-Cl chloride ion to afford chloride bridged linear [Cu2(Li)Cl2]2 tetramers that run along the c-axis of the unit cell. The magnetic susceptibility data for (2.4) reveal the occurrence of weak antiferromagnetic interactions between the copper(II) ions. In contrast, variable temperature magnetic susceptibility measurements for (2.5) reveal more complex magnetic properties with the presence of ferromagnetic exchange between the central dimeric pair of copper atoms and weak antiferromagnetic exchange between the outer pairs of copper atoms. The Schiff-base bis-imine ligand (L2, 2.7) was found to be highly reactive; single crystals grown from dry methanol afforded compound (2.14) for which two methanol molecules had added across the imine double bond. The susceptibility of this ligand to nucleophilic attack at its imine functionality assisted via chelation to Lewis acidic metal ions adds an interesting dimension to its coordination chemistry. In this respect, a Co(II) quaterpyridine-type complex was prepared via a one-pot transformation of ligand L2 in the presence of a Lewis acidic metal salt. The rearranged complex was characterized by X-ray crystallography and a reaction mechanism for its formation has been proposed. Three additional rearranged complexes (2.13), (2.17) and (2.19) were also isolated when ligand (L2, 2.7) was reacted with transition metal ions. The molecular structures of all three complexes have been determined by X-ray crystallography. The second class of compounds that are reported in this thesis, are the two diacetyl pyridine derivatives, 4-pyridyl-2,6-diacetylpyridine (5.5) and 2,2'-6,6'-tetraacetyl-4,4'-bipyridine (5.15). Both of these compounds have been designed as intermediates for the metal templated assembly of a Schiff-base N3O2 macrocycle. From compound (5.15), a covalently tethered dimeric Mn(II) macrocyclic compound of general formula {[Mn^C^XJCl-FkO^Cl-lO.SFbO (5.16) was prepared and characterized. The X-ray analysis of (5.16) reveals that the two manganese ions assume a pentagonal-bipyramidal geometry with the macrocycle occupying the pentagonal plane and the axial positions being filled by a halide ion and a H2O molecule. Magnetic susceptibility data reveal the occurrence of antiferromagnetic interactions between covalently tethered Mn(II)-Mn(II) dimeric units. Following this methodology a Co(II) analogue (5.17) has also been prepared which is isostructural with (5.16).

Molecular Magnetism: The Design, Synthesis, and Characterization of New Building Blocks for Molecule-based Magnetic Materials

Two classes of building blocks have been prepared and characterized and their coordination chemistry explored working towards the preparation of new molecule-based magnetic materials. In the first project, the amine functionality of 3,3'-diamino-2,2'- bipyridine was exploited for the preparation of a new family of ligands (H2L 1)-(H2L 4). The molecular structures of three ligands have been fully characterized by X-ray crystallography. [molecular structure diagram will not copy here, but is available in full pdf.] The coordination chemistry of these ligands with divalent first row transition metal ions was investigated. For ligand (H2L1), the molecular structures of four coordination complexes with stoichiometries [Zn2(Ll)(OAc)(MeO)]2 (I), [Cu2(L1)(OAc)2 (II), [Li(L1)]3 (III), and [Ni(L1)]3 (IV) were determined by X-ray crystallography. For ligand (H2L2), a Cu(II) complex of stoichiometry [Cu3(L2)(OAc)3MeO] (V) was determined by X-ray crystallography. The magnetic properties of complexes (II), (III), and (V) have been fully elucidated. In project two, synthetic strategies for the preparation of porphyrin molecules bearing triol substituents is presented. Following this approach, three new porphyrin derivatives have been prepared and characterized [Zn(HPTPP-CH2C(CH20H)3)] (VI), [P(TPP)(OCH2C(CH2)H)3)2]+CL- (VII), and [P(OEP)(C6H5)(OCH2C(CH2OH)3)]+Cl- (VIII). Attempts to exchange the labile methoxide bridges of a tetraironIIl single molecule magnet of stoichiometry [Fe4(OMe)6(dpm)6] (Hdpm = dipivaloylmethane) with the triol appended porphyrins will be discussed. [molecular structure diagram will not copy here, but is available in full pdf.]

Triarylamminium radical-cation complexes: design and magnetic properties and Base-catalyzed hydrosilylation: mechanism and substrate scope

1. Triarylamminium radical-cation complexes. The detailed study of manganese, copper and nickel metal-radical complexes with triarylamminium ligands was conducted. Stable, neutral and pseudo-octahedral coordination monometallic complexes with simple monodentate 2,2`-bipyridine ligand containing a redox-active N,N`-(4,4`-dimethoxydiphenyl-amino) substituent were synthesized and fully characterized. The one-electron oxidation process and formation of persistent radical-cation complexes was observed by cyclic voltammetry and spectroelectrochemical measurements. Evans method measurements were performed with radical-cation complexes generated by chemical one-electron oxidation with NOPF6 in acetonitrile. The experimental results indicate ferromagnetic coupling between metal and triarylamminium cation in manganese (II) complex and antiferromagnetic coupling in nickel (II) complex. This data is supported by DFT calculations which also lend weight to the  spin polarization mechanism as an operative model for magnetic exchange coupling. Neutral bimetallic complexes with a new ditopic ligand were synthesized and fully characterized, including magnetic and electrochemical studies. Chemical oxidation of these precursor complexes did not generate radical-cations, but dicationic complexes, which was confirmed by UV-vis and EPR-experiments, as well as varied temperature magnetic measurements. DFT calculations for radical-cation complexes are included. A synthetic pathway for polytopic ligand with multiple redox-active triarylamine sites was developed. The structure of the ligand is presumably suitable for -spin polarization exchange model and allows for production of polymetallic complexes having high spin ground states. 2. Base-catalyzed hydrosilylation. A simple reductive base-catalyzed hydrosilation of aldehydes and ketones was adapted to the use of the cheap, safe, and non-toxic polymethylhydrosiloxane (PMHS) instead of the common PhSiH3 and (EtO)3SiH, which present significant cost and safety concerns, respectively. The conversion of silane into pentacoordinate silicate species upon addition of a base was studied in details for the cases of phenyl silane and PMHS and is believed to be essential for the hydrosilylation process. We discovered that nucleophiles (a base or fluoride-anion) induced the rearrangement of PMHS and TMDS into light silanes: MeSiH3 and Me2SiH2, respectively. The reductive properties of PMHS under basic conditions can be attributed to the formation of methyl silane and its conversion into a silicate species. A procedure for the generation of methyl silane and its use in further efficient reductions of aldehydes and ketones has been developed. The protocol was extended to the selective reduction of esters and tertiary amides into alcohols and aldimines into amines with good isolated yields and reduction of heterocyclic compounds was attempted.

Use of cellular impedance to characterize ligand functional selectivity at G protein-coupled receptors

Les récepteurs couplés aux protéines G (RCPGs) représentent la plus grande famille de cibles thérapeutiques pour le traitement d’une panoplie de pathologies humaines. Bien que plusieurs décennies de recherche aient permis de façonner nos connaissances sur ces protéines membranaires, notre compréhension des déterminants moléculaires de leur activité signalétique reste encore limitée. De ces domaines de recherche, une avancée récente a mis à jour un nouveau phénomène, appelé sélectivité fonctionnelle des ligands, qui a bouleversé les paradigmes décrivant leu fonctionnement de ces récepteurs. Ce concept émane d’observations montrant que l’activité pharmacologique de certains ligands n’est pas nécessairement conservée sur tout le répertoire signalétiques connu du récepteur et peu se restreindre à l'activation sélective d’un sous-groupe de voies de signalisation.Ce nouveau modèle pharmacologique de l'activation des RCPG ouvre de nouvelles possibilités pour la découverte de médicaments plus efficace et sûr, ciblant les RCPGs. En effet, il permet la conception de molécules modulant spécifiquement les voies signalétiques d’intérêt thérapeutique, sans engager les autres voies qui pourraient mener à des effets secondaires indésirables ou de la tolérance. Cette thèse décrit l'utilisation d'une nouvelle approche sans marquage, basée sur la mesure du changement l'impédance cellulaire. Par la mesure des changements cellulaires, comme la morphologie, l’adhésion et/ou la redistribution des macromolécules, cette approche permet de mesurer de façon simultanée l'activité de plusieurs voies de signalisation impliqués dans ces réponses. Utilisant le récepteur β2-adrénergique (β2AR) comme modèle, nous avons démontré que les variations dans l’impédance cellulaire étaient directement liées à l’activation de multiples voies de signalisation suite à la stimulation du récepteur par son ligand. L’agoniste type du β2AR, l’isoprotérénol, s’est avéré induire une réponse d’impédance dose-dépendante constituée, dans le temps, de plusieurs caractéristiques distinctes pouvant être bloquées de façon compétitive par l’antagoniste ICI118,551 Par l’utilisation d’inhibiteurs sélectifs, nous avons été en mesure de déterminer la contribution de plusieurs voies signalétiques canoniques, comme les voies dépendantes de Gs et Gi, la production d’AMPc et l’activation de ERK1/2, sur ces changements. De plus, la dissection de la réponse d’impédance a permis d’identifier une nouvelle voie de mobilisation du Ca2+ contribuant à la réponse globale des changements initiés par la stimulation du β2AR. Dans une autre étude, nous avons rapporté que la réponse calcique induite par le β2AR serait attribuable à une transactivation Gs-dépendant du récepteur purinergique P2Y11, lui-même couplé à la protéine Gq. La mesure d’impédance permettant de distinguer et de décrire une pléiade d’activités signalétiques, nous avons émis l’hypothèse que des ligands arborant des profils signalétiques différents généreraient des réponses d’impédance distinctes. Le criblage d’une librairie de ligands spécifiques au β2AR a révélé une grande variété de signatures d’impédance. Grâce au développement d’une approche computationnelle innovatrice, nous avons été en mesure de regrouper ces signatures en cinq classes de composés, un regroupement qui s’est avéré hautement corrélé avec le profil signalétique des différents ligands. Nous avons ensuite combiné le criblage de composés par impédance avec l’utilisation d’inhibiteurs sélectifs de voies signalétiques afin d’augmenter la résolution du regroupement. En évaluant l’impact d’une voie signalétique donnée sur la signature d’impédance, nous avons été en mesure de révéler une plus grande variété de textures parmi les ligands. De plus, cette méthode s’est avérée efficace pour prédire le profil signalétique d’une librairie de composés non caractérisés, ciblant le β2AR. Ces travaux ont mené à l’élaboration d’une méthode permettant d’exprimer visuellement la sélectivité fonctionnelle de ligands et ont révélé de nouvelles classes de composés pour ce récepteur. Ces nouvelles classes de composés ont ensuite été testées sur des cardiomyocytes humains, confirmant que les composés regroupés dans différentes classes produisent des effets distincts sur la contractilité de ces cellules. Globalement, ces travaux démontrent la pertinence de l’utilisation de l’impédance cellulaire pour une évaluation précise des différences fonctionnelles parmi les composés ciblant les RCPGs. En fournissant une représentation pluridimensionnelle de la signalisation émanant des RCPGs à l’aide d’un seul essai ne requérant pas de marquage, les signatures d’impédance représentent une stratégie simple et innovante pour l’évaluation de la fonctionnalité sélective des ligands. Cette méthode pourrait être d’une grande utilité dans le processus de découverte de nouveaux médicaments.

A novel Bayesian approach to drug design with simple and complex enzymes: Use with lens aldehyde dehydrogenase and the glyoxalase pathway

Purpose: Acquiring details of kinetic parameters of enzymes is crucial to biochemical understanding, drug development, and clinical diagnosis in ocular diseases. The correct design of an experiment is critical to collecting data suitable for analysis, modelling and deriving the correct information. As classical design methods are not targeted to the more complex kinetics being frequently studied, attention is needed to estimate parameters of such models with low variance. Methods: We have developed Bayesian utility functions to minimise kinetic parameter variance involving differentiation of model expressions and matrix inversion. These have been applied to the simple kinetics of the enzymes in the glyoxalase pathway (of importance in posttranslational modification of proteins in cataract), and the complex kinetics of lens aldehyde dehydrogenase (also of relevance to cataract). Results: Our successful application of Bayesian statistics has allowed us to identify a set of rules for designing optimum kinetic experiments iteratively. Most importantly, the distribution of points in the range is critical; it is not simply a matter of even or multiple increases. At least 60 % must be below the KM (or plural if more than one dissociation constant) and 40% above. This choice halves the variance found using a simple even spread across the range.With both the glyoxalase system and lens aldehyde dehydrogenase we have significantly improved the variance of kinetic parameter estimation while reducing the number and costs of experiments. Conclusions: We have developed an optimal and iterative method for selecting features of design such as substrate range, number of measurements and choice of intermediate points. Our novel approach minimises parameter error and costs, and maximises experimental efficiency. It is applicable to many areas of ocular drug design, including receptor-ligand binding and immunoglobulin binding, and should be an important tool in ocular drug discovery.

Design and theoretical study of nickel catalysts for syndiotactic polyolefins

A nickel catalyst was modeled with ligand L-2, [ NH = CH-CH = CH-O](-), which should have potential use as a syndiotactic polyolefin catalyst, and the reaction mechanism was studied by theoretical calculations using the density functional method at the B3LYP/ LANL2MB level. The mechanism involves the formation of the intermediate [(NiLMe)-Me-2](+), in which the metal occuples a T-shaped geometry. - This intermediate has two possible structures with the methyl group trans either to the oxygen or to the nitrogen atom of L-2. The results show that both structures can lead to the desired product via similar reaction paths, A and B. Thus, the polymerization could be considered as taking place either with the alkyl group occupying the position trans to the Ni-O or trans to the Ni-N bond in the catalyst. The polymerization process thus favors the catalysis of syndiotactic polyolefins. The syndiotactic synthesis effects could also be enhanced by variations in the ligand substituents. From energy considerations, we can conclude that it is more favorable for the methyl group to occupy the trans-O position to form a complex than to occupy the trans-N position. From bond length considerations, it is also more favoured for ethene to occupy the trans-O position than to occupy the trans-N position.

Molecular modelling studies of binding of DACD derivatives into G-Quadruplex DNA: comparison of force field and quantum polarized ligand docking methods

The DNA G-qadruplexes are one of the targets being actively explored for anti-cancer therapy by inhibiting them through small molecules. This computational study was conducted to predict the binding strengths and orientations of a set of novel dimethyl-amino-ethyl-acridine (DACA) analogues that are designed and synthesized in our laboratory, but did not diffract in Synchrotron light.Thecrystal structure of DNA G-Quadruplex(TGGGGT)4(PDB: 1O0K) was used as target for their binding properties in our studies.We used both the force field (FF) and QM/MM derived atomic charge schemes simultaneously for comparing the predictions of drug binding modes and their energetics. This study evaluates the comparative performance of fixed point charge based Glide XP docking and the quantum polarized ligand docking schemes. These results will provide insights on the effects of including or ignoring the drug-receptor interfacial polarization events in molecular docking simulations, which in turn, will aid the rational selection of computational methods at different levels of theory in future drug design programs. Plenty of molecular modelling tools and methods currently exist for modelling drug-receptor or protein-protein, or DNA-protein interactionssat different levels of complexities.Yet, the capasity of such tools to describevarious physico-chemical propertiesmore accuratelyis the next step ahead in currentresearch.Especially, the usage of most accurate methods in quantum mechanics(QM) is severely restricted by theirtedious nature. Though the usage of massively parallel super computing environments resulted in a tremendous improvement in molecular mechanics (MM) calculations like molecular dynamics,they are still capable of dealing with only a couple of tens to hundreds of atoms for QM methods. One such efficient strategy that utilizes thepowers of both MM and QM are the QM/MM hybrid methods. Lately, attempts have been directed towards the goal of deploying several different QM methods for betterment of force field based simulations, but with practical restrictions in place. One of such methods utilizes the inclusion of charge polarization events at the drug-receptor interface, that is not explicitly present in the MM FF.

In silico identification and characterization of protein-ligand binding sites

Protein–ligand binding site prediction methods aim to predict, from amino acid sequence, protein–ligand interactions, putative ligands, and ligand binding site residues using either sequence information, structural information, or a combination of both. In silico characterization of protein–ligand interactions has become extremely important to help determine a protein’s functionality, as in vivo-based functional elucidation is unable to keep pace with the current growth of sequence databases. Additionally, in vitro biochemical functional elucidation is time-consuming, costly, and may not be feasible for large-scale analysis, such as drug discovery. Thus, in silico prediction of protein–ligand interactions must be utilized to aid in functional elucidation. Here, we briefly discuss protein function prediction, prediction of protein–ligand interactions, the Critical Assessment of Techniques for Protein Structure Prediction (CASP) and the Continuous Automated EvaluatiOn (CAMEO) competitions, along with their role in shaping the field. We also discuss, in detail, our cutting-edge web-server method, FunFOLD for the structurally informed prediction of protein–ligand interactions. Furthermore, we provide a step-by-step guide on using the FunFOLD web server and FunFOLD3 downloadable application, along with some real world examples, where the FunFOLD methods have been used to aid functional elucidation.

Identification of a novel ligand binding motif in the transthyretin channel

The design of therapeutic compounds targeting transthyretin (TTR) is challenging due to the low specificity of interaction in the hormone binding site. Such feature is highlighted by the interactions of TTR with diclofenac, a compound with high affinity for TTR, in two dissimilar modes, as evidenced by crystal structure of the complex. We report here structural analysis of the interactions of TTR with two small molecules, 1-amino-5-naphthalene sulfonate (1,5-AmNS) and 1-anilino-8-naphthalene sulfonate (1,8-ANS). Crystal structure of TTR: 1,8-ANS complex reveals a peculiar interaction, through the stacking of the naphthalene ring between the side-chain of Lys15 and Leu17. The sulfonate moiety provides additional interaction with Lys15` and a water-mediated hydrogen bond with Thr119`. The uniqueness of this mode of ligand recognition is corroborated by the crystal structure of TTR in complex with the weak analogue 1,5-AmNS, the binding of which is driven mainly by hydrophobic partition and one electrostatic interaction between the sulfonate group and the Lys15. The ligand binding motif unraveled by 1,8-ANS may open new possibilities to treat TTR amyloid diseases by the elucidation of novel candidates for a more specific pharmacophoric pattern. (C) 2009 Published by Elsevier Ltd.

Differential effects of TR ligands on hormone dissociation rates: Evidence for multiple ligand entry/exit pathways

Some nuclear receptor (NR) ligands promote dissociation of radiolabeled bound hormone from the buried ligand binding cavity (LBC) more rapidly than excess unlabeled hormone itself This result was interpreted to mean that challenger ligands bind allosteric sites on the LBD to induce hormone dissociation, and recent findings indicate that ligands bind weakly to multiple sites on the LBD surface. Here we show, that a large fraction of thyroid hormone receptor (TR) ligands promote rapid dissociation (T(1/2) < 2 h) of , radiolabeled T(3) vs. T(3) (T(1/2), approximate to 5-7 h). We cannot discern relationships between this effect and ligand size, activity or affinity for TR beta. One ligand, GC-24, binds the TR LBC and (weakly) to the TR beta-LBD surface that mediates dimer/heterodimer interaction, but we cannot link this interaction to rapid T(3) dissociation. Instead, several lines of evidence suggest that the challenger ligand must interact with the buried LBC to promote rapid T(3) release. Since previous molecular dynamics simulations suggest that TR ligands leave the LBC by several routes, we propose that a subset of challenger ligands binds and stabilizes a partially unfolded intermediate state of TR that arises during T(3) release and that this effect enhances hormone dissociation. (C) 2009 Elsevier Ltd. All rights reserved.

Ligand induced interaction of thyroid hormone receptor beta with its coregulators

Thyroid hormones exert most of their physiological effects through two thyroid hormone receptor (TR) subtypes, TR alpha and TR beta, which associate with many transcriptional coregulators to mediate activation or repression of target genes. The search for selective TR beta ligands has been stimulated by the finding that several pharmacological actions mediated by TR beta might be beneficial in medical conditions such as obesity, hypercholesterolemia and diabetes. Here, we present a new methodology which employs surface plasmon resonance to investigate the interactions between TR beta ligand binding domain (LBD) complexes and peptides derived from the nuclear receptor interaction motifs of two of its coregulators, SRC2 and DAX1. The effect of several TR beta ligands, including the TR beta selective agonist GC-I and the TR beta selective antagonist NH-3, were investigated. We also determined the kinetic rate constants for the interaction of TR beta-T3 with both coregulators, and accessed the thermodynamic parameters for the interaction with DAX1. Our findings Suggest that flexibility plays an important role in the interaction between the receptor and its coregulators. and point out important aspects of experimental design that should be addressed when using TR beta LBD and its agonists. Furthermore, the methodology described here may be useful for the identification of new TR beta ligands. (C) 2008 Elsevier Ltd. All rights reserved.

Only subtle protein conformational adaptations are required for ligand binding to thyroid hormone receptors: Simulations using a novel multipoint steered molecular dynamics approach

Thyroid hormone receptors (TR) are hormone-dependent transcription regulators that play a major role in human health, development, and metabolic functions. The thyroid hormone resistance syndrome, diabetes, obesity, and some types of cancer are just a few examples of important diseases that are related to TR malfunctioning, particularly impaired hormone binding. Ligand binding to and dissociation from the receptor ultimately control gene transcription and, thus, detailed knowledge of binding and release mechanisms are fundamental for the comprehension of the receptor`s biological function and development of pharmaceuticals. In this work, we present the first computational study of ligand entry into the ligand binding domain (LBD) of a nuclear receptor. We report molecular dynamics simulations of ligand binding to TRs using a generalization of the steered molecular dynamics technique designed to perform single-molecule pulling simulations along arbitrarily nonlinear driving pathways. We show that only gentle protein movements and conformational adaptations are required for ligand entry into the LBDs and that the magnitude of the forces applied to assist ligand binding are of the order of the forces involved in ligand dissociation. Our simulations suggest an alternative view for the mechanisms ligand binding and dissociation of ligands from nuclear receptors in which ligands can simply diffuse through the protein surface to reach proper positioning within the binding pocket. The proposed picture indicates that the large-amplitude protein motions suggested by the apo- and holo-RXR alpha crystallographic structures are not required, reconciling conformational changes of LBDs required for ligand entry with other nuclear receptors apo-structures that resemble the ligand-bound LBDs.

Ligand dissociation from estrogen receptor is mediated by receptor dimerization: Evidence from molecular dynamics simulations

Estrogen Receptor (ER) is an important target for pharmaceutical design. Like other ligand-dependent transcription factors, hormone binding regulates ER transcriptional activity. Nevertheless, the mechanisms by which ligands enter and leave ERs and other nuclear receptors remain poorly understood. Here, we report results of locally enhanced sampling molecular dynamics simulations to identify dissociation pathways of two ER ligands [the natural hormone 17 beta-estradiol (E-2) and the selective ER modulator raloxifene (RAL)] from the human ER alpha ligand-binding domain in monomeric and dimeric forms. E-2 dissociation occurs via three different pathways in ER monomers. One resembles the mousetrap mechanism (Path I), involving repositioning of helix 12 (H12), others involve the separation of H8 and H11 (Path II), and a variant of this pathway at the bottom of the ligand-binding domain (Path II`). RAL leaves the receptor through Path I and a Path I variant in which the ligand leaves the receptor through the loop region between H11 and H12 (Path I`). Remarkably, ER dimerization strongly suppresses Paths II and II` for E-2 dissociation and modifies RAL escape routes. We propose that differences in ligand release pathways detected in the simulations for ER monomers and dimers provide an explanation for previously observed effects of ER quaternary state on ligand dissociation rates and suggest that dimerization may play an important, and hitherto unexpected, role in regulation of ligand dissociation rates throughout the nuclear receptor family.

Natural products biological screening and ligand-based virtual screening for the discovery of new antileishmanial agents

In the course of our research program to discover novel antileishmanial agents, a biological screening of natural products against Leishmania major promastigotes allowed the identification of a furoquinoline alkaloid (1) and a furanocoumarin (2) as new hits. Subsequently, an integrated ligand-based virtual screening approach was employed to search for new antileishmanial compounds using these naturally occurring molecules as templates. Fourteen out of 40 compounds selected from a database of about 800,000 compounds (extracted from ZINC, a free database for virtual screening) were experimentally confirmed to possess significant in vitro antileishmanial properties. The application of ligand-based virtual screening as a complementary approach to experimental natural product screening was a useful strategy to facilitate the identification of new promising lead candidates.