A comprehensive conformational analysis of tryptophan,its ionic and dimeric forms

Tryptophan is an essential amino acid, and understanding the conformational preferences of monomer and dimer is a subject of outstanding relevance in biological systems. An exhaustive first principles investigation of tryptophan ( W ) and its ionized counterparts cations (WC) , anions (WA) , and zwitterions (WZ) has been carried out. A comprehensive and systematic study of tryptophan dimer (WD) conformations resulted in about 62 distinct minima on the potential energy surface. The hydrogen bonds and a variety of noncovalent interactions such as OH‐π, NH‐π, CH‐π, CH‐O, and π‐π interactions stabilized different forms of tryptophan and its dimers. Over all in monomeric conformers which have NH‐O, hydrogen bonds showed higher stability than other conformers. A cursory analysis reveal that the most stable dimers stabilized by hydrogen bonding interactions while the less stable dimers showed aromatic side chain interactions. Protein Data Bank analysis of tryptophan dimers reveals that at a larger distance greater than 5 Å, T‐shaped orientations (CH‐π interactions) are more prevalent, while stacked orientations (π‐π interactions) are predominant at a smaller distance. © 2013 Wiley Periodicals, Inc.     

DFT Study on the Effect of Hydrogen-bond Formation on the Adsorption of Aminotriazines on Graphene

DFT calculations have been performed to explore the aminotriazine adsorption on graphene surfaces.Relative energies,equilibrium geometries and electronic structures of monomer and dimer of aminotriazine molecules adsorbed at the surface were investigated and analyzed in details.It was found that the hydrogen atoms in the NH2 group of aminotriazine molecules are directed toward the graphene surface,and the adsorption energy increases as the NH2 group is added.The adsorbed aminotriazine molecules facilely form a dimer through the hydrogen bonding interactions,and the two aromatic rings of optimized structure of 2-amino-1,3,5-triazine(B) dimmer(denoted by B2) and melamine(D) dimmer(denoted by D2) are parallel to the graphene sheet.The large deviation of the averaged adsorption energy of B2 and D2 compared to monor adsorption may reflect the increase of π-π repulsion and the effect of hydrogen bond formation.The electronic structure analyses reveal that the formation of hydrogen bonds in melamine dimer has great influence on the adsorption mode at the graphene surface.     

NTO二聚体分子间相互作用的理论研究

在DFT-B3LYP/6-311＋＋G**水平上求得NTO二聚体势能面上六种优化构型和电子结构. 经基组叠加误差(BSSE)和零点能(ZPE)校正, 求得分子间最大相互作用能为－53.66 kJ/mol. 二子体系间的电荷转移很少. 由自然键轨道分析揭示了相互作用的本质. 对优化构型进行振动分析, 并基于统计热力学求得200.0～800.0 K温度范围从单体形成二聚体的热力学性质变化. 发现二聚主要由强氢键所贡献, 但结合能大小并不为氢键所完全决定. 二聚过程在较低温度或常温下能自发进行.     

Force-induced insulin dimer dissociation: a molecular dynamics study

Understanding the forces and dynamics of insulin dissociation is critical for devising formulations for the treatment of insulin-dependent diabetes. In earlier work, we applied AFM-based force spectroscopy to covalently tethered and oriented insulin monomers to assess the effect of molecular orientation on insulin-insulin binding forces. We report here on steered molecular dynamics simulations of the insulin dissociation force spectroscopy experiment. Consistent with our experiments, our simulation results suggest that insulin dimer dissociation occurs near the limit of extensibility of the B-chain. We have also found that the forced dissociation of the insulin dimer is a rate-dependent process, involving significant conformational changes to the monomer(s). The insulin dimer dissociation pathway also depends on the relative strength of the inter-monomer interactions across the antiparallel beta-sheet interface and the intra-monomer interactions of residues A1 and A30 with the insulin B-chain. Our simulation results strongly support the design of bioactive insulin analogues that involves altering hydrogen bonding and hydrophobic interactions across the beta-sheet dimer interface.     

Transformation between α‐helix and β‐sheet structures of one and two polyglutamine peptides in explicit water molecules by replica‐exchange molecular dynamics simulations

Aggregation of polyglutamine peptides with β‐sheet structures is related to some important neurodegenerative diseases such as Huntington's disease. However, it is not clear how polyglutamine peptides form the β‐sheets and aggregate. To understand this problem, we performed all‐atom replica‐exchange molecular dynamics simulations of one and two polyglutamine peptides with 10 glutamine residues in explicit water molecules. Our results show that two polyglutamine peptides mainly formed helix or coil structures when they are separated, as in the system with one‐polyglutamine peptide. As the interpeptide distance decreases, the intrapeptide β‐sheet structure sometimes appear as an intermediate state, and finally the interpeptide β‐sheets are formed. We also find that the polyglutamine dimer tends to form the antiparallel β‐sheet conformations rather than the parallel β‐sheet, which is consistent with previous experiments and a coarse‐grained molecular dynamics simulation. © 2014 Wiley Periodicals, Inc.     

Virtual‐system‐coupled adaptive umbrella sampling to compute free‐energy landscape for flexible molecular docking

A novel enhanced conformational sampling method, virtual‐system‐coupled adaptive umbrella sampling (V‐AUS), was proposed to compute 300‐K free‐energy landscape for flexible molecular docking, where a virtual degrees of freedom was introduced to control the sampling. This degree of freedom interacts with the biomolecular system. V‐AUS was applied to complex formation of two disordered amyloid‐β (Aβ_30–35) peptides in a periodic box filled by an explicit solvent. An interpeptide distance was defined as the reaction coordinate, along which sampling was enhanced. A uniform conformational distribution was obtained covering a wide interpeptide distance ranging from the bound to unbound states. The 300‐K free‐energy landscape was characterized by thermodynamically stable basins of antiparallel and parallel β‐sheet complexes and some other complex forms. Helices were frequently observed, when the two peptides contacted loosely or fluctuated freely without interpeptide contacts. We observed that V‐AUS converged to uniform distribution more effectively than conventional AUS sampling did. © 2015 Wiley Periodicals, Inc.     

Crystal structure of nylon 12

The crystal structure of nylon 12 prepared by polymerization of dodecalactam has been determined by x-ray diffraction. Nylon 12 fiber exhibits only the γ form as its stable crystal structure. The unit cell of nylon 12 was determined with the aid of the x-ray diffraction pattern of a doubly oriented specimen. The unit cell is monoclinic with a = 9.38 Å, b = 32.2 Å (fiber axis), c = 4.87 Å and β = 121.5° and contains four repeating monomer units. The chain is planar zigzag for the most part but is twisted at the position of amide groups, forming hydrogen bonds between neighboring parallel chains. The chain conformation is similar to that of the γ form of nylon 6 proposed by Arimoto. It was deduced from the calculations that there are two chain conformations statistically coexistent according to the direction of twisting. In each conformation, hydrogen bonds are formed between parallel chains to make pleated sheetlike structures. The sheets are nearly parallel to (200) and in the sheet the directions of the neighboring chains are antiparallel, as is the case with nylon 6.     

A density-functional theory investigation of 3-nitro-1,2,4-triazole-5-one dimers and crystal

Density-functional method with different basis sets was applied to the study of the highly efficient and low sensitive explosive 3-nitro-1,2,4-triazole-5-one (NTO) in both gaseous dimer and its bulk state. The binding energies have been corrected for the basis set superposition errors. Six stable dimers (II-VII) were located. The corrected binding energy of the most stable dimer VII is predicted to be -53.66 kJ/mol at the B3LYP/6-311++G(**) level. It was found that the structures of the more stable dimers (V-VII) are through the hydrogen bonding interaction between the carbonyl oxygen and the azole hydrogen of 3-nitro-1,2,4-triazole-5-one. The changes of Gibbs free energies (DeltaG) in the processes from the monomer to the dimers at 298.15 K are 8.51, 0.90, 0.35, -8.74, -10.67, and -11.06 kJ/mol for dimers from II to VII, respectively. Dimers V-VII, possessing cyclic structures, can be spontaneously produced from the isolated monomer at room temperature. The lattice energy is -156.14 kJ/mol, and this value becomes to -150.43 kJ/mol when a 50% correction of the basis set superposition error was adopted. The frontier bands are quite flat. Judged from the value of band gap of 4.0 eV, it may be predicted that 3-nitro-1,2,4-triazole-5-one is an insulator. Most atoms in NTO, with the exception of C(5) atom and the nitro atoms, make up the upper valence bands. In contrast, the lower conduction bands mainly consist of the nitro N and O atoms. The population of the C-NO(2) bond is much less than those of the other bonds and the detonation may be initiated by the breakdown of this bond.     

Structure and matrix isolation infrared spectrum of formyl fluoride dimer: blue-shift of the C-H stretching frequency

Infrared spectroscopy (IR) of formyl fluoride (HCOF) dimer is studied in low-temperature argon and krypton matrixes. New IR absorptions, ca. 17 cm(-1) blue shifted from the monomer C-H stretching fundamental, are assigned to the HCOF dimer. The MP2/6-311++G calculations were utilized to define structures and harmonic frequencies of various HCOF dimers. Among the four optimized structures, the dimer having two C-H...O hydrogen bonds possesses strongest intermolecular bonding. The calculated harmonic frequencies of this dimer structure are shifted from the monomer similarly as observed in the experiment. Thus, we suggest that the experimentally observed blue shifted C-H bands belong to the dimer with two C-H...O hydrogen bonds. This observation includes the HCOF dimer to the class of hydrogen bonded complexes showing blue shift in their vibrational energies.     

How Important Is Molecular Rigidity for the Complex Stability of Artificial Host–Guest Systems? A Theoretical Study on Self‐Assembly of Gas‐Phase Arginine

Arginine forms much less stable dimers than 2-(guanidiniocarbonyl)-1H-pyrrole-5-carboxylate although the principal binding interactions are very similar. The reasons for this difference are addressed in this work by state-of-the-art ab initio computations. The investigation shows that the extraordinary high stability of the 2-(guanidiniocarbonyl)-1H-pyrrole-5-carboxylate dimer results to about 50 % from the rigidity of its monomer. Within this study monomer and dimer conformers of arginine were calculated leading to new low lying structures which have not been reported before as well as new global minima are predicted. In these structures stacking interactions with the guanidinium moiety are especially important. For the monomer we predict the energy minimum to be the canonical form with the lowest lying zwitterionic structure being only 9 kJ mol(-1) less stable. During the course of these calculations we found that DFT did not predict the structures and their relative energy correctly in comparison to perturbation theory (MP2) and some potential reasons for the failure of DFT in these cases are discussed. Vibrational frequencies of the various structures are presented and a suitable wavenumber region for an experimental determination of the global minimum of the arginine monomer is identified. The effect of molecular rigidity on the self-assembly is probed using a local minimum of the arginine monomer which does not possess any intramolecular stabilizing effects. Our results suggest that the deliberate control of the conformational flexibility is a powerful instrument to steer the complex affinity of artificial hosts.     

Theoretical scrutinization of nine benzoic acid dimers: Stability and energy decomposition analysis

Aromatic carboxylic acids are able to form diverse dimers and multimers due to their hydrogen bond donor and acceptor cites, as well as the aromatic rings. In this work, we examine nine benzoic acid dimers stabilized by hydrogen bonding and stacking interactions. Interacting quantum atoms methodology revealed that dominant attractive interactions in all of them, including hydrogen bonded systems, are due to exchange-correlation. Coulomb interactions are significant only in the most stable dimer with a double hydrogen bond, although the corresponding energy term is almost two times lower compared to the nonclassical one. Since interacting quantum atoms approach treats monomers binding by considering electronic energy only, in order to examine dissociation kinetics we performed density functional theory-based molecular dynamics simulations of selected stacked dimers: in 40% of the studied systems at 300 K thermal energy was sufficient to overpower barrier for dissociation within 1 ps, which resulted in the separation of the monomers, whereas 20% of them remained in the stacked position even after 5 ps. These results highlight the importance of noncovalent interactions, particularly weak stacking interactions, on the structure and dynamics of carboxylic acids and their derivatives.     

Geometries and stabilities of l-ascorbic acid dimer and its derivatives: a computational investigation by density functional methods

Geometries, relative stabilities, and hydrogen bonds of l-ascorbic acid (LAA) and d-Erythroascorbate (DEAA) dimers as well as their S- and Se-substituted isomers in gas phase and water solvent are studied using density functional method. Furthermore, the hydrogen bond lengths in LAA and DEAA dimers are generally increased along with the binding dissociation energy of the dimers being decreased as apex O atoms in the five-membered C5 rings of LAA and DEAA dimers are substituted by S and Se atoms in gas phase and water solvent. Interestingly, one LAA dimer and its S- or Se-substituted isomer with four hydrogen bonds in gas and water solvent are the three-centers structures. In addition, the chemical bonding and charge distributions of all the dimers are discussed. A good agreement with available experimental results is reached.     

Dimerization of the keto tautomer of acetohydroxamic acid-infrared matrix isolation and theoretical study

Dimerization of the keto tautomer of acetohydroxamic acid has been studied using FTIR matrix isolation spectroscopy and DFT(B3LYP)/6-31+G(d,p) calculations. Analysis of CH3CONHOH/Ar matrix spectra indicates formation of two dimers in which two intramolecular CO...HON bonds within two interacting acetohydroxamic acid molecules are retained. A chain dimer I is stabilized by the intermolecular CO...HN hydrogen bond, whereas the cyclic dimer II is stabilized by two intermolecular NH...O(H)N bonds. Twelve vibrations were identified for dimer I and six vibrations for dimer II; the observed frequency shifts show a good agreement with the calculated ones for the structures I and II. Both dimers have comparable binding energies (DeltaE(ZPE)(CP)I, II=-7.02, -6.34 kcal mol-1) being less stable than calculated structures III and IV (DeltaE(ZPE)(CP)III, IV=-9.50, -8.87 kcal mol-1) in which one or two intramolecular hydrogen bonds are disrupted. In the most stable 10-membered cyclic dimer III, two intermolecular CO...HON hydrogen bonds are formed at expense of intramolecular hydrogen bonds of the same type. The formation of the less stable (AHA)2 dimers in the studied matrixes indicates that the formation of (AHA)2 is kinetically and not thermodynamically controlled.     

Noncovalent interactions: defining cooperativity. Ligand binding aided by reduced dynamic behavior of receptors. Binding of bacterial cell wall analogues to ristocetin A

Changes in the relative populations of the monomer and asymmetric dimer forms of ristocetin A, upon binding of two molecules of ligand, suggest that ligand binding is negatively cooperative with respect to dimerization. However, strong hydrogen bonds formed in the binding sites of the ligands are reinforced in the dimer relative to the monomer, and the barrier to dissociation of the dimer is increased upon binding of the ligands. It is concluded that the interactions which are common in the binding of both ligands are made with positive cooperativity with respect to those involved in dimerization. The conclusions are relevant to the binding of ligands to proteins, where ligand binding energy can be derived from stabilization of the protein in its ligand-bound form.     

The d'--d--d' vertical triad is less discriminating than the a'--a--a' vertical triad in the antiparallel coiled-coil dimer motif

Elucidating relationships between the amino-acid sequences of proteins and their three-dimensional structures, and uncovering non-covalent interactions that underlie polypeptide folding, are major goals in protein science. One approach toward these goals is to study interactions between selected residues, or among constellations of residues, in small folding motifs. The α-helical coiled coil has served as a platform for such studies because this folding unit is relatively simple in terms of both sequence and structure. Amino acid side chains at the helix-helix interface of a coiled coil participate in so-called "knobs-into-holes" (KIH) packing whereby a side chain (the knob) on one helix inserts into a space (the hole) generated by four side chains on a partner helix. The vast majority of sequence-stability studies on coiled-coil dimers have focused on lateral interactions within these KIH arrangements, for example, between an a position on one helix and an a' position of the partner in a parallel coiled-coil dimer, or between a--d' pairs in an antiparallel dimer. More recently, it has been shown that vertical triads (specifically, a'--a--a' triads) in antiparallel dimers exert a significant impact on pairing preferences. This observation provides impetus for analysis of other complex networks of side-chain interactions at the helix-helix interface. Here, we describe a combination of experimental and bioinformatics studies that show that d'--d--d' triads have much less impact on pairing preference than do a'--a--a' triads in a small, designed antiparallel coiled-coil dimer. However, the influence of the d'--d--d' triad depends on the lateral a'--d interaction. Taken together, these results strengthen the emerging understanding that simple pairwise interactions are not sufficient to describe side-chain interactions and overall stability in antiparallel coiled-coil dimers; higher-order interactions must be considered as well.     

Dissociation energies of deoxyribose nucleotide dimer anions measured using blackbody infrared radiative dissociation

The dissociation kinetics of deprotonated deoxyribose nucleotide dimers were measured using blackbody infrared radiative dissociation. Experiments were performed with noncovalently bound dimers of phosphate, adenosine (dAMP), cytosine (dCMP), guanosine (dGMP), thymidine (dTMP), and the mixed dimers dAMP.dTMP and dGMP.dCMP. The nucleotide dimers fragment through two parallel pathways, resulting in formation of the individual nucleotide or nucleotide + HPO3 ion. Master equation modeling of this kinetic data was used to determine threshold dissociation energies. The dissociation energy of (dGMP.dCMP-H)- is much higher than that for the other nucleotide dimers. This indicates that there is a strong interaction between the nucleobases in this dimer, consistent with the existence of Watson-Crick hydrogen bonding between the base pairs. Molecular mechanics simulations indicate that Watson-Crick hydrogen bonding occurs in the lowest energy structures of (dGMP.dCMP-H)-, but not in (dAMP.dTMP-H)-. The trend in gas phase dissociation energies is similar to the trend in binding energies measured in nonaqueous solutions within experimental error. Finally, the acidity ordering of the nucleotides is determined to be dTMP < dGMP < dCMP < dAMP, where dAMP has the highest acidity (largest delta Gacid).     

Ab initio calculations of structures and interaction energies of toluene dimers including CCSD(T) level electron correlation correction

The intermolecular interaction energy of the toluene dimer has been calculated with the ARS-F model (a model chemistry for the evaluation of intermolecular interaction energy between ARomatic Systems using Feller's method), which was formerly called as the AIMI model III. The CCSD(T) (coupled cluster calculations with single and double substitutions with noniterative triple excitations) interaction energy at the basis set limit has been estimated from the second-order Moller-Plesset perturbation interaction energy at the basis set limit obtained by Feller's method and the CCSD(T) correction term obtained using a medium-size basis set. The cross (C(2)) dimer has the largest (most negative) interaction energy (-4.08 kcal/mol). The antiparallel (C(2h)) and parallel (C(S)) dimers (-3.77 and -3.41 kcal/mol, respectively) are slightly less stable. The dispersion interaction is found to be the major source of attraction in the toluene dimer. The dispersion interaction mainly determines the relative stability of the stacked three dimers. The electrostatic interaction of the stacked three dimers is repulsive. Although the T-shaped and slipped-parallel benzene dimers are nearly isoenergetic, the stacked toluene dimers are substantially more stable than the T-shaped toluene dimer (-2.62 kcal/mol). The large dispersion interaction in the stacked toluene dimers is the cause of their enhanced stability.     

In-silico exploration of noble gas dimer enforced by noncovalent interaction

Density functional calculations have been carried out to investigate the possibility of trapping of noble gas dimers by cyclo[18]carbon dimer. Parallel-displaced conformation of the cyclo[18]carbon dimer is found to be the minimum energy structure. Noncovalent interaction is found to hold the noble gas dimers. The lighter noble gases (He, Ne) posses weaker attractive interactions while the heavier one (Ar, Kr) are held by stronger attractive interactions forming genuine bonds. Each of the noble gas atoms in turn forms noncovalent interaction with the cyclo[18]carbon monomers. The bond dissociation energy of the noble gas dimers dramatically increases inside the cyclo[18]carbon dimer. Energy decomposition analysis reveals that dispersion plays the major role toward the stabilization energy.     

Theoretical Studies on Intermolecular Interactions of 4-Amino-5-nitro-1,2,3-triazole Dimers

1 INTRODUCTION Triazole, a five-membered heterocyclic compound with three nitrogen atoms, is an important interme- diate product of medicine and chemical industry as well as insecticide [1]. Due to its small volume and high nitrogenous density, triazole holds more and more attraction for the material researchers, espe- cially the researchers of high-energy insensitive explosive. It is reported that its nitro and amino deri- vatives are a sort of important high-energy mate- rials[2]. Previ…