Rational modification of the hierarchy of intermolecular interactions in molecular crystal structures by using tunable halogen bonds

A family of 16 isomolecular salts (3-XpyH)(2)[MX'(4)] (3-XpyH=3-halopyridinium; M=Co, Zn; X=(F), Cl, Br, (I); X'=Cl, Br, I) each containing rigid organic cations and tetrahedral halometallate anions has been prepared and characterized by X-ray single crystal and/or powder diffraction. Their crystal structures reflect the competition and cooperation between non-covalent interactions: N-H center dot center dot center dot X'-M hydrogen bonds, C-X center dot center dot center dot X'-M halogen bonds and pi-pi stacking. The latter are essentially unchanged in strength across the series, but both halogen bonds and hydrogen bonds are modified in strength upon changing the halogens involved. Changing the organic halogen (X) from F to I strengthens the C-X center dot center dot center dot X'-M halogen bonds, whereas an analogous change of the inorganic halogen (X') weakens both halogen bonds and N-H center dot center dot center dot X'-M hydrogen bonds. By so tuning the strength of the putative halogen bonds from repulsive to weak to moderately strong attractive interactions, the hierarchy of the interactions has been modified rationally leading to systematic changes in crystal packing. Three classes of crystal structure are obtained. In type A (C-F center dot center dot center dot X'-M) halogen bonds are absent. The structure is directed by N-H center dot center dot center dot X'-M hydrogen bonds and pi-stacking interactions. In type B structures, involving small organic halogens (X) and large inorganic halogens (X'), long (weak) C-X center dot center dot center dot X'-M interactions are observed with type I halogen-halogen interaction geometries (C-X center dot center dot center dot X' approximate to X center dot center dot center dot X'-M approximate to 155 degrees), but hydrogen bonds still dominate. Thus, minor but quite significant perturbations from the type A structure arise. In type C, involving larger organic halogens (X) and smaller inorganic halogens (X'), stronger halogen bonds are formed with a type II halogen-halogen interaction geometry (C-X center dot center dot center dot X' approximate to 180 degrees; X center dot center dot center dot X'-M approximate to 110 degrees) that is electrostatically attractive. The halogen bonds play a major role alongside hydrogen bonds in directing the type C structures, which as a result are quite different from type A and B.

A model for stretch-bend interactions in molecular force fields

Interaction force constants between bond-stretching and angle-bending co-ordinates in polyatomic molecules have been attributed, by some authors, to changes of hybridization due to orbital-following of the bending co-ordinate, and consequent changes of bond length due to the change of hybridization. A method is described for using this model quantitatively to reduce the number of independent force constants in the potential function of a polyatomic molecule, by relating stretch-bend interaction constants to the corresponding diagonal stretching constants. It is proposed to call this model the Hybrid Orbital Force Field. The model is applied to the tetrahedral four co-ordinated carbon atom (as in methane) and to the trigonal planar three coordinated carbon atom (as in formaldehyde).

Phytohormones and plant-herbivore-pathogen interactions: integrating the molecular with the ecological

Current research into indirect phytopathogen–herbivore interactions (i.e., interactions mediated by the host plant) is carried out in two largely independent directions: ecological/mechanistic and molecular. We investigate the origin of these approaches and their strengths and weaknesses. Ecological studies have determined the effect of herbivores and phytopathogens on their host plants and are often correlative: the need for long-term manipulative experiments is pressing. Molecular/cellular studies have concentrated on the role of signaling pathways for systemic induced resistance, mainly involving salicylic acid and jasmonic acid, and more recently the cross-talk between these pathways. This cross-talk demonstrates how interactions between signaling mechanisms and phytohormones could mediate plant–herbivore–pathogen interactions. A bridge between these approaches may be provided by field studies using chemical induction of defense, or investigating whole-organism mechanisms of interactions among the three species. To determine the role of phytohormones in induced resistance in the field, researchers must combine ecological and molecular methods. We discuss how these methods can be integrated and present the concept of “kaleidoscopic defense.” Our recent molecular-level investigations of interactions between the herbivore Gastrophysa viridula and the rust fungus Uromyces rumicis on Rumex obtusifolius, which were well studied at the mechanistic and ecological levels, illustrate the difficulty in combining these different approaches. We suggest that the choice of the right study system (possibly wild relatives of model species) is important, and that molecular studies must consider the environmental conditions under which experiments are performed. The generalization of molecular predictions to ecologically realistic settings will be facilitated by “middle-ground studies” concentrating on the outcomes of the interactions.

A supramolecular polymer based on tweezer-type pi-pi stacking interactions: molecular design for healability and enhanced toughness

An elastomeric, healable, supramolecular polymer blend comprising a chain-folding polyimide and a telechelic polyurethane with pyrenyl end groups is compatibilized by aromatic pi-pi stacking between the pi-electron-deficient diimide groups and the pi-electron-rich pyrenyl units. This interpolymer interaction is the key to forming a tough, healable, elastomeric material. Variable-temperature FTIR analysis of the bulk material also conclusively demonstrates the presence of hydrogen bonding, which complements the pi-pi stacking interactions. Variable-temperature SAXS analysis shows that the healable polymeric blend has a nanophase-separated morphology and that the X-ray contrast between the two types of domain increases with increasing temperature, a feature that is repeatable over several heating and cooling cycles. A fractured sample of this material reproducibly regains more than 95% of the tensile modulus, 91% of the elongation to break, and 77% of the modulus of toughness of the pristine material.

Effect of headgroup size, charge, and solvent structure on polymer−micelle interactions, studied by molecular dynamics simulations

We performed atomistic molecular dynamics simulations of anionic and cationic micelles in the presence of poly(ethylene oxide) (PEO) to understand why nonionic water-soluble polymers such as PEO interact strongly with anionic micelles but only weakly with cationic micelles. Our micelles include sodium n-dodecyl sulfate (SDS), n-dodecyl trimethylammonium chloride (DTAC), n-dodecyl ammonium chloride (DAC), and micelles in which we artificially reverse the sign of partial charges in SDS and DTAC. We observe that the polymer interacts hydrophobically with anionic SDS but only weakly with cationic DTAC and DAC, in agreement with experiment. However, the polymer also interacts with the artificial anionic DTAC but fails to interact hydrophobically with the artificial cationic SDS, illustrating that large headgroup size does not explain the weak polymer interaction with cationic micelles. In addition, we observe through simulation that this preference for interaction with anionic micelles still exists in a dipolar "dumbbell" solvent, indicating that water structure and hydrogen bonding alone cannot explain this preferential interaction. Our simulations suggest that direct electrostatic interactions between the micelle and polymer explain the preference for interaction with anionic micelles, even though the polymer overall carries no net charge. This is possible given the asymmetric distribution of negative charges on smaller atoms and positive charges oil larger units in the polymer chain.

Molecular dynamics simulation of interactions between a sodium dodecyl sulfate micelle and a poly(ethylene oxide) polymer

We have performed atomistic molecular dynamics simulations of an anionic sodium dodecyl sulfate (SDS) micelle and a nonionic poly(ethylene oxide) (PEO) polymer in aqueous solution. The micelle consisted of 60 surfactant molecules, and the polymer chain lengths varied from 20 to 40 monomers. The force field parameters for PEO were adjusted by using 1,2-dimethoxymethane (DME) as a model compound and matching its hydration enthalpy and conformational behavior to experiment. Excellent agreement with previous experimental and simulation work was obtained through these modifications. The simulated scaling behavior of the PEO radius of gyration was also in close agreement with experimental results. The SDS-PEO simulations show that the polymer resides on the micelle surface and at the hydrocarbon-water interface, leading to a selective reduction in the hydrophobic contribution to the solvent-accessible surface area of the micelle. The association is mainly driven by hydrophobic interactions between the polymer and surfactant tails, while the interaction between the polymer and sulfate headgroups on the micelle surface is weak. The 40-monomer chain is mostly wrapped around the micelle, and nearly 90% of the monomers are adsorbed at low PEO concentration. Simulations were also performed with multiple 20-monomer chains, and gradual addition of polymer indicates that about 120 monomers are required to saturate the micelle surface. The stoichiometry of the resulting complex is in close agreement with experimental results, and the commonly accepted "beaded necklace" structure of the SDS-PEO complex is recovered by our simulations.

Self-assembled molecular complexes and coordination polymers of CdII, hexamine, and monocarboxylates: structural analysis and theoretical studies of supramolecular interactions

Four new cadmium(II) complexes [Cd-2(bz)(4)(H2O)(4)(mu 2-hmt)]center dot Hbz center dot H2O (1), [Cd-3(bz)(6)(H2O)(6)(mu 2-hmt)(2)]center dot 6H(2)O (2), [Cd(pa)(2)(H2O)(mu(2)-hmt)](n) (3), and {[Cd-3(ac)(6)(H2O)(3)(mu(3)-hmt)(2)]center dot 6H(2)O}(n) (4) with hexamine (hmt) and monocarboxylate ions, benzoate (bz), phenylacetate (pa), or acetate (ac) have been synthesized and characterized structurally. Structure determinations reveal that 1 is dinuclear, 2 is trinuclear, 3 is a one-dimensional (1D) infinite chain, and 4 is a two-dimensional (2D) polymer with fused hexagonal rings consisting of Cd-II and hmt. All the Cd-II atoms in the four complexes (except one CdII in 2) possess seven-coordinate pentagonal bipyramidal geometry with the various chelating bidentate carboxylate groups in equatorial sites. One of the CdII ions in 2, a complex that contains two monodentate carboxylates is in a distorted octahedral environment. The bridging mode of hmt is mu 2- in complexes 1-3 but is mu 3- in complex 4. In all complexes, there are significant numbers of H-bonds, C-H/pi, and pi-pi interactions which play crucial roles in forming the supramolecular networks. The importance of the noncovalent interactions in terms of energies and geometries has been analyzed using high level ab initio calculations. The effect of the cadmium coordinated to hmt on the energetic features of the C-H/pi interaction is analyzed. Finally, the interplay between C-H/pi and pi-pi interactions observed in the crystal structure of 3 is also studied.

Weak intermolecular interactions in an ionically bound molecular adsorbate: cyclopentadienyl=Cu(111)

The dissociative adsorption of cyclopentadiene (C5H6) on Cu(111) yields a cyclopentadienyl (Cp) species with strongly anionic characteristics. The Cp potential energy surface and frictional coupling to the substrate are determined from measurements of dynamics of the molecule together with density functional calculations. The molecule is shown to occupy degenerate threefold adsorption sites and molecular motion is characterized by a low diffusional energy barrier of 40 +/- 3 meV with strong frictional dissipation. Repulsive dipole-dipole interactions are not detected despite charge transfer from substrate to adsorbate.

Proteomic profiling of the molecular targets of interactions of the mastoparan peptide Protopolybia MP-III at the level of endosomal membranes from rat mast cells

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Vibrational spectroscopy for probing molecular-level interactions in organic films mimicking biointerfaces

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Role of peptide peptide interactions in aggregation: protonectins observed in equilibrium and replica exchange molecular dynamics simulations

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Contribution of Non-Covalent Interactions and Electronic Effects on the Conformational Landscape and Tautomeric Equilibria of Molecules and Molecular Complexes: Structural and Dynamical Data from Rotational Spectroscopy

This thesis concerns the study of complex conformational surfaces and tautomeric equilibria of molecules and molecular complexes by quantum chemical methods and rotational spectroscopy techniques. In particular, the focus of this research is on the effects of substitution and noncovalent interactions in determining the energies and geometries of different conformers, tautomers or molecular complexes. The Free-Jet Absorption Millimeter Wave spectroscopy and the Pulsed-Jet Fourier Transform Microwave spectroscopy have been applied to perform these studies and the obtained results showcase the suitability of these techniques for the study of conformational surfaces and intermolecular interactions. The series of investigations of selected medium-size molecules and complexes have shown how different instrumental setups can be used to obtain a variety of results on molecular properties. The systems studied, include molecules of biological interest such as anethole and molecules of astrophysical interest such as N-methylaminoethanol. Moreover halogenation effects have been investigated on halogen substituted tautomeric systems (5-chlorohydroxypyridine and 6-chlorohydroxypyridine), where it has shown that the position of the inserted halogen atom affects the prototropic equilibrium. As for fluorination effects, interesting results have been achieved investigating some small complexes where a molecule of water is used as a probe to reveal the changes on the electrostatic potential of different fluorinated compounds: 2-fluoropyridine, 3-fluoropyridine and penta-fluoropyridine. While in the case of the molecular complex between water and 2-fluoropyridine and 3-fluoropyridine the geometry of the complex with one water molecule is analogous to that of pyridine with the water molecule linked to the pyridine nitrogen, the case of pentafluoropyridine reveals the effect of perfluorination and the water oxygen points towards the positive center of the pyridine ring. Additional molecular adducts with a molecule of water have been analyzed (benzylamine-water and acrylic acid-water) in order to reveal the stabilizing driving forces that characterize these complexes.

Numerical analysis of the relation between interactions and structure in a molecular fluid

Coarse graining is a popular technique used in physics to speed up the computer simulation of molecular fluids. An essential part of this technique is a method that solves the inverse problem of determining the interaction potential or its parameters from the given structural data. Due to discrepancies between model and reality, the potential is not unique, such that stability of such method and its convergence to a meaningful solution are issues.rnrnIn this work, we investigate empirically whether coarse graining can be improved by applying the theory of inverse problems from applied mathematics. In particular, we use the singular value analysis to reveal the weak interaction parameters, that have a negligible influence on the structure of the fluid and which cause non-uniqueness of the solution. Further, we apply a regularizing Levenberg-Marquardt method, which is stable against the mentioned discrepancies. Then, we compare it to the existing physical methods - the Iterative Boltzmann Inversion and the Inverse Monte Carlo method, which are fast and well adapted to the problem, but sometimes have convergence problems.rnrnFrom analysis of the Iterative Boltzmann Inversion, we elaborate a meaningful approximation of the structure and use it to derive a modification of the Levenberg-Marquardt method. We engage the latter for reconstruction of the interaction parameters from experimental data for liquid argon and nitrogen. We show that the modified method is stable, convergent and fast. Further, the singular value analysis of the structure and its approximation allows to determine the crucial interaction parameters, that is, to simplify the modeling of interactions. Therefore, our results build a rigorous bridge between the inverse problem from physics and the powerful solution tools from mathematics. rn

Specific Molecular Chaperone Interactions and an ATP-dependent Conformational Change Are Required during Posttranslational Protein Translocation into the Yeast ER

The posttranslational translocation of proteins across the endoplasmic reticulum (ER) membrane in yeast requires ATP hydrolysis and the action of hsc70s (DnaK homologues) and DnaJ homologues in both the cytosol and ER lumen. Although the cytosolic hsc70 (Ssa1p) and the ER lumenal hsc70 (BiP) are homologous, they cannot substitute for one another, possibly because they interact with specific DnaJ homologues on each side of the ER membrane. To investigate this possibility, we purified Ssa1p, BiP, Ydj1p (a cytosolic DnaJ homologue), and a GST–63Jp fusion protein containing the lumenal DnaJ region of Sec63p. We observed that BiP, but not Ssa1p, is able to associate with GST–63Jp and that Ydj1p stimulates the ATPase activity of Ssa1p up to 10-fold but increases the ATPase activity of BiP by <2-fold. In addition, Ydj1p and ATP trigger the release of an unfolded polypeptide from Ssa1p but not from BiP. To understand further how BiP drives protein translocation, we purified four dominant lethal mutants of BiP. We discovered that each mutant is defective for ATP hydrolysis, fails to undergo an ATP-dependent conformational change, and cannot interact with GST–63Jp. Measurements of protein translocation into reconstituted proteoliposomes indicate that the mutants inhibit translocation even in the presence of wild-type BiP. We conclude that a conformation- and ATP-dependent interaction of BiP with the J domain of Sec63p is essential for protein translocation and that the specificity of hsc70 action is dictated by their DnaJ partners.