Relaxation and jump dynamics of water at the mica interface

The orientational relaxation dynamics of water confined between mica surfaces is investigated using molecular dynamics simulations. The study illustrates the wide heterogeneity that exists in the dynamics of water adjacent to a strongly hydrophilic surface such as mica. Analysis of the survival probabilities in different layers is carried out by normalizing the corresponding relaxation times with bulk water layers of similar thickness. A 10-fold increase in the survival times is observed for water directly in contact with the mica surface and a non-monotonic variation in the survival times is observed moving away from the mica surface to the bulk-like interior. The orientational relaxation time is highest for water in the contact layer, decreasing monotonically away from the surface. In all cases the ratio of the relaxation times of the 1st and 2nd rank Legendre polynomials of the HH bond vector is found to lie between 1.5 and 1.9 indicating that the reorientational relaxation in the different water layers is governed by jump dynamics. The orientational dynamics of water in the contact layer is particularly novel and is found to undergo distinct two-dimensional hydrogen bond jump reorientational dynamics with an average waiting time of 4.97 ps. The waiting time distribution is found to possess a long tail extending beyond 15 ps. Unlike previously observed jump dynamics in bulk water and other surfaces, jump events in the mica contact layer occur between hydrogen bonds formed by the water molecule and acceptor oxygens on the mica surface. Despite slowing down of the water orientational relaxation near the surface, life-times of water in the hydration shell of the K ion are comparable to that observed in bulk salt solutions. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4717710]

Nanoindentation as a probe for mechanically-induced molecular migration in layered organic donor-acceptor complexes

Nanoindentation and scratch experiments on 1:1 donor-acceptor complexes, 1 and 2, of 1,2,4,5-tetracyanobenzene with pyrene and phenanthrene, respectively, reveal long-range molecular layer gliding and large interaction anisotropy. Due to the layered arrangements in these crystals, these experiments that apply stress in particular directions result in the breaking of interlayer interactions, thus allowing molecular sheets to glide over one another with ease. Complex 1 has a layered crystal packing wherein the layers are 68° skew under the (002) face and the interlayer space is stabilized by van der Waals interactions. Upon indenting this surface with a Berkovich tip, pile-up of material was observed on just one side of the indenter due to the close angular alignment of the layers with the half angle of the indenter tip (65.35°). The interfacial differences in the elastic modulus (21 ) and hardness (16 ) demonstrate the anisotropic nature of crystal packing. In 2, the molecular stacks are arranged in a staggered manner; there is no layer arrangement, and the interlayer stabilization involves C-H�N hydrogen bonds and ��� interactions. This results in a higher modulus (20 ) for (020) as compared to (001), although the anisotropy in hardness is minimal (4 ). The anisotropy within a face was analyzed using AFM image scans and the coefficient of friction of four orthogonal nanoscratches on the cleavage planes of 1 and 2. A higher friction coefficient was obtained for 2 as compared to 1 even in the cleavage direction due to the presence of hydrogen bonds in the interlayer region making the tip movement more hindered. Copyright © 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Topological electron density analysis and electrostatic properties of aspirin: an experimental and theoretical study

The topological and the electrostatic properties of the aspirin drug molecule were determined from high-resolution X-ray diffraction data at 90 K, and the corresponding results are compared with the theoretical calculations. The electron density at the bond critical point of all chemical bonds induding the intermolecular interactions of aspirin has been quantitatively described using Bader's quantum theory of ``Atoms in Molecules''. The electrostatic potential of the molecule emphasizes the preferable binding sites of the drug and the interaction features of the molecule, which are crucial for drug-receptor recognition. The topological analysis of hydrogen bonds reveals the strength of intermolecular interactions.

Tracing a common ``origin'' of phase transformation, polymorphism, disorder, isosterism, and isostructuralism in fluorobenzoylcarvacryl thiourea

The terms phase transformation, polymorphism, disorder, isosterism, and isostructuralism are often the keywords used in the design and engineering of molecular crystals. Three benzoylcarvacryl thiourea derivatives with -NH-C(S)-NH-C(O)-] cores generate molecular crystals, which provide the basis for exploring a common link between the structures related by aforementioned terms. The apparent ``origin'' of all these structural modifications has been traced to the formation of a planar molecular dimeric chain built with homomeric R-2(2)(12) and R-2(2)(8) synthons occurring in tandem, one formed with N-H center dot center dot center dot O and the other with N-H center dot center dot center dot S hydrogen bonds.

Polymorphs and hydrates of Etoricoxib, a selective COX-2 inhibitor

The crystal structures of two polymorphs and two polymorphic hemihydrates of Etoricoxib are reported. Etoricoxib is a non-steroidal anti-inflammatory drug (NSAID) that is a selective inhibitor of COX-2. It is used in the treatment of various types of inflammation, pain and fever. Clas et al. have reported four polymorphs (labeled I through IV) and two solvates (hemi-and sesquihydrate) of the API in US patent 6,441,002 (Clas et al, US patent 6,441,002, 2002). However, no crystal structures have been reported for any of these forms. A comparison was made between the PXRD patterns reported in patent `002 and the powder spectra simulated from single crystal data. The two polymorphs characterized here correspond to form I and form IV of the patent. Form II of the patent could not be obtained by us with a variety of experimental conditions. Form III of the patent corresponds to hemihydrate II of this study. Form III is therefore not a polymorph of form I and form IV. What we have termed hemihydrate I in this study is obtained under a wide variety of conditions and it is also the only hemihydrate reported as such in the patent. Because the Etoricoxib molecule contains no conventional hydrogen bond donors, there cannot be any strong hydrogen bonds in the crystal structures of forms I and IV. The packing is accordingly characterized by weak hydrogen bonds of the C-H center dot center dot center dot O=S and C-H center dot center dot center dot N type. Thermal data were collected for form I, form IV and hemihydrate I to shed some light on relative stabilities. PXRD diffractograms show the transformation of form IV to form I at elevated temperature, indicating that form I is more stable than form IV. However, this transformation occurs only in samples of form IV that contain some form I; it does not occur in pure form IV. The formation of the two hemihydrates could follow from the known tendency of an acceptor-rich molecule to crystallize as a hydrate.

From Molecular to Supramolecular: An Exploration into the Modes of Self-assembly in Conformationally Locked Polycyclitols

This brief account highlights the notable findings of our investigation into the supramolecular chemistry of conformationally locked polycyclitols in the solid state. The study was aimed at analyzing the crystal packing and unraveling the modalities of non-covalent interactions (particularly, intramolecular vis-a-vis intermolecular OH center dot center dot center dot O hydrogen bonds) in polyols. The know-how obtained thereof, was successfully utilized to engineer self-assemblies of designer polycyclitols, having hydrogen bond donors and acceptors fettered onto a trans-decalin scaffold. The results seek to draw particular attention to the intrinsic attribute of this rigid carbocyclic framework to lock functional groups into spatially invariant positions and bring potential intramolecular hydrogen bonding partners into favorable interaction geometry to engender predictability in the self-assembly patterns.

Crystal structures and binding studies of atovaquone and its derivatives with cytochrome bc(1): a molecular basis for drug design

Crystal structure of trans-atovaquone (antimalarial drug), its polymorph and its stereoisomer (cis) along with five other derivatives with different functional groups have been analyzed. Based on the conformational features of these compounds and the characteristics of the nature of intermolecular interactions, valuable insights into the atomistic details of protein-inhibitor interactions have been derived by docking studies. Atovaquone and its derivatives pack in the crystal lattice using intermolecular O-H center dot center dot center dot O hydrogen bond dimer motifs supported by surrogate weak interactions including C-H center dot center dot center dot O and C-H center dot center dot center dot Cl hydrogen bonds. The docking results of these molecules with cytochrome bc(1) show preferences to form N-H center dot center dot center dot O, O-H center dot center dot center dot O and O-H center dot center dot center dot Cl hydrogen bonds. The involvement of halogen atoms in the binding pocket appears to be significant and is contrary to the theoretically predicted mechanism of protein-ligand docking reported earlier based on mimicking experimental binding results of stigmatellin with cytochrome bc(1). The significance of subtle energy factors controlled by weak intermolecular interactions appears to play a major role in drug binding.

Analysis of designed beta-hairpin peptides: molecular conformation and packing in crystals

The crystal structures of several designed peptide hairpins have been determined in order to establish features of molecular conformations and modes of aggregation in the crystals. Hairpin formation has been induced using a centrally positioned (D)Pro-Xxx segment (Xxx = (L)Pro, Aib, Ac(6)c, Ala; Aib = alpha-aminoisobutyric acid; Ac(6)c = 1-aminocyclohexane-1-carboxylic acid). Structures of the peptides Boc-Leu-Phe-Val-(D)Pro-(L)Pro-Leu-Phe-Val-OMe (1), Boc-Leu-Tyr-Val-(D)Pro-(L)Pro-Leu-Phe-Val-OMe (2, polymorphic forms labeled as 2a and 2b), Boc-Leu-Val-Val-(D)Pro-(L)Pro-Leu-Val-Val-OMe (3), Boc-Leu-Phe-Val-(D)Pro-Aib-Leu-Phe-Val-OMe (4, polymorphic forms labeled as 4a and 4b), Boc-Leu-Phe-Val-(D)Pro-Ac(6)c-Leu-Phe-Val-OMe (5) and Boc-Leu-Phe-Val-(D)Pro-Ala-Leu-Phe-Val-OMe (6) are described. All the octapeptides adopt type II' beta-turn nucleated hairpins, stabilized by three or four cross-strand intramolecular hydrogen bonds. The angle of twist between the two antiparallel strands lies in the range of -9.8 degrees to -26.7 degrees. A detailed analysis of packing motifs in peptide hairpin crystals is presented, revealing three broad modes of association: parallel packing, antiparallel packing and orthogonal packing. An attempt to correlate aggregation modes in solution with observed packing motifs in crystals has been made by indexing of crystal faces in the case of three of the peptide hairpins. The observed modes of hairpin aggregation may be of relevance in modeling multiple modes of association, which may provide insights into the structure of insoluble polypeptide aggregates.

Heat and SDS insensitive NDK dimers are largely stabilised by hydrophobic interaction to form functional hexamer in Mycobacterium smegmatis

The primary structure and function of nucleoside diphosphate kinase (NDK), a substrate non-specific enzyme involved in the maintenance of nucleotide pools is also implicated to play pivotal roles in many other cellular processes. NDK is conserved from bacteria to human and forms a homotetramer or hexamer to exhibit its biological activity. However, the nature of the functional oligomeric form of the enzyme differs among different organisms. The functional form of NDKs from many bacterial systems, including that of the human pathogen, Mycobacterium tuberculosis (MtuNDK), is a hexamer, although some bacterial NDKs are tetrameric in nature. The present study addresses the oligomeric property of MsmNDK and how a dimer, the basic subunit of a functional hexamer, is stabilized by hydrogen bonds and hydrophobic interactions. Homology modeling was generated using the three-dimensional structure of MtuNDK as a template; the residues interacting at the monomer-monomer interface of MsmNDK were mapped. Using recombinant enzymes of wild type, catalytically inactive mutant, and monomer-monomer interactive mutants of MsmNDK, the stability of the dimer was verified under heat, SDS, low pH, and methanol. The predicted residues (Gln17, Ser24 and Glu27) were engaged in dimer formation, however the mutated proteins retained the ATPase and GTPase activity even after introducing single (MsmNDK- Q17A, MsmNDK-E27A, and MsmNDK-E27Q) and double (MsmNDK-E27A/Q17A) mutation. However, the monomer monomer interaction could be abolished using methanol, indicating the stabilization of the monomer-monomer interaction by hydrophobic interaction.

Simulations reveal that the HIV-1 gp120-CD4 complex dissociates via complex pathways and is a potential target of the polyamidoamine (PAMAM) dendrimer

The polyamidoamine (PAMAM) dendrimer prevents HIV-1 entry into target cells in vitro. Its mechanism of action, however, remains unclear and precludes the design of potent dendrimers targeting HIV-1 entry. We employed steered molecular dynamics simulations to examine whether the HIV-1 gp120-CD4 complex is a target of PAMAM. Our simulations mimicked single molecule force spectroscopy studies of the unbinding of the gp120-CD4 complex under the influence of a controlled external force. We found that the complex dissociates via complex pathways and defies the standard classification of adhesion molecules as catch and slip bonds. When the force loading rate was large, the complex behaved as a slip bond, weakening gradually. When the loading rate was small, the complex initially strengthened, akin to a catch bond, but eventually dissociated over shorter separations than with large loading rates. PAMAM docked to gp120 and destabilized the gp120-CD4 complex. The rupture force of the complex was lowered by PAMAM. PAMAM disrupted salt bridges and hydrogen bonds across the gp120-CD4 interface and altered the hydration pattern of the hydrophobic cavity in the interface. In addition, intriguingly, PAMAM suppressed the distinction in the dissociation pathways of the complex between the small and large loading rate regimes. Taken together, our simulations reveal that PAMAM targets the gp120-CD4 complex at two levels: it weakens the complex and also alters its dissociation pathway, potentially inhibiting HIV-1 entry.

Antithyroid drugs and their analogues: synthesis, structure, and mechanism of action

Thyroid hormones are essential for the development and differentiation of all cells of the human body. They regulate protein, fat, and carbohydrate metabolism. In this Account, we discuss the synthesis, structure, and mechanism of action of thyroid hormones and their analogues. The prohormone thyroxine (14) is synthesized on thyroglobulin by thyroid peroxidase (TPO), a heme enzyme that uses iodide and hydrogen peroxide to perform iodination and phenolic coupling reactions. The monodeiodination of T4 to 3,3',5-triiodothyronine (13) by selenium-containing deiodinases (ID-1, ID-2) is a key step in the activation of thyroid hormones. The type 3 deiodinase (ID-3) catalyzes the deactivation of thyroid hormone in a process that removes iodine selectively from the tyrosyl ring of T4 to produce 3,3',5'-triiodothyronine (rT3). Several physiological and pathological stimuli influence thyroid hormone synthesis. The overproduction of thyroid hormones leads to hyperthyroidism, which is treated by antithyroid drugs that either inhibit the thyroid hormone biosynthesis and/or decrease the conversion of T4 to T3. Antithyroid drugs are thiourea-based compounds, which indude propylthiouracil (PTU), methimazole (MM I), and carbimazole (CBZ). The thyroid gland actively concentrates these heterocyclic compounds against a concentration gradient Recently, the selenium analogues of PTU, MMI, and CBZ attracted significant attention because the selenium moiety in these compounds has a higher nucleophilicity than that of the sulfur moiety. Researchers have developed new methods for the synthesis of the selenium compounds. Several experimental and theoretical investigations revealed that the selone (C=Se) in the selenium analogues is more polarized than the thione (C=S) in the sulfur compounds, and the selones exist predominantly in their zwitterionic forms. Although the thionamide-based antithyroid drugs have been used for almost 70 years, the mechanism of their action is not completely understood. Most investigations have revealed that MMI and PTU irreversibly inhibit TPO. PTU, MTU, and their selenium analogues also inhibit ID-1, most likely by reacting with the selenenyl iodide intermediate. The good ID-1 inhibitory activity of Pill and its analogues can be ascribed to the presence of the -N(H)-C(=O)- functionality that can form hydrogen bonds with nearby amino add residues in the selenenyl sulfide state. In addition to the TPO and ID-1 inhibition, the selenium analogues are very good antioxidants. In the presence of cellular reducing agents such as GSH, these compounds catalytically reduce hydrogen peroxide. They can also efficiently scavenge peroxynitrite, a potent biological oxidant and nitrating agent.

Molecular Dynamics Perspective on the Protein Thermal Stability: A Case Study Using SAICAR Synthetase

The enzyme SAICAR synthetase ligates aspartate with CAIR (5'-phosphoribosyl-4-carboxy-5-aminoimidazole) forming SAICAR (5-amino-4-imidazole-N-succinocarboxamide ribonucleotide) in the presence of ATP. In continuation with our previous study on the thermostability of this enzyme in hyper-/thermophiles based on the structural aspects, here, we present the dynamic aspects that differentiate the mesophilic (E. coli, E. chaffeensis), thermophilic (G. kaustophilus), and hyperthermophilic (M. jannaschii, P. horikoshii) SAICAR synthetases by carrying out a total of 11 simulations. The five functional dimers from the above organisms were simulated using molecular dynamics for a period of 50 ns each at 300 K, 363 K, and an additional simulation at 333 K for the thermophilic protein. The basic features like root-mean-square deviations, root-mean-square fluctuations, surface accessibility, and radius of gyration revealed the instability of mesophiles at 363 K. Mean square displacements establish the reduced flexibility of hyper-/thermophiles at all temperatures. At the simulations time scale considered here, the long-distance networks are considerably affected in mesophilic structures at 363 K. In mesophiles, a comparatively higher number of short-lived (having less percent existence time) C alpha, hydrogen bonds, hydrophobic interactions are formed, and long-lived (with higher percentage existence time) contacts are lost. The number of time-averaged salt-bridges is at least 2-fold higher in hyperthermophiles at 363 K. The change in surface accessibility of salt-bridges at 363 K from 300 K is nearly doubled in mesophilic protein compared to proteins from other temperature classes.

Sensitivity of Water Dynamics to Biologically Significant Surfaces of Monomeric Insulin: Role of Topology and Electrostatic Interactions

In addition to the biologically active monomer of the protein insulin circulating in human blood, the molecule also exists in dimeric and hexameric forms that are used as storage. The insulin monomer contains two distinct surfaces, namely, the dimer forming surface (DFS) and the hexamer forming surface (HFS), that are specifically designed to facilitate the formation of the dimer and the hexamer, respectively. In order to characterize the structural and dynamical behavior of interfacial water molecules near these two surfaces (DFS and HFS), we performed atomistic molecular dynamics simulations of insulin with explicit water. Dynamical characterization reveals that the structural relaxation of the hydrogen bonds formed between the residues of DFS and the interfacial water molecules is faster than those formed between water and that of the HFS. Furthermore, the residence times of water molecules in the protein hydration layer for both the DFS and HFS are found to be significantly higher than those for some of the other proteins studied so far, such as HP-36 and lysozyme. In particular, we find that more structured water molecules, with higher residence times (similar to 300-500 ps), are present near HFS than those near DFS. A significant slowing down is observed in the decay of associated rotational auto time correlation functions of O-H bond vector of water in the vicinity of HFS. The surface topography and the arrangement of amino acid residues work together to organize the water molecules in the hydration layer in order to provide them with a preferred orientation. HFS having a large polar solvent accessible surface area and a convex extensive nonpolar region, drives the surrounding water molecules to acquire predominantly an outward H-atoms directed, clathrate-like structure. In contrast, near the DFS, the surrounding water molecules acquire an inward H-atoms directed orientation owing to the flat curvature of hydrophobic surface and the interrupted hydrophilic residual alignment. We have followed escape trajectory of several such quasi-bound water molecules from both the surfaces that reveal the significant differences between the two hydration layers.

Intermolecular interactions, charge-density distribution and the electrostatic properties of pyrazinamide anti-TB drug molecule: an experimental and theoretical charge-density study

An experimental charge-density analysis of pyrazinamide (a first line antitubercular drug) was performed using high-resolution X-ray diffraction data (sin theta/lambda)(max) = 1.1 angstrom(-1)] measured at 100 (2) K. The structure was solved by direct methods using SHELXS97 and refined by SHELXL97. The total electron density of the pyrazinamide molecule was modeled using the Hansen-Coppens multipole formalism implemented in the XD software. The topological properties of electron density determined from the experiment were compared with the theoretical results obtained from CRYSTAL09 at the B3LYP/6-31G** level of theory. The crystal structure was stabilized by N-H center dot center dot center dot N and N-H center dot center dot center dot O hydrogen bonds, in which the N3-H3B center dot center dot center dot N1 and N3-H3A center dot center dot center dot O1 interactions form two types of dimers in the crystal. Hirshfeld surface analysis was carried out to analyze the intermolecular interactions. The fingerprint plot reveals that the N center dot center dot center dot H and O center dot center dot center dot H hydrogen-bonding interactions contribute 26.1 and 18.4%, respectively, of the total Hirshfeld surface. The lattice energy of the molecule was calculated using density functional theory (B3LYP) methods with the 6-31G** basis set. The molecular electrostatic potential of the pyrazinamide molecule exhibits extended electronegative regions around O1, N1 and N2. The existence of a negative electrostatic potential (ESP) region just above the upper and lower surfaces of the pyrazine ring confirm the pi-electron cloud.

PROPENSITIES OF AMINO ACID RESIDUES IN PROTEINS FOR DIFFERENT REGIONS OF THE RAMACHANDRAN MAP

Using a dataset of 1164 crystal structures of largely non-homologous proteins defined at a resolution of 1.5 angstrom or better, we have investigated the (phi,psi) preferences of 20 residue types by considering the residues which occur in loops. Propensities of residue types to occur in the loops with (phi,psi) values in the aa region of the Ramachandran map has a poor correlation coefficient of 0.48 to the Chou-Fasman propensities of the residue types to occur in the a-helical segments. However the correlation coefficient between propensities of residues in loops to adopt beta conformations and those in beta-sheet is much higher (0.95). These observations suggest that a-helix formation is well influenced by the local amino acid sequence while intrinsic preference of residue types for beta-sheet plays a major role in the formation of beta-sheet. The main chain polar groups of residues in loops, that can affect the (phi,psi) values, can be involved in intra-molecular hydrogen bonding. Therefore we investigated further by considering subset of residues in loops with low (0 to 2) number of intra-molecular hydrogen bonds per residue involving main chain polar atoms. For this subset, the correlation coefficients between propensities for alpha-helix and alpha(R) region and between beta-sheet and beta-region are 0.26 and 0.64 respectively. This reiterates higher intrinsic tendency of beta-region favouring residues to adopt beta-sheet than alpha(R) region favouring residues to adopt alpha-helical structure.