Calculations on crystal packing of a flexible molecule, Leu-enkephalin

Crystal packing calculations have been carried out on a substantial number of conformations of Leu-enkephalin; namely, those obtained both from crystal structures and from energy minimizations on isolated molecules, and with and without waters of crystallization. The known crystal structures represent the most energetically stable packings found. The conformations of the enkephalin molecules in the crystal are not the most stable for an isolated molecule; i.e. intermolecular interactions force the isolated molecule to change conformation in order to achieve a small packing volume and an optimal packing energy in the crystal. It is found that the packing energy of an enkephalin molecule is a reasonably smooth function of its molecular volume in the unit cell, if structures with intermolecular hydrogen bonding are excluded, and is substantially independent of other details of the molecular conformation or of the crystal packing. Hydrogen bonding provides additional stabilization of the crystal structure, and would likely permit crystallization of the system if it is sufficiently dense. Solvent molecules further stabilize the structure when they can also provide intermolecular hydrogen bonds.     

Hydrogen-bonding interactions in adrenaline–water complexes: DFT and QTAIM studies of structures, properties, and topologies

ωB97XD/6-311++G(d,p) calculations were carried out to investigate the hydrogen-bonding interactions between adrenaline (Ad) and water. Six Ad-H(2)O complexes possessing various types of hydrogen bonds (H-bonds) were characterized in terms of their geometries, energies, vibrational frequencies, and electron-density topology. Natural bond orbital (NBO) and quantum theory of atoms in molecules (QTAIM) analyses were performed to elucidate the nature of the hydrogen-bonding interactions in these complexes. The intramolecular H-bond between the amino and carboxyl oxygen atom of Ad was retained in most of the complexes, and cooperativity between the intra- and intermolecular H-bonds was present in some of the complexes. H-bonds in which hydroxyls of Ad/water acted as proton donors were stronger than other H-bonds. Both hydrogen-bonding interactions and structural deformation play important roles in the relative stabilities of the complexes. The intramolecular H-bond was broken during the formation of the most stable complex, which indicates that Ad tends to break the intramolecular H-bond and form two new intermolecular H-bonds with the first water molecule.     

Evaluation of free energy landscape for base-amino acid interactions using ab initio force field and extensive sampling

Structural data of protein-DNA complex show redundancy and flexibility in base-amino acid interactions. To understand the origin of the specificity in protein-DNA recognition, we calculated the interaction free energy, enthalpy, entropy, and minimum energy maps for AT-Asn, GC-Asn, AT-Ser, and GC-Ser by means of a set of ab initio force field with extensive conformational sampling. We found that the most preferable interactions in these pairs are stabilized by hydrogen bonding, and are mainly enthalpy driven. However, minima in the free energy maps are not necessarily the same as those in the minimum energy map or enthalpy maps, due to the entropic effect. The effect of entropy is particularly important in the case of GC-Asn. Experimentally observed structures of base-amino acid interactions are within preferable regions in the calculated free energy maps, where there are many different interaction configurations with similar energy. The full geometry optimization procedure using ab initio molecular orbital method was applied to get the optimal interaction geometries for AT-Asn, GC-Asn, AT-Ser, and GC-Ser. We found that there are various base-amino acid combinations with similar interaction energies. These results suggest that the redundancy and conformational flexibility in the base-amino acid interactions play an important role in the protein-DNA recognition.     

The Calculation of 3D Solubility Parameters Using Molecular Models

Abstract

In this paper we describe the use of molecular mechanics models to examine detailed intermolecular interactions within the liquid state of a common nonionic surfactant system, nonyl phenol ethoxylate (NPE). Using constant energy molecular dynamics simulations we have studied the relative strengths of dispersive interactions versus polar interactions and have estimated three dimensional solubility parameters for NPE systems as a function of temperature and ethylene oxide content. The predictions at 300 K are in good agreement with three dimensional solubility parameters predicted using group contribution tables. Models of the amorphous liquid state were represented by single molecular structures of NPE in a periodic cell. The solubility parameters predicted with these models were in good agreement with those values derived from models having eight NPE molecules packed into a cell with the exception of the electrostatic interactions, which are the most sensitive to system size effects.     

Design and conformational analysis of peptoids containing N-hydroxy amides reveals a unique sheet-like secondary structure

N-hydroxy amides can be found in many naturally occurring and synthetic compounds and are known to act as both strong proton donors and chelators of metal cations. We have initiated studies of peptoids, or N-substituted glycines which contain N-hydroxy amide side chains to investigate the potential effects of these functional groups on peptoid backbone amide rotamer equilibria and local conformations. We reasoned that the propensity of these functional groups to participate in hydrogen bonding could be exploited to enforce intramolecular or intermolecular interactions that yield new peptoid structures. Here, we report the design, synthesis, and detailed conformational analysis of a series of model N-hydroxy peptoids. These peptoids were readily synthesized, and their structures were analyzed in solution by 1D and 2D NMR and in the solid-state by X-ray crystallography. The N-hydroxy amides were found to strongly favor trans conformations with respect to the peptoid backbone in chloroform. More notably, unique sheet-like structures held together via intermolecular hydrogen bonds were observed in the X-ray crystal structures of an N-hydroxy amide peptoid dimer, which to our knowledge represent the first structure of this type reported for peptoids. These results suggest that the N-hydroxy amide can be utilized to control both local backbone geometries and longer-range intermolecular interactions in peptoids, and represents a new functional group in the peptoid design toolbox.     

Effects of pD on ring-stacking interactions between dinucleoside phosphates and tryptamine

Complex formation between tryptamine and mononucleotides and dinucleoside phosphates containing adenine and/or cytosine has been studied at five pD's ranging from 1.1 to 7.4 by proton magnetic resonance spectroscopy. Chemical shifts of base ring protons and the ribose anomeric proton in the nucleotides and indole ring protons in tryptamine were monitored and their changes with pD and intermolecular interactions interpreted qualitatively. Stacked complexes were found to exist at all pD's in the range studied. Complex geometries differ depending on pD. An electrostatic interaction between the tryptamine amino group and the nucleotide phosphate group contributes to complex formation above pD 4 but is not strong enough to shift the dinucleoside phosphate equilibrium towards the unstacked conformer.     

Towards an understanding of the arginine-aspartate interaction.

We have made a comparison of the geometries of intra- and intermolecular arginine-aspartate interactions by extracting orientation information from protein co-ordinate data. The results show a pronounced difference, with both types of interaction preferring to form twin N-H . . . O = C hydrogen bonds, but involving different nitrogen atoms. In intramolecular interactions, the aspartate favours a "side on" geometry, forming hydrogen bonds with N epsilon and N eta 2; in the intermolecular case, however, "end on" contacts involving N eta 1 and N eta 2 of the arginine are preferred. We have used Distributed Multipole Analysis of the methylguanidinium-acetate system to model the electrostatic component of the arginine-aspartate ion pair interaction in vacuo. We find, in agreement with the experimental arginine-aspartate distribution, that side on and end on doubly N-H . . . O = C hydrogen-bonded configurations are clearly the most favourable, with the side on being marginally lower in energy. Thus, despite the many competing side-chain interactions in proteins, many arginine-aspartate pairs adopt one of the minimum electrostatic energy conformations, or one close to a minimum. Within each of the two regions (side on and end on) we find only a small energy gap between the "symmetric" doubly hydrogen-bonded and slightly displaced "staggered" structures, again in agreement with the crystal structure data. Further calculations of the total ab initio interaction energy show that this follows the electrostatic term in its orientational variation, this phenomenon of "electrostatic domination" being well known in hydrogen-bonded systems. The end on arginine nitrogen atoms are observed to be more surface-exposed than N epsilon, as demonstrated by their greater accessibilities over a large sample of proteins. This helps explain the side on and end on preferences of intra- and intermolecular interactions, respectively. We also note the effect of short sequence intervals, particularly i in equilibrium with i + 2 relationships, in forcing many intramolecular contacts to be side on.     

DNA polymorphism and local variation in base-pair orientation: a theoretical rationale.

Basepair stacking calculations have been carried out to understand the conformational polymorphism of DNA and its sequence dependence. The recently developed self-consistent parameter set, which is specially suitable for describing irregular DNA structures, has been used to describe the geometry of a basepair doublet. While for basepairs without any propeller, the favourable stacking patterns do not appear to have very strong features, much more noticeable sequence dependent stacking patterns emerge once a propeller is applied to the basepairs. The absolute minima for most sequences occurs for a doublet geometry close to the B-DNA fibre models. Hence in the B-DNA region, no strong sequence dependent features are found, but the range of doublet geometries observed in the crystal structures generally lie within the low energy contours, obtained from stacking energy calculations. The doublet geometry corresponding to the A-DNA fibre model is not energetically favourable for the purine-pyrimidine sequences, which prefer small roll angle values when the slide has a large negative value as in A-DNA. However positive roll with large negative slide is allowed for GG, GA, AG and the pyrimidine-purine steps. This is consistent with the observed geometries of various steps in A-DNA crystals. Thus the general features of the basepair doublets predicted from these theoretical studies agree very well with the results from crystal structure analysis. However, since most sequences show an overall preference for B-type doublet geometry, the B----A transition for random sequence DNA cannot be explained on the basis of basepair stacking interactions.     

C-H. . .pi interactions in the metal-porphyrin complexes with chelate ring as the H acceptor

Specific C-H. . .pi interactions with the pi-system of porphyrinato chelate ring were found in crystal structures of transition metal complexes from the Cambridge Structural Database and statistical analysis of geometrical parameters for intramolecular and intermolecular interactions was done. By density functional theory calculations on a model system it was evaluated that an interaction energy is above 1.5 kcal/mol and that the strongest interaction occurs when the distance between hydrogen atom and the center of the chelate ring is 2.6 A. This prediction is in good agreement with the distances for intermolecular interactions found in the crystal structures. In many cases the intramolecular interaction distances are much shorter than 2.6 A, and these short distances are caused by geometrical constrains. The C-H. . .pi interactions with chelate ring of porphyrinato ligand can influence the structure, contribute to its stability, and play some role in the function of biomolecules with metalo porphyrins.     

Quantum-mechanical computations on the electronic structure of trans-resveratrol and trans-piceatannol: a theoretical study of the stacking interactions in trans-resveratrol dimers

Accurate quantum-chemical calculations based on the second-order M?ller-Plesset perturbation method (MP2) and density functional theory (DFT) were performed for the first time to investigate the electronic structures of trans-resveratrol and trans-piceatannol, as well as to study the stacking interaction between trans-resveratrol molecules. Ab initio MP2 calculations performed with using standard split-valence Pople basis sets led us to conclude that these compounds have structures that deviate strongly from planarity, whereas the DFT computations for the same basis sets revealed that the equilibrium geometries of these bioactive polyphenols are planar. Furthermore, the results obtained at the MP2(full)/aug-cc-pVTZ and B3LYP/aug-cc-pVTZ levels indicated that the geometries of trans-resveratrol and trans-piceatannol are practically planar at their absolute energy minima. The relative energies of the equilibrium geometries of trans-resveratrol on its potential energy surface were computed at the MP2(full)/aug-cc-pVTZ level. According to the results obtained, a T-shaped (edge-to-phase) conformer of trans-resveratrol dimer is the most stable in vacuum. This T-shaped conformer is mainly stabilized by strong hydrogen bonding and weak C-H...π interactions. Stacked structures with parallel-displaced trans-stilbene skeletons were also found to be energetically stable. The vertical separation and twist angle dependencies of the stacking energy were investigated at the MP2(full)/aug-cc-pVTZ, B3LYP/aug-cc-pVTZ, and HF/aug-cc-pVTZ levels. The standard B3LYP functional and the Hartree-Fock method neglect long-range attractive dispersion interactions. The MP2 computations revealed that the London dispersion energy cannot be neglected at long or short distances. The stacked model considered here may be useful for predicting the quantum nature of the interactions in π-stacked systems of other naturally occurring stilbenoids, and can help to enhance our understanding of the antioxidant and anticancer activities of trans-resveratrol.     

Geometries of stacking interactions between phenanthroline ligands in crystal structures of square-planar metal complexes

Stacking interactions of phenanthroline square-planar complexes in crystal structures were studied by analyzing data from the Cambridge Structural Database. In most of the crystal structures, two phenanthroline complexes were oriented “head to tail.” Phenanthroline complexes show a wide range of overlap geometries in stacking interactions, while short metal–metal distances were not observed. Stacking chains with alternating overlaps were the predominant type of packing in the crystal structures.     

Heteronuclear NMR and soft docking: an experimental approach for a structural model of the cytochrome c553-ferredoxin complex

The combination of docking algorithms with NMR data has been developed extensively for the studies of protein-ligand interactions. However, to extend this development for the studies of protein-protein interactions, the intermolecular NOE constraints, which are needed, are more difficult to access. In the present work, we describe a new approach that combines an ab initio docking calculation and the mapping of an interaction site using chemical shift variation analysis. The cytochrome c553-ferredoxin complex is used as a model of numerous electron-transfer complexes. The 15N-labeling of both molecules has been obtained, and the mapping of the interacting site on each partner, respectively, has been done using HSQC experiments. 1H and 15N chemical shift analysis defines the area of both molecules involved in the recognition interface. Models of the complex were generated by an ab initio docking software, the BiGGER program (bimolecular complex generation with global evaluation and ranking). This program generates a population of protein-protein docked geometries ranked by a scoring function, combining relevant stabilization parameters such as geometric complementarity surfaces, electrostatic interactions, desolvation energy, and pairwise affinities of amino acid side chains. We have implemented a new module that includes experimental input (here, NMR mapping of the interacting site) as a filter to select the accurate models. Final structures were energy minimized using the X-PLOR software and then analyzed. The best solution has an interface area (1037.4 A2) falling close to the range of generally observed recognition interfaces, with a distance of 10.0 A between the redox centers.     

Comparison of energy-minimized crystal structures of 2,3-dilauroyl-D-glycerol, 3-palmitoyl-DL-glycerol-1-phosphorylethanolamine and 1,2-dilauroyl-DL-phosphatidylethanolamine:acetic acid

Computational procedures have been developed by which the total energy of a lipid multibilayer can be calculated and minimized. The energy is expressed as a sum of non-bonded, electrostatic, hydrogen bonded and torsional energy terms and includes intramolecular and intermolecular components. Calculations were carried out on three lipid crystals for which structural data are available from X-ray diffraction analysis. For each crystal, the energy was minimized as a function of all bond rotations, molecular rotations and translations and the lattice constants. The minimized structures differed by only small amounts from the experimental structures, which confirms the validity of the current set of energy functions and parameters for use with lipids. The intermolecular energy of each crystal is analyzed in terms of lateral interactions, interactions between the two monolayers of the same bilayer and interactions between bilayers. The intermolecular non-bonded energy per CH2 or CH3 group in the acyl chains is also given.     

Computation of mixed phosphatidylcholine-cholesterol bilayer structures by energy minimization.

The energetically preferred structures of dimyristoylphosphatidylcholine (DMPC)-cholesterol bilayers were determined at a 1:1 mole ratio. Crystallographic symmetry operations were used to generate planar bilayers of cholesterol and DMPC. Energy minimization was carried out with respect to bond rotations, rigid body motions, and the two-dimensional lattice constants. The lowest energy structures had a hydrogen bond between the cholesterol hydroxyl and the carbonyl oxygen of the sn-2 acyl chain, but the largest contribution to the intermolecular energy was from the nonbonded interactions between the flat alpha surface of cholesterol and the acyl chains of DMPC. Two modes of packing in the bilayer were found; in structure A (the global minimum), unlike molecules are nearest neighbors, whereas in structure B (second lowest energy) like-like intermolecular interactions predominate. Crystallographic close packing of the molecules in the bilayer was achieved, as judged from the molecular areas and the bilayer thickness. These energy-minimized structures are consistent with the available experimental data on mixed bilayers of lecithin and cholesterol, and may be used as starting points for molecular dynamics or other calculations on bilayers.     

Interaction between bound water molecules and local protein structures: A statistical analysis of the hydrogen bond structures around bound water molecules

Water molecules play an important role in protein folding and protein interactions through their structural association with proteins. Examples of such structural association can be found in protein crystal structures, and can often explain protein functionality in the context of structure. We herein report the systematic analysis of the local structures of proteins interacting with water molecules, and the characterization of their geometric features. We first examined the interaction of water molecules with a large local interaction environment by comparing the preference of water molecules in three regions, namely, the protein–protein interaction (PPI) interfaces, the crystal contact (CC) interfaces, and the non‐interfacial regions. High preference of water molecules to the PPI and CC interfaces was found. In addition, the bound water on the PPI interface was more favorably associated with the complex interaction structure, implying that such water‐mediated structures may participate in the shaping of the PPI interface. The pairwise water‐mediated interaction was then investigated, and the water‐mediated residue–residue interaction potential was derived. Subsequently, the types of polar atoms surrounding the water molecules were analyzed, and the preference of the hydrogen bond acceptor was observed. Furthermore, the geometries of the structures interacting with water were analyzed, and it was found that the major structure on the protein surface exhibited planar geometry rather than tetrahedral geometry. Several previously undiscovered characteristics of water–protein interactions were unfolded in this study, and are expected to lead to a better understanding of protein structure and function. Proteins 2016; 84:43–51. © 2015 Wiley Periodicals, Inc.     

Thermochemistry of dissolving glycine, glycyl-glycine, and diglycyl-glycine in a mixed water-dimethylsulfoxide solvent at 298.15 K

The enthalpies of dissolving glycine, glycyl-glycine and diglycyl-glycine (deltaH(soln)0) in a mixed water-dimethylsulfoxide (DMSO) solvent was determined by the calorimetric method in the range of concentrations of the organic component 0 < X2 < 0.4 m.d. at 298.15 K. The enthalpies of solvation ((deltaH(solv)0) and transfer ((deltaH(tr)0) of these compounds from water to a mixed solvent were calculated. The dependencies deltaH(tr)0 =f(X2) were found to be extreme, indicating complex intermolecular interactions between the solution components. The influence of the structure and the properties of the substances dissolved and the composition of the mixture and the nature of organic solvent on their thermochemical characteristics was studied. The coefficients of enthalpy for pair interactions of glycine and its oligomers with DMSO molecules were calculated. These have positive values and increase in the order: glycyl-glycine < glycine < diglycyl-glycine. The changes in the thermochemical characteristics of dissolving, transfer, and solvation of glycine and its olygomers were shown to be determined by the energy of the mixed solvent formation, the nature of the organic solvents, and the structure of amino acids and peptides.     

Hydrogen bond connectivity patterns and hydrophobic interactions in crystal structures of small, acyclic peptides.

Crystal structures of all available unblocked linear peptides with two to five residues were retrieved from the Cambridge Structural Database and their intermolecular contacts and packing modes studied using molecular graphics. This survey reveals that interactions between hydrophobic portions of the molecules are critically important in determining the overall features of their crystal packing patterns. Distinct hydrophobic columns or layers are observed in almost all crystal structures. Analyses of the relationships between these interactions and crystal growth properties of small peptides are given. It is suggested that needle growth is promoted by hydrophobic packing, usually along a short crystallographic axis (4.6-6.0 angstroms). Also contributing to these morphologic characteristics are entropic factors associated with hydrophobic aggregation as well as tightly bound water molecules on hydrophobic faces. The paper also provides a comprehensive overview of hydrogen bond patterns in acyclic peptide crystals. It is demonstrated that one of their primary roles is to provide a scaffolding within which hydrophobic groups can aggregate. Even though there is a high density of hydrogen bonds in the crystals, often with complex patterns and networks, certain motifs are found to recur in a number of structures indicating specific hydrogen bond preferences. Water, for example, is an integral part of the hydrogen bond networks in these crystals, usually acting as the primary donor for main-chain carboxylate groups in peptide hydrates.     

Competition and interplay between the lithium bonding and hydrogen bonding: R3C···HY···LiY and R3C···LiY···HY triads as a working model (R=H, CH3; Y=CN, NC)

UMP2 calculations with aug-cc-pVDZ basis set were used to analyze intermolecular interactions in R₃C···HY···LiY and R₃C···LiY···HY triads (R=H, CH₃; Y=CN, NC), which are connected via lithium and hydrogen bonds. To better understand the properties of these systems, the corresponding dyads were also studied. Molecular geometries and binding energies of dyads, and triads were investigated at the UMP2/aug-cc-pVDZ computational level. Particular attention was paid to parameters such as cooperative energies, and many-body interaction energies. All studied complexes, with the simultaneous presence of a lithium bond and a hydrogen bond, showed cooperativity with energy values ranging between ?1.71 and ?9.03 kJ mol^?1. The electronic properties of the complexes were analyzed using parameters derived from atoms in molecules (AIM) methodology. Energy decomposition analysis revealed that the electrostatic interactions are the major source of the attraction in the title complexes.     

Insights into scFv:drug binding using the molecular dynamics simulation and free energy calculation

Molecular dynamics simulations and free energy calculation have been performed to study how the single-chain variable fragment (scFv) binds methamphetamine (METH) and amphetamine (AMP). The structures of the scFv:METH and the scFv:AMP complexes are analyzed by examining the time-dependence of their RMSDs, by analyzing the distance between some key atoms of the selected residues, and by comparing the averaged structures with their corresponding crystallographic structures. It is observed that binding an AMP to the scFv does not cause significant changes to the binding pocket of the scFv:ligand complex. The binding free energy of scFv:AMP without introducing an extra water into the binding pocket is much stronger than scFv:METH. This is against the first of the two scenarios postulated in the experimental work of Celikel et al. (Protein Science 18, 2336 (2009)). However, adding a water to the AMP (at the position of the methyl group of METH), the binding free energy of the scFv:AMP-H2O complex, is found to be significantly weaker than scFv:METH. This is consistent with the second of the two scenarios given by Celikel et al. Decomposition of the binding energy into ligand-residue pair interactions shows that two residues (Tyr175 and Tyr177) have nearly-zero interactions with AMP in the scFv:AMP-H2O complex, whereas their interactions with METH in the scFv:METH complex are as large as -0.8 and -0.74 kcal mol^-1. The insights gained from this study may be helpful in designing more potent antibodies in treating METH abuse.