Ion-molecule reactions catalyzed by a single gold atom

We report that Au atoms within van der Waals complexes serve as catalysts for the first time. This was observed in ionization-induced chemistry of 1,6-hexanediol–Au and 1,8-octanediol–Au complexes formed in superfluid helium nanodroplets, where the addition of Au atom(s) made C₂H₄⁺ the sole prominent product in dissociative reactions. Density functional theory (DFT) calculations showed that the Au atom significantly strengthens all of the C–C bonds and weakens the C–O bonds in the meantime, making the C–C bonds stronger than the two C–O bonds in the ionized complexes. This leads to a preferential cleavage of the C–O bonds and thus a strong catalytic effect of the Au atoms in the reactions.

Single Au atoms within van der Waals complexes are found to serve as catalysts in ionisation-induced chemistry for the first time.     

Active control of strong plasmon–exciton coupling in biomimetic pigment–polymer antenna complexes grown by surface-initiated polymerisation from gold nanostructures

Plexcitonic antenna complexes, inspired by photosynthetic light-harvesting complexes, are formed by attachment of chlorophylls (Chl) to poly(cysteine methacrylate) (PCysMA) scaffolds grown by atom-transfer radical polymerisation from gold nanostructure arrays. In these pigment–polymer antenna complexes, localised surface plasmon resonances on gold nanostructures are strongly coupled to Chl excitons, yielding hybrid light–matter states (plexcitons) that are manifested in splitting of the plasmon band. Modelling of the extinction spectra of these systems using a simple coupled oscillator model indicates that their coupling energies are up to twice as large as those measured for LHCs from plants and bacteria. Coupling energies are correlated with the exciton density in the grafted polymer layer, consistent with the collective nature of strong plasmon–exciton coupling. Steric hindrance in fully-dense PCysMA brushes limits binding of bulky chlorophylls, but the chlorophyll concentration can be increased to ∼2 M, exceeding that in biological light-harvesting complexes, by controlling the grafting density and polymerisation time. Moreover, synthetic plexcitonic antenna complexes display pH- and temperature-responsiveness, facilitating active control of plasmon–exciton coupling. Because of the wide range of compatible polymer chemistries and the mild reaction conditions, plexcitonic antenna complexes may offer a versatile route to programmable molecular photonic materials.

Excitons in pigment–polymer antenna complexes formed by attachment of chlorophyll to surface grafted polymers are coupled strongly to plasmon modes, with coupling energies twice those for biological light-harvesting complexes and active control of plasmon–exciton coupling.     

Structural basis for recognition of bacterial cell wall teichoic acid by pseudo-symmetric SH3b-like repeats of a viral peptidoglycan hydrolase

Endolysins are bacteriophage-encoded peptidoglycan hydrolases targeting the cell wall of host bacteria via their cell wall-binding domains (CBDs). The molecular basis for selective recognition of surface carbohydrate ligands by CBDs remains elusive. Here, we describe, in atomic detail, the interaction between the Listeria phage endolysin domain CBD500 and its cell wall teichoic acid (WTA) ligands. We show that 3′O-acetylated GlcNAc residues integrated into the WTA polymer chain are the key epitope recognized by a CBD binding cavity located at the interface of tandem copies of beta-barrel, pseudo-symmetric SH3b-like repeats. This cavity consists of multiple aromatic residues making extensive interactions with two GlcNAc acetyl groups via hydrogen bonds and van der Waals contacts, while permitting the docking of the diastereomorphic ligands. Our multidisciplinary approach tackled an extremely challenging protein–glycopolymer complex and delineated a previously unknown recognition mechanism by which a phage endolysin specifically recognizes and targets WTA, suggesting an adaptable model for regulation of endolysin specificity.

Combining genetic, biochemical and computational approaches, we elucidated the molecular mechanisms underlying the recognition of Listeria wall teichoic acid by bacteriophage-encoded SH3b repeats.     

Chlorophyll excitation energies and structural stability of the CP47 antenna of photosystem II: a case study in the first-principles simulation of light-harvesting complexes

Natural photosynthesis relies on light harvesting and excitation energy transfer by specialized pigment–protein complexes. Their structure and the electronic properties of the embedded chromophores define the mechanisms of energy transfer. An important example of a pigment–protein complex is CP47, one of the integral antennae of the oxygen-evolving photosystem II (PSII) that is responsible for efficient excitation energy transfer to the PSII reaction center. The charge-transfer excitation induced among coupled reaction center chromophores resolves into charge separation that initiates the electron transfer cascade driving oxygenic photosynthesis. Mapping the distribution of site energies among the 16 chlorophyll molecules of CP47 is essential for understanding excitation energy transfer and overall antenna function. In this work, we demonstrate a multiscale quantum mechanics/molecular mechanics (QM/MM) approach utilizing full time-dependent density functional theory with modern range-separated functionals to compute for the first time the excitation energies of all CP47 chlorophylls in a complete membrane-embedded cyanobacterial PSII dimer. The results quantify the electrostatic effect of the protein on the site energies of CP47 chlorophylls, providing a high-level quantum chemical excitation profile of CP47 within a complete computational model of “near-native” cyanobacterial PSII. The ranking of site energies and the identity of the most red-shifted chlorophylls (B3, followed by B1) differ from previous hypotheses in the literature and provide an alternative basis for evaluating past approaches and semiempirically fitted sets. Given that a lot of experimental studies on CP47 and other light-harvesting complexes utilize extracted samples, we employ molecular dynamics simulations of isolated CP47 to identify which parts of the polypeptide are most destabilized and which pigments are most perturbed when the antenna complex is extracted from PSII. We demonstrate that large parts of the isolated complex rapidly refold to non-native conformations and that certain pigments (such as chlorophyll B1 and β-carotene h1) are so destabilized that they are probably lost upon extraction of CP47 from PSII. The results suggest that the properties of isolated CP47 are not representative of the native complexed antenna. The insights obtained from CP47 are generalizable, with important implications for the information content of experimental studies on biological light-harvesting antenna systems.

Advanced QM/MM simulations explore the excited states of a photosynthetic light-harvesting antenna in its physiologically complexed state and model the consequences of extraction on conformational and electronic properties.     

A Carbohydrate–Anion Recognition System in Aprotic Solvents

A carbohydrate–anion recognition system in nonpolar solvents is reported, in which complexes form at the B‐faces of β‐D ‐pyranosides with H1‐, H3‐, and H5‐cis patterns similar to carbohydrate–π interactions. The complexation effect was evaluated for a range of carbohydrate structures; it resulted in either 1:1 carbohydrate–anion complexes, or 1:2 complex formation depending on the protection pattern of the carbohydrate. The interaction was also evaluated with different anions and solvents. In both cases it resulted in significant binding differences. The results indicate that complexation originates from van der Waals interactions or weak CH ??? A^? hydrogen bonds between the binding partners and is related to electron‐withdrawing groups of the carbohydrates as well as increased hydrogen‐bond‐accepting capability of the anions.     

Anisotropy of the van der Waals Configuration of Atoms in Complex,Condensed, and Gas-Phase Molecules

Anisotropic van der Waals (vdW) radii of 5b–7bsubgroup elements were determined from structures of gas-phase van der Waals complexes and crystal molecular compounds. The anisotropy of the van der Waals configuration of atoms was shown to decrease when going from isolated molecules to the condensed state. Variations in intermolecular distances, which are usually explained in terms of the formation of hydrogen bonds, are substantially governed by the anisotropic effect.     

A polarizable dipole–dipole interaction model for evaluation of the interaction energies for NH···OC and CH···OC hydrogen‐bonded complexes

In this article, a polarizable dipole–dipole interaction model is established to estimate the equilibrium hydrogen bond distances and the interaction energies for hydrogen‐bonded complexes containing peptide amides and nucleic acid bases. We regard the chemical bonds N? H, C?O, and C? H as bond dipoles. The magnitude of the bond dipole moment varies according to its environment. We apply this polarizable dipole–dipole interaction model to a series of hydrogen‐bonded complexes containing the N? H···O?C and C? H···O?C hydrogen bonds, such as simple amide‐amide dimers, base‐base dimers, peptide‐base dimers, and β‐sheet models. We find that a simple two‐term function, only containing the permanent dipole–dipole interactions and the van der Waals interactions, can produce the equilibrium hydrogen bond distances compared favorably with those produced by the MP2/6‐31G(d) method, whereas the high‐quality counterpoise‐corrected (CP‐corrected) MP2/aug‐cc‐pVTZ interaction energies for the hydrogen‐bonded complexes can be well‐reproduced by a four‐term function which involves the permanent dipole–dipole interactions, the van der Waals interactions, the polarization contributions, and a corrected term. Based on the calculation results obtained from this polarizable dipole–dipole interaction model, the natures of the hydrogen bonding interactions in these hydrogen‐bonded complexes are further discussed. © 2013 Wiley Periodicals, Inc.     

Playing with the weakest supramolecular interactions in a 3D crystalline hexakis[60]fullerene induces control over hydrogenation selectivity

Weak forces can play an essential role in chemical reactions. Controlling such subtle forces in reorganization processes by applying thermal or chemical stimuli represents a novel synthetic strategy and one of the main targets in supramolecular chemistry. Actually, to separate the different supramolecular contributions to the stability of the 3D assemblies is still a major challenge. Therefore, a clear differentiation of these contributions would help in understanding the intrinsic nature as well as the chemical reactivity of supramolecular ensembles. In the present work, a controlled reorganization of an hexakis[60]fullerene-based molecular compound purely governed by the weakest van der Waals interactions known, i.e. the dihydrogen interaction – usually called sticky fingers – is illustrated. This pre-reorganization of the hexakis[60]fullerene under mild conditions allows a further selective hydrogenation of the crystalline material via hydrazine vapors exposure. This unique two-step transformation process is monitored by single-crystal to single-crystal diffraction (SCSC) which allows the direct observation of the molecular movements in the lattice and the subsequent solid–gas hydrogenation reaction.

Weak forces play an essential role in chemistry. Controlling these supramolecular interactions will contribute to the creation dynamic absorbent materials with a variety of technological applications as chemosensors and environmental remediation.     

Interaction identification of Zif268 and TATAZF proteins with GC‐/AT‐rich DNA sequence: A theoretical study

Molecular dynamics (MD) simulations for Zif268 (a zinc‐finger‐protein binding specifically to the GC‐rich DNA)‐d(A₁G₂C₃G₄T₅G₆G₇G₈C₉A₁₀C₁₁)₂ and TATA_ZF (a zinc‐finger‐protein recognizing the AT‐rich DNA)‐d(A₁C₂G₃C₄T₅A₆T₇A₈A₉A₁₀A₁₁G₁₂G₁₃)₂ complexes have been performed for investigating the DNA binding affinities and specific recognitions of zinc fingers to GC‐rich and AT‐rich DNA sequences. The binding free energies for the two systems have been further analyzed by using the molecular mechanics Poisson‐Boltzmann surface area (MM‐PBSA) method. The calculations of the binding free energies reveal that the affinity energy of Zif268‐DNA complex is larger than that of TATA_ZF‐DNA one. The affinity between the zinc‐finger‐protein and DNA is mainly driven by more favorable van‐der‐Waals and nonpolar/solvation interactions in both complexes. However, the affinity energy difference of the two binding systems is mainly caused by the difference of van‐der‐Waals interactions and entropy components. The decomposition analysis of MM‐PBSA free energies on each residue of the proteins predicts that the interactions between the residues with the positive charges and DNA favor the binding process; while the interactions between the residues with the negative charges and DNA behave in the opposite way. The interhydrogen‐bonds at the protein‐DNA interface and the induced intrafinger hydrogen bonds between the residues of protein for the Zif268‐DNA complex have been identified at some key contact sites. However, only the interhydrogen‐bonds between the residues of protein and DNA for TATA_ZF‐DNA complex have been found. The interactions of hydrogen‐bonds, electrostatistics and van‐der‐Waals type at some new contact sites have been identified. Moreover, the recognition characteristics of the two studied zinc‐finger‐proteins have also been discussed. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

The origin of unidirectional charge separation in photosynthetic reaction centers: nonadiabatic quantum dynamics of exciton and charge in pigment–protein complexes

Exciton charge separation in photosynthetic reaction centers from purple bacteria (PbRC) and photosystem II (PSII) occurs exclusively along one of the two pseudo-symmetric branches (active branch) of pigment–protein complexes. The microscopic origin of unidirectional charge separation in photosynthesis remains controversial. Here we elucidate the essential factors leading to unidirectional charge separation in PbRC and PSII, using nonadiabatic quantum dynamics calculations in conjunction with time-dependent density functional theory (TDDFT) with the quantum mechanics/molecular mechanics/polarizable continuum model (QM/MM/PCM) method. This approach accounts for energetics, electronic coupling, and vibronic coupling of the pigment excited states under electrostatic interactions and polarization of whole protein environments. The calculated time constants of charge separation along the active branches of PbRC and PSII are similar to those observed in time-resolved spectroscopic experiments. In PbRC, Tyr-M210 near the accessary bacteriochlorophyll reduces the energy of the intermediate state and drastically accelerates charge separation overcoming the electron–hole interaction. Remarkably, even though both the active and inactive branches in PSII can accept excitons from light-harvesting complexes, charge separation in the inactive branch is prevented by a weak electronic coupling due to symmetry-breaking of the chlorophyll configurations. The exciton in the inactive branch in PSII can be transferred to the active branch via direct and indirect pathways. Subsequently, the ultrafast electron transfer to pheophytin in the active branch prevents exciton back transfer to the inactive branch, thereby achieving unidirectional charge separation.

Essential factors leading to unidirectional charge separation in photosynthetic reaction centers are clarified via nonadiabatic quantum dynamics calculations.     

Anisotropy of Atomic Van der Waals Radii in the Gas-Phase and Condensed Molecules\sp{1}

The anisotropic components of van der Waals radii for 5b, 6b and 7b Group elements were determined from the structures of gaseous van der Waals complexes and of crystalline molecular solids. The transition from the gaseous to the solid state reduces the anisotropy of a van der Waals atomic shape. This anisotropy is largely responsible for the changes of intermolecular distances, which are often misinterpreted as an effect of hydrogen bonds.     

Room temperature conductance switching in a molecular iron(iii) spin crossover junction

Herein, we report the first room temperature switchable Fe(iii) molecular spin crossover (SCO) tunnel junction. The junction is constructed from [Fe^III(qsal-I)₂]NTf₂ (qsal-I = 4-iodo-2-[(8-quinolylimino)methyl]phenolate) molecules self-assembled on graphene surfaces with conductance switching of one order of magnitude associated with the high and low spin states of the SCO complex. Normalized conductance analysis of the current–voltage characteristics as a function of temperature reveals that charge transport across the SCO molecule is dominated by coherent tunnelling. Temperature-dependent X-ray absorption spectroscopy and density functional theory confirm the SCO complex retains its SCO functionality on the surface implying that van der Waals molecule—electrode interfaces provide a good trade-off between junction stability while retaining SCO switching capability. These results provide new insights and may aid in the design of other types of molecular devices based on SCO compounds.

Herein, we report the first room temperature switchable Fe(iii) molecular spin crossover (SCO) tunnel junction.     

The orientation of N-H...O=C and N-H...N hydrogen bonds in biological systems: How good is a point charge as a model for a hydrogen bonding atom?

In order to design new ligands for protein-binding sites of unknownstructure, it would be useful to predict the likely sites of hydrogenbonding of an unknown protein fragment to a known molecule. The positions ofmaxima and minima in the electrostatic potential at appropriate distancesfrom the van der Waals surface were calculated for various small molecules,nucleic acid bases, peptide units and amino acid side chains containinggroups which can form the biologically important N-H...O=C andN-H...N hydrogen bonds. Their ability to predict the positions of H andO/N in hydrogen bonded complexes, as predicted by optimising theelectrostatic interactions of pairs of such molecules constrained by themolecular shapes, was assessed. It is shown that extrema in theelectrostatic potential around the isolated molecules give worthwhilepredictions for the locations of hydrogen bonding partners. For moleculesbound by a single N-H...O=C hydrogen bond, the electrostatic maximumassociated with the H is usually less than 1 Å from an acceptor atom,while a C=O electrostatic minimum is generally less than 1.5 Å fromthe hydrogen bond proton. However, a significant number of hydrogen bondsform to the opposite lone pair from the electrostatic minimum, in which casethe separation is up to 3.3 Å. This reflects the broad electrostaticpotential well around a carbonyl oxygen between the lone pair directions.The model predicts when neighbouring atoms drastically change the hydrogenbonding characteristics of an N-H or C=O group. Although the geometries ofhydrogen bonded complexes are influenced by the other van der Waals contactsbetween the molecules, particularly multiple hydrogen bonds, theseinfluences are constant when considering hydrogen bonding to a specificuncharacterised binding site. Hence, the consideration of stericallyaccessible electrostatic extrema will be useful in the design of newligands.     

Structures and Chemical Bonding in Antimony(III) Bromide Complexes with Pyridine

Weakly or “partially” bonded molecules are an important link between the chemical and van der Waals interactions. Molecular structures of six new SbBr₃-Py complexes in the solid state have been determined by single-crystal X-ray diffraction analysis. In all complexes all Sb atoms adopt a pseudo-octahedral coordination geometry which is completed by additional Sb⋅⋅⋅Br contacts shorter than the sum of the van der Waals radii, with Br−Sb⋅⋅⋅Br angles close to 180°. To reveal the nature of Sb–Br and Sb–N interactions, the DFT calculations were performed followed by the analysis of the electrostatic potentials, the orbital interactions and the topological analysis. Based on Natural Bond Orbital (NBO) analysis, the Sb–Br interactions range from the covalent bonds to the pnictogen bonds. A simple structural parameter, non-covalence criterion (NCC) is defined as a ratio of the atom-atom distance to the linear combination of sums of covalent and van der Waals radii. NCC correlates with E(2) values for Sb−N, Sb−Cl and Sb−Br bonds, and appears to be useful criterion for a preliminary evaluation of the bonding situation.     

Steroid binding by antibodies and artificial receptors: Exploration of theoretical methods to determine the origins of binding affinities and specificities

Binding mode calculations for complexes between an artificial paracyclophane receptor and digoxins, cholic acids as well as cortisone steroids show encapsulation of different ring combinations. Docking experiments were performed between the 26-10 antibody and digoxins. Coordination affinity arises from hydrophobic desolvation and van der Waals interactions rather than from hydrogen bonds. The specificity and affinity arises mainly from shape complementarity. Computed binding free energies and Kohonen neural network computations both point to physicochemical and structural similarities of natural antibodies and artificial receptors.     

1‐Benzoyl‐3‐[(2‐benzylsulfanyl)phenyl]thiourea

In the title compound, C₂₁H₁₈N₂OS₂, a strong intramolecular N—H...O hydrogen bond [N...O = 2.642 (3) Å] between the amide N atom and the benzoyl O atom forms an almost planar six‐membered ring in the central part of the molecule. In the crystal, molecules are packed through weak N—H...S interactions. Intra‐ and intermolecular hydrogen bonds and van der Waals interactions are the stabilizing forces for the crystal structure.     

An electric field-based approach for quantifying effective volumes and radii of chemically affected space

Chemical shape and size play a critical role in chemistry. The van der Waals (vdW) radius, a familiar manifold used to quantify size by assuming overlapping spheres, provides rapid estimates of size in atoms, molecules, and materials. However, the vdW method may be too rigid to describe highly polarized systems and chemical species that stray from spherical atomistic environments. To deal with these exotic chemistries, numerous alternate methods based on electron density have been presented. While each boasts inherent generality, all define the size of a chemical system, in one way or another, by its electron density. Herein, we revisit the longstanding problem of assessing sizes of atoms and molecules, instead through examination of the local electric field produced by them. While conceptually different than nuclei-centered methods like that of van der Waals, the field assesses chemically affected volumes. This approach implicitly accounts for long-range fields in highly polar systems and predicts that cations should affect more space than neutral counterparts.

Computing atomic and molecular volumes from DFT and ab initio-based electric fields.     

Towards photoswitchable quadruple hydrogen bonds via a reversible “photolocking” strategy for photocontrolled self-assembly

Developing new photoswitchable noncovalent interaction motifs with controllable bonding affinity is crucial for the construction of photoresponsive supramolecular systems and materials. Here we describe a unique “photolocking” strategy for realizing photoswitchable control of quadruple hydrogen-bonding interactions on the basis of modifying the ureidopyrimidinone (UPy) module with an ortho-ester substituted azobenzene unit as the “photo-lock”. Upon light irradiation, the obtained Azo-UPy motif is capable of unlocking/locking the partial H-bonding sites of the UPy unit, leading to photoswitching between homo- and heteroquadruple hydrogen-bonded dimers, which has been further applied for the fabrication of novel tunable hydrogen bonded supramolecular systems. This “photolocking” strategy appears to be broadly applicable in the rational design and construction of other H-bonding motifs with sufficiently photoswitchable noncovalent interactions.

A photolocking strategy is described to achieve the construction of effectively photoswitchable quadruple hydrogen bonds featuring with photoregulable H-bonding affinities, which is further applied in the photocontrollable H-bonded self-assemblies.     

Self-healing and shape memory functions exhibited by supramolecular liquid-crystalline networks formed by combination of hydrogen bonding interactions and coordination bonding

We here report a new approach to develop self-healing shape memory supramolecular liquid-crystalline (LC) networks through self-assembly of molecular building blocks via combination of hydrogen bonding and coordination bonding. We have designed and synthesized supramolecular LC polymers and networks based on the complexation of a forklike mesogenic ligand with Ag⁺ ions and carboxylic acids. Unidirectionally aligned fibers and free-standing films forming layered LC nanostructures have been obtained for the supramolecular LC networks. We have found that hybrid supramolecular LC networks formed through metal–ligand interactions and hydrogen bonding exhibit both self-healing properties and shape memory functions, while hydrogen-bonded LC networks only show self-healing properties. The combination of hydrogen bonds and metal–ligand interactions allows the tuning of intermolecular interactions and self-assembled structures, leading to the formation of the dynamic supramolecular LC materials. The new material design presented here has potential for the development of smart LC materials and functional LC membranes with tunable responsiveness.

New supramolecular hybrid liquid-crystalline networks exhibiting self-healing and shape memory properties are developed by self-assembly of small components through hydrogen bonding interactions and coordination bonding.     

表皮生长因子受体和抑制剂之间分子对接的研究

研究了表皮生长因子受体(EGFR)和4-苯胺喹唑啉类抑制剂之间的相互作用模式，表皮生长因子受体的三维结构通过同源蛋白模建的方法得到，而抑制剂和靶酶结合复合物结构则通过分子力学和分子动力学结合的方法计算得到。从模拟结果得到的抑制剂和靶酶之间的相互作用模式表明范德华相互作用、疏水相互作用以及氢键相互作用对抑制剂的活性都有重要的影响，抑制剂的苯胺部分位于活性口袋的底部，能够与受体残基的非极性侧链产生很强的范德华和疏水相互作用，抑制剂双环上的取代基团也能和活性口袋外部的部分残基形成一定的范德华和疏水性相互作用，而抑制剂喹唑啉环上的氮原子能和周围的残基形成较强的氢键相互作用，对抑制剂的活性有较大的影响，计算得到抑制剂和靶酶之间的非键相互作用能以及抑制剂和靶酶之间的相互作用信息能够很好地解释抑制剂活性和结构的关系，为全新抑制剂的设计提供了重要的结构信息。