Structural studies of analgesics and their interactions. Part 4. Crystal structures of phenylbutazone and a 2 : 1 complex between phenylbutazone and piperazine

The X-ray crystal structures of 4-butyl-1,2-diphenylpyrazolidine-3,5-dione (phenylbutazone)(I). and its 2 : 1 complex (II) with piperazine have been determined by direct methods and the structures refined to R 0.096 (2 300 observed reflections measured by diffractometer) and 0.074 (2 494 observed reflections visuallyestimated). Crystals are monoclinic, space group P21/c; for (I)a= 21.695(4), b= 5.823(2), c= 27.881(4)Å, = 108.06 (10)°, Z= 8, and for (II)a= 8.048(4), b= 15.081(4), c= 15.583(7)Å, = 95.9(3)°, Z= 2. The two crystallographically independant molecules in the structure of (I) are similar except for the conformation of the butyl group, which is disordered in one of the molecules. In the pyrazolidinedione group, the two C–C bonds are single and the two C–O bonds double. The two nitrogen atoms in the five-membered ring are pyramidal with the attached phenyl groups lying on the opposite sides of the mean plane of the ring. The phenylbutazone molecule in (II) exists as a negative ion owing to deprotonation of C-4. C-4 is therefore trigonal and the orientation of the Bu group with respect to the pyrazolidinedione group is considerably different from that in (I); there is also considerable electron delocalization along the C–O and C–C bonds. These changes in geometry and electronic structure may relate to biological activity. The doubly charged cationic piperazine molecule exists in the chair form with the nitrogen atoms at the apices. The crystal structure of (II) is stabilized by ionic interactions and N–H O hydrogen bonds.

Diffusion of localized optical excitations and trap emission in ruby

Using the critical percolation conductance method the energy-dependent diffusion coefficient associated with thermally assisted transfer of the R1 line excitation between single Cr3+ ions with strain-induced randomness has been calculated in the 4A2 to E(2E) transition energies. For localized states sufficiently far away from the mobility edge the energy transfer is dominated by dipolar interactions, while very close to the mobility edge it is determined by short-range exchange interactions. Using the above energy-dependent diffusion coefficient a macroscopic diffusion equation is solved for the rate of light emission by Cr3+ ion-pair traps to which single-ion excitations are transferred. The dipolar mechanism leads to good agreement with recent measurements of the pair emission rate by Koo et al. (Phys. Rev. Lett., vol.35, p.1669 (1975)) right up to the mobility edge.

Platelet endothelial cell adhesion molecule 1 (PECAM-1) and its interactions with glycosaminoglycans: 2. Biochemical analyses

Platelet endothelial cell adhesion molecule 1 (PECAM-1) (CD31), a member of the immunoglobulin (Ig) superfamily of cell adhesion molecules with six Ig-like domains, has a range of functions, notably its contributions to leukocyte extravasation during inflammation and in maintaining vascular endothelial integrity. Although PECAM-1 is known to mediate cell adhesion by homophilic binding via domain 1, a number of PECAM-1 heterophilic ligands have been proposed. Here, the possibility that heparin and heparan sulfate (HS) are ligands for PECAM-1 was reinvestigated. The extracellular domain of PECAM-1 was expressed first as a fusion protein with the Fc region of human IgG1 fused to domain 6 and second with an N-terminal Flag tag on domain 1 (Flag-PECAM-1). Both proteins bound heparin immobilized on a biosensor chip in surface plasmon resonance (SPR) binding experiments. Binding was pH-sensitive but is easily measured at slightly acidic pH. A series of PECAM-1 domain deletions, prepared in both expression systems, were tested for heparin binding. This revealed that the main heparin-binding site required both domains 2 and 3. Flag-PECAM-1 and a Flag protein containing domains 1-3 bound HS on melanoma cell surfaces, but a Flag protein containing domains 1-2 did not. Heparin oligosaccharides inhibited Flag-PECAM-1 from binding immobilized heparin, with certain structures having greater inhibitory activity than others. Molecular modeling similarly identified the junction of domains 2 and 3 as the heparin-binding site and further revealed the importance of the iduronic acid conformation for binding. PECAM-1 does bind heparin/HS but by a site that is distinct from that required for homophilic binding.

Platelet endothelial cell adhesion molecule 1 (PECAM-1) and its interactions with glycosaminoglycans: 1. Molecular modeling studies

Platelet endothelial cell adhesion molecule 1 (PECAM-1) has many functions, including its roles in leukocyte extravasation as part of the inflammatory response and in the maintenance of vascular integrity through its contribution to endothelial cell−cell adhesion. PECAM-1 has been shown to mediate cell−cell adhesion through homophilic binding events that involve interactions between domain 1 of PECAM-1 molecules on adjacent cells. However, various heterophilic ligands of PECAM-1 have also been proposed. The possible interaction of PECAM-1 with glycosaminoglycans (GAGs) is the focus of this study. The three-dimensional structure of the extracellular immunoglobulin (Ig) domains of PECAM-1 were constructed using homology modeling and threading methods. Potential heparin/heparan sulfate-binding sites were predicted on the basis of their amino acid consensus sequences and a comparison with known structures of sulfate-binding proteins. Heparin and other GAG fragments have been docked to investigate the structural determinants of their protein-binding specificity and selectivity. The modeling has predicted two regions in PECAM-1 that appear to bind heparin oligosaccharides. A high-affinity binding site was located in Ig domains 2 and 3, and evidence for a low-affinity site in Ig domains 5 and 6 was obtained. These GAG-binding regions were distinct from regions involved in PECAM-1 homophilic interactions.

Computational analyses of the catalytic and heparin-binding sites and their interactions with glycosaminoglycans in glycoside hydrolase family 79 endo-d-glucuronidase (heparanase)

Mammalian heparanase is an endo-β-glucuronidase associated with cell invasion in cancer metastasis, angiogenesis and inflammation. Heparanase cleaves heparan sulfate proteoglycans in the extracellular matrix and basement membrane, releasing heparin/heparan sulfate oligosaccharides of appreciable size. This in turn causes the release of growth factors, which accelerate tumor growth and metastasis. Heparanase has two glycosaminoglycan-binding domains; however, no three-dimensional structure information is available for human heparanase that can provide insights into how the two domains interact to degrade heparin fragments. We have constructed a new homology model of heparanase that takes into account the most recent structural and bioinformatics data available. Heparin analogs and glycosaminoglycan mimetics were computationally docked into the active site with energetically stable ring conformations and their interaction energies were compared. The resulting docked structures were used to propose a model for substrates and conformer selectivity based on the dimensions of the active site. The docking of substrates and inhibitors indicates the existence of a large binding site extending at least two saccharide units beyond the cleavage site (toward the nonreducing end) and at least three saccharides toward the reducing end (toward heparin-binding site 2). The docking of substrates suggests that heparanase recognizes the N-sulfated and O-sulfated glucosamines at subsite +1 and glucuronic acid at the cleavage site, whereas in the absence of 6-O-sulfation in glucosamine, glucuronic acid is docked at subsite +2. These findings will help us to focus on the rational design of heparanase-inhibiting molecules for anticancer drug development by targeting the two heparin/heparan sulfate recognition domains.

The structure of glycosaminoglycans and their interactions with proteins

Glycosaminoglycans (GAGs) are important complex carbohydrates that participate in many biological processes through the regulation of their various protein partners. Biochemical, structural biology and molecular modelling approaches have assisted in understanding the molecular basis of such interactions, creating an opportunity to capitalize on the large structural diversity of GAGs in the discovery of new drugs. The complexity of GAG–protein interactions is in part due to the conformational flexibility and underlying sulphation patterns of GAGs, the role of metal ions and the effect of pH on the affinity of binding. Current understanding of the structure of GAGs and their interactions with proteins is here reviewed: the basic structures and functions of GAGs and their proteoglycans, their clinical significance, the three-dimensional features of GAGs, their interactions with proteins and the molecular modelling of heparin binding sites and GAG–protein interactions. This review focuses on some key aspects of GAG structure–function relationships using classical examples that illustrate the specificity of GAG–protein interactions, such as growth factors, anti-thrombin, cytokines and cell adhesion molecules. New approaches to the development of GAG mimetics as possible new glycotherapeutics are also briefly covered.

Computational methods for the prediction of the structure and interactions of coiled-coil peptides

The past several years have seen significant advances in the development of computational methods for the prediction of the structure and interactions of coiled-coil peptides. These methods are generally based on pairwise correlations of amino acids, helical propensity, thermal melts and the energetics of sidechain interactions, as well as statistical patterns based on Hidden Markov Model (HMM) and Support Vector Machine (SVM) techniques. These methods are complemented by a number of public databases that contain sequences, motifs, domains and other details of coiled-coil structures identified by various algorithms. Some of these computational methods have been developed to make predictions of coiled-coil structure on the basis of sequence information; however, structural predictions of the oligomerisation state of these peptides still remains largely an open question due to the dynamic behaviour of these molecules. This review focuses on existing in silico methods for the prediction of coiled-coil peptides of functional importance using sequence and/or three-dimensional structural data.

Free energy calculations of glycosaminoglycan - Protein interactions

Glycosaminoglycans (GAGs) are complex highly charged linear polysaccharides that have a variety of roles in biological processes. We report the first use of molecular dynamics (MD) free energy calculations using the MM/PBSA method to investigate the binding of GAGs to protein molecules, namely the platelet endothelial cell adhesion molecule 1 (PECAM-1) and annexin A2. Calculations of the free energy of the binding of heparin fragments of different sizes reveal the existence of a region of low GAG-binding affinity in domains 5-6 of PECAM-1 and a region of high affinity in domains 2-3, consistent with experimental data and ligand-protein docking studies. A conformational hinge movement between domains 2 and 3 was observed, which allows the binding of heparin fragments of increasing size (pentasaccharides to octasaccharides) with an increasingly higher binding affinity. Similar simulations of the binding of a heparin fragment to annexin A2 reveal the optimization of electrostatic and hydrogen bonding interactions with the protein and protein-bound calcium ions. In general, these free energy calculations reveal that the binding of heparin to protein surfaces is dominated by strong electrostatic interactions for longer fragments, with equally important contributions from van der Waals interactions and vibrational entropy changes, against a large unfavorable desolvation penalty due to the high charge density of these molecules.

Molecular dynamics simulations of CXCL-8 and its interactions with a receptor peptide, heparin fragments, and sulfated linked cyclitols

CXCL-8 (Interleukin 8) is a CXC chemokine with a central role in the human immune response. We have undertaken extensive in silico analyses to elucidate the interactions of CXCL-8 with its various binding partners, which are crucial for its biological function. Sequence and structure analyses showed that residues in the thirdq β-sheet and basic residues in the heparin binding site are highly variable, while residues in the second β-sheet are highly conserved. Molecular dynamics simulations in aqueous solution of dimeric CXCL-8 have been performed with starting geometries from both X-ray and NMR structures showed shearing movements between the two antiparallel C-terminal helices. Dynamic conservation analyses of these simulations agreed with experimental data indicating that structural differences between the two structures at quaternary level arise from changes in the secondary structure of the N-terminal loop, the 310-helix, the 30s, 40s, and 50s loops and the third β-sheet, resulting in a different interhelical separation. Nevertheless, the observation of these different states indicates that CXCL-8 has the potential to undergo conformational changes, and it seems likely that this feature is relevant to the mode of binding of glycosaminoglycan (GAG) mimetics such as cyclitols. Simulations of the receptor peptide fragment−CXCL-8 complex identified several specific interactions of the receptor peptide with CXCL-8 that could be exploited in the structure-based design of competitive peptides and nonpeptidic molecules targeting CXCL-8 for combating inflammatory diseases. Simulations of the CXCL-8 dimer complexed with a 24-mer heparin fragment and of the CXCL-8−receptor peptide complex revealed that Arg60, Lys64, and Arg68 in the dimer bind to cyclitols in a horseshoe pattern, defining a region which is spatially distinct from the receptor binding site. There appears to be an optimum number of sulfates and an optimum length of alkyl spacers required for the interaction of cyclitol inhibitors with the dimeric form of CXCL-8. Calculation of the binding affinities of cyclitol inhibitors reflected satisfactorily the ranking of experimentally determined inhibitory potencies. The findings of these molecular modeling studies will help in the search for inhibitors which can modulate various CXCL-8 biological activities and serve as an excellent model system to study CXC-inhibitor interactions.

Free energy calculations of the interactions of c-Jun-based synthetic peptides with the c-Fos protein

The c-Fos–c-Jun complex forms the activator protein 1 transcription factor, a therapeutic target in the treatment of cancer. Various synthetic peptides have been designed to try to selectively disrupt the interaction between c-Fos and c-Jun at its leucine zipper domain. To evaluate the binding affinity between these synthetic peptides and c-Fos, polarizable and nonpolarizable molecular dynamics (MD) simulations were conducted, and the resulting conformations were analyzed using the molecular mechanics generalized Born surface area (MM/GBSA) method to compute free energies of binding. In contrast to empirical and semiempirical approaches, the estimation of free energies of binding using a combination of MD simulations and the MM/GBSA approach takes into account dynamical properties such as conformational changes, as well as solvation effects and hydrophobic and hydrophilic interactions. The predicted binding affinities of the series of c-Jun-based peptides targeting the c-Fos peptide show good correlation with experimental melting temperatures. This provides the basis for the rational design of peptides based on internal, van der Waals, and electrostatic interactions.

Solid phase fermentation of leaf biomass to biogas

This paper describes a simple technique for the fermentation of untreated or partly-treated leafy biomass in a digester of novel design without incurring the normal problems of feeding, floating and scum formation of feed, etc. The solid phase fermentation studied consists of a bed of biomass frequently sprinkled with an aqueous bacterial inoculum and recycling the leachate to conserve moisture and improve the bacterial dispersion in the bed. The decomposition of the leaf biomass and water hyacinth substrates used in this study was rapid, taking 45 and 30 days for the production of 250 and 235 l biogas per kg total solids (TS) respectively, for the above mentioned substrates at a daily sprinkled volume of 26 ml cm−2 of bed per day sprinkled at 12 h intervals. Very little volatile fatty acid (VFA) intermediates accumulated in the liquid sprinkled, suggesting acidogenesis to be rate-limiting in this process. From the pattern of VFA and gas produced it is concluded that most of the biogas produced is from the biomass bed, thus making the operation of a separate methanogenic reactor unnecessary.

Stereospecific photodimerization of coumarins in crystalline inclusion complexes. Molecular and crystal structure of 1:2 complex of (S,S)-(-)-1,6-bis(o-chlorophenyl)-1,6-diphenyl-hexa-2,4-diyne-1,6-diol and coumarin

Proximity of molecules is a crucial factor in many solid- state photochemical processes.'S2 The biomolecular photodimerization reactions in the solid state depend on the relative geometry of reactant molecules in the crystal lattice with center-to-center distance of nearest neighbor double bonds of the order of ca. 4 A. This fact emanates from the incisive studies of Schmidt and Cohen.2 One of the two approaches to achieve this distance requirement is the so-called "Crystal-Engineering" of structures, which essentially involves the introduction of certain functional groups that display in-plane interstacking interactions (Cl...Cl, C-He-0, etc.) in the crystal The chloro group is by far the most successful in promoting the /3- packing m ~ d e ,th~o,u~gh recent studies have shown its limitations? Another approach involves the use of constrained media in which the reactants could hopefully be aligned.