A fast algorithm for macromolecular packing calculation

An algorithm to improve the computation time of packing calculations for macromolecules is presented. This is achieved by reducing the three-dimensional search to a small set of two-dimensional searches.

Thermodynamic and structural studies of cavity formation in proteins suggest that loss of packing interactions rather than the hydrophobic effect dominates the observed energetics

The hydrophobic effect is widely believed to be an important determinant of protein stability. However, it is difficult to obtain unambiguous experimental estimates of the contribution of the hydrophobic driving force to the overall free energy of folding. Thermodynamic and structural studies of large to small substitutions in proteins are the most direct method of measuring this contribution. We have substituted the buried residue Phe8 in RNase S with alanine, methionine, and norleucine, Binding thermodynamics and structures were characterized by titration calorimetry and crystallography, respectively. The crystal structures of the RNase S F8A, F8M, and F8Nle mutants indicate that the protein tolerates the changes without any main chain adjustments, The correlation of structural and thermodynamic parameters associated with large to small substitutions was analyzed for nine mutants of RNase S as well as 32 additional cavity-containing mutants of T4 lysozyme, human lysozyme, and barnase. Such substitutions were typically found to result in negligible changes in Delta C-p and positive values of both Delta Delta H degrees and aas of folding. Enthalpic effects were dominant, and the sign of Delta Delta S is the opposite of that expected from the hydrophobic effect. Values of Delta Delta G degrees and Delta Delta H degrees correlated better with changes in packing parameters such as residue depth or occluded surface than with the change in accessible surface area upon folding. These results suggest that the loss of packing interactions rather than the hydrophobic effect is a dominant contributor to the observed energetics for large to small substitutions. Hence, estimates of the magnitude of the hydrophobic driving force derived from earlier mutational studies are likely to be significantly in excess of the actual value.

Topological analysis off charge density distribution in concomitant polymorphs of 3-acetylcoumarin, a case off packing polymorphism

Detailed investigation of the charge density distribution in concomitant polymorphs of 3-acetylcoumarin in terms of experimental and theoretical densities shows significant differences in the intermolecular features when analyzed based on the topological properties via the quantum theory of atoms in molecules. The two forms, triclinic and monoclinic (Form A and Form B), pack in the crystal lattice via weak C-H---O and C-H---pi interactions. Form A results in a head-to-head molecular stack, while Form B generates a head-to-tail stack. Form A crystallizes in PI (Z' = 2) and Form B crystallizes in P2(1)/n (Z = 1). The electron density maps of the polymorphs demonstrate the differences in the nature of the charge density distribution in general. The charges derived from experimental and theoretical analysis show significant differences with respect to the polymorphic forms. The molecular dipole moments differ significantly for the two forms. The lattice energies evaluated at the HF and DFT (B3LYP) methods with 6-31G** basis set for the two forms clearly suggest that Form A is the thermodynamically stable form as compared to Form B. Mapping of electrostatic potential over the molecular surface shows dominant variations in the electronegative region, which bring out the differences between the two forms.

Structural features of B-DNA dodecamer crystal structures: Influence of crystal packing versus base sequence

We have analyzed the set of inter and intra base pair parameters for each dinucleotide step in single crystal structures of dodecamers, solved at high and medium resolution and all crystallized in P2(1)2(1)2(1) space group. The objective was to identify whether all the structures which have either the Drew-Dickerson (DD) sequence d[CGCGAATTCGCG] with some base modification or related sequence (non-DD), would display the same sequence dependent structural variability about its palindromic sequence, despite the molecule being bent at one end because of similar crystal lattice packing effect. Most of the local doublet parameters for base pairs steps G2-C3 and G10-C11 positions, symmetrically situated about the lateral twofold, were significantly correlated between themselves. In non-DD sequences, significant correlations between these positional parameters were absent. The different range of local step parameter values at each sequence position contributed to the gross feature of smooth helix axis bending in all structures. The base pair parameters in some of the positions, for medium resolution DD sequence, were quite unlike the high-resolution set and encompassed a higher range of values. Twist and slide are the two main parameters that show wider conformational range for the middle region of non-DD sequence structures in comparison to DD sequence structures. On the contrary, the minor and major groove features bear good resemblance between DD and non-DD sequence crystal structure datasets. The sugar-phosphate backbone torsion angles are similar in all structures, in sharp contrast to base pair parameter variation for high and low resolution DD and non-DD sequence structures, consisting of unusual (epsilon =g(-), xi =t) B-II conformation at the 10(th) position of the dodecamer sequence. Thus examining DD and non-DD sequence structures packed in the same crystal lattice arrangement, we infer that inter and intra base pair parameters are as symmetrically equivalent in its value as the symmetry related step for the palindromic DD sequence about lateral two-fold axis. This feature would lead us to agree with the conclusion that DNA conformation is not substantially affected by end-to-end or lateral inter-molecular interaction due to crystal lattice packing effect. Non-DD sequence structures acquire step parameter values which reflect the altered sequence at each of the dodecamer sequence position in the orthorhombic lattice while showing similar gross features of DD sequence structures

Insights into conformational and packing features in a series of aryl substituted ethyl-6-methyl-4-phenyl-2-oxo-1,2,3,4-tetrahydropyrimidine-5-carboxylate

The supramolecular structures of eight aryl protected ethyl-6-methyl-4-phenyl-2-oxo-1,2,3,4-tetrahydropyrimidine- 5-carboxylates have been analyzed to determine the role of different functional groups on the molecular geometry, conformational characteristics and the packing of these molecules in the crystal lattice. Out of these the para fluoro substituted compound on the aryl ring exhibits conformational polymorphism, due to the different conformation of the ester moiety. This behaviour has been characterized using both powder and single-crystal X-ray diffraction, optical microscopy and differential scanning calorimetry performed on both these polymorphs. The compounds pack via the cooperative interplay of strong N-H center dot center dot center dot O=C intermolecular dimers and chains forming a sheet like structure. In addition, weak C-H center dot center dot center dot O=C and C-H center dot center dot center dot pi interactions impart additional stability to the crystal packing.

Studies in crystal engineering: Effect of fluorine substitution in crystal packing and topological photodimerization of styryl coumarins in the solid state

Styryl coumarins generally yield centrosymmetric (alpha-mode, anti-HT) photodimers when subjected to irradiation in the solid state, However, the substitution of fluorine dramatically alters the packing mode and steers the molecules 4-(4-fluorostyryl)coumarin 1 and 4-(2-fluorostyryl)coumarin 2 to form a stereospecific photodimer, beta-mode, syn-HH across the styrenic double bond (yield 78-85%). The stereochemistry of the photodimer 2a has been established by X-ray crystallography. There is no evidence for the presence of C-H ... F interactions. The true nature of the weak atom-atom interactions called into play when fluorine is substituted is not clear, It is observed that the fluoro substituted compounds have greater crystal density than the corresponding unsubstituted ones.

Solid-state photobehaviour and crystal packing of o-chlorobenzylidene-DL-piperitone: influence of molecular topology on photobehaviour

C17H19ClO, M(r) = 274.7, triclinic, P1BAR, a = 11.154 (3), b = 12.685 (2), c = 12.713 (2) angstrom, alpha = 100.68 (1), beta = 113.58 (1), gamma = 104.50 (2)-degrees, V = 1511.1 (6) angstrom3, Z = 4, D(m) = 1.22, D(x) = 1.215 Mg m-3, Cu K-alpha, lambda = 1.5418 angstrom, mu = 2.16 mm-1, F(000) = 584, T = 293 K, R = 0.057 for 3481 observed reflections. The title compound is photostable in the crystalline state and lattice-energy calculations have been employed to rationalize the photobehaviour. The well-known beta-steering ability of the chloro group is not operative in this system as there are no Cl...Cl interactions in the crystal lattice. All five benzylidene-DL-piperitone structures so far studied are alpha-packed and the molecular topology appears to be a deciding factor even in the presence of steering groups.

Studies in crystal engineering: Crystal packing, topological photodimerization and structure-reactivity correlations in fluoro-substituted styrylcoumarins

In continuation of our studies on the influence of fluoro substitution on the solid state photobehaviour and packing pattern of styrylcoumarins, the results obtained for 4-(3-fluorostyryl)coumarin 1, 4-styryl-6-fluorocoumarin 2 and 4-styryl-7-fluorocoumarin 3 are presented. The configuration of the dimers was established on the basis of crystal packing of 1 and 2 (alpha-packed). A rationale for the significantly lower dimer yield in the crystal for 2 is proposed. In the observed centrosymmetric arrangement of the reactants the C=O ...pi (phenyl) contacts seem to provide additional attractive interactions. C-H ... O and C-H ... F hydrogen bonding seems to provide stability in these structures.

Effect of packing density and adsorption parameters on the throughput of a thermal compressor

Vapour adsorption refrigeration systems (VAdS) have the advantage of scalability over a wide range of capacities ranging from a few watts to several kilowatts. In the first instance, the design of a system requires the characteristics of the adsorbate-adsorbent pair. Invariably, the void volume in the adsorbent reduces the throughput of the thermal compressor in a manner similar to the clearance volume in a reciprocating compressor. This paper presents a study of the activated carbon +HFC-134a (1,1,1,2-tetrafluoroethane) system as a possible pair for a typical refrigeration application. The aim of this study is to unfold the nexus between the adsorption parameters, achievable packing densities of charcoal and throughput of a thermal compressor. It is shown that for a thermal compressor, the adsorbent should not only have a high surface area, but should also be able to provide a high packing density. Given the adsorption characteristics of an adsorbent-adsorbate pair and the operating conditions, this paper discloses a method for the calculation of the minimum packing density necessary for an effective throughput of a thermal compressor. (C) 2002 Elsevier Science Ltd. All rights reserved.

Crystal packing and melting temperatures of small oxalate esters: the role of C-H center dot center dot center dot O hydrogen bonding

The simple dialkyl oxalates are generally liquids at room temperature except for dimethyl and di-tert-butyl oxalate which melt at 327 and 343 K. The crystal structures of diethyl, di-iso-propyl, di-n-butyl, di-tert-butyl and methyl ethyl oxalates were determined. The liquid esters were crystallized using the cryocrystallization technique. A comparison of the intermolecular interactions and packing features in these crystal structures was carried out. The crystal structure of dimethyl oxalate was redetermined at various temperatures. The other compounds were also studied at several temperatures in order to assess the attractive nature of the hydrogen bonds therein. A number of moderate to well defined C-H center dot center dot center dot O interactions account for the higher melting points of the two solid esters. Additionally, a diminished entropic contribution Delta S(m) in di-tert-butyl oxalate possibly increases the melting point of this compound further.

Evaluation of the role of disordered organic fluorine in crystal packing: insights from halogen substituted benzanilides

The determination of the crystal and molecular structures of a large number of compounds containing the C(sp(2))-F bond has been investigated in detail in halogenated benzanilides and also in liquids, namely the fluorinated amines. It has been observed that when the fluorine atom is present in the ortho or meta position with respect to the amide functionality in benzanilides or the amino group in fluorinated amines which are liquids at room temperature, the fluorine atom exhibits positional disorder. This is associated with changes in patterns of intermolecular interactions which affect crystal packing. Furthermore, the presence of a fluorine atom on the benzanilide framework, in the presence of a heavier halogen (chloro, bromo and iodo), meta or ortho to the amide group does not eliminate the disorder associated with these molecules. In this article, we highlight the salient features present in halogenated compounds exhibiting disorder in the position of organic fluorine with concomitant changes in crystal packing. This feature is also compared with related compounds exhibiting similarity in electronic features, namely positional disorder.