Molecular organization of n-cyanobiphenyl liquid crystals on a molybdenite surface

The alignement and anchoring of liquid crystals on solid surfaces is a key problem for modern device technology that until now has been treated empirically, but that can now be tackled by atomistic computer simulations. Molecular dynamics (MD) simulations were used in this thesis work to study two films of 7 and 8 n-alkyl-4’cyanobiphenyl (7CB and 8CB) liquid crystals , with a thickness of 15 nm, confined between two (001) surfaces of MoS2 (molybdenite). The isotropic and nematic phases of both liquid crystals were simulated, and the resulting structures characterized structurally. A new force field was designed to model the interactions between the liquid crystal (LC) molecules and the surface of molybdenite, while an accurate force field developed previously was used to model the 7CB and 8CB molecules. The results show that the (001) molybdenite surface induces a planar orientation in both the liquid crystals. For the nematic phase of 8CB, one of the two solid/LC interfaces is composed of a first layer of molecules aligned parallel to the surface, followed by a second layer of molecules aligned perpendicular to the surface (also called, homeotropic). The effect of the surface appears to be local in nature as it is confined to the first 15 Angström of the LC film. Conversely, for the nematic phase of 7CB, a planar ordering is established into the LC film. The LC molecules at the interface with the molybdenite appear to align preferentially their alkyl chains toward the solid substrate. The resulting tilt angle of molecules was found to be in good agreement with experimental measurements available in literature. Despite the fact that the MD simulations spanned a time range of more than 100 ns, the nematic phases of both 7CB and 8CB were found not to be completely formed. In order to confirm the findings presented in this thesis, we propose to extend the current study.

Spectroscopic properties from hybrid QM, MM molecular dynamics simulation

The central aim of this thesis work is the application and further development of a hybrid quantum mechanical/molecular mechanics (QM/MM) based approach to compute spectroscopic properties of molecules in complex chemical environments from electronic structure theory. In the framework of this thesis, an existing density functional theory implementation of the QM/MM approach is first used to calculate the nuclear magnetic resonance (NMR) solvent shifts of an adenine molecule in aqueous solution. The findings show that the aqueous solvation with its strongly fluctuating hydrogen bond network leads to specific changes in the NMR resonance lines. Besides the absolute values, also the ordering of the NMR lines changes under the influence of the solvating water molecules. Without the QM/MM scheme, a quantum chemical calculation could have led to an incorrect assignment of these lines. The second part of this thesis describes a methodological improvement of the QM/MM method that is designed for cases in which a covalent chemical bond crosses the QM/MM boundary. The development consists in an automatized protocol to optimize a so-called capping potential that saturates the electronic subsystem in the QM region. The optimization scheme is capable of tuning the parameters in such a way that the deviations of the electronic orbitals between the regular and the truncated (and "capped") molecule are minimized. This in turn results in a considerable improvement of the structural and spectroscopic parameters when computed with the new optimized capping potential within the QM/MM technique. This optimization scheme is applied and benchmarked on the example of truncated carbon-carbon bonds in a set of small test molecules. It turns out that the optimized capping potentials yield an excellent agreement of NMR chemical shifts and protonation energies with respect to the corresponding full molecules. These results are very promising, so that the application to larger biological complexes will significantly improve the reliability of the prediction of the related spectroscopic properties.

Solid-state NMR on columnar supramolecular structures

The present work deals with the characterisation of three columnar self-assembled systems, that is, benzene-1,3,5-tricarboxamides, a peripherally thioalkyl-substituted phthalocyanine, and several oligo-(p-phenylenevinylene)s. In order to probe the supramolecular organisation solid-state NMR has been used as the main technique, supported by X-ray measurements, theoretical methods, and thermal analysis. rnrnBenzene-1,3,5-tricarboxamides (BTAs) turned out to be well suited model compounds to study various fundamental supramolecular interactions, such as π-π-interactions, hydrogen bonding, as well as dynamic and steric effects of attached side chains. Six BTAs have been investigated in total, five with a CO-centred amide group bearing different side chains and one with an inverted N-centred amide group. The physical properties of these BTAs have been investigated as a function of temperature. The results indicated that in case of the CO-centred BTAs the stability of the columnar mesophase depends strongly on the nature of the side chains. Further experiments revealed a coplanar orientation of adjacent BTA molecules in the columnar assembly of CO-centred BTAs, whereas the N-centred BTA, showed a deviating not fully coplanar arrangement. These differences were ascribed to distinct hydrogen bonding schemes, involving a parallel alignment of hydrogen bonds in case of CO-centred BTAs and an antiparallel alignment in case of the N-centred counterpart.rnrn The fundamental insights of the supramolecular organisation of BTAs could be partially adapted to an octa-substituted phthalocyanine with thiododecyl moieties. Solid-state NMR in combination with chemical shift calculations determined a tilted herringbone arrangement of phthalocyanine rings in the crystalline phase as well as in the mesophase. Moreover, 1H NMR measurements in the mesophase of this compound suggested an axial rotation of molecules, which is inhibited in the crystalline phase.rnrnAs a third task, the supramolecular assembly of oligo-(p-phenylenevinylene)s of varying length and with different polar head groups have been investigated by a combined X-ray and solid-state NMR study. The results revealed a columnar structure formation of these compounds, being promoted by phase separation of alkyl side chains and aromatic rigid rods. In this system solid-state NMR yielded meaningful insight into the isotropisation process of butoxy and 2-S-methylbutoxy substituted oligo-(p-phenylenevinylene) rods.rn

Thermodynamics of Forming Water Clusters at Various Temperatures and Pressures by Gaussian-2, Gaussian-3, Complete Basis Set-QB3, and Complete Basis Set-APNO Model Chemistries; Implications for Atmospheric Chemistry

The Gaussian-2, Gaussian-3, complete basis set- (CBS-) QB3, and CBS-APNO methods have been used to calculate ΔH° and ΔG° values for neutral clusters of water, (H2O)n, where n = 2−6. The structures are similar to those determined from experiment and from previous high-level calculations. The thermodynamic calculations by the G2, G3, and CBS-APNO methods compare well against the estimated MP2(CBS) limit. The cyclic pentamer and hexamer structures release the most heat per hydrogen bond formed of any of the clusters. While the cage and prism forms of the hexamer are the lowest energy structures at very low temperatures, as temperature is increased the cyclic structure is favored. The free energies of cluster formation at different temperatures reveal interesting insights, the most striking being that the cyclic trimer, cyclic tetramer, and cyclic pentamer, like the dimer, should be detectable in the lower troposphere. We predict water dimer concentrations of 9 × 1014 molecules/cm3, water trimer concentrations of 2.6 × 1012 molecules/cm3, tetramer concentrations of approximately 5.8 × 1011 molecules/cm3, and pentamer concentrations of approximately 3.5 × 1010 molecules/cm3 in saturated air at 298 K. These results have important implications for understanding the gas-phase chemistry of the lower troposphere.

N-H center dot center dot center dot pi hydrogen-bonding and large-amplitude tipping vibrations in jet-cooled pyrrole-benzene

The N-H center dot center dot center dot pi hydrogen bond is an important intermolecular interaction in many biological systems. We have investigated the infrared (IR) and ultraviolet (UV) spectra of the supersonic-jet cooled complex of pyrrole with benzene and benzene-d(6) (Pyr center dot Bz, Pyr center dot Bz-d(6)). DFT-D density functional, SCS-MP2 and SCS-CC2 calculations predict a T-shaped and (almost) C(s) symmetric structure with an N-H center dot center dot center dot pi hydrogen bond to the benzene ring. The pyrrole is tipped by omega(S(0)) = +/- 13 degrees relative to the surface normal of Bz. The N center dot center dot center dot ring distance is 3.13 angstrom. In the S(1) excited state, SCS-CC2 calculations predict an increased tipping angle omega(S(1)) = +/- 21 degrees. The IR depletion spectra support the T-shaped geometry: The NH stretch is redshifted by -59 cm(-1), relative to the "free" NH stretch of pyrrole at 3531 cm(-1), indicating a moderately strong N-H center dot center dot center dot pi interaction. The interaction is weaker than in the (Pyr)(2) dimer, where the NH donor shift is -87 cm(-1) [Dauster et al., Phys. Chem. Chem. Phys., 2008, 10, 2827]. The IR C-H stretch frequencies and intensities of the Bz subunit are very similar to those of the acceptor in the (Bz)(2) dimer, confirming that Bz acts as the acceptor. While the S(1) <- S(0) electronic origin of Bz is forbidden and is not observable in the gas-phase, the UV spectrum of Pyr center dot Bz in the same region exhibits a weak 0(0)(0) band that is red-shifted by 58 cm(-1) relative to that of Bz (38 086 cm(-1)). The origin appears due to symmetry-breaking of the p-electron system of Bz by the asymmetric pyrrole NH center dot center dot center dot pi hydrogen bond. This contrasts with (Bz)(2), which does not exhibit a 0(0)(0) band. The Bz moiety in Pyr center dot Bz exhibits a 6a(0)(1) band at 0(0)(0) + 518 cm(-1) that is about 20x more intense than the origin band. The symmetry breaking by the NH center dot center dot center dot pi hydrogen bond splits the degeneracy of the v(6)(e(2g)) vibration, giving rise to 6a' and 6b' sub-bands that are spaced by similar to 6 cm(-1). Both the 0(0)(0) and 6(0)(1) bands of Pyr center dot Bz carry a progression in the low-frequency (10 cm(-1)) excited-state tipping vibration omega', in agreement with the change of the omega tipping angle predicted by SCS-MP2 and SCS-CC2 calculations.

Crystal structure of the Leishmania major phosphodiesterase LmjPDEB1 and insight into the design of the parasite-selective inhibitors

Human leishmaniasis is a major public health problem in many countries, but chemotherapy is in an unsatisfactory state. Leishmania major phosphodiesterases (LmjPDEs) have been shown to play important roles in cell proliferation and apoptosis of the parasite. Thus LmjPDE inhibitors may potentially represent a novel class of drugs for the treatment of leishmaniasis. Reported here are the kinetic characterization of the LmjPDEB1 catalytic domain and its crystal structure as a complex with 3-isobutyl-1-methylxanthine (IBMX) at 1.55 A resolution. The structure of LmjPDEB1 is similar to that of human PDEs. IBMX stacks against the conserved phenylalanine and forms a hydrogen bond with the invariant glutamine, in a pattern common to most inhibitors bound to human PDEs. However, an extensive structural comparison reveals subtle, but significant differences between the active sites of LmjPDEB1 and human PDEs. In addition, a pocket next to the inhibitor binding site is found to be unique to LmjPDEB1. This pocket is isolated by two gating residues in human PDE families, but constitutes a natural expansion of the inhibitor binding pocket in LmjPDEB1. The structure particularity might be useful for the development of parasite-selective inhibitors for the treatment of leishmaniasis.

Synthesis, structure and magnetic properties of cobalt(II) and copper(II) coordination polymers assembled by phthalate and 4-methylimidazole

New coordination polymers [M(Pht)(4-MeIm)2(H2O)]n (M=Co (1), Cu (2); Pht2−=dianion of o-phthalic acid; 4-MeIm=4-methylimidazole) have been synthesized and characterized by IR spectroscopy, X-ray crystallography, thermogravimetric analysis and magnetic measurements. The crystal structures of 1 and 2 are isostructural and consist of [M(4-MeIm)2(H2O)] building units linked in infinite 1D helical chains by 1,6-bridging phthalate ions which also act as chelating ligands through two O atoms from one carboxylate group in the case of 1. In complex 1, each Co(II) atom adopts a distorted octahedral N2O4 geometry being coordinated by two N atoms from two 4-MeIm, three O atoms of two phthalate residues and one O atom of a water molecule, whereas the square-pyramidal N2O3 coordination of the Cu(II) atom in 2 includes two N atoms of N-containing ligands, two O atoms of two carboxylate groups from different Pht, and a water molecule. An additional strong O–H⋯O hydrogen bond between a carboxylate group of the phthalate ligand and a coordinated water molecule join the 1D helical chains to form a 2D network in both compounds. The thermal dependences of the magnetic susceptibilities of the polymeric helical Co(II) chain compound 1 were simulated within the temperature range 20–300 K as a single ion case, whereas for the Cu(II) compound 2, the simulations between 25 and 300 K, were made for a linear chain using the Bonner–Fisher approximation. Modelling the experimental data of compound 1 with MAGPACK resulted in: g=2.6, |D|=62 cm−1. Calculations using the Bonner–Fisher approximation gave the following result for compound 2: g=2.18, J=–0.4 cm−1.

Hydrogen-bonding complexes of functionalized BEDT-TTF derivatives with chloranilic acid

The two crystalline donor-acceptor complexes showing hydrogen-bondings between bis(ethylenedithio) tetrathiofulvalene (BEDT-TTF) derivatives containing pyridine and pyrazine groups and 2,5-dichloro-3,6-dihydroxyl-1,4-benzoquinone (chloranilic acid) were prepared. X-ray structure analyses revealed that functional groups play an important role in constructing the unique crystal structures.

X-Ray Diffraction, Calorimetric and Dielectric Relaxation Study of the Amorphous and Smectic States of a Main Chain Liquid Crystalline Polymer

Los polímeros cristales líquidos (LCP) son sistemas complejos que forman mesofases que presentan orden orientacional y polímeros amorfos. Con frecuencia, el estado amorfo isotrópico no puede ser estudiado debido a la rápida formación de mesofases. En este trabajo se ha sintetizado y estudiado un nuevo LCP: poli(trietilenglicol metil p, p '-bibenzoato), PTEMeB. Este polímero presenta una formación de mesofase bastante lenta haciendo posible estudiar de forma independiente tanto los estados amorfo y de cristal líquidos. La estructura y las transiciones de fase del PTEMeB han sido investigados por calorimetría (DSC), con MAXS / WAXS con temperatura variable que emplean radiación de sincrotrón y con difracción de rayos X. Estos estudios han mostrado la existencia de dos transiciones vítreas, relacionadas con las fases amorfa y cristal líquido. Se ha realizado un estudio de relajación dieléctrica en amplios intervalos de temperatura y presión. Se ha encontrado que la transición vítrea dinámica de la fase amorfa es más lenta que la del cristal líquido. El estudio de la relajación ? nos ha permitido seguir la formación isoterma de la mesofase a presión atmosférica. Además, con el estudio el comportamiento dinámico a alta presión se ha encontrado que se produce la formación rápida de la mesofase inducida por cambios bruscos de presión. Liquid crystalline polymers (LCPs) are complex systems that include features of both orientationally ordered mesophases and amorphous polymers. Frequently, the isotropic amorphous state cannot be studied due to the rapid mesophase formation. Here, a new main chain LCP, poly(triethyleneglycol methyl p,p'-bibenzoate), PTEMeB, has been synthesized. It shows a rather slow mesophase formation making possible to study independently both the amorphous and the liquid crystalline states. The structure and phase transitions of PTEMeB have been investigated by calorimetry, variable-temperature MAXS/WAXS employing synchrotron radiation, and X-ray diffraction in oriented fibers. These experiments have pointed out the presence of two glass transitions, related to the amorphous or to the liquid crystal phases. Additionally, the mesophase seems to be a coexistence of orthogonal and tilted smectic phases. A dielectric relaxation study of PTEMeB over broad ranges of temperature and pressure has been performed. The dynamic glass transition turns out to be slower for the amorphous state than for the liquid crystal. Monitoring of the α relaxation has allowed us to follow the isothermal mesophase formation at atmospheric pressure. Additionally, the dynamical behavior at high pressures has pointed out the fast formation of the mesophase induced by sudden pressure changes.

A dimeric crystal structure for the N-terminal two domains of intercellular adhesion molecule-1

The 3.0-Å structure of a 190-residue fragment of intercellular adhesion molecule-1 (ICAM-1, CD54) reveals two tandem Ig-superfamily (IgSF) domains. Each of two independent molecules dimerizes identically with a symmetry-related molecule over a hydrophobic interface on the BED sheet of domain 1, in agreement with dimerization of ICAM-1 on the cell surface. The residues that bind to the integrin LFA-1 are well oriented for bivalent binding in the dimer, with the critical Glu-34 residues pointing away from each other on the periphery. Residues that bind to rhinovirus are in the flexible BC and FG loops at the tip of domain 1, and these and the upper half of domain 1 are well exposed in the dimer for docking to virus. By contrast, a residue important for binding to Plasmodium falciparum-infected erythrocytes is in the dimer interface. The presence of A′ strands in both domains 1 and 2, conserved hydrogen bonds at domain junctions, and elaborate hydrogen bond networks around the key integrin binding residues in domain 1 make these domains suited to resist tensile forces during adhesive interactions. A subdivision of the intermediate (I) set of IgSF domains is proposed in which domain 1 of ICAM-1 and previously described I set domains belong to the I1 set and domain 2 of ICAM-1, ICAM-2, and vascular cell adhesion molecule-1 belong to the I2 set.

NMR reveals hydrogen bonds between oxygen and distal histidines in oxyhemoglobin

Compared with free heme, the proteins hemoglobin (Hb) and myoglobin (Mb) exhibit greatly enhanced affinity for oxygen relative to carbon monoxide. This physiologically vital property has been attributed to either steric hindrance of CO or stabilization of O2 binding by a hydrogen bond with the distal histidine. We report here the first direct evidence of such a hydrogen bond in both α- and β-chains of oxyhemoglobin, as revealed by heteronuclear NMR spectra of chain-selectively labeled samples. Using these spectra, we have assigned the imidazole ring 1H and 15N chemical shifts of the proximal and distal histidines in both carbonmonoxy- and oxy-Hb. Because of their proximity to the heme, these chemical shifts are extremely sensitive to the heme pocket conformation. Comparison of the measured chemical shifts with values predicted from x-ray structures suggests differences between the solution and crystal structures of oxy-Hb. The chemical shift discrepancies could be accounted for by very small displacements of the proximal and distal histidines. This suggests that NMR could be used to obtain very high-resolution heme pocket structures of Hb, Mb, and other heme proteins.

Toward a quantum-mechanical description of metal-assisted phosphoryl transfer in pyrophosphatase

The wealth of kinetic and structural information makes inorganic pyrophosphatases (PPases) a good model system to study the details of enzymatic phosphoryl transfer. The enzyme accelerates metal-complexed phosphoryl transfer 1010-fold: but how? Our structures of the yeast PPase product complex at 1.15 Å and fluoride-inhibited complex at 1.9 Å visualize the active site in three different states: substrate-bound, immediate product bound, and relaxed product bound. These span the steps around chemical catalysis and provide strong evidence that a water molecule (Onu) directly attacks PPi with a pKa vastly lowered by coordination to two metal ions and D117. They also suggest that a low-barrier hydrogen bond (LBHB) forms between D117 and Onu, in part because of steric crowding by W100 and N116. Direct visualization of the double bonds on the phosphates appears possible. The flexible side chains at the top of the active site absorb the motion involved in the reaction, which may help accelerate catalysis. Relaxation of the product allows a new nucleophile to be generated and creates symmetry in the elementary catalytic steps on the enzyme. We are thus moving closer to understanding phosphoryl transfer in PPases at the quantum mechanical level. Ultra-high resolution structures can thus tease out overlapping complexes and so are as relevant to discussion of enzyme mechanism as structures produced by time-resolved crystallography.

Crystal structure of the complex of a catalytic antibody Fab fragment with a transition state analog: structural similarities in esterase-like catalytic antibodies.

The x-ray structure of the complex of a catalytic antibody Fab fragment with a phosphonate transition-state analog has been determined. The antibody (CNJ206) catalyzes the hydrolysis of p-nitrophenyl esters with significant rate enhancement and substrate specificity. Comparison of this structure with that of the uncomplexed Fab fragment suggests hapten-induced conformational changes: the shape of the combining site changes from a shallow groove in the uncomplexed Fab to a deep pocket where the hapten is buried. Three hydrogen-bond donors appear to stabilize the charged phosphonate group of the hapten: two NH groups of the heavy (H) chain complementarity-determining region 3 (H3 CDR) polypeptide chain and the side-chain of histidine-H35 in the H chain (His-H35) in the H1 CDR. The combining site shows striking structural similarities to that of antibody 17E8, which also has esterase activity. Both catalytic antibody ("abzyme") structures suggest that oxyanion stabilization plays a significant role in their rate acceleration. Additional catalytic groups that improve efficiency are not necessarily induced by the eliciting hapten; these groups may occur because of the variability in the combining sites of different monoclonal antibodies that bind to the same hapten.

Disulfide bond mutagenesis and the structure and function of the head-to-tail macrocyclic trypsin inhibitor SFTI-1

SFTI-1 is a novel 14 amino acid peptide comprised of a circular backbone constrained by three proline residues, a hydrogen-bond network, and a single disulfide bond. It is the smallest and most potent known Bowman-Birk trypsin inhibitor and the only one with a cyclic peptidic backbone. The solution structure of [ABA(3,11)]SFTI-1, a disulfide-deficient analogue of SFTI-1, has been determined by H-1 NMR spectroscopy. The lowest energy structures of native SFTI-1 and [ABA(3,11)]SFTI-1 are similar and superimpose with a root-mean-square deviation over the backbone and heavy atoms of 0.26 +/- 0.09 and 1.10 +/- 0.22 Angstrom, respectively. The disulfide bridge in SFTI-1 was found to be a minor determinant for the overall structure, but its removal resulted in a slightly weakened hydrogen-bonding network. To further investigate the role of the disulfide bridge, NMR chemical shifts for the backbone H-alpha protons of two disulfide-deficient linear analogues of SFTI-1, [ABA(3,11)]SFTI-1[6,5] and [ABA(3,11)]SFTI-1[1,14] were measured. These correspond to analogues of the cleavage product of SFTI-1 and a putative biosynthetic precursor, respectively. In contrast with the cyclic peptide, it was found that the disulfide bridge is essential for maintaining the structure of these open-chain analogues. Overall, the hydrogen-bond network appears to be a crucial determinant of the structure of SFTI-1 analogues.

Comparison of the binding of 3-fluoromethyl-7-sulfonyl-1,2,3,4-tetrahydroisoquinolines with their isosteric sulfonamides to the active site of phenylethanolamine N-methyltransferase

3-Fluoromethyl-7-(N-substituted aminosulfonyl)-1,2,3,4-tetrahydroisoquinolines (14, 16, and 18-22) are highly potent and selective inhibitors of phenylethanolamine N-methyltransferase (PNMT). Molecular modeling studies with 3-fluoromethyl-7-(N-alkyl aminosulfonyl)-1,2,3,4-tetrahydroisoquinolines, such as 16, suggested that the sulfonamide -NH-could form a hydrogen bond with the side chain of Lys57. However, SAR studies and analysis of the crystal structure of human PNMT (hPNMT) in complex with 7 indicated that the sulfonamide oxygens, and not the sulfonamide -NH-, formed favorable interactions with the enzyme. Thus, we hypothesized that replacement of the sulfonamide -NH-with a methylene group could result in compounds that would retain potency at PNMT and that would have increased lipophilicity, thus increasing the likelihood they will cross the blood brain barrier. A series of 3-fluoromethyl-7-sulfonyl-1,2,3,4-tetrahydroisoquinolines (23-30) were synthesized and evaluated for their PNMT inhibitory potency and affinity for the R2-adrenoceptor. A comparison of these compounds with their isosteric sulfonamides (14, 16, and 18-22) showed that the sulfones were more lipophilic but less potent than their corresponding sulfonamides. Sulfone 24 (hPNMT K-i = 1.3 mu M) is the most potent compound in this series and is quite selective for PNMT versus the R2-adrenoceptor, but 24 is less potent than the corresponding sulfonamide, 16 (hPNMT K-i = 0.13 mu M). We also report the crystal structure of hPNMT in complex with sulfonamide 15, from which a potential hydrogen bond acceptor within the hPNMT active site has been identified, the main chain carbonyl oxygen of Asn39. The interaction of this residue with the sulfonamide -NH-is likely responsible for much of the enhanced inhibitory potency of the sulfonamides versus the sulfones.