Halogen Bonds in Crystal Engineering: Like Hydrogen Bonds yet Different

CONSPECTUS: The halogen bond is an attractive interaction in which an electrophilic halogen atom approaches a negatively polarized species. Short halogen atom contacts in crystals have been known for around 50 years. Such contacts are found in two varieties: type I, which is symmetrical, and type II, which is bent. Both are influenced by geometric and chemical considerations. Our research group has been using halogen atom interactions as design elements in crystal engineering, for nearly 30 years. These interactions include halogen center dot center dot center dot halogen interactions (X center dot center dot center dot X) and halogen center dot center dot center dot heteroatom interactions (X center dot center dot center dot B). Many X center dot center dot center dot X and almost all X center dot center dot center dot B contacts can be classified as halogen bonds. In this Account, we illustrate examples of crystal engineering where one can build up from previous knowledge with a focus that is provided by the modern definition of the halogen bond. We also comment on the similarities and differences between halogen bonds and hydrogen bonds. These interactions are similar because the protagonist atoms halogen and hydrogen are both electrophilic in nature. The interactions are distinctive because the size of a halogen atom is of consequence when compared with the atomic sizes of, for example, C, N, and O, unlike that of a hydrogen atom. Conclusions may be drawn pertaining to the nature of X center dot center dot center dot X interactions from the Cambridge Structural Database (CSD). There is a clear geometric and chemical distinction between type I and type II, with only type II being halogen bonds. Cl/Br isostructurality is explained based on a geometric model. In parallel, experimental studies on 3,4-dichlorophenol and its congeners shed light on the nature of halogen center dot center dot center dot halogen interactions and reveal the chemical difference between Cl and Br. Variable temperature studies also show differences between type I and type II contacts. In terms of crystal design, halogen bonds offer a unique opportunity in the strength, atom size and interaction gradation; this may be used in the design of ternary cocrystals. Structural modularity in which an entire crystal structure is defined as a combination of modules is rationalized on the basis of the intermediate strength of a halogen bond. The specific directionality of the halogen bond makes it a good tool to achieve orthogonality in molecular crystals. Mechanical properties can be tuned systematically by varying these orthogonally oriented halogen center dot center dot center dot halogen interactions. In a further development, halogen bonds are shown to play a systematic role in organization of LSAMs (long range synthon aufbau module), which are bigger structural units containing multiple synthons. With a formal definition in place, this may be the right time to look at differences between halogen bonds and hydrogen bonds and exploit them in more subtle ways in crystal engineering.

Effect of contact stresses on shape recovery of NiTiCu thin films

Permanent plastic deformation induced by mechanical contacts affects the shape recovery of shape memory alloys. To understand the shape recovery of NiTiCu thin films subjected to local contact stresses, systematic investigations are carried out by inducing varying levels of contact stresses using nanoindentation. The resulting indents are located precisely for imaging using a predetermined array consisting of different sized indents. Morphology and topography of these indents before and after shape recovery are characterized using Scanning Electron Microscope and Atomic Force Microscope quantitatively. Shape recovery is found to be dependent on the contact stresses at the low loads while the recovery ratio remains constant at 0.13 for higher loads. Shape recovery is found to occur mainly in depth direction of the indent, while far field residual stresses play very little role in the recovery. (C) 2014 Elsevier B.V. All rights reserved.

Re-analysis of cryoEM data on HCV IRES bound to 40S subunit of human ribosome integrated with recent structural information suggests new contact regions between ribosomal proteins and HCV RNA

In this study, we combine available high resolution structural information on eukaryotic ribosomes with low resolution cryo-EM data on the Hepatitis C Viral RNA (IRES) human ribosome complex. Aided further by the prediction of RNA-protein interactions and restrained docking studies, we gain insights on their interaction at the residue level. We identified the components involved at the major and minor contact regions, and propose that there are energetically favorable local interactions between 40S ribosomal proteins and IRES domains. Domain II of the IRES interacts with ribosomal proteins S5 and S25 while the pseudoknot and the downstream domain IV region bind to ribosomal proteins S26, S28 and S5. We also provide support using UV cross-linking studies to validate our proposition of interaction between the S5 and IRES domains II and IV. We found that domain IIIe makes contact with the ribosomal protein S3a (S1e). Our model also suggests that the ribosomal protein S27 interacts with domain IIIc while S7 has a weak contact with a single base RNA bulge between junction IIIabc and IIId. The interacting residues are highly conserved among mammalian homologs while IRES RNA bases involved in contact do not show strict conservation. IRES RNA binding sites for S25 and S3a show the best conservation among related viral IRESs. The new contacts identified between ribosomal proteins and RNA are consistent with previous independent studies on RNA-binding properties of ribosomal proteins reported in literature, though information at the residue level is not available in previous studies.

Hydrophobic hydration driven self-assembly of curcumin in water: Similarities to nucleation and growth under large metastability, and an analysis of water dynamics at heterogeneous surfaces

As the beneficial effects of curcumin have often been reported to be limited to its small concentrations, we have undertaken a study to find the aggregation properties of curcumin in water by varying the number of monomers. Our molecular dynamics simulation results show that the equilibrated structure is always an aggregated state with remarkable structural rearrangements as we vary the number of curcumin monomers from 4 to 16 monomers. We find that the curcumin monomers form clusters in a very definite pattern where they tend to aggregate both in parallel and anti-parallel orientation of the phenyl rings, often seen in the formation of beta-sheet in proteins. A considerable enhancement in the population of parallel alignments is observed with increasing the system size from 12 to 16 curcumin monomers. Due to the prevalence of such parallel alignment for large system size, a more closely packed cluster is formed with maximum number of hydrophobic contacts. We also follow the pathway of cluster growth, in particular the transition from the initial segregated to the final aggregated state. We find the existence of a metastable structural intermediate involving a number of intermediate-sized clusters dispersed in the solution. We have constructed a free energy landscape of aggregation where the metatsable state has been identified. The course of aggregation bears similarity to nucleation and growth in highly metastable state. The final aggregated form remains stable with the total exclusion of water from its sequestered hydrophobic core. We also investigate water structure near the cluster surface along with their orientation. We find that water molecules form a distorted tetrahedral geometry in the 1st solvation layer of the cluster, interacting rather strongly with the hydrophilic groups at the surface of the curcumin. The dynamics of such quasi-bound water molecules near the surface of curcumin cluster is considerably slower than the bulk signifying a restricted motion as often found in protein hydration layer. (C) 2014 AIP Publishing LLC.

Computational Study of the Formation of Short Centrosymmetric N-center dot center dot center dot S Supramolecular Synthon and Related Weak Interactions in Crystalline 1,2,4-Triazoles

A comprehensive analysis of the crystal packing and the energetic features of a series of four biologically active molecules belonging to the family of substituted 4-(benzylideneamino)-3-(4-fluoro-3-phenoxyphenyl)-1H-1,2,4-triazole-5-(4 H)-thione derivatives have been performed based on the molecular conformation and the supramolecular packing. This involves the formation of a short centrosymmetric R-2(2)(8) NH...S supramolecular synthon in the solid state, including the presence of CH...S, CH...O, CH...N, CH...F, CH...Cl, CF...FC, CCl...ClC, and CH...pi intermolecular interactions along with pp stacking to evaluate the role of noncovalent interactions in the crystal. The presence of such synthons has a substantial contribution toward the interaction energy (-18 to -20 kcal/mol) as obtained from the PIXEL calculation, wherein the Coulombic and polarization contribution are more significant than the dispersion contribution. The geometrical characteristics of such synthons favor short distance, and the population of related molecules having these geometries is rare as has been obtained from the Cambridge Structural Database (CSD). Furthermore, their interaction energies have been compared with those present in our molecules in the solid state. The topological characteristics of the NH...S supramolecular synthon, in addition to related weak interactions, CH...N, CH...Cl, CF...FC, and CCl...ClC, have been estimated using the quantum theory of atoms in molecules (QTAIM). In addition, an analysis of the Hirshfeld surface and associated fingerprint plots of these four molecules also have provided a platform for the evaluation of the contribution of different atom...atom contacts, which contribute toward the packing of the molecules in solids.

Mapping Post-translational Modifications of Mammalian Testicular Specific Histone Variant TH2B in Tetraploid and Haploid Germ Cells and Their Implications on the Dynamics of Nucleosome Structure

Histones regulate a variety of chromatin templated events by their post-translational modifications (PTMs). Although there are extensive reports on the PTMs of canonical histones, the information on the histone variants remains very scanty. Here, we report the identification of different PTMs, such as acetylation, methylation, and phosphorylation of a major mammalian histone variant TH2B. Our mass spectrometric analysis has led to the identification of both conserved and unique modifications across tetraploid spermatocytes and haploid spermatids. We have also computationally derived the 3-dimensional model of a TH2B containing nucleosome in order to study the spatial orientation of the PTMs identified and their effect on nucleosome stability and DNA binding potential. From our nucleosome model, it is evident that substititution of specific amino acid residues in TH2B results in both differential histone-DNA and histone-histone contacts. Furthermore, we have also observed that acetylation on the N-terminal tail of TH2B weakens the interactions with the DNA. These results provide direct evidence that, similar to somatic H2B, the testis specific histone TH2B also undergoes multiple PTMs, suggesting the possibility of chromatin regulation by such covalent modifications in mammalian male germ cells.

The key role of polymer grafted nanoparticles in the phase miscibility of an LCST mixture

Blends of bromo-terminated polystyrene (PS-Br) and poly(vinyl methylether) (PVME) exhibit lower critical solution temperatures. In this study, PS-Br was designed by atom transfer radical polymerization and was converted to thiol-capped polystyrene (PS-SH) by reacting with thiourea. The silver nanoparticles (nAg) were then decorated with covalently bound PS-SH macromolecules to improve the phase miscibility in the PS-Br-PVME blends. Thermally induced demixing in this model blend was followed in the presence of polystyrene immobilized silver nanoparticles (PS-g-nAg). The graft density of the PS macromolecules was estimated to be ca. 0.78 chains per nm(2). Although the matrix and the grafted molecular weights were similar, PS-g-nAg particles were expelled from the PS phase and were localized in the PVME phase of the blends. This was addressed with respect to intermediate graft density and favourable PS-PVME contacts from microscopic interactions point of view. Interestingly, blends with 0.5 wt% PS-g-nAg delayed the spinodal decomposition temperature in the blends by ca. 18 degrees C with respect to the control blends. The scale of cooperativity, as determined by differential scanning calorimetry, increased only marginally in the case of PS-g-nAg; however, it increased significantly in the presence of bare nAg particles.

Use of a structural alphabet to find compatible folds for amino acid sequences

The structural annotation of proteins with no detectable homologs of known 3D structure identified using sequence-search methods is a major challenge today. We propose an original method that computes the conditional probabilities for the amino-acid sequence of a protein to fit to known protein 3D structures using a structural alphabet, known as Protein Blocks (PBs). PBs constitute a library of 16 local structural prototypes that approximate every part of protein backbone structures. It is used to encode 3D protein structures into 1D PB sequences and to capture sequence to structure relationships. Our method relies on amino acid occurrence matrices, one for each PB, to score global and local threading of query amino acid sequences to protein folds encoded into PB sequences. It does not use any information from residue contacts or sequence-search methods or explicit incorporation of hydrophobic effect. The performance of the method was assessed with independent test datasets derived from SCOP 1.75A. With a Z-score cutoff that achieved 95% specificity (i.e., less than 5% false positives), global and local threading showed sensitivity of 64.1% and 34.2%, respectively. We further tested its performance on 57 difficult CASP10 targets that had no known homologs in PDB: 38 compatible templates were identified by our approach and 66% of these hits yielded correctly predicted structures. This method scales-up well and offers promising perspectives for structural annotations at genomic level. It has been implemented in the form of a web-server that is freely available at http://www.bo-protscience.fr/forsa.

Molecular effects of encapsulation of glucose oxidase dimer by graphene

Knowing the nature of the enzyme-graphene interface is critical for a design of graphene-based biosensors. Extensive contacts between graphene and enzyme could be obtained by employing a suitable encapsulation which does not impede its enzymatic reaction. We have performed molecular dynamics simulations to obtain an insight on many forms of contact between glucose oxidase dimer and the single-layer graphene nano-sheets. The unconnected graphene sheets tended to form a flat stack regardless of their initial positions around the enzyme, whereas the same graphene sheets linked together formed a flower-like shape engendering different forms of wrapping of the enzyme. During the encapsulation no core hydrophobic residues of the enzyme were exposed. Since the polar and charged amino acids populated the enzyme's surface we also estimated, using DFT calculations, the interaction energies of individual polar and charged amino acid residues with graphene. It was found that the negatively charged residues can bind to graphene unexpectedly strongly; however, the main effect of encapsulation comes from the overlap of adjacent edges of graphene sheets.

Characterization of friction at three contact pairs by photoelastic Isotropic point (IP)

Friction coefficient between a circular-disk periphery and V-block surface was determined by introducing the concept of isotropic point (IP) in isochromatic field of the disk under three-point symmetric loading. IP position on the symmetry axis depends on active coefficient of friction during experiment. We extend this work to asymmetric loading of circular disk in which case two frictional contact pairs out of three loading contacts, independently control the unconstrained IP location. Photoelastic experiment is conducted on particular case of asymmetric three-point loading of circular disk. Basics of digital image processing are used to extract few essential parameters from experimental image, particularly IP location. Analytical solution by Flamant for half plane with a concentrated load, is utilized to derive stress components for required loading configurations of the disk. IP is observed, in analytical simulations of three-point asymmetric normal loading, to move from vertical axis to the boundary along an ellipse-like curve. When friction is included in the analysis, IP approaches the center with increase in loading friction and it goes away with increase in support friction. With all these insights, using experimental IP information, friction angles at three contact pairs of circular disk under asymmetric loading, are determined.

BLOC-2 targets recycling endosomal tubules to melanosomes for cargo delivery

Hermansky-Pudlak syndrome (HPS) is a group of disorders characterized by the malformation of lysosome-related organelles, such as pigment cell melanosomes. Three of nine characterized HPS subtypes result from mutations in subunits of BLOC-2, a protein complex with no known molecular function. In this paper, we exploit melanocytes from mouse HPS models to place BLOC-2 within a cargo transport pathway from recycling endosomal domains to maturing melanosomes. In BLOC-2-deficient melanocytes, the melanosomal protein TYRP1 was largely depleted from pigment granules and underwent accelerated recycling from endosomes to the plasma membrane and to the Golgi. By live-cell imaging, recycling endosomal tubules of wild-type melanocytes made frequent and prolonged contacts with maturing melanosomes; in contrast, tubules from BLOC-2-deficient cells were shorter in length and made fewer, more transient contacts with melanosomes. These results support a model in which BLOC-2 functions to direct recycling endosomal tubular transport intermediates to maturing melanosomes and thereby promote cargo delivery and optimal pigmentation.

The crystal structure of a lectin from Butea monosperma: Insight into its glycosylation and binding of ligands

Crystal structure of a lectin purified from Butea monosperma seeds was determined by Molecular Replacement method. Its primary structure was determined by Tandem Mass Spectroscopy and electron density maps from X-ray diffraction data. Its quaternary structure was tetrameric, formed of two monomers, alpha and beta, beta appearing as truncated alpha. The occurrence of two tetramers in the asymmetric unit of the crystal might be a consequence of asymmetric contacts due to difference in glycosylation and variable loops structures, to form an `octamer-structure'. The crystal structure showed binding pockets for gamma Abu, having a proposed role in plant defense, at the interface of canonical dimer-partners. Hemagglutination studies, enzyme kinetics, isothermal titration calorimetry and molecular dynamics showed that the lectin is specific to N-acetyl D-galactosamine, galactose and lactose in decreasing order, and alpha-amylase inhibitor. (C) 2014 Elsevier B.V. All rights reserved.

Suspended core-shell Pt-PtOx nanostructure for ultrasensitive hydrogen gas sensor

High sensitivity gas sensors are typically realized using metal catalysts and nanostructured materials, utilizing non-conventional synthesis and processing techniques, incompatible with on-chip integration of sensor arrays. In this work, we report a new device architecture, suspended core-shell Pt-PtOx nanostructure that is fully CMOS-compatible. The device consists of a metal gate core, embedded within a partially suspended semiconductor shell with source and drain contacts in the anchored region. The reduced work function in suspended region, coupled with builtin electric field of metal-semiconductor junction, enables the modulation of drain current, due to room temperature Redox reactions on exposure to gas. The device architecture is validated using Pt-PtO2 suspended nanostructure for sensing H-2 down to 200 ppb under room temperature. By exploiting catalytic activity of PtO2, in conjunction with its p-type semiconducting behavior, we demonstrate about two orders of magnitude improvement in sensitivity and limit of detection, compared to the sensors reported in recent literature. Pt thin film, deposited on SiO2, is lithographically patterned and converted into suspended Pt-PtO2 sensor, in a single step isotropic SiO2 etching. An optimum design space for the sensor is elucidated with the initial Pt film thickness ranging between 10 nm and 30 nm, for low power (< 5 mu W), room temperature operation. (C) 2015 AIP Publishing LLC.

Halogen bonds in some dihalogenated phenols: applications to crystal engineering

3,4-Dichlorophenol (1) crystallizes in the tetragonal space group I4(1)/a with a short axis of 3.7926 (9) angstrom. The structure is unique in that both type I and type II Cl.....Cl interactions are present, these contact types being distinguished by the angle ranges of the respective C-Cl....Cl angles. The present study shows that these two types of contacts are utterly different. The crystal structures of 4-bromo-3-chlorophenol (2) and 3-bromo-4-chlorophenol (3) have been determined. The crystal structure of (2) is isomorphous to that of (1) with the Br atom in the 4-position participating in a type II interaction. However, the monoclinic P2(1)/c packing of compound (3) is different; while the structure still has O-H....O hydrogen bonds, the tetramer O-H.....O synthon seen in (1) and (2) is not seen. Rather than a type I Br....Br interaction which would have been mandated if (3) were isomorphous to (1) and (2), Br forms a Br....O contact wherein its electrophilic character is clearly evident. Crystal structures of the related compounds 4-chloro-3-iodophenol (4) and 3,5-dibromophenol (5) were also determined. A computational survey of the structural landscape was undertaken for (1), (2) and (3), using a crystal structure prediction protocol in space groups P2(1)/c and I4(1)/a with the COMPASS26 force field. While both tetragonal and monoclinic structures are energetically reasonable for all compounds, the fact that (3) takes the latter structure indicates that Br prefers type II over type I contacts. In order to differentiate further between type I and type II halogen contacts, which being chemically distinct are expected to have different distance fall-off properties, a variable-temperature crystallography study was performed on compounds (1), (2) and (4). Length variations with temperature are greater for type II contacts compared with type I. The type II Br....Br interaction in (2) is stronger than the corresponding type II Cl....Cl interaction in (1), leading to elastic bending of the former upon application of mechanical stress, which contrasts with the plastic deformation of (1). The observation of elastic deformation in (2) is noteworthy; in that it finds an explanation based on the strengths of the respective halogen bonds, it could also be taken as a good starting model for future property design. Cl/Br isostructurality is studied with the Cambridge Structural Database and it is indicated that this isostructurality is based on shape and size similarity of Cl and Br, rather than arising from any chemical resemblance.

Accelerating Information Diffusion in Social Networks Under the Susceptible-Infected-Susceptible Epidemic Model

Standard Susceptible-Infected-Susceptible (SIS) epidemic models assume that a message spreads from the infected to the susceptible nodes due to only susceptible-infected epidemic contact. We modify the standard SIS epidemic model to include direct recruitment of susceptible individuals to the infected class at a constant rate (independent of epidemic contacts), to accelerate information spreading in a social network. Such recruitment can be carried out by placing advertisements in the media. We provide a closed form analytical solution for system evolution in the proposed model and use it to study campaigning in two different scenarios. In the first, the net cost function is a linear combination of the reward due to extent of information diffusion and the cost due to application of control. In the second, the campaign budget is fixed. Results reveal the effectiveness of the proposed system in accelerating and improving the extent of information diffusion. Our work is useful for devising effective strategies for product marketing and political/social-awareness/crowd-funding campaigns that target individuals in a social network.