Evidence for nonhydrogen bonded compound II in cyclic reaction of hemoglobin I from Lucina pectinata with hydrogen peroxide

Studies that elucidate the behavior of the hemoglobins (Hbs) and myoglobins upon reaction with hydrogen peroxide are essential to the development of oxygen carrier substitutes. Stopped-flow kinetics and resonance Raman data show that the reaction between hydrogen peroxide and oxygenated and deoxygenated ferric Hb I (oxy- and deoxy-HbI) from Lucina pectinata produce compound I and compound II ferryl species. The rate constants ratio (k23/k41) between the formation of compound II from compound I (k23) and the oxidation of the ferrous HbI (k41, i.e., 25 M(-1) s(-1)) of 12 x 10(-4) M suggests that HbI has a peroxidative capacity for removing H2O2 from solution. Resonance Raman presents the formation of both, met-aquo-HbI and compound II ferryl species in the cyclic reaction of HbI with H2O2. The ferric HbI species is maintained by the presence of H2O2; it can produce HbI compound I, or it can be reduced to a deoxy-HbI derivative to form HbI compound II upon reaction with H2O2. The compound II ferryl vibration frequency appears at 805 and 769 cm(-1) for HbIFe(IV)=(16)O and HbIFe(IV)=(18)O species, respectively. This ferryl mode indicates the absence of hydrogen bonding between the carbonyl group of the distal Q64 and the HbIFe(IV)=O ferryl moiety. The observation suggests that both the trans-ligand effect and the polarizabilty of the HbI heme pocket are responsible for the observed ferryl oxo vibrational energy. The vibrational mode also suggests that the carbonyl group of the distal Q64 is oriented toward the iron of the heme group, increasing the distal pocket electron density.     

Expression and purification of recombinant hemoglobin I from Lucina pectinata

Hemoglobin I (HbI) from Lucina pectinata reacts with hydrogen sulfide to form the ferric sulfide complex needed to transport H₂S to the bacterial endosymbiont. To further study HbI, expression studies of this protein were performed in Escherichia coli. This is the first time that the recombinant HbI was produced using a recombinant DNA expression system. Hemoglobin I cDNA was amplified and cloned into the TOPO-P_BAD expression vector, which contains a fusion tag of six histidine residues (6XHis tag). Plasmid clone sequence analysis was carried out in order to ensure that the insert was in the correct reading frame for proper protein expression in E. coli. The expression of recombinant HbI was optimal when induced for 5 hr with 0.002% of l-arabinose as detected by Western blot analysis. The proto-porphyrin group was inserted into the recombinant HbI. Purification of the heme-bound recombinant protein was performed under native conditions by affinity chromatography using Ni-NTA and Probond resins. The sodium dithionite-reduced recombinant protein presented a shift from the Soret band at 413-435 nm, indicating the presence of the heme group in the adequate amino acid environment of HbI. These results indicate that recombinant HbI from Lucina pectinata can be successfully expressed in a prokaryotic system retaining its activity toward reduction, oxidation, and ligand binding.     

Functional Characterization of the Purified Holo Form of Hemoglobin I from Lucina pectinata Overexpressed in Escherichia coli

The clam Lucina pectinata inhabits the sulfide-rich west coast of the island of Puerto Rico. It contains three different hemoglobins. Hemoglobin I (HbI), which is monomeric at all concentrations, carries H2S in its ferric state. Overexpression of recombinant HbI from Lucina pectinata in BL21STAR(DE3) Escherichia coli cells was performed in the presence of delta-aminolevulinic acid (delta-ALA). Purification of the protein was achieved using FPLC anion exchange and size exclusion chromatography. Functional characterization of the recombinant holo-protein was assessed by detection of the protein heme O2, CO, and H2S derivatives by UV-Vis spectroscopy, with Soret maxima at 416, 421, and 425 nm, respectively. The results indicated that the recombinant HbI binds H2S and forms a heme sulfide complex like the HbI wild-type hemoglobin. Kinetic measurements were performed to determine the H2S affinity of the recombinant HbI. The H2S dissociation and association rate constants were 0.055 x 10(-3)s(-1) and 0.068 x 10(5) M(-1)s(-1), respectively. The H2S affinity constant of the recombinant HbI (0.124 x 10(9) M(-1)) is eightfold lower than that of the native clam HbI reported earlier. This effect is attributed mostly to the first of two missense mutations [Met 61 (E4)-->Val 61 and Ile101 (FG4)-->Val 101] and additional amino acids present in our construct as demonstrated by measurements of the association rate with a new construct lacking most of the additional residues and the missense mutations. The elimination of these residues restores the similarity between the expressed and wild-type hemoglobins, as evidenced by H2S association kinetics. A pH dependence on the H2S association rate was also contributing to the overall affinity constant and was taken into account in the measurements of the functional properties of the new HbI construct.     

Conjugated double bonds in lipid bilayers: a molecular dynamics simulation study

Conjugated linoleic acids (CLA) are found naturally in dairy products. Two isomers of CLA, that differ only in the location of cis and trans double bonds, are found to have distinct and different biological effects. The cis 9 trans 11 (C9T11) isomer is attributed to have the anti-carcinogenic effects, while the trans 10 cis 12 (T10C12) isomer is believed to be responsible for the anti-obesity effects. Since dietary CLA are incorporated into membrane phospholipids, we have used Molecular Dynamics (MD) simulations to investigate the comparative effects of the two isomers on lipid bilayer structure. Specifically, simulations of phosphatidylcholine lipid bilayers in which the sn-2 chains contained one of the two isomers of CLA were performed. Force field parameters for the torsional potential of double bonds were obtained from ab initio calculations. From the MD trajectories we calculated and compared structural properties of the two lipid bilayers, including areas per molecule, density profiles, thickness of bilayers, tilt angle of tail chains, order parameters profiles, radial distribution function (RDF) and lateral pressure profiles. The main differences found between bilayers of the two CLA isomers, are (1) the order parameter profile for C9T11 has a dip in the middle of sn-2 chain while the profile for T10C12 has a deeper dip close to terminal of sn-2 chain, and (2) the lateral pressure profiles show differences between the two isomers. Our simulation results reveal localized physical structural differences between bilayers of the two CLA isomers that may contribute to different biological effects through differential interactions with membrane proteins or cholesterol.     

Aechmea pectinata: a hummingbird-dependent bromeliad with inconspicuous flowers from the rainforest in south-eastern Brazil

The pollination biology of Aechmea pectinata (Bromeliaceae) was studied in a submontane rainforest in south-eastern Brazil. This species has a mainly clumped distribution and its aggregated individuals are likely to be clones. From October to January, during the flowering period, the distal third of its leaves becomes red. The inflorescence produces 1-15 flowers per day over a period of 20-25 d. The flowers are inconspicuous, greenish-white coloured, tubular shaped with a narrow opening, and the stigma is situated just above the anthers. Anthesis begins at 0400 h and flowers last for about 13 h. The highest nectar volume and sugar concentration occur between 0600 and 1000 h, and decrease throughout the day. Aechmea pectinata is self-incompatible and therefore pollinator-dependent. Hummingbirds are its main pollinators (about 90 % of the visits), visiting flowers mainly in the morning. There is a positive correlation between the number of hummingbird visits per inflorescence and the production of nectar, suggesting that the availability of this resource is important in attracting and maintaining visitors. The arrangement of the floral structures favours pollen deposition on the bill of the hummingbirds. Flowers in clumps promote hummingbird territoriality, and a consequence is self-pollination in a broader sense (geitonogamy) as individuals in assemblages are genetically close. However, trap-lining and intruding hummingbirds promote cross-pollination. These observations suggest that successful fruit set of A. pectinata depends on both the spatial distribution of its individuals and the interactions among hummingbirds.     

Influence of point mutations on the flexibility of cytochrome b5: molecular dynamics simulations of holoproteins

Two membrane-bound isoforms of cytochrome b5 have been identified in mammals, one associated with the outer mitochondrial membrane (OM b5) and the other with the endoplasmic reticulum (microsomal, or Mc b5). The soluble heme binding domains of OM and Mc b5 have highly similar three-dimensional structures but differ significantly in physical properties, with OM b5 exhibiting higher stability due to stronger heme association. In this study, we present results of 8.5-ns length molecular dynamics simulations for rat Mc b5, bovine Mc b5, and rat OM b5, as well as for two rat OM b5 mutants that were anticipated to exhibit properties intermediate between those of rat OM b5 and the two Mc proteins: the A18S/I32L/L47R triple mutant (OM3M) and the A18S/I25L/I32L/L47R/L71S quintuple mutant (OM5M). Analysis of the structure, fluctuations, and interactions showed that the five b5 variants used in this study differed in organization of their molecular surfaces and heme binding cores in a way that could be used to explain certain experimentally observed physical differences. Overall, our simulations provided qualitative microscopic explanations of many of the differences in physical properties between OM and Mc b5 and two mutants in terms of localized changes in structure and flexibility. They also reveal that opening of a surface cleft between hydrophobic cores 1 and 2 in bovine Mc b5, observed in two previously reported simulations (E. M. Storch and V. Daggett, Biochemistry, 1995, Vol. 34, pp. 9682-9693; A. Altuve, Biochemistry, 2001, Vol. 40, pp. 9469-9483), probably resulted from removal of crystal contacts and likely does not occur on the nanosecond time scale. Finally, the MD simulations of OM5M b5 verify that stability and dynamic properties of cytochrome b5 are remarkably resistant to mutations that dramatically alter the stability and structure of the apoprotein.     

重组旋毛虫plancitoxin-1-like活性位点突变体蛋白的核酸酶活性

旋毛虫plancitoxin-1-like(Ts-Pt)是旋毛虫125种DNaseⅡ家族蛋白中唯一具有典型DNaseⅡ活性区域HKD基序的核酸酶,且普遍认为,组氨酸位点是DNaseⅡ的活性氨基酸位点。为研究Ts-Pt活性位点突变体蛋白的核酸酶活性,利用重叠PCR方法获得Ts-Pt活性位点突变体片段,以p ET-28a(+)为载体构建重组表达质粒并在大肠杆菌中诱导表达。重组Ts-Pt突变体蛋白经亲和层析纯化后进行SDS-PAGE分析。利用琼脂糖凝胶电泳法和核酸酶酶谱分析重组Ts-Pt突变体蛋白的核酸酶活性。成功构建含Ts-Pt突变体重组质粒的基因工程菌,SDS-PAGE和亲和层析纯化结果显示,重组Ts-Pt突变体蛋白呈包涵体表达。重组蛋白经复性后并没有表现出核酸酶活性,但核酸酶酶谱分析结果显示,包涵体表达的重组Ts-Pt突变体蛋白表现出降解DNA的能力。同时,N端和C端活性位点H及HCK和DHSK突变并不影响Ts-Pt的核酸酶活性,研究结果为进一步研究庞大的DNaseⅡ家族蛋白在旋毛虫发育和感染方面的作用提供一定的参考。     

Impact of disease-causing mutations on inter-domain interactions in cMyBP-C: a steered molecular dynamics study

The molecular interactions of the sarcomeric proteins are essential in the regulation of various cardiac functions. Mutations in the gene MYBPC3 coding for cardiac myosin-binding protein-C (cMyBP-C), a multi-domain protein, are the most common cause of hypertrophic cardiomyopathy (HCM). The N-terminal complex, C1-motif-C2 is a central region in cMyBP-C for the regulation of cardiac muscle contraction. However, the mechanism of binding/unbinding of this complex during health and disease is unknown. Here, we study possible mechanisms of unbinding using steered molecular dynamics simulations for the complex in the wild type, in single mutations (E258K in C1, E441K in C2), as well as in a double mutation (E258K in C1 + E441K in C2), which are associated with severe HCM. The observed molecular events and the calculation of force utilized for the unbinding suggest the following: (i) double mutation can encourage the formation of rigid complex that required large amount of force and long-time to unbind, (ii) C1 appears to start to unbind ahead of C2 regardless of the mutation, and (iii) unbinding of C2 requires larger amount of force than C1. This molecular insight suggests that key HCM-causing mutations might significantly modify the native affinity required for the assembly of the domains in cMyBP-C, which is essential for normal cardiac function.     

Insilico modeling and molecular dynamic simulation of claudin-1 point mutations in HCV infection

Claudin-1 (CLDN1) in association with envelope glycoprotein (CD81) mediates the fusion of HCV into the cytosol. Recent studies have indicated that point mutations in CLDN1 are important for the entry of hepatitis C virus (HCV). To validate these findings, we employed a computational platform to investigate the structural effect of two point mutations (I32M and E48K). Initially, three-dimensional co-ordinates for CLDN1 receptor sequence were generated. Then, three mutant models were built using the point mutation including a double mutant (I32M/E48K) model from the native model structure. Finally, all the four model structures including the native and three mutant models were subjected to molecular dynamics (MD) simulation for a period of 25?ns to appreciate their dynamic behavior. The MD trajectory files were analyzed using cluster and principal component method. The analysis suggested that either of the single mutation has negligible effect on the overall structure of CLDN1 compared to the double mutant form. However, the double mutant model of CLDN1 shows significant negative impact through the impairment of H-bonds and the simultaneous increase in solvent accessible surface area. Our simulation results are visibly consistent with the experimental report suggesting that the CLDN1 receptor distortion is prominent due to the double mutation with large surface accessibility. This increase in accessible surface area due to the coexistence of double mutation may be presumed as one of the key factor that results in permissive action of HCV attachment and infection.     

Role of the amino acid invariants in the active site of MurG as evaluated by site-directed mutagenesis

To evaluate their role in the active site of the MurG enzyme from Escherichia coli, 13 residues conserved in the sequences of 73 MurG orthologues were submitted to site-directed mutagenesis. All these residues lay within, or close to, the active site of MurG as defined by its tridimensional structure [Ha et al., Prot. Sci. 9 (2000) 1045-1052, and Hu et al., Proc. Natl. Acad. Sci. USA 100 (2003) 845-849]. Thirteen mutants proteins, in which residues T15, H18, Y105, H124, E125, N127, N134, S191, N198, R260, E268, Q288 or N291 have been replaced by alanine, were obtained as the C-terminal His-tagged forms. The effects of the mutations on the activity were checked: (i) by functional complementation of an E. coli murG mutant strain by the mutated genes; and (ii) by the determination of the steady-state kinetic parameters of the purified proteins. Most mutations resulted in an important loss of activity and, in the case of N134A, in the production of a highly unstable protein. The results correlated with the assigned or putative functions of the residues based on the tridimensional structure.     

Chondrotinase ABC I thermal stability is enhanced by site-directed mutagenesis: a molecular dynamic simulations approach

Chondroitin sulfate proteoglycans (CSPGs) are potent inhibitors of growth in the adult central nervous system. Use of the enzyme chondroitinase ABC I (ChABC I) as a strategy to reduce CSPG inhibition in experimental models of spinal cord injury has led to observations of its remarkable capacity for repair. More importantly, ChABC therapy has been demonstrated to promote significant recovery of function to spinal injured animals. Despite this incomparable function of ChABC I, its clinical application has been limited because of its thermal instability as reported in the literature. In a recent study by Nazari-Robati et al., thermal stability of ChABC I was improved by protein engineering using site-directed mutagenesis method. Here, in this study, molecular dynamics simulations were used to take a closer look into the phenomenon leading to the experimentally observed thermal stability improvement followed by the corresponding site-directed mutagenesis. We concluded that the mutations induce local flexibility along with a re-conformation into the native structure which consequently increase the protein thermal stability.     

Rapid binding of a cationic active site inhibitor to wild type and mutant mouse acetylcholinesterase: Brownian dynamics simulation including diffusion in the active site gorge

It is known that anionic surface residues play a role in the long-range electrostatic attraction between acetylcholinesterase and cationic ligands. In our current investigation, we show that anionic residues also play an important role in the behavior of the ligand within the active site gorge of acetylcholinesterase. Negatively charged residues near the gorge opening not only attract positively charged ligands from solution to the enzyme, but can also restrict the motion of the ligand once it is inside of the gorge. We use Brownian dynamics techniques to calculate the rate constant k_on for wild type and mutant acetylcholinesterase with a positively charged ligand. These calculations are performed by allowing the ligand to diffuse within the active site gorge. This is an extension of previously reported work in which a ligand was allowed to diffuse only to the enzyme surface. By setting the reaction criteria for the ligand closer to the active site, better agreement with experimental data is obtained. Although a number of residues influence the movement of the ligand within the gorge, Asp74 is shown to play a particularly important role in this function. Asp74 traps the ligand within the gorge, and in this way helps to ensure a reaction. © 1998 John Wiley & Sons, Inc. Biopoly 46: 465–474, 1998     

Redox control in heme proteins: electrostatic substitution in the active site of leghemoglobin

The effects of electrostatic substitutions on the spectroscopic, ligand binding, and redox properties of the heme in leghemoglobin have been examined by replacement of the proximal leucine 88 residue with an aspartic acid residue (Leu88Asp). Electronic and resonance Raman spectra of the ferric derivative of Leu88Asp indicate a mixture of 6-coordinate, high-spin and 6-coordinate, low-spin hemes, analogous to that observed in the recombinant wild-type protein (rLb). At alkaline pH, formation of hydroxide-bound heme is indicated for Leu88Asp; the pK(a) for this transition (8.7 +/- 0.2, micro = 0.10 M, 25.0 degrees C) is 0.4 pH units higher than for rLb. Equilibrium dissociation constants (sodium phosphate, pH 7.0, micro = 0.10 M, 25.0 +/- 0.1 degrees C) for binding of anionic ligands (N(-)(3), nicotinate) to Leu88Asp are higher (K(d,nicotinate) = 6.8 +/- 0.2 microM; K(d,azide) = 33 +/- 0.6 microM) than the corresponding values for rLb (K(d,nicotinate) = 1.4 +/- 0.3 microM (pH 5.5, micro = 0.10 M, 25.0 +/- 0.1 degrees C); K(d,azide) = 4.8 +/- 0.2 microM). Resonance Raman spectra (sodium phosphate, pH 7.0, micro = 0.10 M) for the ferrous derivatives of Leu88Asp and rLb exhibit a strong nu(Fe-His) stretching frequency at 223 cm(-1) in both cases, indicating that the hydrogen bonding structure on the proximal side is not substantially altered in the variant. The reduction potential of Leu88Asp is -14 +/- 2 mV vs standard hydrogen electrode (SHE) (25.0 degrees C, micro = 0.10 M, pH 7.0), a decrease of 35 mV over the corresponding value for the wild-type protein under the same conditions (21 +/- 3 mV vs SHE). An assessment of these data in terms of electrostatic and hydrogen bonding considerations is presented.     

The dominant interaction between peptide and urea is electrostatic in nature: a molecular dynamics simulation study

The conformational equilibrium of a blocked valine peptide in water and aqueous urea solution is studied using molecular dynamics simulations. Pair correlation functions indicate enhanced concentration of urea near the peptide. Stronger hydrogen bonding of urea-peptide compared to water-peptide is observed with preference for helical conformation. The potential of mean force, computed using umbrella sampling, shows only small differences between urea and water solvation that are difficult to quantify. The changes in solvent structure around the peptide are explained by favorable electrostatic interactions (hydrogen bonds) of urea with the peptide backbone. There is no evidence for significant changes in hydrophobic interactions in the two conformations of the peptide in urea solution. Our simulations suggest that urea denatures proteins by preferentially forming hydrogen bonds to the peptide backbone, reducing the barrier for exposing protein residues to the solvent, and reaching the unfolded state.     

Effects of environment on the structure of Pyrococcus furiosus rubredoxin: a molecular dynamics study

Molecular dynamics simulations based on a 0.95-A resolution crystal structure of Pyrococcus furiosus have been performed to elucidate the effects of the environment on the structure of rubredoxin, and proteins in general. Three 1-ns simulations are reported here: two crystalline state simulations at 123 and 300 K, and a solution state simulation at 300 K. These simulations show that temperature has a greater impact on the protein structure than the close molecular contacts of the crystal matrix in rubredoxin, although both have an effect on its dynamic properties. These results indicate that differences between NMR solution structures and X-ray crystal structures will be relatively minor if they are done at similar temperatures. In addition, the crystal simulations appears to mimic previous crystallographic experiments on the effects of cryo-temperature on temperature factors, and might provide a useful tool in the structural analysis of protein structures solved at cryo-temperatures.     

Identification of a single HNH active site in type IIS restriction endonuclease Eco31I

Type IIS restriction endonuclease Eco31I is a "short-distance cutter", which cleaves DNA strands close to its recognition sequence, 5'-GGTCTC(1/5). Previously, it has been proposed that related endonucleases recognizing a common sequence core GTCTC possess two active sites for cleavage of both strands in the DNA substrate. Here, we present bioinformatic identification and experimental evidence for a single nuclease active site. We identified a short region of homology between Eco31I and HNH nucleases, constructed a three-dimensional model of the putative catalytic domain and validated our predictions by random and site-specific mutagenesis. The restriction mechanism of Eco31I is suggested by analogy to the mechanisms of phage T4 endonuclease VII and homing endonuclease I-PpoI. We propose that residues D311 and N334 coordinate the cofactor. H312 acts as a general base-activating water molecule for the nucleophilic attack. K337 together with R340 and D345 are located in close proximity to the active center and are essential for correct folding of catalytic motif, while D345 together with R264 and D273 could be directly involved in DNA binding. We also predict that the Eco31I catalytic domain contains a putative Zn-binding site, which is essential for its structural integrity. Our results suggest that the HNH-like active site is involved in the cleavage of both strands in the DNA substrate. On the other hand, analysis of site-specific mutants in the region, previously suggested to harbor the second active site, revealed its irrelevance to the nuclease activity. Thus, our data argue against the earlier prediction and indicate the presence of a single conserved active site in type IIS restriction endonucleases that recognize common sequence core GTCTC.