Water-mediated contacts in the trp-repressor operator complex recognition process

Water-mediated contacts are known as an important recognition tool in trp-repressor operator systems. One of these contacts involves two conserved base pairs (G(6).C(-6) and A(5). T(-5)) and three amino acids (Lys 72, Ile 79, and Ala 80). To investigate the nature of these contacts, we analyzed the X-ray structure (PDB code: 1TRO) of the trp-repressor operator complex by means of molecular dynamics simulations. This X-ray structure contains two dimers that exhibit structural differences. From these two different starting structures, two 10 ns molecular dynamics simulations have been performed. Both of our simulations show an increase of water molecules in the major groove at one side of the dimer, while the other side remains unchanged compared to the X-ray structure. Though the maximum residence time of the concerned water molecules decreases with an increase of solvent at the interface, these water molecules continue to play an important role in mediating DNA-protein contacts. This is shown by new stable amino acids-DNA distances and a long water residence time compared to free DNA simulation. To maintain stability of the new contacts, the preferential water binding site on O6(G6) is extended. This extension agrees with mutation experiment data on A5 and G6, which shows different relative affinity due to mutation on these bases [A. Joachimiak, T. E. Haran, P. B. Sigler, EMBO Journal 1994, Vol. 13, No. (2) pp. 367-372]. Due to the rearrangements in the system, the phosphate of the base G6 is able to interconvert to the B(II) substate, which is not observed on the other half side of the complex. The decrease of the number of hydrogen bonds between protein and DNA backbone could be the initial step of the dissociation process of the complex, or in other words an intermediate complex conformation of the association process. Thus, we surmise that these features show the importance of water-mediated contacts in the trp-repressor operator recognition process.     

Amino acid–base interactions: a three-dimensional analysis of protein–DNA interactions at an atomic level

To assess whether there are universal rules that govern amino acid–base recognition, we investigate hydrogen bonds, van der Waals contacts and water-mediated bonds in 129 protein–DNA complex structures. DNA–backbone interactions are the most numerous, providing stability rather than specificity. For base interactions, there are significant base–amino acid type correlations, which can be rationalised by considering the stereochemistry of protein side chains and the base edges exposed in the DNA structure. Nearly two-thirds of the direct read-out of DNA sequences involves complex networks of hydrogen bonds, which enhance specificity. Two-thirds of all protein–DNA interactions comprise van der Waals contacts, compared to about one-sixth each of hydrogen and water-mediated bonds. This highlights the central importance of these contacts for complex formation, which have previously been relegated to a secondary role. Although common, water-mediated bonds are usually non-specific, acting as space-fillers at the protein–DNA interface. In conclusion, the majority of amino acid–base interactions observed follow general principles that apply across all protein–DNA complexes, although there are individual exceptions. Therefore, we distinguish between interactions whose specificities are ‘universal’ and ‘context-dependent’. An interactive Web-based atlas of side chain–base contacts provides access to the collected data, including analyses and visualisation of the three-dimensional geometry of the interactions.     

Molecular dynamics analysis of the engrailed homeodomain-DNA recognition

Molecular dynamics (MD) simulations were performed for investigating the role of Gln50 in the engrailed homeodomain-DNA recognition. Employing the crystal structure of free engrailed homeodomain and homeodomain-DNA complex as a starting structure, we carried out MD simulations of: (i) the complex between engrailed homeodomain and a 20 base-pair DNA containing TAATTA core sequence; (ii) the free engrailed homeodomain. The simulations show that homeodomain flexibility does not depend on its ligation state. The engrailed homeodomain shows similar flexibility, and the recognition helix-3 shows very similar characteristic of high rigidity and limited conformational space in two complexation states. At the same time, DNA structure has also no obvious conformational fluctuations. These results preclude the possibility of the side chain of Gln50 forming direct hydrogen bonds to the core DNA bases. MD simulations confirm a few well-conserved sites for water-mediated hydrogen bonds from protein to DNA are occupied by water molecules, and Gln50 interacts with corresponding core DNA bases through water-mediated hydrogen bonds. So Gln50 plays a relatively modest role in determining the affinity and specificity of the engrailed homeodomain. In addition, the electrostatic interaction between homeodomain and phosphate backbone of the DNA is a main factor for N- and C-terminal arm becoming ordered upon DNA binding.     

Recognition of DNA sequences by the repressor of bacteriophage 434

The structure of a complex between the DNA-binding domain of phage 434 repressor and a 14 base-pair synthetic DNA operator reveals the molecular interactions important for sequence-specific recognition. A set of contacts with DNA backbone, notably involving hydrogen bonds between peptide-NH groups and DNA phosphates, position the repressor and fix the DNA configuration. Direct interactions between amino acid side chains and DNA bases involve nonpolar van der Waals contacts as well as hydrogen bonds. The structures of the repressor domain and of the 434 cro protein are extremely similar. There appear to be no major conformational changes in the proteins when they bind to DNA.     

Mapping the active site polarity in structures of endothelial nitric oxide synthase heme domain complexed with isothioureas

Analyzing the active site topology and plasticity of nitric oxide synthase (NOS) and understanding enzyme-drug interactions are crucial for the development of potent, isoform-selective NOS inhibitors. A small hydrophobic pocket in the active site is identified in the bovine eNOS heme domain structures complexed with potent isothiourea inhibitors: seleno analogue of S-ethyl-isothiourea, S-isopropyl-isothiourea, and 2-aminothiazoline, respectively. These structures reveal the importance of nonpolar van der Waals contacts in addition to the well-known hydrogen bonding interactions between inhibitor and enzyme. The scaffold of a potent NOS inhibitor should be capable of donating hydrogen bonds to as well as making nonpolar contacts with amino acids in the NOS active site.     

Computational analysis of BACE1-ligand complex crystal structures and linear discriminant analysis for identification of BACE1 inhibitors with anti P-glycoprotein binding property

More than 100 years of research on Alzheimer’s disease didn’t yield a potential cure for this dreadful disease. Poor Blood Brain Barrier (BBB) permeability and P-glycoprotein binding of BACE1 inhibitors are the major causes for the failure of these molecules during clinical trials. The design of BACE1 inhibitors with a balance of sufficient affinity to the binding site and little or no interaction with P-glycoproteins is indispensable. Identification and understanding of protein–ligand interactions are essential for ligand optimization process. Structure-based drug design (SBDD) efforts led to a steady accumulation of BACE1-ligand crystal complexes in the PDB. This study focuses on analyses of 153 BACE1-ligand complexes for the direct contacts (hydrogen bonds and weak interactions) observed between protein and ligand and indirect contacts (water-mediated hydrogen bonds), observed in BACE1-ligand complex crystal structures. Intraligand hydrogen bonds were analyzed, with focus on ligand P-glycoprotein efflux. The interactions are dissected specific to subsites in the active site and discussed. The observed protein-ligand and intraligand interactions were used to develop the linear discriminant model for the identification of BACE1 inhibitors with less or no P-glycoprotein binding property. Excellent statistical results and model’s ability to correctly predict a new data-set with an accuracy of 92% is achieved. The results are retrospectively analyzed to give input for the design of potential BACE1 inhibitors.     

Crystal structure of paromomycin docked into the eubacterial ribosomal decoding A site

BACKGROUND: Aminoglycoside antibiotics interfere with translation in both gram-positive and gram-negative bacteria by binding to the tRNA decoding A site of the 16S ribosomal RNA. RESULTS: Crystals of complexes between oligoribonucleotides incorporating the sequence of the ribosomal A site of Escherichia coli and the aminoglycoside paromomycin have been solved at 2.5 A resolution. Each RNA fragment contains two A sites inserted between Watson-Crick pairs. The paromomycin molecules interact in an enlarged deep groove created by two bulging and one unpaired adenines. In both sites, hydroxyl and ammonium side chains of the antibiotic form 13 direct hydrogen bonds to bases and backbone atoms of the A site. In the best-defined site, 8 water molecules mediate 12 other hydrogen bonds between the RNA and the antibiotics. Ring I of paromomycin stacks over base G1491 and forms pseudo-Watson-Crick contacts with A1408. Both the hydroxyl group and one ammonium group of ring II form direct and water-mediated hydrogen bonds to the U1495oU1406 pair. The bulging conformation of the two adenines A1492 and A1493 is stabilized by hydrogen bonds between phosphate oxygens and atoms of rings I and II. The hydrophilic sites of the bulging A1492 and A1493 contact the shallow groove of G=C pairs in a symmetrical complex. CONCLUSIONS: Water molecules participate in the binding specificity by exploiting the antibiotic hydration shell and the typical RNA water hydration patterns. The observed contacts rationalize the protection, mutation, and resistance data. The crystal packing mimics the intermolecular contacts induced by aminoglycoside binding in the ribosome.     

Importance of hydrophobic and polar residues in ligand binding in the family 15 carbohydrate-binding module from Cellvibrio japonicus Xyn10C

Modular glycoside hydrolases that degrade the plant cell wall often contain noncatalytic carbohydrate-binding modules (CBMs) that interact with specific polysaccharides within this complex macromolecule. CBMs, by bringing the appended catalytic module into intimate and prolonged association with the substrate, increase the rate at which these enzymes are able to hydrolyze glycosidic bonds. Recently, the crystal structure of the family 15 CBM (CBM15) from Cellvibrio japonicus (formerly Pseudomonas cellulosa) Xyn10C was determined in complex with the ligand xylopentaose. In this report we have used a rational design approach, informed by the crystal structure of the CBM15-ligand complex, to probe the importance of hydrophobic stacking interactions and both direct and water-mediated hydrogen bonds in the binding of this protein to xylan and xylohexaose. The data show that replacing either Trp 171 or Trp 186, which stack against xylose residues n and n + 2 in xylopentaose, with alanine abolished ligand binding. Similarly, replacing Asn 106, Gln 171, and Gln 217, which make direct hydrogen bonds with xylopentaose, with alanine greatly reduced the affinity of the protein for its saccharide ligands. By contrast, disrupting water-mediated hydrogen bonds between CBM15 and xylopentaose by introducing the mutations S108A, Q167A, Q221A, and K223A had little effect on the affinity of the protein for xylan or xylohexaose. These data indicate that CBM15 binds xylan and xylooligosaccharides via the same interactions and provide clear evidence that direct hydrogen bonds are a key determinant of affinity in a type B CBM. The generic importance of these data is discussed.     

Double coned inverse sandwich complexes [M-(η4-C4H4)-M′] of Gr-IA and Gr-IIA Metals: theoretical study of electronic of structure and second hyperpolarizability

Molecular interaction between dioxane and methanol involves certain polar and nonpolar bonding to form a one to one complex. Interatomic distances between hydrogen and oxygen within 3 Å have been considered as hydrogen bonding. Optimizations of the structures of dioxane-methanol complexes were carried out considering any spatial orientation of a methanol molecule around a chair/boat/twisted-boat conformation of dioxane. From 45 different orientations of dioxane and water, 23 different structures with different local minima were obtained and the structural characteristics like interatomic distances, bond angles, dihedral angles, dipole moment of each complex were discussed. The most stable structure, i.e., with minimum heat of formation is found to have a chair form dioxane, one O-H…O, and two C-H…O hydrogen bonds. In general, the O-H…O hydrogen bonds have an average distance of 1.8 Å while C-H…O bonds have 2.6 Å. The binding energy of the dioxane-methanol complex is found to be a linear function of number of O-H…O and C-H…O bonds, and hydrogen bond length. Graphical Abstract

Sixteen orientations of methanol around dioxane converge to six local minima including the global minima with one H-O…H and two C-H…O hydrogen bonds.     

Crystal structure of the p-hydroxybenzoate hydroxylase-substrate complex refined at 1.9 A resolution. Analysis of the enzyme-substrate and enzyme-product complexes

Using synchrotron radiation, the X-ray diffraction intensities of crystals of p-hydroxy-benzoate hydroxylase, complexed with the substrate p-hydroxybenzoate, were measured to a resolution of 1.9 A. Restrained least-squares refinement alternated with rebuilding in electron density maps yielded an atom model of the enzyme-substrate complex with a crystallographic R-factor of 15.6% for 31,148 reflections between 6.0 and 1.9 A. A total of 330 solvent molecules was located. In the final model, only three residues have deviating phi-psi angle combinations. One of them, the active site residue Arg44, has a well-defined electron density and may be strained to adopt this conformation for efficient catalysis. The mode of binding of FAD is distinctly different for the different components of the coenzyme. The adenine ring is engaged in three water-mediated hydrogen bonds with the protein, while making only one direct hydrogen bond with the enzyme. The pyrophosphate moiety makes five water-mediated versus three direct hydrogen bonds. The ribityl and ribose moieties make only direct hydrogen bonds, in all cases, except one, with side-chain atoms. The isoalloxazine ring also makes only direct hydrogen bonds, but virtually only with main-chain atoms. The conformation of FAD in p-hydroxybenzoate hydroxylase is strikingly similar to that in glutathione reductase, while the riboflavin-binding parts of these two enzymes have no structural similarity at all. The refined 1.9 A structure of the p-hydroxybenzoate hydroxylase-substrate complex was the basis of further refinement of the 2.3 A structure of the enzyme-product complex. The result was a final R-factor of 16.7% for 14,339 reflections between 6.0 and 2.3 A and an improved geometry. Comparison between the complexes indicated only small differences in the active site region, where the product molecule is rotated by 14 degrees compared with the substrate in the enzyme-substrate complex. During the refinements of the enzyme-substrate and enzyme-product complexes, the flavin ring was allowed to bend or twist by imposing planarity restraints on the benzene and pyrimidine ring, but not on the flavin ring as a whole. The observed angle between the benzene ring and the pyrimidine ring was 10 degrees for the enzyme-substrate complex and 19 degrees for the enzyme-product complex. Because of the high temperature factors of the flavin ring in the enzyme-product complex, the latter value should be treated with caution. Six out of eight peptide residues near the flavin ring are oriented with their nitrogen atom pointing towards the ring.(ABSTRACT TRUNCATED AT 400 WORDS)     

Solution structure of the nogalamycin-DNA complex

The nogalamycin-d(A-G-C-A-T-G-C-T) complex (two drugs per duplex) has been generated in aqueous solution and its structure characterized by a combined application of two-dimensional NMR experiments and molecular dynamics calculations. Two equivalents of nogalamycin binds to the self-complementary octanucleotide duplex with retention of 2-fold symmetry in solution. We have assigned the proton resonances of nogalamycin and the d(A1-G2-C3-A4-T5-G6-C7-T8) duplex in the complex and identified the intermolecular proton-proton NOEs that define the alignment of the antitumor agent at its binding site on duplex DNA. The analysis was greatly aided by a large number of intermolecular NOEs involving exchangeable protons on both the nogalamycin and the DNA in the complex. The molecular dynamics calculations were guided by 274 intramolecular nucleic acid distance constraints, 90 intramolecular nogalamycin distance constraints, and 104 intermolecular distance constraints between nogalamycin and the nucleic acid protons in the complex. The aglycon chromophore intercalates at (C-A).(T-G) steps with the long axis of the aglycon approximately perpendicular to the long axis of the flanking C3.G6 and A4.T5 base pairs. The aglycon selectively stacks over T5 and G6 on the T5-G6-containing strand with the aglycon edge containing OH-4 and OH-6 substituents directed toward the C3-A4-containing strand. The C3.G6 and A4.T5 base pairs are intact but buckled at the intercalation site with a wedge-shaped alignment of C3 and A4 on the C3-A4 strand compared to the parallel alignment of T5 and G6 on the T5-G6 strand in the complex. The nogalose sugar in a chair conformation, the aglycon ring A in a half-chair conformation, and the COOCH3-10 side chain form a continuous domain that is sandwiched within the walls of the minor groove and spans the three base pair (G2-C3-A4).(T5-G6-C7) segment. The nogalose ring is positioned in the minor groove such that its nonpolar face is directed toward the G6-C7 sugar-phosphate backbone while its polar face containing OCH3 groups is directed toward the G2-C3 sugar-phosphate backbone in the complex. The intermolecular contacts include a nonpolar patch of aglycon (CH3-9) and nogalose (CH3-3') methyl groups forming van der Waals contacts with the base-sugar residues in the minor groove and intermolecular hydrogen bonds involving the amino groups of G2 and G6 with the ether oxygens OCH3-3' and O7, respectively, on the nogalose sugar.(ABSTRACT TRUNCATED AT 400 WORDS)     

Discovering the interaction propensities of amino acids and nucleotides from protein-RNA complexes

With the availability of many genome sequences, the mining of biological data is attracting much attention, most of it limited to the sequences of macromolecules. Sequence data are easy to analyze as they can be treated as strings of characters, whereas the structure of a macromolecule is much more complex. We developed a set of algorithms to analyze the structures of protein-RNA complexes at the atomic level and used them to analyze protein-RNA interactions using structural data on 51 protein-RNA complexes. The analysis revealed, among other things, that: (1) polar and charged amino acids have a strong tendency to interact with nucleotides, (2) arginine and asparagine tend to hydrogen bond with uracil, and (3) histidine favors uracil in water-mediated bonding with RNA. We analyzed a large set of structural data of protein-RNA complexes involving water-mediated hydrogen bonds as well as direct hydrogen bonds. The interaction patterns discovered from the analysis provide useful information for predicting the structure of RNA that binds proteins, and of proteins that bind RNA.     

An asymmetric complex of restriction endonuclease MspI on its palindromic DNA recognition site

Most well-known restriction endonucleases recognize palindromic DNA sequences and are classified as Type IIP. Due to the recognition and cleavage symmetry, Type IIP enzymes are usually found to act as homodimers in forming 2-fold symmetric enzyme-DNA complexes. Here we report an asymmetric complex of the Type IIP restriction enzyme MspI in complex with its cognate recognition sequence. Unlike any other Type IIP enzyme reported to date, an MspI monomer and not a dimer binds to a palindromic DNA sequence. The enzyme makes specific contacts with all 4 base pairs in the recognition sequence, by six direct and five water-mediated hydrogen bonds and numerous van der Waal contacts. This MspI-DNA structure represents the first example of asymmetric recognition of a palindromic DNA sequence by two different structural motifs in one polypeptide. A few possible pathways are discussed for MspI to cut both strands of DNA, either as a monomer or dimer.     

Structural evidence for a dominant role of nonpolar interactions in the binding of a transport/chemosensory receptor to its highly polar ligands.

The receptor, a maltose/maltooligosaccharide-binding protein, has been found to be an excellent system for the study of molecular recognition because its polar and nonpolar binding functions are segregated into two globular domains. The X-ray structures of the "closed" and "open" forms of the protein complexed with maltose and maltotetraitol have been determined. These sugars have approximately 3 times more accessible polar surface (from OH groups) than nonpolar surface (from small clusters of sugar ring CH bonds). In the closed structures, the oligosaccharides are buried in the groove between the two domains of the protein and bound by extensive hydrogen bonding interactions of the OH groups with the polar residues confined mostly in one domain and by nonpolar interactions of the CH clusters with four aromatic residues lodged in the other domain. Substantial contacts between the sugar hydroxyls and aromatic residues are also formed. In the open structures, the oligosaccharides are bound almost exclusively in the domain rich in aromatic residues. This finding, along with the analysis of buried surface area due to complex formations in the open and closed structures, supports a major role for nonpolar interactions in initial ligand binding even when the ligands have significantly greater potential for highly specific polar interactions.     

紫胶虫的生物学研究

紫胶虫Laccifer lacca(Kerr)Targ.在云南自然分布区一年发生两代,各世代有涌散、固定、泌胶、泌蜡和排泄蜜露等活动。由于幼虫在饥饿状态下的存活期较短,迁移能力较差,必须及时实行人工放养和科学管理,才能获得紫胶高产。本文报道了紫胶虫各世代的泌胶量、生殖力和性比,并提出发展紫胶生产的建议。     

Conformational analysis by NMR and distance geometry techniques of a peptide mimetic of the third helix of the Antennapedia homeodomain.

The Antennapedia homeodomain structure consists of four helices. The helices II and III are connected by a tripeptide that forms a turn, and constitute the well-known helix-turn-helix motif. The recognition helix penetrates the DNA major groove, gives specific protein-DNA contacts and forms direct, or water-mediated, intermolecular hydrogen bonds. It was suggested that helix III (and perhaps also helix IV) might represent the recognition helix of Antennapedia homeodomain, which makes contact with the surface of the major groove of the DNA. In an attempt to clarify the helix III capabilities of assuming an helical conformation when separated from the rest of the protein, we carried out the structural determination of the recognition helix III in different solvent media. The conformational study of fragments 42-53, where residues W48 and F49, not involved in the protein-DNA interaction, were substituted by two alanines, was conducted in sodium dodecyl sulfate (SDS), trifluoroethanol (TFE) and TFE/water, using circular dichroism, nuclear magnetic resonance (NMR) and distance geometry (DG) techniques. The fragment assumes a well-defined secondary structure in TFE and in TFE/water (90/10, v/v) with an alpha-helix encompassing residues 4-9, while in TFE/water (70/30, v/v) a less regular structure was found. The DG results in the micellar system evidence the presence of a distorted alpha-helical conformation involving residues 4-8. Our results reveal that the isolated Antennapedia recognition helix III tend to preserve in solution the alpha-helical conformation even if separated from the rest of the molecule.