Structure of a bis-amidinium derivative of hoechst 33258 complexed to dodecanucleotide d(CGCGAATTCGCG)2: the role of hydrogen bonding in minor groove drug-DNA recognition.

The crystal structure is reported of a complex between the dodecanucleotide sequence d(CGCGAATTCGCG)2and an analogue of the DNA binding drug Hoechst 33258, in which the piperazine ring has been replaced by an amidinium group and the phenol ring by a phenylamidinium group. The structure has been refined to an R factor of 19.5% at 2.2 A resolution. The drug is held in the minor groove by five strong hydrogen bonds, together with bridging water molecules at both ends. There are few other contacts with the floor of the groove, indicating a lack of isohelicity with the groove and suggesting (i) that the observed high DNA affinity of this drug is primarily due to the array of hydrogen bonds and (ii) that these more than compensate for its poor isohelicity.     

Mapping the lifetimes of local opening events in a native state protein.

The rate constants for the processes that lead to local opening and closing of the structures around hydrogen bonds in native proteins have been determined for most of the secondary structure hydrogen bonds in the four-helix protein acyl coenzyme A binding protein. In an analysis that combines these results with the energies of activation of the opening processes and the stability of the local structures, three groups of residues in the protein structure have been identified. In one group, the structures around the hydrogen bonds have frequent openings, every 600 to 1,500 s, and long lifetimes in the open state, around 1 s. In another group of local structures, the local opening is a very rare event that takes place only every 15 to 60 h. For these the lifetime in the open state is also around 1 s. The majority of local structures have lifetimes between 2,000 and 20,000 s and relatively short lifetimes of the open state in the range between 30 and 400 ms. Mapping of these groups of amides to the tertiary structure shows that the openings of the local structures are not cooperative at native conditions, and they rarely if ever lead to global unfolding. The results suggest a mechanism of hydrogen exchange by progressive local openings.     

The role of specific 2'-hydroxyl groups in the stabilization of the folded conformation of kink-turn RNA

The role of 2'-hydroxyl groups in stabilizing the tightly kinked geometry of the kink-turn (K-turn) has been investigated. Individual 2'-OH groups have been removed by chemical synthesis, and the kinking of the RNA has been studied by gel electrophoresis and fluorescence resonance energy transfer. The results have been analyzed by reference to a database of 11 different crystallographic structures of K-turns. The potential hydrogen bonds fall into several classes. The most important are those in the core of the turn and ribose-phosphate interactions around the bulge. Of these the single most important hydrogen bond is one donated from the 2'-OH of the 5' nucleotide of the bulge to the N1 of the adenine of the kink-proximal A*G pair. This is present in all known K-turn structures, and removal of the 2'-OH completely prevents metal ion-induced folding. Hydrogen bonds formed in the minor grooves of the helical stems are less important, and removal of the participating 2'-OH groups leads to reduced impairment of folding. These interactions are generally more polymorphic, and hydrogen bonds probably form where possible, as permitted by the global structure.     

Structural basis of DNA flexibility

It has been recently shown by us, on the basis of crystal structure database that the flexibility of B-DNA double helices depends significantly on their base sequence. Our model building studies further indicated that the existence of bifurcated cross-strand hydrogen bonds between successive base pairs is possibly the main factor behind the sequence directed DNA flexibility. These cross-strand hydrogen bonds are, of course, weaker than the usual Watson-Crick hydrogen bonds and their bond geometry is characterized by relatively larger bond lengths and smaller bond angles. We have tried to improve our model structures by incorporating non-planarity of the amino groups in DNA bases due to the presence of lone pair electrons at the nitrogen atoms. Energy minimization studies have been carried out by using different quantum chemical methods, whereby it is found that in all cases of N-H....O type cross-strand hydrogen bonds, the bond geometry improves significantly. In the cases of N-H....N type hydrogen bonds, however, no such consistent improvements can be noticed. Perhaps the true picture would emerge only if all the other interactions present in the DNA macromolecule could be appropriately taken into account.     

Monensin A acid complexes as a model of electrogenic transport of sodium cation

New Monensin A acid complexes with water molecule, sodium chloride and sodium perchlorate were obtained and studied by X-ray and (1)H, (13)C NMR and FT-IR methods as well as ab initio calculations. The crystal structure of the complexes indicates the complexation of the water molecule and Na(+) cation in the pseudo-cycle conformation of the Monensin acid molecule stabilised by intramolecular hydrogen bonds. Important for stabilisation of this structure is also the intermolecular hydrogen bonds with water molecule or the coordination bonds with Na(+) cation. It is demonstrated that the counterions forming intermolecular hydrogen bonds with OH groups influence the strength of the intramolecular hydrogen bonds, but they have no influence on the formation of pseudo-cyclic structure. Spectroscopic studies of the complexes in dichloromethane solution have shown that the pseudo-cyclic structure of the compounds is conserved. As follows from the ab initio calculations, the interactions between the Na(+) cation and the electronegative oxygen atoms of Monensin acid totally change the molecular electrostatic potential around the supramolecular Monensin acid-Na(+) cationic complex relative to that of the neutral Monensin acid molecule.     

A method for determining the positions of polar hydrogens added to a protein structure that maximizes protein hydrogen bonding.

An automated method for the optimal placement of polar hydrogens in a protein structure is described. This method treats the polar, side chain hydrogens of lysine, serine, threonine, and tyrosine and the amino terminus of a protein. The program, called NETWORK, divides the potential hydrogen-bonding pairs of a protein into groups of interacting donors and acceptors. A search is conducted on each of the local groups to find an arrangement which forms the most hydrogen bonds. If two or more arrangements have the same number of hydrogen bonds, the arrangement with the shortest set of hydrogen bonds is selected. The polar hydrogens of the histidyl side chain are specifically treated, and the ionization state of this residue is allowed to change, if this change results in additional hydrogen bonds for the local group. The program will accept Protein Data Bank as well as Biosym-format coordinate files. Input and output routines can be easily modified to accept other coordinate file formats. The predictions from this method are compared to known hydrogen positions for bovine pancreatic trypsin inhibitor, insulin, RNase-A, and trypsin for which the neutron diffraction structures have been determined. The usefulness of this program is further demonstrated by a comparison of molecular dynamics simulations for the enzyme cytochrome P-450cam with and without using NETWORK.     

Active-site mobility in human immunodeficiency virus, type 1, protease as demonstrated by crystal structure of A28S mutant.

The mutation Ala28 to serine in human immunodeficiency virus, type 1, (HIV-1) protease introduces putative hydrogen bonds to each active-site carboxyl group. These hydrogen bonds are ubiquitous in pepsin-like eukaryotic aspartic proteases. In order to understand the significance of this difference between HIV-1 protease and homologous, eukaryotic aspartic proteases, we solved the three-dimensional structure of A28S mutant HIV-1 protease in complex with a peptidic inhibitor U-89360E. The structure has been determined to 2.0 A resolution with an R factor of 0.194. Comparison of the mutant enzyme structure with that of the wild-type HIV-1 protease bound to the same inhibitor (Hong L, Treharne A, Hartsuck JA, Foundling S, Tang J, 1996, Biochemistry 35:10627-10633) revealed double occupancy for the Ser28 hydroxyl group, which forms a hydrogen bond either to one of the oxygen atoms of the active-site carboxyl or to the carbonyl oxygen of Asp30. We also observed marked changes in orientation of the Asp25 catalytic carboxyl groups, presumably caused by the new hydrogen bonds. These observations suggest that catalytic aspartyl groups of HIV-1 protease have significant conformational flexibility unseen in eukaryotic aspartic proteases. This difference may provide an explanation for some unique catalytic properties of HIV-1 protease.     

Interesting coordination abilities of antiulcer drug famotidine and antimicrobial activity of drug and its cobalt(III) complex

Crystal structure of a novel cobalt(III) complex with antiulcer drug famotidine and ethylenediamine was determined. This is the second structure of a transition metal complex with famotidine resolved by a single crystal X-ray analysis, in which famotidine shows different mode of coordination than that observed in the other cases. Drug molecule is coordinated to metal ion as a tetradentate ligand through guanidine N6, thiazole N4, thioether S2 and terminal N3 atom. Two NH(2) groups (N3H(2) and N6H(2)) are deprotonated and drug coordinates as dianion. In the asymmetric unit, one chloride anion and one water molecule were found to complete the complex stoichiometry. The structure of the complex is abundant in atoms, which can be involved in hydrogen bond formation either as hydrogen acceptors or hydrogen donors. Because of that, a great number of hydrogen bonds dominates the crystal packing. Beside the hydrogen bonds, there are two interesting noncovalent interactions: CH(...)pi and NH(...)pi within the famotidine anion, which stabilize the complex structure. The pi(...)pi stacking interactions between neighboring complex cations are also observed. Antibacterial and antifungal activity of famotidine and its newly synthesized complex against representative bacteria: Escherichia coli, Staphilococcus aureus and Micrococcus lysodeikticus and fungi: Aspergillus niger and Candida albicans were examined. The results indicate a higher selectivity of the famotidine-Co(III) complex, as well as better growth inhibitory activity (lower MIC values (MIC, minimal inhibitory concentration)) in comparison with the drug alone.     

NMR evaluation of adipocyte fatty acid binding protein (aP2) with R- and S-ibuprofen

We have examined global chemical shift perturbations for aP2 ligand complexes and compared these with amide temperature coefficients. Hydrogen bond potential was monitored by amide chemical shift's temperature coefficient. Based on this information, we propose that the binding energy contribution can be spread out to multiple distant residues. For aP2, the ability of the receptor protein to change its hydrogen bond interactions in the beta-strands to accommodate different ligand scaffolds seems to make this receptor difficult for structure based drug design. While stabilization energy differential on hydrogen bonds is likely to be small for individual residues, the accumulative effect on multiple hydrogen bonds may have a dramatic impact on ligand affinity.     

C-H...O hydrogen bonds in the nuclear receptor RARgamma--a potential tool for drug selectivity

Hydrogen bonds between polarized atoms play a crucial role in protein interactions and are often used in drug design, which usually neglects the potential of C-H...O hydrogen bonds. The 1.4 A resolution crystal structure of the ligand binding domain of the retinoic acid receptor RARgamma complexed with the retinoid SR11254 reveals several types of C-H...O hydrogen bonds. A striking example is the hydroxyl group of the ligand that acts as an H bond donor and acceptor, leading to a synergy between classical and C-H...O hydrogen bonds. This interaction introduces both specificity and affinity within the hydrophobic ligand pocket. The similarity of intraprotein and protein-ligand C-H...O interactions suggests that such bonds should be considered in rational drug design approaches.     

Hydrophobic tendencies of polar groups as a major force in molecular recognition

Proteins and nucleic acids are able to adopt their native conformation and perform their biological role only in the presence of water with which they actively interact in a mutually modifying way. Traditionally, hydrophobic effect has been considered to be the major factor stabilizing biopolymeric structures. However, solvent reorganization around polar groups is an event thermodynamically more unfavorable than solvent reorganization around nonpolar groups. Consequently, burial of polar groups with formation of complementary solute-solute hydrogen bonds out of contact with water is an energetically favorable process that also provides a major force driving macromolecular association and folding. In contrast to nonpolar groups, polar groups may form their complementary intra- or intersolute hydrogen bonds out of contact with water only provided that an appropriate solute structure has been formed with properly positioned hydrogen bond donors and acceptors. Formation of such structures is disfavored entropically and may not be possible due to steric reasons. However, the interior of a folded protein, alpha-helices and beta-sheets, double helical nucleic acid structures, and protein-ligand interfaces all provide rigid matrices where polar groups may form their complementary hydrogen bonds. For these structures, the inward drive of polar groups represents a considerable stabilizing factor.     

Molecular and crystal structures of N-arylglycopyranosylamines formed by reaction between sulfanilamide and D-ribose, D-arabinose and D-mannose

The X-ray crystal structures of three monosaccharide derivatives prepared by the reaction of sulfanilamide with D-ribose, D-arabinose, and D-mannose have been determined. The derivatives are N-(p-sulfamoylphenyl)-alpha-D-ribopyranosylamine (1), N-(p-sulfamoylphenyl)-alpha-D-arabinopyranosylamine (2), and N-(p-sulfamoylphenyl)-beta-D-mannopyranosylamine monohydrate (3). The monosaccharide ring of 1 and 2 has the 1C4 conformation, stabilized in 1 by an intramolecular hydrogen bond from 0-2 to 0-4. Compound 3 has the 4C1 conformation at the monosaccharide ring and the gt conformation at the C-6-O-6 side chain. Occupancy of the water molecule in the crystal of 3 actually examined was 22%. The degree of interaction between sulfamoyl groups and monosaccharide moieties varies from structure to structure. The packing arrangement of 2 involves hydrogen bonding between sulfamoyl groups and monosaccharide hydroxyl groups, but interactions of this type are fewer in 1, and in 3 the hydrogen bonds are either strictly between monosaccharide hydroxyl groups or strictly between sulfamoyl groups. Pairs of hydrogen bonds (two-point contacts) link neighboring molecules in all three structures, between screw-axially related molecules in 1 and 2 and between translationally related molecules in 3. The contact in 3 defined by the O-3-H...O-5 and O-6-H...O-4 hydrogen bonds is found in several other N-aryl-beta-D-mannopyranosylamine crystal structures and is apparently an especially favorable mode of intermolecular interaction in these compounds.     

Role of Hydrogen Bonds in Protein-DNA Recognition: A Comparison of Generalized Born and Finite Difference Poisson-Boltzmann Solvation Treatments

Abstract

Hydrogen bonds have been accredited with a major role historically, in the formation and stabilization of biomolecular structures. The formation of hydrogen bonds at protein-DNA interfaces in aqueous medium involves not only favorable interactions of the donor and acceptor functional groups but also a loss of interactions between these groups with the solvent water. We have investigated the energetics of about 500 potential hydrogen bonds occuring at protein-DNA interfaces incorporating some recent improvements in biomolecular force fields and solvation treatments. We present here results of our assessment of hydrogen bond contributions to the overall standard free energy of formation of protein-DNA complexes obtained with the generalized Born model and finite difference Poisson- Boltzmann methodology for solvation in conjunction with AMBER force field. Our results support the emerging view on the role of electrostatics in general and that of hydrogen bonds in particular which is that hydrogen bonds do not drive protein-DNA complex formation by virtue of the unfavourable cost of the electrostatics of desolvation. They however, act to stabilize the complex once it is formed.     

Energetics of hydrogen bonding in proteins: a model compound study.

Differences in the energetics of amide-amide and amide-hydroxyl hydrogen bonds in proteins have been explored from the effect of hydroxyl groups on the structure and dissolution energetics of a series of crystalline cyclic dipeptides. The calorimetrically determined energetics are interpreted in light of the crystal structures of the studied compounds. Our results indicate that the amide-amide and amide-hydroxyl hydrogen bonds both provide considerable enthalpic stability, but that the amide-amide hydrogen bond is about twice that of the amide-hydroxyl. Additionally, the interaction of the hydroxyl group with water is seen most readily in its contributions to entropy and heat capacity changes. Surprisingly, the hydroxyl group shows weakly hydrophobic behavior in terms of these contributions. These results can be used to understand the effects of mutations on the stability of globular proteins.     

Crystal structure of the Y52F/Y73F double mutant of phospholipase A2: increased hydrophobic interactions of the phenyl groups compensate for the disrupted hydrogen bonds of the tyrosines.

The enzyme phospholipase A2 (PLA2) catalyzes the hydrolysis of the sn-2 ester bond of membrane phospholipids. The highly conserved Tyr residues 52 and 73 in the enzyme form hydrogen bonds to the carboxylate group of the catalytic Asp-99. These hydrogen bonds were initially regarded as essential for the interfacial recognition and the stability of the overall catalytic network. The elimination of the hydrogen bonds involving the phenolic hydroxyl groups of the Tyr-52 and -73 by changing them to Phe lowered the stability but did not significantly affect the catalytic activity of the enzyme. The X-ray crystal structure of the double mutant Y52F/Y73F has been determined at 1.93 A resolution to study the effect of the mutation on the structure. The crystals are trigonal, space group P3(1)21, with cell parameters a = b = 46.3 A and c = 102.95 A. Intensity data were collected on a Siemens area detector, 8,024 reflections were unique with an R(sym) of 4.5% out of a total of 27,203. The structure was refined using all the unique reflections by XPLOR to a final R-factor of 18.6% for 955 protein atoms, 91 water molecules, and 1 calcium ion. The root mean square deviation for the alpha-carbon atoms between the double mutant and wild type was 0.56 A. The crystal structure revealed that four hydrogen bonds were lost in the catalytic network; three involving the tyrosines and one involving Pro-68. However, the hydrogen bonds of the catalytic triad, His-48, Asp-99, and the catalytic water, are retained. There is no additional solvent molecule at the active site to replace the missing hydroxyl groups; instead, the replacement of the phenolic OH groups by H atoms draws the Phe residues closer to the neighboring residues compared to wild type; Phe-52 moves toward His-48 and Asp-99 of the catalytic diad, and Phe-73 moves toward Met-8, both by about 0.5 A. The closing of the voids left by the OH groups increases the hydrophobic interactions compensating for the lost hydrogen bonds. The conservation of the triad hydrogen bonds and the stabilization of the active site by the increased hydrophobic interactions could explain why the double mutant has activity similar to wild type. The results indicate that the aspartyl carboxylate group of the catalytic triad can function alone without additional support from the hydrogen bonds of the two Tyr residues.     

Spectroscopic, mass spectrometry, and semiempirical investigation of a new ester of Monensin A with ethylene glycol and its complexes with monovalent metal cations

A new ester of Monensin A with ethylene glycol (MON2) has been synthesized by a new method and its ability to form complexes with Li+, Na+, and K+ cations has been studied by ESI MS, 1H and 13C NMR, FT-IR, and PM5 semiempirical methods. It is demonstrated that MON2 forms stable complexes of 1:1 stoichiometry with monovalent metal cations. The structures of the complexes are stabilized by intramolecular hydrogen bonds in which the OH groups are always involved. In the structure of MON2 the oxygen atom of the C=O ester group is involved in very weak bifurcated intramolecular hydrogen bonds with two hydroxyl groups, whereas in the complexes of MON2 with monovalent metal cations the C=O ester group is not engaged in any intramolecular hydrogen bonds. The structures of the MON2 and its complexes with Li+, Na+, and K+ cations are visualized and discussed in detail.     

The Ribbon of Hydrogen Bonds in the Three-Dimensional Structure of Globular Proteins

We report quantum mechanical computations and experimental evidence which suggest that the backbone conformation of globular proteins depends generally on the conservation of that part of the hydrogen bond network or ribbon which is joined, in general, directly to the backbone and is largely independent of the remainder of this whole network of hydrogen bonds. The familiar hydrogen bonds of the helix and the sheet form about one-half of this ribbon of hydrogen bonds. Both water molecules and hydrogen bonding side chain groups are involved in the formation of the ribbon.This view of the three-dimensional structure of globular proteins in terms of the `molecule' allows us to deal with the non-secondary structure as well as with the familiar secondary structure. It also suggests that the ribbon contains approximately the same number of hydrogen bonds within all three structures – the helix, the sheet and the coil – and that this is the reason for the ease of interconversion of these three structures.The quantum mechanical computations on hydrogen bonding suggest that delocalised water molecules which have substantial mobility are an essential part of the ribbon. This situation arises because the hydrogen bonding groups of the protein molecule are not free to move to optimise the hydrogen bonding geometries as are the oxygen atoms in the waters and ices. Such delocalised water molecules either have high B values or are invisible in the X-ray data and yet are able to form a structure which is as strong as a normal hydrogen bond.The experimental data on the point mutations of the THRI57 residue of the T4 phage lysome provides an initial test of this model. Both the local backbone conformation and the ribbon of hydrogen bonds are conserved throughout all the mutations of residue 157,providing that the delocalised water molecules are accepted as a genuine part of the structure. These mutations include the introduction of hydrocarbon side chains at position 157 when water molecules or other side chain groups take over the formation of the hydrogen bonds.We suggest that, provided steric effects are not important, many point mutations succeed because they leave the ribbon of hydrogen bonds (and so the backbone conformation) largely unchanged.     

The Influence of Simple Phenols on Collagen Type I Fibrillogenesis in vitro

The influence of phenolic compounds with different numbers of hydroxy groups (phenol, pyrocatechol, resorcinol, and pyrogallol) on the kinetics of in vitro fibrillogenesis of collagen and on fibril structure has been studied. It has been shown that these phenols accelerate fibril formation mainly by shortening the lag phase, presumably facilitating the formation of collagen dimers and their subsequent association to linear aggregates. The accelerating activity of phenols is proportional to the number of hydroxy groups in the molecule. It increases in the series: phenol < resorcinol < pyrogallol. Therefore, the ability of phenols to accelerate fibril formation is likely to stem from the formation of hydrogen bonds with amino-acid residues in collagen chains. The hydrogen bonds may stabilize the structure of the intermediates, facilitating their interaction during fibrillogenesis.     

Refined structure of alpha-lytic protease at 1.7 A resolution. Analysis of hydrogen bonding and solvent structure

The structure of alpha-lytic protease, a serine protease produced by the bacterium Lysobacter enzymogenes, has been refined at 1.7 A resolution. The conventional R-factor is 0.131 for the 14,996 reflections between 8 and 1.7 A resolution with I greater than or equal to 2 sigma (I). The model consists of 1391 protein atoms, two sulfate ions and 156 water molecules. The overall root-meansquare error is estimated to be about 0.14 A. The refined structure was compared with homologous enzymes alpha-chymotrypsin and Streptomyces griseus protease A and B. A new sequence numbering was derived based on the alignment of these structures. The comparison showed that the greatest structural homology is around the active site residues Asp102, His57 and Ser195, and that basic folding pathways are maintained despite chemical changes in the hydrophobic cores. The hydrogen bonds in the structure were tabulated and the distances and angles of interaction are similar to those found in small molecules. The analysis also revealed the presence of close intraresidue interactions. There are only a few direct intermolecular hydrogen bonds. Most intermolecular interactions involve bridging solvent molecules. The structural importance of hydrogen bonds involving the side-chain of Asx residues is discussed. All the negatively charged groups have a counterion nearby, while the excess positively charged groups are exposed to the solvent. One of the sulfate ions is located near the active site, whereas the other is close to the N terminus. Of the 156 water molecules, only seven are not involved in a hydrogen bond. Six of these have polar groups nearby, while the remaining one is in very weak density. There are nine internal water molecules, consisting of two monomers, two dimers and one trimer. No significant second shell of solvent is observed.     

Computational full electron structure study of biological activity in Cyclophilin A

Cyclosporine (CsA) is widely used in organ transplant patients to help prevent the patient’s body from rejecting the organ. CsA has been shown to be a safe and highly effective immunosuppressive drug that binds with the protein Cyclophilin A (CypA) at active sites. However, the exact mechanism of this binding at the molecular level remains unknown. In this project, we elucidate the binding of CsA to CypA at the molecular level by computing their electron structures and revealing their interactions. We employ a novel technique called electron Computer-Aided Drug Design (eCADD) on the protein’s full electron structure along with its hydrophobic pocket and the perturbation theory of the interaction between two wave functions. We have identified the wave function of CypA, the biological active residues and active atoms of CypA and CsA, the interaction site between CypA and CsA, and the hydrogen bonds in the ligand CsA binding site. All these calculated active residues, active atoms, and hydrogen bonds are in good agreement with recorded laboratory experiments and provide guidelines for designing new ligands of CypA. We believe that our eCADD framework can provide researchers with a cost-efficient new method of drug design based on the full electron structure of proteins.