Lone pair ··· π interactions between water oxygens and aromatic residues: Quantum chemical studies based on high‐resolution protein structures and model compounds

The π electron cloud of aromatic centers is known to be involved in several noncovalent interactions such as C—H···π, O—H···π, and π···π interactions in biomolecules. Lone-pair (lp) ··· π interactions have gained attention recently and their role in biomolecular structures is being recognized. In this article, we have carried out systematic analysis of high-resolution protein structures and identified more than 400 examples in which water oxygen atoms are in close contact (distance < 3.5 Å) with the aromatic centers of aromatic residues. Three different methods were used to build hydrogen atoms and we used a consensus approach to find out potential candidates for lp···π interactions between water oxygen and aromatic residues. Quantum mechanical calculations at MP2/6-311++G(d,p) level on model systems based on protein structures indicate that majority of the identified examples have energetically favorable interactions. The influence of water hydrogen atoms was investigated by sampling water orientations as a function of two parameters: distance from the aromatic center and the angle between the aromatic plane and the plane formed by the three water atoms. Intermolecular potential surfaces were constructed using six model compounds representing the four aromatic amino acids and 510 different water orientations for each model compound. Ab initio molecular orbital calculations at MP2/6-311++G(d,p) level show that the interaction energy is favorable even when hydrogen atoms are farthest from the aromatic plane while water oxygen is pointing toward the aromatic center. The strength of such interaction depends upon the distance of water hydrogen atoms from the aromatic substituents. Our calculations clearly show that the lp···π interactions due to the close approach of water oxygen and aromatic center are influenced by the positions of water hydrogen atoms and the aromatic substituents.     

Secondary structure of Aβ(1–16) complexes with zinc: A study in the gas phase using deuterium/hydrogen exchange and ultra-high-resolution mass spectrometry

Complexes of peptide fragment 1–16 of beta-amyloid with transition metals play an important role in the development of a broad class of neurodegenerative diseases, which determines the interest in investigating the structures of these complexes. In this work, we have applied the method of the deuterium/hydrogen exchange in combination with ultra-high-resolution mass spectrometry to study conformational changes in (1–16) beta-amyloid peptide induced by binding of zinc(II) atoms. The efficiency of the deuterium/hydrogen exchange depended on the number of zinc atoms bound to the peptide and on the temperature of the ionization source region. Deuterium/hydrogen exchange reactions have been performed directly in the ionization source. The number of exchanges decreased considerably with an increasing numbers of zinc atoms. The relationship has been described with a damped exponential curve, which indicated that the binding of zinc atoms altered the conformation of the peptide ion by making it less open, which limits the access to inner areas of the molecule.     

Complete oxidation of catechol by the strictly anaerobic sulfate-reducing Desulfobacterium catecholicum sp. nov.

From an anaerobic enrichment culture with vanillate as substrate, a catechol-degrading lemon-shaped nonsporing sulfate-reducing bacterium, strain NZva20, was isolated in pure culture. Growth occurred in defined, bicarbonate-buffered, sulfide-reduced freshwater medium with catechol as sole electron donor and carbon source. Catechol was completely oxidized to CO₂ with an average growth yield of 31 g cell dry mass per mol of catechol, corresponding to 9.5 g cell dry mass per mol of sulfate reduced. Further substrates utilized as electron donors and carbon sources were resorcinol, hydroquinone, benzoate and several other aromatic compounds, hydrogen plus carbon dioxide, formate, lactate, pyruvate, alcohols including methanol, dicarboxylic acids, acetate, propionate and higher fatty acids up to 18 carbon atoms. Instead of sulfate, sulfite, thiosulfate, dithionite or nitrate served as electron acceptors. Nitrate was reduced to ammonium. Strain NZva20 is the first bacterium in which the complete oxidation of organic substrates is linked to the ammonification of nitrate. Elemental sulfur was not utilized as electron acceptor. In the absence of an electron acceptor slow growth occurred on pyruvate or fumarate. The G+C content of the DNA of strain NZva20 was 52.4 mol%. Cytochromes were present. Desulfoviridin could not be detected. Strain NZva20 is described as type strain of a new species, Desulfobacterium catecholicum sp. nov.Affectionately dedicated to Professor Ralph S. Wolfe on the occassion of his 65th birthday     

锌金属蛋白酶家族介绍及结构机理研究进展

锌金属蛋白酶是一类分布广泛,种类繁多的水解酶家族,是近年来人们研究的热点．一般来说,HEXXH保守序列一直都作为锌金属蛋白酶家族分类的依据．除此之外,不同种类的酶还有一些其它的序列特征．谷氨酸锌蛋白（gluzincin）在锌离子配体处共有2组保守区域;甲硫氨酸锌蛋白（metzincin）则拥有1个延长的保守序列HEBXHXBGBXH,里面包含了第3个His配体．大量的金属蛋白酶晶体结构被解析出来,从中可以发现该类酶的活性中心都包含有1个锌离子．锌金属蛋白酶的催化机理通常认为是锌离子与1个水分子结合活化而成．但是,最近的发现证实了并不是所有的锌金属蛋白酶催化都需要水分子参与活化．在本文中,综述了这些已发表的锌金属蛋白酶家族的分类、结构特征、底物特异性识别和催化机理．     

Synergism between coenzyme and alcohol binding to liver alcohol dehydrogenase

Heterotropic cooperativity effects in the binding of alcohols and NAD+ or NADH to liver alcohol dehydrogenase have been examined by equilibrium measurements and stopped-flow kinetic studies. Equilibrium data are reported for benzyl alcohol, 2-chloroethanol, 2,2-dichloroethanol, and trifluoroethanol binding to free enzyme over the pH range 6-10. Binary-complex formation between enzyme and alcohols leads to inner-sphere coordination of the alcohol to catalytic zinc and shows a pH dependence reflecting the ionization states of zinc-bound water and the zinc-bound alcohol. The affinity of the binding protonation state of the enzyme for unionized alcohols increases approximately by a factor of 10 on complex formation between enzyme and NAD+ or NADH. The rate and kinetic cooperativity with coenzyme binding of the alcohol association step indicates that enzyme-bound alcohols participate in hydrogen bonding interactions which affect the rates of alcohol and coenzyme equilibration with the enzyme without providing any pronounced contribution to the net energetics of alcohol binding. The pKa values determined for alcohol deprotonation at the binary-complex level are linearly dependent on those of the free alcohols, and can be readily reconciled with the pKa values attributed to ionization of zinc-bound water. Alcohol coordination to catalytic zinc provides a major contribution to the pKa shift which ensures that the substrate is bound predominantly as an alcoholate ion in the catalytically productive ternary complex at physiological pH. The additional pKa shift contributed by NAD+ binding is less pronounced, but may be of particular mechanistic interest since it increases the acidity of zinc-bound alcohols relatively to that of zinc-bound water.     

Structure of the zinc-binding site in the crystal structure of a zinc endoprotease from Streptomyces caespitosus at 1 A resolution

A zinc endoprotease produced by Streptomyces caespitosus (ScNP) contains a H83E84TGH87VLG90LPD93-Met103 sequence. Except for D93, this amino acid sequence is the same as a characteristic consensus HEXXHXXGXXH-M motif found in one class of zinc endoprotease called 'metzincins'. We analyzed the structural and functional role of the consensus sequence located around a catalytically essential zinc ion based on the crystal structure of ScNP. The structure was determined at the highest level on resolution (1 A resolution) and accuracy among crystal structures of zinc endoproteases ever determined. The zinc ion of ScNP is tetrahedrally coordinated by three amino acid side-chains (H83, H87 and D93) and a water molecule. The distances between the zinc ion and the coordinating atoms are 2.01, 2.01 and 1.95 A for H83N epsilon, H87N epsilon and D93O delta, respectively. These distances agree very well with those normally found in crystal structures of zinc-containing small molecules in the Cambridge Structural Database. On the other hand, the distance between the zinc ion and the coordinating water molecule (1.93 A) is slightly shorter than the typical value (2.01 A) found in the database. In addition, E84O epsilon makes a short hydrogen bond to this water molecule with the distance of 2.54 A. Two hydrogen bonds (H83N delta-L102O, H87N delta-L91O) and van der Waals interactions between the side-chain of M103 and the two imidazole rings of H83 and H87 are also observed. These interactions are probably important for the imidazole rings to construct the tetrahedral coordination arrangement toward the zinc ion.     

Do water molecules mediate protein-DNA recognition?

A comprehensive analysis of interfacial water molecules in the structures of 109 unique protein-DNA complexes is presented together with a new view on their role in protein-DNA recognition. Location of interfacial water molecules as reported in the crystal structures and as emerging from a series of molecular dynamics studies on protein-DNA complexes with explicit solvent and counterions, was analyzed based on their acceptor, donor hydrogen bond relationships with the atoms and residues of the macromolecules, electrostatic field calculations and packing density considerations. Water molecules for the purpose of this study have been categorized into four classes: viz. (I) those that contact both the protein and the DNA simultaneously and thus mediate recognition directly; (II) those that contact either the protein or the DNA exclusively via hydrogen bonds solvating each solute separately; (III) those that contact the hydrophobic groups in either the protein or the DNA; and, lastly (IV) those that contact another water molecule. Of the 17,963 crystallographic water molecules under examination, about 6% belong to class I and 76% belong to class II. About three-fourths of class I and class II water molecules are exclusively associated with hydrogen bond acceptor atoms of both protein and DNA. Noting that DNA is polyanionic, it is significant that a majority of the crystallographically observed water molecules as well as those from molecular dynamics simulations should be involved in facilitating binding by screening unfavorable electrostatics. Less than 2% of the reported water molecules occur between hydrogen bond donor atoms of protein and acceptor atoms of DNA. These represent cases where protein atoms cannot reach out to DNA to make favorable hydrogen bond interactions due to packing/structural restrictions and interfacial water molecules provide an extension to side-chains to accomplish hydrogen bonding.     

The structure of 2Zn pig insulin crystals at 1.5 A resolution

The paper describes the arrangement of the atoms within rhombohedral crystals of 2Zn pig insulin as seen in electron density maps calculated from X-ray data extending to 1.5 A (1 A = 10(-10) m = 10(-1) nm) at room temperature and refined to R = 0.153. The unit cell contains 2 zinc ions, 6 insulin molecules and about 3 x 283 water molecules. The atoms in the protein molecules appear well defined, 7 of the 102 side chains in the asymmetric unit have been assigned alternative disordered positions. The electron density over the water molecules has been interpreted in terms of 349 sites, 217 weighted 1.0, 126 weighted 0.5, 5 at 0.33 and 1 at 0.25 giving ca. 282 molecules. The positions and contacts of all the residues belonging to the two A and B chains of the asymmetric unit are shown first and then details of their arrangement in the two insulin molecules, 1 and 2, which are different. The formation from these molecules of a compact dimer and the further aggregation of three dimers to form a hexamer around two zinc ions, follows. It appears that in the packing of the hexamers in the crystal there are conflicting influences; too-close contacts between histidine B5 residues in neighbouring hexamers are probably responsible for movements of atoms at the beginning of the A chain of one of the two molecules of the dimer that initiate movements in other parts, particularly near the end of the B chain. At every stage of the building of the protein structure, residues to chains of definite conformation, molecules, dimers, hexamers and crystals, we can trace the effect of the packing of like groups to like, aliphatic groups together, aromatic groups together, hydrogen-bonded structures, positive and negative ions. Between the protein molecules, the water is distributed in cavities and channels that are continuous throughout the crystals. More than half the water molecules appear directly hydrogen bonded to protein atoms. These are generally in contact with other water molecules in chains and rings of increasing disorder, corresponding with their movement through the crystals. Within the established crystal structure we survey next the distribution of hydrogen bonds within the protein molecules and between water and protein and water and water; all but eight of the active atoms in the protein form at least one hydrogen bond.(ABSTRACT TRUNCATED AT 400 WORDS)     

Effect of pH on the liver alcohol dehydrogenase reaction.

New transient kinetic methods, which allow kinetics to be carried out under conditions of excess substrate, have been employed to investigate the kinetics of hydride transfer from NADH to aromatic aldehydes and from aromatic alcohols to NAD+ as a function of pH. The hydride transfer rate from 4-deuterio-NADH to beta-naphthaldehyde is nearly pH independent from pH 6.0 to pH 9.9; the isotope effect is also pH independent with kappa-H/kappaD congruent to 2.3. Likewise, the rate of oxidation of benzyl alcohol by NAD+ changes little with pH between pH 8.75 and pH 5.9; the isotope effect for this process is between 3.0 and 4.4. Earlier substituent effect studies on the reduction of aromatic aldehydes were consistent with electrophilic catalysis by either zinc or a protonic acid. The pH independence of hydride transfer is consistent with electrophilic catalysis by zinc since such catalysis by protonic acid (with a pK between 6.0 and 10.0) would show strong pH dependence. However, protonic acid catalysis cannot be excluded if the pKa of the acid catalyst in the ternary NADH-E-RCOH complex were smaller than 6.0 or smaller than 10.0. The two kinetic parameters changing significantly with pH are the kinetic binding constant for ternary complex formation with aromatic alcohol and the rate of dissociation of aromatic alcohols from enzyme. This is consistent with base-catalyzed removal of a proton from alcohol substrated and consequent acid catalysis of protonation of a zinc-alcoholate complex. The equilibrium constant for hydride transfer from benzaldehyde to benzyl alcohol at pH 8.75 is K-eq equals kappa-H/kappa-H equals 42; this constant has important consequences concerning subunit interactions during liver alcohol dehydrogenase catalysis.     

Oxidation of Alcohols by Botrytis cinerea

Crude cell-free preparations of Botrytis cinerea were found to oxidize straight-chain primary alcohols (except methanol), aromatic primary alcohols, and unsaturated primary alcohols. The resulting products were the corresponding aldehydes and an equal molar quantity of hydrogen peroxide.     

Effect of alcohol compounds found in hemicellulose hydrolysate on the growth and fermentation of ethanologenic Escherichia coli

Lignocellulose can be readily hydrolyzed into a mixture of sugars using dilute mineral acids. During hydrolysis, a variety of inhibitors are also produced which include aromatic alcohols from lignin and furfuryl alcohol from pentose destruction. Seven compounds were investigated individually and in binary combinations (catechol, coniferyl alcohol, furfuryl alcohol, guaiacol, hydroquinone, methylcatechol, and vanillyl alcohol). Aromatic alcohols and furfuryl alcohol inhibited ethanol production from xylose in batch fermentations primarily by inhibiting the growth of Escherichia coli LY01, the biocatalyst. The toxicities of these compounds were directly related to their hydrophobicity. Methylcatechol was the most toxic compound tested (MIC = 1.5 g/L). In binary combination, the extent of growth inhibition was roughly additive for most compounds tested. However, combinations with furfuryl alcohol and furfural (furaldehyde) appear synergistic in toxicity. When compared individually, alcohol components which are formed during hemicellulose hydrolysis are less toxic for growth than the aldehydes and organic acids either on a weight basis or a molar basis.     

Apo and Holo structures of an NADPH-dependent cinnamyl alcohol dehydrogenase from Saccharomyces cerevisiae

The crystal structure of Saccharomyces cerevisiae ScAdh6p has been solved using the anomalous signal from the two zinc atoms found per subunit, and it constitutes the first structure determined from a member of the cinnamyl alcohol dehydrogenase family. ScAdh6p subunits exhibit the general fold of the medium-chain dehydrogenases/reductases (MDR) but with distinct specific characteristics. In the three crystal structures solved (two trigonal and one monoclinic), ScAdh6p molecules appear to be structural heterodimers composed of one subunit in the apo and the second subunit in the holo conformation. Between the two conformations, the relative disposition of domains remains unchanged, while two loops, Cys250-Asn260 and Ile277-Lys292, experience large movements. The apo-apo structure is disfavoured because of steric impairment involving the loop Ile277-Lys292, while in the holo-holo conformation some of the hydrogen bonds between subunits would break apart. These suggest that the first NADPH molecule would bind to the enzyme much more tightly than the second. In addition, fluorimetric analysis of NADPH binding demonstrates that only one cofactor molecule binds per dimer. Therefore, ScAdh6p appears to function according to a half-of-the-sites reactivity mechanism, resulting from a pre-existing (prior to cofactor binding) tendency for the structural asymmetry in the dimer. The specificity of ScAdh6p towards NADPH is mainly due to the tripod-like interactions of the terminal phosphate group with Ser210, Arg211 and Lys215. The size and the shape of the substrate-binding pocket correlate well with the substrate specificity of ScAdh6p towards cinnamaldehyde and other aromatic compounds. The structural relationships of ScAdh6p with other MDR structures are analysed.     

Cloning and Sequencing of a Novel meta-Cleavage Dioxygenase Gene Whose Product Is Involved in Degradation of γ-Hexachlorocyclohexane in Sphingomonas paucimobilis

Sphingomonas (formerly Pseudomonas) paucimobilis UT26 utilizes gamma-hexachlorocyclohexane (gamma-HCH), a halogenated organic insecticide, as a sole source of carbon and energy. In a previous study, we showed that gamma-HCH is degraded to chlorohydroquinone (CHQ) and then to hydroquinone (HQ), although the rate of reaction from CHQ to HQ was slow (K. Miyauchi, S. K. Suh, Y. Nagata, and M. Takagi, J. Bacteriol. 180:1354-1359, 1998). In this study, we cloned and characterized a gene, designated linE, which is located upstream of linD and is directly involved in the degradation of CHQ. The LinE protein consists of 321 amino acids, and all of the amino acids which are reported to be essential for the activity of meta-cleavage dioxygenases are conserved in LinE. Escherichia coli overproducing LinE could convert both CHQ and HQ, producing gamma-hydroxymuconic semialdehyde and maleylacetate, respectively, with consumption of O(2) but could not convert catechol, which is one of the major substrates for meta-cleavage dioxygenases. LinE seems to be resistant to the acylchloride, which is the ring cleavage product of CHQ and which seems to react with water to be converted to maleylacetate. These results indicated that LinE is a novel type of meta-cleavage dioxygenase, designated (chloro)hydroquinone 1, 2-dioxygenase, which cleaves aromatic rings with two hydroxyl groups at para positions preferably. This study represents a direct demonstration of a new type of ring cleavage pathway for aromatic compounds, the hydroquinone pathway.     

Substrate-induced deactivation of penicillinases. Studies of beta-lactamase I by hydrogen exchange.

The conformational motility of beta-lactamase I from Bacillus cereus was studied by hydrogen exhange. The time course of the isotopic replacement of peptide hydrogen atoms was followed by 'exchange-in' or 'exchange-out' experiments. Many of the substrates for this enzyme that have o-substituted aromatic or heterocyclic side chains (e.g. methicillin or cloxacillin) are known to effect a decrease in enzymic activity ('substrate-induced deactivation'). There was a marked discontinuity in the exchange-out curve when methicillin or cloxacillin was diffused into the enzyme solution. About one-half of the hydrogen atoms that were probed were affected by the presence of these substrates, and the change in the reactivity of the hydrogen atoms was also large. Substrates that do not bring about deactivation (benzylpenicillin and cephalosporin C) do not affect the hydrogen exchange, nor do reversible competitive inhibitors such as the penicilloic acid or penilloic acid. On the other hand, Zn2+ ions do affect the hydrogen exchange; their effect is similar to that of methicillin or cloxacillin.     

ACOMPARATIVE STUDY ON THE ELECTROANTENOGRAM RESPONSES OF THREE APHID SPECIES TO PLANT VOLATILES

Abstract  Electroantennogram responses (EAGs) to 39 plant volatiles from Aphis glycines, A. gossypii and Myzus persicae were recorded. The olfactory system of these three aphid species differentially responded to green leaf, aromatic, mono-and sequiterpene chemicals. The EAG depolarizatin to green leaf chemicals and their isomers were more pronounced than to terpenes, alcohols and aldehydes, 6 carbon atom compounds being the most effective. It was shown that there was certain difference in the olfactory spectra of these three aphid species to plant volatiles. Higher responses in A. glycines and A. gossypii were elicited by the leaf aldehydes and esters than by the corresponding alcohols, which are in contrast with those in M. persicue . The activity of aphid olfactory system was related to molecular structures of the compounds. The responses to saturated alcohols and aldehydes were higher than those to unsaturated. The responses of spring emigrants of A. glycines to terpene derivatives (alcohols, aldehydes and esters) were higher than those of A. gossypii , but were lower to terpene hydrocarbons. The dissimilarity in olfactory spectra of the aphid species reflected the evolution of their host plant selection specificity. In the evolution, their chemosensory system adapted to perceive host plants. The shapes of the EAGs evoked by the various components were consistently different, with the slowest recovery for geraniol, 1-decanol and (±)-β-cit-ronellol, and ware related to number of carbon atoms in straight chain saturated alcohols.     

Analysis of zinc binding sites in protein crystal structures

The geometrical properties of zinc binding sites in a dataset of high quality protein crystal structures deposited in the Protein Data Bank have been examined to identify important differences between zinc sites that are directly involved in catalysis and those that play a structural role. Coordination angles in the zinc primary coordination sphere are compared with ideal values for each coordination geometry, and zinc coordination distances are compared with those in small zinc complexes from the Cambridge Structural Database as a guide of expected trends. We find that distances and angles in the primary coordination sphere are in general close to the expected (or ideal) values. Deviations occur primarily for oxygen coordinating atoms and are found to be mainly due to H-bonding of the oxygen coordinating ligand to protein residues, bidentate binding arrangements, and multi-zinc sites. We find that H-bonding of oxygen containing residues (or water) to zinc bound histidines is almost universal in our dataset and defines the elec-His-Zn motif. Analysis of the stereochemistry shows that carboxyl elec-His-Zn motifs are geometrically rigid, while water elec-His-Zn motifs show the most geometrical variation. As catalytic motifs have a higher proportion of carboxyl elec atoms than structural motifs, they provide a more rigid framework for zinc binding. This is understood biologically, as a small distortion in the zinc position in an enzyme can have serious consequences on the enzymatic reaction. We also analyze the sequence pattern of the zinc ligands and residues that provide elecs, and identify conserved hydrophobic residues in the endopeptidases that also appear to contribute to stabilizing the catalytic zinc site. A zinc binding template in protein crystal structures is derived from these observations.     

Antigelling and antisickling bisphenyl oligopeptides and peptide analogues have similar structural features

Single-crystal X-ray diffraction was used to determine the three-dimensional structures of two antigelling oligopeptides, L-lysyl-L-phenylalanyl-L-phenylalanine and L-phenylalanylglycylglycyl-D-phenylalanine, and two antisickling peptide analogues, L-phenylalanine benzyl ester and N-phenylacetyl-L-phenylalanine. Although these bisphenyl compounds are chemically quite different from one another, they demonstrate unusual structural similarities: The molecules have compact conformations in which the two phenyl rings are positioned approximately 5 A apart with interplanar angles approaching 90 degrees, thereby making intramolecular edge-to-face interactions. In addition, the polar atoms, nitrogen and oxygen, are in close proximity without forming intramolecular hydrogen bonds. The relative spatial distribution of polar and nonpolar atoms renders the structures compact and amphipathic. The intramolecular edge-to-face interaction between two aromatic rings, which brings a hydrogen atom with relative positive charge near the pi-electron cloud with relative negative charge, is enthalpically favorable and maintains the molecules in a compact and amphipathic conformation. Nonbonded potential energy calculations were used to characterize the energetics of the aromatic-aromatic interaction, and they showed that the observed geometry is stabilized enthalpically by a favorable interaction on the order of -1 to -2 kcal/mol. Structural differences between the two antisickling and the two antigelling agents suggest that molecular volume limits red cell membrane passage. These data provide a molecular structural framework from which to design and synthesize amphipathic bisphenyl compounds that both bind to deoxy sickle cell hemoglobin and cross the erythrocyte membrane.