Pramod Patil
Kunturu
a,
Özlem
Kap
b,
Kai
Sotthewes
b,
Pierre
Cazade
c,
Harold J. W.
Zandvliet
*b,
Damien
Thompson
*c,
Ofer
Reany
*d and
Jurriaan
Huskens
*a
aMolecular Nanofabrication Group, MESA + Institute for Nanotechnology, Faculty of Science and Technology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands. E-mail: j.huskens@utwente.nl
bPhysics of Interfaces and Nanomaterials, MESA+ Institute for Nanotechnology, University of Twente, P. O. Box, 217, 7500AE Enschede, The Netherlands. E-mail: h.j.w.zandvliet@utwente.nl
cDepartment of Physics, Bernal Institute, University of Limerick, Limerick V94 T9PX, Ireland. E-mail: damien.thompson@ul.ie
dDepartment of Natural Sciences, The Open University of Israel, 1 University Road, Ra'anana 43107, Israel. E-mail: oferre@openu.ac.il
First published on 17th December 2019
Semithio-bambusurils are a unique family of anion-binding host macrocycles that form self-assembled monolayers (SAMs) on Au(111). SAMs of semithio-bambus[n]uril homologs with different cage sizes (1: n = 4; 2: n = 6) have been investigated using electrochemistry, X-ray photoelectron spectroscopy (XPS), scanning tunneling microscopy (STM), and molecular dynamics (MD) simulations. Electrochemical measurements showed that electron transfer occurs via tunneling through the SAMs, and the low resistivity indicated an open layer architecture. XPS confirmed that thiocarbonyl sulfur atoms are chemisorbed to the Au(111) surface, and STM revealed the formation of ordered domains in a rectangular lattice for 1 and a highly ordered triangular/hexagonal lattice for 2. MD simulations substantiated the STM data by quantifying the balance between molecule–surface bonding, molecular conformations, and supramolecular packing that drive the formation of SAMs that maximize their surface coverage within the limits of conformational strain.
Design, System, ApplicationThe attachment of receptors to a transducer surface provides a general strategy to design sensor devices. Yet, these receptor molecules may undergo large conformational changes that may remain obscure when not inspecting the molecular structure of the receptor in its surface-immobilized state. The current manuscript provides an example of a class of receptor molecules with thiocarbonyl groups purposefully integrated in the molecular design in order to provide anion recognition capabilities. Upon binding of the receptors to a gold surface, which can be used as an electrode for electrochemical sensing, the molecules appear to undergo large conformational changes that optimize their layer formation properties but at the expense of their anion binding. The lessons learned in this detailed surface-analytical and molecular dynamics study may provide insights into the effect of surface binding to the molecular conformations of the receptors and may assist in developing new strategies to promote receptor binding with retention of recognition. |
Bambus[n]urils (BUs, n = 4,6) are a family of cavitands composed of 2,4-disubstituted glycoluril subunits interconnected by methylene bridges in an alternate arrangement.24,25 The glycoluril subunits in the BU structure present their convex side towards the interior with the methine hydrogens pointing inward, forming a highly electron-deficient cavity with both polarization and electrostatic interactions contributing to anion binding. Following the predictions that replacement of oxygen atoms in bambusurils by other heteroatoms could significantly modify their anion binding properties, the calculated molecular electrostatic potentials of hetero-BUs were found to follow a general trend of anion-affinity: S > O > NH.26 These predictions were also confirmed experimentally indicating that, although all analogs are excellent anion binders, semithio-BU[6] is an effective transmembrane transporter capable of polarizing lipid membranes through the uniport mechanism of chloride anions.27 Similarly, semiaza-BU[6]s exhibit the simultaneous accommodation of three anions, linearly positioned along the main symmetry axis, which is reminiscent of natural chloride channels in E. coli.28
Here, we report the SAM formation ability of both semithio-BU homologs (n = 4, 6), 1 and 2, respectively (Fig. 1, Fig. S1†)29 on a gold surface. The semithio-BU SAMs were characterized by their contact angles, X-ray photoelectron spectroscopy (XPS), and electrochemistry. Anion binding was studied using electrochemical impedance spectroscopy (EIS) and cyclic voltammetry (CV). We used scanning tunneling microscopy (STM) combined with molecular dynamics (MD) simulations to determine the geometry of the semithio-BU SAMs in order to understand their anchoring on gold, intermolecular packing and anion recognition properties.
Fig. 1 Molecular structures (left) of semithio-BU[4], 1 (top) and semithio-BU[6], 2 (bottom), and the corresponding side and top views of the X-ray structures.29 Hydrogen atoms and solvent molecules were omitted for clarity (unit cells with solvent shown in Fig. S1†). Color code: oxygen, red; sulfur, yellow; nitrogen, blue; carbon, grey. |
The contact angles show almost no dependence on preparation temperature in the range of 25 to 80 °C (Table S1†), indicating that the temperature has no noticeable effect on the SAM formation. The cleaned bare gold layer is more hydrophobic than the SAMs, which indicates that both semithio-BU homologs have been successfully adsorbed onto the Au surface. The thiocarbonyl (CS) moieties within both portals of 1 and 2 are expected to have a high affinity for Au, which drives SAM formation on the Au surface. Regarding SAM stability, monitoring the contact angle of the SAMs over the course of a week (58° and 62° for 1 and 2, respectively) did not show significant changes.
Infrared reflection absorption spectroscopy (IRRAS) showed three characteristic peaks of 1 and 2 (Fig. S2,† top and bottom, respectively): at 1260 cm−1, 1060 cm−1 and 950 cm−1 corresponding to (C)–(N), (C)(S) and (C)–(N–CS) stretching vibrations, respectively.
XPS measurements (Fig. 2a and b) clearly identified carbon, sulfur, nitrogen, and oxygen in the SAMs (Fig. 2 and S3†). The C1s signal can be deconvoluted into three components at 284.7, 286.8 and 289.1 eV, which are assigned to C–C, CS and CO bonds, respectively (Fig. 2c and d). There is qualitative agreement between the XPS spectra and elemental composition for the SAM and powder samples of 1 and 2 (Table S2†), but a more detailed comparison is complicated by the presence of solvent and/or salts in the powders (see Fig. S4† and the discussion following it). XPS was also used to determine the average number of sulfur atoms that are involved in binding the macrocycles to the gold surface. For sulfur-containing SAMs, e.g., thiols and sulfides, it is known that adsorption on gold results in a negative shift of about 1.5 eV in the XPS signal of the S2p electrons. Hence, with bound sulfur 2p3/2 to the gold surface the XPS signal is observed at 161.9 eV compared to unbound sulfur 2p3/2 observed at 163.4 eV.10,12 Both XPS spectra of SAMs of 1 and 2 in the S2p region showed broad peaks with a binding energy (BE) of sulfur (S2p3/2) located around ∼162 eV (Fig. 2e and f). An additional set of signals, with much lower intensity, is visible at around 168 eV as well, indicative of oxidized S. In particular, no detectable intensity is present in the BE region at ∼164 eV for any of these two SAMs, which means that unbound sulfur atoms (S2p1/2 BE of ca. 165 eV) are not detected by XPS. The most intense doublet at 162 eV is attributed to surface-bound Au–S species. The smaller doublets, at 168 eV, are tentatively attributed to oxidized sulfur in 1, while we note that SAMs of 2 are fully unoxidized. Potentially, some of the sulfur atoms may get oxidized by Au oxide, present after piranha treatment, or by ambient oxygen,30 but evidence of oxidation is also seen in the powder spectra (Fig. S4†). The peak separation in the doublets is about 1.2 eV and the intensity ratio is 2:1, corresponding to the S 2p3/2 and S 2p1/2 lines of semithio-BUs. These data confirm that the majority or all of the sulfur atoms are connected to the gold surface upon SAM formation, with a minor fraction (approximately 15%) in SAMs of 1 oxidized which prohibits their precise assignment as free or bound species.
Electrochemical measurements, employing the redox behavior of the redox couple Fe(CN)63−/4– supplied to solution, are commonly used to probe the thickness and compactness of a monolayer on gold. Previous studies have shown that the diffusion-limited current, measured during cyclic voltammetry experiments at a monolayer-covered electrode, depends on the presence of pinholes in the monolayer.32,33 The efficiency of the monolayer in blocking redox reactions can be estimated by the kinetics of the redox probe, as expressed by the peak-to-peak separation ΔEp. CV in a 0.5 M K2SO4 electrolyte was used to determine the capacitance C of the monolayers (Table 1). The electrode adsorbs ions, resulting in an electrochemical double layer, which provides a capacitance to the surface. The capacitance of the monolayer electrical double layer (Cdl) is related to the effective thickness and order of the SAM on the gold substrate10 and can be calculated from the measured current density.
CV measurements (Fig. 3a) showed that the capacitance decreases upon SAM formation, with capacitances reduced from 11.6 μF cm−2 for bare gold to 7.67 μF cm−2 for 1 and 6.07 μF cm−2 for 2 (Table 1). All tested samples (1 and 2) showed a slightly enhanced resistance compared to bare gold samples. Overall, the relatively low resistance (Fig. S5†) of the SAMs of 1 and 2 and the concomitantly relatively high capacitance in the presence of redox couple (Fig. 3b) indicate a rather thin and open SAM structure, possibly with a high number of pinhole defects present in the monolayer. This electrochemically open structure is expected from the macrocyclic structure of these hosts and the absence of well-packing long alkyl chains in these molecules.
In line with earlier data obtained in solution,26,29 we attempted to use the SAMs as anion binders, and therefore we tried to probe the electrochemical performance at various concentrations of KBr solutions (Fig. S6†). Surprisingly, although semithio-BU[6] is a known strong bromide anion binder, addition of KBr solution (in the range of 1 mM to 30 mM) did not lead to detectable changes in the electrochemical parameters of the SAM of 2. Therefore, we assumed, also based on the XPS results which indicate involvement of most of the sulfur atoms in binding, that the attachment and assembly of 1 and 2 on gold was accompanied by large changes in their conformations relative to their conformational state in solution. Recently, Reany and co-workers reported that perthio-BUs exhibit a yet unknown compact conformation, which explains their diminished anion binding properties but still bind to thiophilic metals.31 Hence, we investigated the changes in the conformation of the current semithio-BUs on gold by scanning tunneling microscopy (STM) and molecular modeling.
Fig. 4a depicts an STM image of a SAM of 1 on Au(111). Several vacancy islands are present at the Au(111) surface. The vacancy islands have a depth of 2.5 Å, i.e., one Au atomic layer, indicating that these vacancy islands are vacancies in the Au surface and not in the semithio-bambusuril SAM. The formation of vacancy islands in the Au(111) surface upon the adsorption of S-containing molecules is a well-known phenomenon that relieves surface stress caused by the formation of the strong Au–S bonds.33,34
Fig. 4b shows a zoomed image of the ordered arrangement of the SAM of 1. A fast Fourier transform (FFT) (Fig. 4c) of the image displayed in Fig. 4b reveals a rectangular symmetry in line with the four-fold symmetry of the semithio-bambus[4]uril 1 (Fig. 1a, top). The two periodicities extracted from the apparent height profiles (Fig. 4d–e) correspond to 1 and 1.15 nm, respectively. Based on these findings we propose that the molecules self-assemble in a (4 × 2 ) unit cell (red box in Fig. 4f) where (1 × 1) refers to the unreconstructed Au(111) surface with a nearest neighbor distance of 0.288 nm. The dark depressions are most likely the cavities present within the semithio-bambusuril monolayers,17,35 whereas the brighter areas are the sulfur bonding groups.36,37 From the XPS data (see Fig. 2 and Table 1) it is clear that the sulfur atoms in 1 form strong bonds with the Au(111) surface. Unfortunately, the resolution of the images is insufficient to draw any firm conclusions regarding the exact number of binding sites and geometry of the molecules.
Similarly, Fig. 5a shows an STM image of a SAM of 2 on Au(111). Again, vacancy islands with a depth of 2.5 Å were observed indicating that the islands are in the Au (111) layer, albeit more and larger vacancy islands were found as compared to 1. Fig. 5b shows a zoomed image of 2 monolayer. The FFT of the image displayed in Fig. 5b (Fig. 5c) reveals a hexagonal symmetry, which is in line with the six-fold symmetry of the semithio-BU[6], 2, molecule (Fig. 1a, bottom).The observed periodicity is 0.86 nm, which corresponds to the distance between the sulfur atoms of 2 and thus, the bright spots are most likely the sulfur atoms.36,38 This is also supported by the relatively larger fluctuations observed in the apparent height profiles (Fig. 5d and e) of the monolayer of 2 on gold in comparison to that of 1.
A closer look at the FFT in Fig. 5c reveals a second, but substantially weaker, (3 × 3)R30° periodicity in a hexagonal symmetry (green circles). This larger periodicity is 30° rotated with respect to the sulfur (3 × 3) unit cell and corresponds to the arrangement of the semithio-BU[6], 2, molecules shown in Fig. 5e. The unit cell is a (3 × 3) hexagon and contains 9 sulfur atoms and three molecules. Based on these findings we propose that all six thiocarbonyl sulfur atoms in the semithio-BU[6], 2, interact with the Au(111) surface. XPS results suggest that all 6 of the sulfur atoms bind to the gold (Table 1 and Fig. 2f) while MD simulations (vide infra, Fig. 8c) predict that 5 to 6 sulfur atoms interact with the surface. We propose that all six sulfur atoms point down towards the surface, however, three of the sulfur atoms make a direct and strong covalent bond with the gold surface, whereas the other three sulfur atoms interact differently with the gold surface, i.e. via a Au adatom.39 This type of bonding could be a result of an oxidized sulfur atom37 of high conformational strain, which leads to weakening of the Au–S bond.40 When Au adatom bonding occurs, the Au–Au contact is weakened, which makes the Au–S complex more mobile and leads to a reduced intensity in the STM images. However, more research is needed to find the exact origin of the different bonding configurations. The latter falls outside scope of our current work. To further resolve the atomic scale features of the exact binding modes we supplemented the STM data with MD computer simulations (Fig. 6).
In the MD simulations (full details in Methods), compound 1 was initially placed in a DMSO solution exhibiting two major conformations. In one conformation, the glycoluril subunits are arranged in alternate conformation with a pair of thiocarbonyl groups pointing to one direction and another pair of thiocarbonyls groups pointing to the opposite direction (2/2). Another conformation includes three of the four thiocarbonyl groups pointing to one direction and the remaining thiocarbonyl group is flipped to the opposite direction (3/1) (see Fig. 6). In this fully solvated phase, the molecule did not spontaneously sample orientations with all four thiocarbonyl groups pointing in the same direction (4/0) during 100 nanoseconds of free, unconstrained molecular dynamics in DMSO, indicating that a putative binding mode with all sulfurs pointing down involves a significant penalty due to intramolecular conformational strain.
Molecule-surface binding modes (for one molecule adhering to a gold substrate) were generated from naturally occurring sulfur–gold interactions: for example, the monovalent complex was prepared from the solution conformation in which a sulfur spontaneously approached within bonding distance to a gold atom on the surface. The two-site binding complex was prepared from a suitably oriented single-site binding mode, and so on. The starting conformation and final conformations of 1 in water are shown in Fig. S7 and S8.†
Fig. 7 shows the calculated changes in energy for all four possibilities of binding states: single site via Au–S bond (1), and so on, compared to free 1 (0) with no Au–S bonds. All binding modes from mono- to tetravalent are stabilized relative to the unbound state (Fig. 7a), with increasing steric strain offset by stabilisation of −164 kJ mol−1 per additional Au–S bond.41 The conformational energy is quite similar for all binding modes up to trivalent (Fig. 7b) with the penalty for tetravalent due to the arched 4/0 conformation that does not occur spontaneously in simulations of 1 in solution. Small solvation penalties also come with each additional Au–S contact (Fig. 7c). The computed energies indicate that formation of the last Au–S bond creates significant strain in 1; the required high-energy conformation does not form spontaneously in the more loosely bound states. XPS (Table 1 and Fig. 2) showed that about 15% of the sulfur atoms are oxidized which corresponds to one sulfur atom per 1 molecule. The presence of an oxygen atom to bridge the fourth sulfur atom could potentially reduce the conformational strain by increasing the effective length of the molecule and reduce the steric penalty for formation of the fourth bond, as calculated previously for multivalent binding of PAMAM vs. PPI dendrimers.42
Computed binding modes of compound 2 to Au(111) are given in Fig. 8. For the two three-binding modes, regardless of the exact order of thiocarbonyl binding pattern, we observe that, after 100 ns of simulation, the remaining unbound thiocarbonyl groups also point towards the gold surface. This observation clearly suggests that the flexibility of 2 allows larger conformational changes in comparison with 1. Regardless of the initial sulfur bonding pattern, 5 to 6 sulfur atoms were found to interact with the surface within 10 ns, reflecting the greater flexibility of 2 compared with 1 which reduces the steric constraints of the maximally bound conformation. Despite the fact that these calculations were performed for individual molecules, we do not expect SAM packing to have a large effect on the structures given the shortness of the anchoring groups.
For IRRAS measurements, the cryostat was inserted into the optical path of a PEM module of a Nicolet FT-IR spectrometer, making sure that the incidence of the IR beam on the shroud windows was normal. The angle of incidence on the gold surface was 82° with a resolution of 4 cm−1, and 200 scans were collected.
XPS measurements were performed on a Quantera SXM (scanning XPS microprobe) from Physical Electronics, equipped with a monochromatic Al Kα X-ray source monochromatic at 1486.6 eV, producing approximately 50 W of X-ray power. Spectra were referenced to the main C 1s peak set at 284.0 eV. A surface area of 600 × 300 μm2 was scanned with an X-ray beam size about 200 μm (mappings are often done with beam sizes down to 9 μm).
All electrochemical measurements were performed on a CH instruments bipotentiostat 760D using a three-electrode setup using the SAM-covered gold plate as the working electrode (area = 0.43 cm2), a platinum disc as the counter electrode, and a Ag/AgCl reference electrode (Eo Ag/AgCl) 0.197 V vs. reversible hydrogen electrode (RHE) in aqueous 0.1 M K2SO4. Electrochemical impedance measurements were carried out from 1 × 106 to 0.1 Hz with amplitude of 10 mV at the formal potential of the measured species after a preconditioning step of 10 s.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9me00149b |
This journal is © The Royal Society of Chemistry 2020 |