Density functional theory study of mixed-phase TiO2: heterostructures and electronic properties

In this work, density functional theory calculations have been performed to study the geometric, electronic, and energetic properties of two-phase TiO2 composites built by joining two single-phase TiO2 slabs, aiming at verifying possible improvement of the photo-activities of the composites through phase separation of excitons. We find that such desired electronic properties can be determined by several factors. When both the HOMO and LUMO levels of one of the two single-phase TiO2 slabs are higher than the corresponding ones of the other, the composite may have native electronic structures with phase-separated HOMO-LUMO states, especially when the two slabs exhibit highly matched surface lattices. For those pairs of TiO2 slabs with the HOMO and LUMO levels of one phase being within the range of those of the other, though the energetically favored composite give HOMO-LUMO states within one phase, one may still be able to separate them and move the HOMO state to the interface region by destabilizing the interactions between the two slabs.

Activation of human TLR4/MD-2 by hypoacylated lipopolysaccharide from a clinical isolate of Burkholderia cenocepacia

Lung infection by Burkholderia species, in particular B. cenocepacia, accelerates tissue damage and increase post-lung transplant mortality in cystic fibrosis patients. Host- microbes interplay largely depends on interactions between pathogen specific molecules and innate immune receptors such as the Toll-like receptor 4 (TLR4), which recognizes the lipid A moiety of the bacterial lipopolysaccharide (LPS). The human TLR4/MD-2 LPS receptor complex is strongly activated by hexa-acylated lipid A and poorly activated by underacylated lipid A. Here, we report that B. cenocepacia LPS strongly activates human TLR4/MD-2 despite its lipid A having only five acyl chains. Further, we show that aminoarabinose residues in lipid A contribute to TLR4-lipid A interactions, and experiments in a mouse model of LPS-induced endotoxic shock confirmed the pro- inflammatory potential of B. cenocepacia penta-acylated lipid A. Molecular modeling, combined with mutagenesis of TLR4-MD2 interactive surfaces, suggests that longer acyl chains and the aminoarabinose residues in the B. cenocepacia lipid A allow exposure of the fifth acyl chain on the surface of MD-2 enabling interactions with TLR4 and its dimerization. Our results provide a molecular model for activation of the human TLR4/MD- 2 complex by penta-acylated lipid A, explaining the ability of hypoacylated B. cenocepacia LPS to promote pro- inflammatory responses associated to the severe pathogenicity of this opportunistic bacterium.

Piezoelectricity and ferroelectricity in amino acid glycine

Bioorganic ferroelectrics and piezoelectrics are becoming increasingly important in view of their intrinsic compatibility with biological environment and biofunctionality combined with strong piezoelectric effect and switchable polarization at room temperature. Here we study piezoelectricity and ferroelectricity in the smallest amino acid glycine, representing a broad class of non-centrosymmetric amino acids. Glycine is one of the basic and important elements in biology, as it serves as a building block for proteins. Three polymorphic forms with different physical properties are possible in glycine (α, β and γ), Of special interest for various applications are non-centrosymmetric polymorphs: β-glycine and γ-glycine. The most useful β-polymorph being ferroelectric took much less attention than the other due to its instability under ambient conditions. In this work, we could grow stable microcrystals of β-glycine by the evaporation of aqueous solution on a (111)Pt/Ti/SiO2/Si substrate as a template. The effects of the solution concentration and Pt-assisted nucleation on the crystal growth and phase evolution were characterized by X-ray diffraction analysis and Raman spectroscopy. In addition, spin-coating technique was used for the fabrication of highly aligned nano-islands of β-glycine with regular orientation of the crystallographic axes relative the underlying substrate (Pt). Further we study both as-grown and tip-induced domain structures and polarization switching in the β-glycine molecular systems by Piezoresponse Force Microscopy (PFM) and compare the results with molecular modeling and computer simulations. We show that β-glycine is indeed a room-temperature ferroelectric and polarization can be switched by applying a bias to non-polar cuts via a conducting tip of atomic force microscope (AFM). Dynamics of these in-plane domains is studied as a function of applied voltage and pulse duration. The domain shape is dictated by both internal and external polarization screening mediated by defects and topographic features. Thermodynamic theory is applied to explain the domain propagation induced by the AFM tip. Our findings suggest that β-glycine is a uniaxial ferroelectric with the properties controlled by the charged domain walls which in turn can be manipulated by external bias. Besides, nonlinear optical properties of β-glycine were investigated by a second harmonic generation (SHG) method. SHG method confirmed that the 2-fold symmetry is preserved in as-grown crystals, thus reflecting the expected P21 symmetry of the β-phase. Spontaneous polarization direction is found to be parallel to the monoclinic [010] axis and directed along the crystal length. These data are confirmed by computational molecular modeling. Optical measurements revealed also relatively high values of the nonlinear optical susceptibility (50% greater than in the z-cut quartz). The potential of using stable β-glycine crystals in various applications are discussed in this work.

A minimal version of a Protein Tyrosine Phosphatase

Phosphatase and tensin homologue (PTEN) protein belongs to the family of protein tyrosine phos-phatase. Mutations on the phosphatase and tensin homologue (PTEN) protein are highly observed in diverse types of human tumors, being mostly identified on the phosphatase domain of the protein. Although PTEN is a modular protein composed by a phosphatase domain and a C2 domain for mem-brane anchoring, this work aimed at developing a minimal version of PTEN´s phosphatase domain. The minimal version (Small Domain) comprises a 28 residue peptide, with the PTEN 8-mer catalytic peptide accommodated between a α-helix and β-turn as observed in PTEN native structure. Firstly, a de novo prediction of the Small Domain´s secondary structure was carried out by molecular modeling tools. The stability of the predicted structures were then evaluated by Molecular Dynamics. Automated molecular docking of PTEN natural substrate PIP3, its analogue (Inositol) and a PTEN inhibitor (L-tar-tare) were performed with the modeled structure, and PTEN used as a positive control. The gene en-coding for Small Domain was designed and cloned into an expression vector at N-terminal of Green Fluorescence Protein (GFP) encoding gene. The fusion protein was then expressed in Escherichia coli cells. Different expression conditions have been explored for the production of the fusion protein to minimize the formation of inclusion bodies.

Structure of the extracellular domains of human and Xenopus Fn14: implications in the evolution of TWEAK and Fn14 interactions.

TWEAK (TNF homologue with weak apoptosis-inducing activity) and Fn14 (fibroblast growth factor-inducible protein 14) are members of the tumor necrosis factor (TNF) ligand and receptor super-families. Having observed that Xenopus Fn14 cross-reacts with human TWEAK, despite its relatively low sequence homology to human Fn14, we examined the conservation in tertiary fold and binding interfaces between the two species. Our results, combining NMR solution structure determination, binding assays, extensive site-directed mutagenesis and molecular modeling, reveal that, in addition to the known and previously characterized β-hairpin motif, the helix-loop-helix motif makes an essential contribution to the receptor/ligand binding interface. We further discuss the insight provided by the structural analyses regarding how the cysteine-rich domains of the TNF receptor super-family may have evolved over time. DATABASE: Structural data are available in the Protein Data Bank/BioMagResBank databases under the accession codes 2KMZ, 2KN0 and 2KN1 and 17237, 17247 and 17252. STRUCTURED DIGITAL ABSTRACT: TWEAK binds to hFn14 by surface plasmon resonance (View interaction) xeFn14 binds to TWEAK by enzyme linked immunosorbent assay (View interaction) TWEAK binds to xeFn14 by surface plasmon resonance (View interaction) hFn14 binds to TWEAK by enzyme linked immunosorbent assay (View interaction).

Intramolecular carbenoid insertions : the reactions of [alpha]-diazoketones derived from furanyl, thienyl, benzofuranyl and benzothienyl acetic acids with rhodium (II) acetate /

A number of synthetically useful ring systems can be prepared via the intramolecular insertion of a metal-stabilized carbenoid into a heteroaromatic systems. The chemical outcome of these reactions are dependent not only on the nature of the heteroatom but also on the length of the aliphatic tether linking the carbenoid moiety with the aromatic fragment. Our work with furanyl and thienyl systems containing a single methylene tether have allowed for some rather atypical chemistry. For example, treatment of l-diazo-3-(2-thienyl)-2-propanone (6) with catalytic rhodium (II) acetate yields 5,6- dihydro-4^-cyclopenta[Z>]thiophen-5-one (3) while, the isomeric l-diazo-3-(3-thienyl)-2- propanone(15) gives a spiro-disulphide (20). Novel chemistry was also exhibited in the analogous furanyl systems. While treatment of l-diazo-3-(3-furanyl)-2-propanone (52) with Rh2(OAc)4 resulted in the expected 2-(4-Oxo-2-cyclopentenyliden)acetaldehyde (54), isomeric l-diazo-3-(2- furanyl)-2-propanone (8) undergoes vinylogous Wolff rearrangement to give a mixture of 6a-methyl-2,3,3a,6a-tetrahydrofuro[2,i-^>]furan-2-one (44) and 2-(2-methyl-3-furyl)acetic acid (43). Rhodium acetate catalyzed decomposition of l-diazo-3-(3-benzofuranyl)-2- propanone (84) and l-diazo-3-(2-benzofuranyl)-2-propanone (69)also allows for vinylogous Wolff rearrangement, a chemistry unseen in benzofuranyl systems with longer tethers. A number of interesting products were isolated from the trapping of intermediate ketenes. Decomposition of l-diazo-3-(3-benzothienyl)-2-propanone (100) resulted in the formation of 2,3-dihydro-l//-benzo[^]cyclopenta[^thiophen-2-one (102). However, in addition to (102), a dimer was also generated from the decomposition of l-diazo-3-(2- benzothienyl)-2-propanone (109). The insight into the mechanistic underpinnings of the above reactions are provided by molecular modeling at a PM3 level.

A comparative study of protoheme and heme d catalases : role of the heme and the heme pocket in catalysis and ligand binding

Catalase dismutes H20 2 to O2 and H20. In successive twoelectron reactions H20 2 induces both oxidation and reduction at the heme group. In the first step the protoheme prosthetic group of beef liver catalase forms compound I, in which the heme has been oxidized from Fe3+ to Fe4+=0 and a porphyrin radical has been created. Compound II is formed by the oneelectron reduction of comp I. It retains Fe4+=0 but lacks the porphyrin radical and is catalytically inert. Molecular structures are available for Escherichia coli Hydroperoxidase II, Micrococcus Iysodeiktus, Penicillium vitale and beef liver enzymes, which contain different hemes and heme pockets. In the present work, the pockets and substrate access channels of protoheme (beef liver & Micrococcus) and heme d (HPII of E. coli and Penicillium) catalases have been analysed using Quanta™ and CharmMTM molecular modeling packages on the Silicon Graphics Iris Indigo 2 computer. Experimental studies have been carried out with two catalases, HPII (and its mutants) and beef liver. Fluoride and formate' are inhibitors of both enzymes, and their binding is modulated by the heme and by distal residues N201 & H128. Both HPII and beef liver enzymes form compound I with H202 or peracetate. The reduction of beef liver enzyme compound I to II and the decay of compound II are accelerated by fluoride. The decay of compound II is also accelerated by formate, and this reagent acts as a 2-electron donor towards compound I of both enzymes. It is concluded that heme d enzymes (Penicillium and HPII of E. coli) are formed by autocatalytic transformation of protoheme in a modified pocket which contains a characteristic serine residue as well as a partially occluded heme channel. They are less active than protoheme enzymes but also do not form the inactive compound II species. Binding of peroxide as well as fluoride and formate is prevented by mutation of H128 and modulated by mutation of N201.

Stereoselective synthesis of substituted hexahydro-3a,4a-diazacyclopentaphenanthren-4-ones and aminoferrocenes

This thesis explored the development of several methodologies for the stereoselective construction of ligand frameworks and some of their applications. The first segment concerns the application of an enantioselective lithiation at an Sp3_ hybridized position adjacent to nitrogen by means of the widely used and typically highly effective enantioselective lithiation with ( -)-sparteine. This investigation was intended to develop a method to install chirality into a system that would be converted into a family of diaminoylidenes for use as phosphine mimics in transition metal catalysis or as nucleophilic reagents. Molecular modeling of the system revealed some key interactions between the substrate and (-)-sparteine that provided general insight into the diamine's mode of action and should lend some predictive value to its future applications. The second portion focuses on the development of methods to access 1,2- disubstituted aminoferrocenes, an underexplored class of metallocenes possessing planar chirality. Two routes were examined involving a diastereoselective and an enantioselective pathway, where the latter method made use of the first BF3-mediated lithiation-substitution to install planar chirality. Key derivatives such as 1,2- aminophosphines, made readily accessible by the new route, were evaluated as ligands for Pd(II), Pt(II) and Ir(I). These complexes show activity in a number of transformations with both achiral and prochiral substrates. Optimization experiments were conducted to prepare enantiomerically enriched 2-substituted-I-aminoferrocenes by direct asymmetric lithiation of BF3-coordinated tertiary aminoferrocenes. A predictive computational model describing the transition state of this reaction was developed in collaboration with Professor Travis Dudding's group (Department of Chemistry, Brock University). The predicted stereochemistry of the process was confirmed by single-crystal X-ray analysis of a 2-phosphino-l-dimethylaminoferrocene derivative. Enantiomerically pure samples of the aminophosphine ligands derived from this new process have given promising preliminary results in the enantioselective hydrogenation of prochiral alkenes and warrant further stUdy in metal-mediated catalysis.

Directed evolution of human dihydrofolate reductase: towards a better understanding of binding at the active site

La dihydrofolate réductase humaine (DHFRh) est une enzyme essentielle à la prolifération cellulaire, ce qui en fait une cible de choix pour le traitement de différents cancers. À cet effet, plusieurs inhibiteurs spécifiques de la DHFRh, les antifolates, ont été mis au point : le méthotrexate (MTX) et le pemetrexed (PMTX) en sont de bons exemples. Malgré l’efficacité clinique certaine de ces antifolates, le développement de nouveaux traitements s’avère nécessaire afin de réduire les effets secondaires liés à leur utilisation. Enfin, dans l’optique d’orienter la synthèse de nouveaux composés inhibiteurs des DHFRh, une meilleure connaissance des interactions entre les antifolates et leur enzyme cible est primordiale. À l’aide de l’évolution dirigée, il a été possible d’identifier des mutants de la DHFRh pour lesquels l’affinité envers des antifolates cliniquement actifs se voyait modifiée. La mutagenèse dite ¬¬de saturation a été utilisée afin de générer des banques de mutants présentant une diversité génétique au niveau des résidus du site actif de l’enzyme d’intérêt. De plus, une nouvelle méthode de criblage a été mise au point, laquelle s’est avérée efficace pour départager les mutations ayant entrainé une résistance aux antifolates et/ou un maintient de l’activité enzymatique envers son substrat natif, soient les phénotypes d’activité. La méthode de criblage consiste dans un premier temps en une sélection bactérienne à haut débit, puis dans un second temps en un criblage sur plaques permettant d’identifier les meilleurs candidats. Plusieurs mutants actifs de la DHFRh, résistants aux antifolates, ont ainsi pu être identifiés et caractérisés lors d’études de cinétique enzymatique (kcat et IC50). Sur la base de ces résultats cinétiques, de la modélisation moléculaire et des données structurales de la littérature, une étude structure-activité a été effectuée. En regardant quelles mutations ont les effets les plus significatif sur la liaison, nous avons commencé à construire un carte moléculaire des contacts impliqués dans la liaison des ligands. Enfin, des connaissances supplémentaires sur les propriétés spécifiques de liaison ont put être acquises en variant l’inhibiteur testé, permettant ainsi une meilleure compréhension du phénomène de discrimination du ligand.

Mécanismes moléculaires d’activation du récepteur A des peptides natriurétiques

Le récepteur A des peptides natriurétiques (NPRA) fait partie de la famille des guanylates cyclases membranaires. L’activation du NPRA par ses agonistes naturels, ANP et BNP, induit une production de GMPc qui est responsable de leur rôle dans l’homéostasie cardiovasculaire, l’inhibition de l’hypertrophie et de la fibrose cardiaques et la régulation de la lipolyse. Le NPRA est un homodimère non covalent composé d’un domaine extracellulaire de liaison du ligand (ECD), d’un unique domaine transmembranaire (TM), d’un domaine d’homologie aux kinases et d’un domaine guanylate cyclase. Bien que le NPRA ait un rôle physiologique important, les mécanismes moléculaires régissant son processus d’activation restent inconnus. Nous avons donc analysé les premières étapes du processus d’activation du NPRA. Nous avons d'abord étudié le rôle de la dimérisation des ECD dans l’activation du récepteur. Nous avons utilisé les techniques de liaison de radioligand, de FRET et de modélisation moléculaire, pour caractériser la liaison à l’ECD des agonistes naturels, d’un superagoniste et d’un antagoniste. L’ANP se lie à un dimère d’ECD préformé et la dimérisation spontanée est l’étape limitante du processus de liaison. De plus, comme le démontrent nos études de FRET, tous les peptides, incluant l’antagoniste, stabilisent le récepteur sous sa forme dimérique. Cependant, l’antagoniste A71915 stabilise le dimère d’ECD dans une conformation différente de celle induite par l’ANP. La dimérisation du NPRA semble donc nécessaire, mais non suffisante à l’activation du récepteur. L’état d’activation du NPRA dépend plutôt de l’orientation des sous unités dans le dimère. Nous avons ensuite étudié le mécanisme moléculaire de transduction du signal à travers la membrane. Plusieurs études ont suggéré que l’activation du NPRA implique un changement de conformation du domaine juxtamembranaire (JM). Cependant, les études de cristallographie de l’ECD soluble de NPRA n’ont pas permis de documenter la structure du JM et le changement de conformation impliqué dans la transduction du signal reste inconnu. Pour analyser ce changement de conformation, nous avons d’abord séquentiellement substitué les neuf acides aminés du JM par une cystéine. En étudiant la capacité des mutants à former des dimères covalents de façon constitutive ou induite par l’ANP, nous avons pu évaluer la proximité relative des résidus du JM, avant et après activation du NPRA. Ces résultats ont démontré la proximité élevée de certains résidus spécifiques et sont en contradiction avec les données cristallographiques. Nous avons également démontré que le domaine intracellulaire impose une contrainte conformationnelle au JM à l’état de base, qui est levée après liaison de l’ANP. En introduisant de 1 à 5 alanines dans l’hélice-α transmembranaire, nous avons montré qu’une rotation des TM de 40° induit une activation constitutive du NPRA. Le signal d’activation pourrait donc être transmis à travers la membrane par un mécanisme de rotation des TM. En utilisant nos données expérimentales, nous avons généré le premier modèle moléculaire illustrant la conformation active du NPRA, où les domaines JM et TM sont représentés. Dans son ensemble, cette étude apporte une meilleure compréhension des mécanismes moléculaires régissant les premières étapes du processus complexe d’activation du NPRA.

Étude des interactions supramoléculaires par modélisation moléculaire

L’avancée des infrastructures informatiques a permis l’émergence de la modélisation moléculaire. À cet effet, une multitude de modèles mathématiques sont aujourd’hui disponibles pour simuler différents systèmes chimiques. À l’aide de la modélisation moléculaire, différents types d’interactions chimiques ont été observés. À partir des systèmes les plus simples permettant l’utilisation de modèles quantiques rigoureux, une série d’approximations a été considérée pour rendre envisageable la simulation de systèmes moléculaires de plus en plus complexes. En premier lieu, la théorie de la fonctionnelle de densité dépendante du temps a été utilisée pour simuler les énergies d’excitation de molécules photoactives. De manière similaire, la DFT indépendante du temps a permis la simulation du pont hydrogène intramoléculaire de structures analogues au 1,3,5-triazapentadiène et la rationalisation de la stabilité des états de transition. Par la suite, la dynamique moléculaire et la mécanique moléculaire ont permis de simuler les interactions d’un trimère d’acide cholique et d’un pyrène dans différents solvants. Cette même méthodologie a été utilisée pour simuler les interactions d’un rotaxane-parapluie à l’interface d’un système biphasique. Finalement, l’arrimage moléculaire et les fonctions de score ont été utilisés pour simuler les interactions intermoléculaires entre une protéine et des milliers de candidats moléculaires. Les résultats ont permis de mettre en place une stratégie de développement d’un nouvel inhibiteur enzymatique.

Étude structure/fonction des cotransporteurs Na+/glucose

Cette thèse porte sur l’étude de la relation entre la structure et la fonction chez les cotransporteurs Na+/glucose (SGLTs). Les SGLTs sont des protéines membranaires qui se servent du gradient électrochimique transmembranaire du Na+ afin d’accumuler leurs substrats dans la cellule. Une mise en contexte présentera d’abord un bref résumé des connaissances actuelles dans le domaine, suivi par un survol des différentes techniques expérimentales utilisées dans le cadre de mes travaux. Ces travaux peuvent être divisés en trois projets. Un premier projet a porté sur les bases structurelles de la perméation de l’eau au travers des SGLTs. En utilisant à la fois des techniques de modélisation moléculaire, mais aussi la volumétrie en voltage imposé, nous avons identifié les bases structurelles de cette perméation. Ainsi, nous avons pu identifier in silico la présence d’une voie de perméation passive à l’eau traversant le cotransporteur, pour ensuite corroborer ces résultats à l’aide de mesures faites sur le cotransporteur Na/glucose humain (hSGLT1) exprimé dans les ovocytes. Un second projet a permis d’élucider certaines caractéristiques structurelles de hSGLT1 de par l’utilisation de la dipicrylamine (DPA), un accepteur de fluorescence dont la répartition dans la membrane lipidique dépend du potentiel membranaire. L’utilisation de la DPA, conjuguée aux techniques de fluorescence en voltage imposé et de FRET (fluorescence resonance energy transfer), a permis de démontrer la position extracellulaire d’une partie de la boucle 12-13 et le fait que hSGLT1 forme des dimères dont les sous-unités sont unies par un pont disulfure. Un dernier projet a eu pour but de caractériser les courants stationnaires et pré-stationaires d’un membre de la famille des SGLTs, soit le cotransporteur Na+/myo-inositol humain hSMIT2 afin de proposer un modèle cinétique qui décrit son fonctionnement. Nous avons démontré que la phlorizine inhibe mal les courants préstationnaires suite à une dépolarisation, et la présence de courants de fuite qui varient en fonction du temps, du potentiel membranaire et des substrats. Un algorithme de recuit simulé a été mis au point afin de permettre la détermination objective de la connectivité et des différents paramètres associés à la modélisation cinétique.

Algorithm for generating defective graphene sheets

An algorithm is presented for the generation of molecular models of defective graphene fragments, containing a majority of 6-membered rings with a small number of 5- and 7-membered rings as defects. The structures are generated from an initial random array of points in 2D space, which are then subject to Delaunay triangulation. The dual of the triangulation forms a Voronoi tessellation of polygons with a range of ring sizes. An iterative cycle of refinement, involving deletion and addition of points followed by further triangulation, is performed until the user-defined criteria for the number of defects are met. The array of points and connectivities are then converted to a molecular structure and subject to geometry optimization using a standard molecular modeling package to generate final atomic coordinates. On the basis of molecular mechanics with minimization, this automated method can generate structures, which conform to user-supplied criteria and avoid the potential bias associated with the manual building of structures. One application of the algorithm is the generation of structures for the evaluation of the reactivity of different defect sites. Ab initio electronic structure calculations on a representative structure indicate preferential fluorination close to 5-ring defects.

Nanosized polymetallic resorcinarene-based host assemblies that strongly bind fullerenes

Polymetallic nanodimensional assemblies have been prepared via metal directed assembly of dithiocarbamate functionalized cavitand structural frameworks with late transition metals (Ni, Pd, Cu, Au, Zn, and Cd). The coordination geometry about the metal centers is shown to dictate the architecture adopted. X-ray crystallographic studies confirm that square planar coordination geometries result in "cagelike" octanuclear complexes, whereas square-based pyramidal metal geometries favor hexanuclear "molecular loop" structures. Both classes of complex are sterically and electronically complementary to the fullerenes (C-60 and C-70). The strong binding of these guests occurred via favorable interactions with the sulfur atoms of multiple dithiocarbamate moieties of the hosts. In the case of the tetrameric copper(II) complexes, the lability of the copper(II)-dithiocarbamate bond enabled the fullerene guests to be encapsulated in the electron-rich cavity of the host, over time. The examination of the binding of fullerenes has been undertaken using spectroscopic and electrochemical methods, electrospray mass spectrometry, and molecular modeling.

Bis(calix[4]diquinone) receptors: cesium- and rubidium-selective redox-active ionophores

A new class of redox-active ionophore comprised of two calix[4]diquinone moieties connected through either alkylene or pyridylene linkages has been developed. Spectroscopic and electrochemical investigations, X-ray crystal structure analyses, and molecular modeling studies show butylene- and propylene-linked members of this family of redox-active receptors exhibit remarkable selectivity preferences and substantial electrochemical recognition effects toward cesium and rubidium cations.