Investigation of inter-ion interactions in N,N,N',N'-tetramethylethylenediammonium dithiocyanate via experimental and theoretical charge density studies

The crystal structure of the N,N,N',N'-tetramethylethylenediammonium dithiocyanate salt has been examined by experimental charge density studies from high-resolution X-ray diffraction data. The corresponding results are compared with multipole refinements, using theoretical structure factors obtained from a periodic density functional theory calculation at the B3LYP level with a 6-31G** basis set. The salt crystallizes in space group P (1) over bar and contains only a single ion pair with an inversion center in the cation. The salt has thus one unique classical N+-H center dot center dot center dot(NCS)(-) hydrogen bond but also has six other weaker interactions: four C-H center dot center dot center dot S, one C-H center dot center dot center dot N, and one C-H center dot center dot center dot C-pi. The nature of all these interactions has been examined topologically using Bader's quantum theory of "atoms in molecules" and all eight of the Koch-Popelier criteria. The experimental and theoretical approaches agree well and both show that the inter-ion interactions, even in this simplest of systems, play an integrated and complex role in the packing of the ions in the crystal. Electrostatic potential maps are derived from experimental charge densities. This is the first time such a system has been examined in detail by these methods.

Electromechanical interactions in a carbon nanotube based thin film field emitting diode

Carbon nanotubes (CNTs) have emerged as promising candidates for biomedical x-ray devices and other applications of field emission. CNTs grown/deposited in a thin film are used as cathodes for field emission. In spite of the good performance of such cathodes, the procedure to estimate the device current is not straightforward and the required insight towards design optimization is not well developed. In this paper, we report an analysis aided by a computational model and experiments by which the process of evolution and self-assembly (reorientation) of CNTs is characterized and the device current is estimated. The modeling approach involves two steps: (i) a phenomenological description of the degradation and fragmentation of CNTs and (ii) a mechanics based modeling of electromechanical interaction among CNTs during field emission. A computational scheme is developed by which the states of CNTs are updated in a time incremental manner. Finally, the device current is obtained by using the Fowler–Nordheim equation for field emission and by integrating the current density over computational cells. A detailed analysis of the results reveals the deflected shapes of the CNTs in an ensemble and the extent to which the initial state of geometry and orientation angles affect the device current. Experimental results confirm these effects.

Functional interaction of diphenols with polyphenol oxidase - Molecular determinants of substrate/inhibitor specificity SO FEBS JOURNAL

Polyphenol oxidase (PPO) catalyzes the oxidation of o-diphenols to their respective quinones. The quinones autopolymerize to form dark pigments, an undesired effect. PPO is therefore the target for the development of antibrowning and antimelanization agents. A series of phenolic compounds experimentally evaluated for their binding affinity and inhibition constants were computationally docked to the active site of catechol oxidase. Docking studies suggested two distinct modes of binding, dividing the docked ligands into two groups. Remarkably, the first group corresponds to ligands determined to be substrates and the second group corresponds to reversible inhibitors. Analyses of the complexes provide structural explanations for correlating subtle changes in the position and nature of the substitutions on o-diphenols to their functional properties as substrates and inhibitors. Higher reaction rates and binding are reckoned by additional interactions of the substrates with key residues that line the hydrophobic cavity. The docking results suggest that inhibition of oxidation stems from an interaction between the aromatic carboxylic acid group and the apical His 109 of the four coordinates of the trigonal pyramidal coordination polyhedron of CuA. The spatial orientation of the hydroxyl in relation to the carboxylic group either allows a perfect fit in the substrate cavity, leading to inhibition, or because of a steric clash flips the molecule vertically, facilitating oxidation. This is the first study to explain, at the molecular level, the determinants Of substrate and inhibitor specificity of a catechol oxidase, thereby providing a platform for the design of selective inhibitors useful to both the food and pharmaceutical industries.

Hydrophobic Peptide Channels and Encapsulated Water Wires

Peptide nanotubes with filled and empty pores and close-packed structures are formed in closely related pentapeptides. Enantiomorphic sequences, Boc-(D)Pro-Aib-Xxx-Aib-Val-OMe (Xxx = Leu, 1; Val, 2; Ala, 3; Phe, 4) and Boc-Pro-Aib-(D)Xxx-Aib-(D)Val-OMe ((XXX)-X-D = (D)Leu, 5; (D)Val, 6; (D)Ala, 7; (D)Phe, 8), yield molecular structures with a very similar backbone conformation but varied packing patterns in crystals. Peptides 1, 2, 5, and 6 show tubular structures with the molecules self-assembling along the crystallographic six-fold axis (c-axis) and revealing a honeycomb arrangement laterally (ab plane). Two forms of entrapped water wires have been characterized in 2: 2a with d(O center dot center dot center dot O) = 2.6 angstrom and 2b with d(O center dot center dot center dot O) = 3.5 angstrom. The latter is observed in 6 (6a) also. A polymorphic form of 6 (6b), grown from a solution of methanol-water, was observed to crystallize in a monoclinic system as a close-packed structure. Single-file water wire arrangements encapsulated inside hydrophobic channels formed by peptide nanotubes could be established by modeling the published structures in the cases of a cyclic peptide and a dipeptide. In all the entrapped water wires, each water molecule is involved in a hydrogen bond with a previous and succeeding water molecule. The O-H group of the water not involved in any hydrogen bond does not seem to be involved in an energetically significant interaction with the nanotube interior, a general feature of the one-dimensional water wires encapsulated in hydrophobic environements. Water wires in hydrophobic channels are contrasted with the single-file arrangements in amphipathic channels formed by aquaporins.

How Effective Is The Data On Co-Occurrence Of Domains In Multi-Domain Proteins In Prediction Of Protein-Protein Interactions?

We explore the fuse of information on co-occurrence of domains in multi-domain proteins in predicting protein-protein interactions. The basic premise of our work is the assumption that domains co-occurring in a polypeptide chain undergo either structural or functional interactions among themselves. In this study we use a template dataset of domains in multidomain proteins and predict protein-protein interactions in a target organism. We note that maximum number of correct predictions of interacting protein domain families (158) is made in S. cerevisiae when the dataset of closely related organisms is used as the template followed by the more diverse dataset of bacterial proteins (48) and a dataset of randomly chosen proteins (23). We conclude that use of multi-domain information from organisms closely-related to the target can aid prediction of interacting protein families.

Enhanced Tetrahedral Ordering of Water Molecules in Minor Grooves of DNA: Relative Role of DNA Rigidity, Nanoconfinement, and Surface Specific Interactions

Confinement and Surface specific interactions call induce Structures otherwise unstable at that temperature and pressure. Here we Study the groove specific water dynamics ill the nucleic acid sequences, poly-AT and poly-GC, in long B-DNA duplex chains by large scale atomistic molecular dynamics simulations, accompanied by thermodynamic analysis. While water dynamics in the major groove remains insensitive to the sequence differences, exactly the opposite is true for the minor groove water. Much slower water dynamics observed in the minor grooves (especially in the AT minor) call be attributed to all enhanced tetrahedral ordering (< t(h)>) of water. The largest value of < t(h)> in the AT minor groove is related to the spine of hydration found in X-ray Structure. The calculated configurational entropy (S-C) of the water molecules is found to be correlated with the self-diffusion coefficient of water in different region via Adam-Gibbs relation D = A exp(-B/TSC), and also with < t(h)>.

Supramolecular Organization in Tetra Aqua (mu-8-Hydroxyquinoline-5-sulfonate) Barium (II) and Ag center dot center dot center dot I Interactions in a Pseudopolymorphic Form of (7-Iodo-8-hydroxyquinoline-5-sulfonate) Silver (I) Monohydrate

The complexes, Ba (HQS) (H2O)(4) (HQS = 8-hydroxyquinoline-5-sulfonic acid) (1) and Ag (HIQS) (H2O) (Ferron = 7-iodo-8-hydroxyquinoline-5-sulfonic acid) (2) have been synthesized and characterized by X-ray diffraction analysis and spectroscopic studies. In compound 1, Ba2+ ion has a nine-coordinate monocapped antiprismatic geometry. In compound 2, Ag+ has distorted tetrahedral coordination and Ag center dot center dot center dot I interactions generate the supramolecular architectures. The complexes have been characterized by FT-IR and UV-Visible measurements. In both the structures, the inversion-related organic ligands are stacked over one another leading to three-dimensional networks.

A network representation of protein structures: Implications for protein stability

This study views each protein structure as a network of noncovalent connections between amino acid side chains. Each amino acid in a protein structure is a node, and the strength of the noncovalent interactions between two amino acids is evaluated for edge determination. The protein structure graphs (PSGs) for 232 proteins have been constructed as a function of the cutoff of the amino acid interaction strength at a few carefully chosen values. Analysis of such PSGs constructed on the basis of edge weights has shown the following: 1), The PSGs exhibit a complex topological network behavior, which is dependent on the interaction cutoff chosen for PSG construction. 2), A transition is observed at a critical interaction cutoff, in all the proteins, as monitored by the size of the largest cluster (giant component) in the graph. Amazingly, this transition occurs within a narrow range of interaction cutoff for all the proteins, irrespective of the size or the fold topology. And 3), the amino acid preferences to be highly connected (hub frequency) have been evaluated as a function of the interaction cutoff. We observe that the aromatic residues along with arginine, histidine, and methionine act as strong hubs at high interaction cutoffs, whereas the hydrophobic leucine and isoleucine residues get added to these hubs at low interaction cutoffs, forming weak hubs. The hubs identified are found to play a role in bringing together different secondary structural elements in the tertiary structure of the proteins. They are also found to contribute to the additional stability of the thermophilic proteins when compared to their mesophilic counterparts and hence could be crucial for the folding and stability of the unique three-dimensional structure of proteins. Based on these results, we also predict a few residues in the thermophilic and mesophilic proteins that can be mutated to alter their thermal stability.

Exploring polymorphism by solvent mediation in potentially active herbicide Metribuzin: A subtle interplay of weak intermolecular interactions

Metribuzin, 4-amino-6-tert-butyl-3-methylthio- 1,2,4-triazin-5-one, exhibits polymorphic behaviour, crystallizing as plates and needles, driven by variation in solvent polarity, a delicate balance of weak intermolecular forces generating different molecular assemblies.

Peptide-lanthanide interactions. Crystal structure of a europium(III)-triglycine complex

Crystals of Eu-(Gly-Gly-Gly).(H2O)5.(ClO4)3 are triclinic, spacegroup P1BAR with a = 9.123 (2), b = 11.185 (5), c = 11.426 (2) angstrom; alpha = 90.79 (2), beta = 98.08 (1), gamma = 98.57 (2)-degrees; Z = 2. The europium cation is surrounded by four oxygens from three different peptide units and four oxygens from water molecules. The geometry around the metal is a distorted bi-capped trigonal prism. The peptide backbone conformation in this complex is compared with those in the free peptide and in various metal complexes. Considerable differences are observed between Eu(III) and Ca(II) complexes of triglycine. (C) Munksgaard 1994.

Interactions of linear dicationic molecules with lipid A: structural requisites for optimal binding affinity

The structural determinants of the binding affinity of linear dicationic molecules toward lipid A have been examined with respect to the distance between the terminal cationic functions, the basicity, and the type of cationic moieties using a series of spermidine derivatives and pentamidine analogs by fluorescence spectroscopic methods, The presence of two terminal cationic groups corresponds to enhanced affinity, A distinct sigmoidal relationship between the intercationic distance and affinity was observed with a sharp increase at 11 Angstrom, levelling off at about 13 Angstrom. The basicity (pK) and nature of the cationic functions are poor correlates of binding potency, since molecules bearing primary amino, imidazolino, or guanido termini are equipotent, The interaction of pentamidine, a bisamidine drug, with lipid A, characterized in considerable detail employing the putative intermolecular excimerization of the drug, suggests a stoichiometry of 1:1 in the resultant complex, The binding is driven almost exclusively by electrostatic forces, and is dependent on the ionization states of both lipid A and the drug, Under conditions when lipid A is highly disaggregated, pentamidine binds specifically to bis-phosphoryl- but not to monophosphoryl-lipid A indicating that both phosphate groups of lipid A are necessary for electrostatic interactions by the terminal amidininium groups of the drug, Based on these data, a structural model is proposed for the pentamidine-lipid A complex, which may be of value in designing endotoxin antagonists from first principles.

Charge density based classification of intermolecular interactions in molecular crystals

Evaluation of intermolecular interactions in terms of both experimental and theoretical charge density analyses has produced a unified picture with which to classify strong and weak hydrogen bonds, along with van der Waals interactions, into three regions.

The role of polar and facial amphipathic character in determining lipopolysaccharide-binding properties in synthetic cationic peptides

Two series of peptides, designated K and NK were synthesized and tested for lipid A binding and neutralizing properties. K-2, which has an 11-residue amphiphilic core, and a branched N-terminus bearing two branched lysinyl residues does not bind lipid A, while NK2, also with an 11-residue amphiphilic core comprised entirely of non-ionizable residues, and a similarly branched, cationic N-terminus, binds lipid A very weakly. Both peptides do not inhibit lipopolysaccharide (LPS) activity in the Limulus assay, nor do they inhibit LPS-induced TNF-alpha and NO production in 5774 cells. These results are entirely unlike a homologous peptide with an exclusively hydrophobic core whose LPS-binding and neutralizing properties are very similar to that of polymyxin B [David SA, Awasthi SK, Wiese A et al. Characterization of the interactions of a polycationic, amphiphilic, terminally branched oligopeptide with lipid A and lipopolysaccharide from the deep rough mutant of Salmonella minnesota. J Endotoxin Res 1996; 3: 369-379]. These data suggest that a clear segregation of charged and apolar domains is crucial in molecules designed for purposes of LPS sequestration and that head-tail (polar) orientation of the cationic/hydrophobic regions is preferable to molecules with mixed or facial cationic/amphipathic character.

Evaluation of solute-solvent interactions from solvent blue-shifts of n → π* transitions of C=O, C=S, NO2 and N=N groups: hydrogen bond energies of various donor-acceptor systems

Effects of non-polar, polar and proton-donating solvents on the n → π* transitions of C=O, C=S, NO2 and N=N groups have been investigated. The shifts of the absorption maxima in non-polar and polar solvents have been related to the electrostatic interactions between solute and solvent molecules, by employing the theory of McRAE. In solvents which can donate protons the solvent shifts are mainly determined by solute-solvent hydrogen bonding. Isobestic points have been found in the n → π* bonds of ethylenetrithio-carbonate in heptane-alcohol and heptane-chloroform solvent systems, indicating the existence of equilibria between the hydrogen bonded and the free species of the solute. Among the different proton-donating solvents studied water produces the largest blue-shifts. The blue-shifts in alcohols decrease in the order 2,2,2-trifluoroethanol, methanol, ethanol, isopropanol and t-butanol, the blue-shift in trifluoroethanol being nearly equal to that in water. This trend is exactly opposite to that for the self-association of alcohols. It is suggested that electron-withdrawing groups not merely decrease the extent of self-association of alcohols, but also increase the ability to donate hydrogen bonds. The approximate hydrogen-bond energies for several donor-acceptor systems have been estimated. In a series of aliphatio ketones and nitro compounds studied, the blue-shifts and consequently the hydrogen bond energies decrease with the decrease in the electron-withdrawing power of the alkyl groups. It is felt that electron-withdrawing groups render the chromophores better proton acceptors, and the alcohols better donors. A linear relationship between n → π* transition frequency and the infrared frequency of ethylenetrithiocarbonate has been found. It is concluded that stabilization of the electronic ground states of solute molecules by electrostatic and/or hydrogen-bond interactions determines the solvent shifts.

Potential functions for hydrogen bond interactions - III. Empirical Potential Function for the Peptide N-H...O=C Hydrogen Bond

A careful comparison of the distribution in the (R, θ)-plane of all NH ... O hydrogen bonds with that for bonds between neutral NH and neutral C=O groups indicated that the latter has a larger mean R and a wider range of θ and that the distribution was also broader than for the average case. Therefore, the potential function developed earlier for an average NH ... O hydrogen bond was modified to suit the peptide case. A three-parameter expression of the form {Mathematical expression}, with △ = R - Rmin, was found to be satisfactory. By comparing the theoretically expected distribution in R and θ with observed data (although limited), the best values were found to be p1 = 25, p3 = - 2 and q1 = 1 × 10-3, with Rmin = 2·95 Å and Vmin = - 4·5 kcal/mole. The procedure for obtaining a smooth transition from Vhb to the non-bonded potential Vnb for large R and θ is described, along with a flow chart useful for programming the formulae. Calculated values of ΔH, the enthalpy of formation of the hydrogen bond, using this function are in reasonable agreement with observation. When the atoms involved in the hydrogen bond occur in a five-membered ring as in the sequence[Figure not available: see fulltext.] a different formula for the potential function is needed, which is of the form Vhb = Vmin +p1△2 +q1x2 where x = θ - 50° for θ ≥ 50°, with p1 = 15, q1 = 0·002, Rmin = 2· Å and Vmin = - 2·5 kcal/mole. © 1971 Indian Academy of Sciences.