Comparative sequence and structure analysis reveal features of cold adaptation of an enzyme in the thermolysin family

Knowledge about the structural features underlying cold adaptation is important for designing enzymes of different industrial relevance. Vibriolysin from Antarctic bacterium strain 643 (VAB) is at present the only enzyme of the thermolysin family from an organism that thrive in extremely cold climate. In this study comparative sequence-structure analysis and molecular dynamics (MD) simulations were used to reveal the molecular features of cold adaptation of VAB. Amino acid sequence analysis of 44 thermolysin enzymes showed that VAB compared to the other enzymes has: (1) fewer arginines, (2) a lower Arg/(Lys + Arg) ratio, (3) a lower fraction of large aliphatic side chains, expressed by the (Ile + Leu)/(Ile + Leu + Val) ratio, (4) more methionines, (5) more serines, and (6) more of the thermolabile amino acid asparagine. A model of the catalytic domain of VAB was constructed based on homology with pseudolysin. MD simulations for 3 ns of VAB, pseudolysin, and thermolysin supported the assumption that cold-adapted enzymes have a more flexible three-dimensional (3D) structure than their thermophilic and mesophilic counterparts, especially in some loop regions. The structural analysis indicated that VAB has fewer intramolecular cation-pi electron interactions and fewer hydrogen bonds than its mesophilic (pseudolysin) and thermophilic (thermolysin) counterparts. Lysine is the dominating cationic amino acids involved in salt bridges in VAB, while arginine is dominating in thermolysin and pseudolysin. VAB has a greater volume of inaccessible cavities than pseudolysin and thermolysin. The electrostatic potentials on the surface of the catalytic domain were also more negative for VAB than for thermolysin and pseudolysin. Thus, the MD simulations, the structural patterns, and the amino acid composition of VAB relative to other enzymes of the thermolysin family suggest that VAB possesses the biophysical properties generally following adaptation to cold climate.     

Potential immunoregulatory role of heme oxygenase-1 in human milk: a combined biochemical and molecular modeling approach

Human milk contains biological factors that are involved in a newborn's growth and immune system regulation. By integrating standard biochemical experimental protocols with computational methods, the present study investigates the presence of heme oxygenase-1 (HO-1), a cytoprotective enzyme, in human milk at different levels of maturation and in milk formulae. Furthermore, we evaluated cytokine and glutathione S-transferase (GSH) levels. Samples were collected from colostrum (on Day 1 after birth), from transition milk (on Postdelivery Days 7 and 14) and from mature milk (on Day 30 after delivery) in 14 healthy women. HO-1 protein, GSH and cytokines levels were measured using enzyme-linked immunosorbent assay and flow cytometry. HO-1 protein levels were significantly higher in colostrum (1.33 ng/ml; 5th centile 0.92; 95th centile 2.38) and in transition milk at 14 days (0.97 ng/ml; 5th centile 0.87; 95th centile 1.45) than in mature milk (0.9 ng/ml; 5th centile 0.8; 95th centile 1.38). Levels of HO-1 in milk formulae were similar to those in colostrum. No significant differences in GSH content were observed in mature milk, transition milk and colostrum, whereas significantly higher GSH levels were observed in milk formulae. No significant levels of cytokines, with the exception of interleukin-8, were found. Computational studies on the possible interactions between HO-1 and CD91 were carried out by a battery of softwares, namely, GRAMM (version 1.03), DALI, CLUSTALW (version 2.0), PatchDock and FireDock, mutually counterchecking and validating each other. The computational results, the strong convergence (to the same “solution”) of which finally leads to an “experimental-like” character, showed that HO-1 may bind to CD91, thus suggesting its major role as a new chaperokine in immune response regulation. These findings, which connect and integrate biochemical data and computational data interpretation, represent a synergistic and powerful means of conducting biological research.     

The protein-protein interactions between SMPI and thermolysin studied by molecular dynamics and MM/PBSA calculations

Thermolysin is a zinc-metalloendopeptidase secreted by the gram-positive thermophilic bacterium Bacillus thermoproteolyticus. Thermolysin belongs to the gluzinicin family of enzymes, which is selectively inhibited by Steptomyces metalloproteinase inhibitor (SMPI). Very little is known about the interaction between SMPI and thermolysin. Knowledge about the protein-protein interactions is very important for designing new thermolysin inhibitors with possible industrial or pharmaceutical applications. In the present study, two binding modes between SMPI and thermolysin were studied by 2300 picoseconds (ps) of comparative molecular dynamics (MD) simulations and calculation of the free energy of binding using the molecular mechanics-Poisson-Boltmann surface area (MM/PBSA) method. One of the positions, the 'horizontal arrow head docking' (HAHD) was similar to the previously proposed binding mode by Tate et al. (Tate, S., Ohno, A., Seeram, S. S., Hiraga, K., Oda, K., and Kainosho, M. J. Mol. Biol. 282, 435-446 (1998)). The other position, the 'vertical arrow head docking' (VAHD) was obtained by a manual docking guided by the shape and charge distribution of SMPI and the binding pocket of thermolysin. The calculations showed that SMPI had stronger interactions with thermolysin in the VAHD than in the HAHD complex, and the VAHD complex was considered more realistic than the HAHD complex. SMPI interacted with thermolysin not only at the active site but had auxiliary binding sites contributing to proper interactions. The VAHD complex can be used for designing small molecule inhibitors mimicking the SMPI-thermolysin binding interfaces.     

Stable Seamless Interfaces and Rapid Ionic Conductivity of Ca–CeO2/LiTFSI/PEO Composite Electrolyte for High‐Rate and High‐Voltage All‐Solid‐State Battery

Stable and seamless interfaces among solid components in all‐solid‐state batteries (ASSBs) are crucial for high ionic conductivity and high rate performance. This can be achieved by the combination of functional inorganic material and flexible polymer solid electrolyte. In this work, a flexible all‐solid‐state composite electrolyte is synthesized based on oxygen‐vacancy‐rich Ca‐doped CeO₂ (Ca–CeO₂) nanotube, lithium bis(trifluoromethanesulfonyl)imide (LiTFSI), and poly(ethylene oxide) (PEO), namely Ca–CeO₂/LiTFSI/PEO. Ca–CeO₂ nanotubes play a key role in enhancing the ionic conductivity and mechanical strength while the PEO offers flexibility and assures the stable seamless contact between the solid electrolyte and the electrodes in ASSBs. The as‐prepared electrolyte exhibits high ionic conductivity of 1.3 × 10^?4 S cm^?1 at 60 °C, a high lithium ion transference number of 0.453, and high‐voltage stability. More importantly, various electrochemical characterizations and density functional theory (DFT) calculations reveal that Ca–CeO₂ helps dissociate LiTFSI, produce free Li ions, and therefore enhance ionic conductivity. The ASSBs based on the as‐prepared Ca–CeO₂/LiTFSI/PEO composite electrolyte deliver high‐rate capability and high‐voltage stability.     

Identification of single nucleotide polymorphisms in the agouti signaling protein (ASIP) gene in some goat breeds in tropical and temperate climates

The agouti-signaling protein (ASIP) plays a major role in mammalian pigmentation as an antagonist to melanocortin-1 receptor gene to stimulate pheomelanin synthesis, a major pigment conferring mammalian coat color. We sequenced a 352 bp fragment of ASIP gene spanning part of exon 2 and part of intron 2 in 215 animals representing six goat breeds from Nigeria and the United States: West African Dwarf, predominantly black; Red Sokoto, mostly red; and Sahel, mostly white from Nigeria; black and white Alpine, brown and white Spanish and white Saanen from the US. Twenty haplotypes from nine mutations representing three intronic, one silent and five missense (p.S19R, p.N35K, p.L36V, p.M42L and p.L45W) mutations were identified in Nigerian goats. Approximately 89 % of Nigerian goats carry haplotype 1 (TGCCATCCG) which seems to be the wild type configuration of mutations in this region of the gene. Although we found no association between these polymorphisms in the ASIP gene and coat color in Nigerian goats, in-silico functional analysis predicts putative deleterious functional impact of the p.L45W mutation on the basic amino-terminal domain of ASIP. In the American goats, two intronic mutations, g.293G>A and g.327C>A, were identified in the Alpine breed, although the g.293G>A mutation is common to American and Nigerian goat populations. All Sannen and Sahel goats in this study belong to haplotypes 1 of both populations which seem to be the wild-type composite ASIP haplotype. Overall, there was no clear association of this portion of the ASIP gene interrogated in this study with coat color variation. Therefore, additional genomic analyses of promoter sequence, the entire coding and non-coding regions of the ASIP gene will be required to obtain a definite conclusion.     

Theoretical calculations of the catalytic triad in short-chain alcohol dehydrogenases/reductases

Three highly conserved active site residues (Ser, Tyr, and Lys) of the family of short-chain alcohol dehydrogenases/reductases (SDRs) were demonstrated to be essential for catalytic activity and have been denoted the catalytic triad of SDRs. In this study computational methods were adopted to study the ionization properties of these amino acids in SDRs from Drosophila melanogaster and Drosophila lebanonensis. Three enzyme models, with different ionization scenarios of the catalytic triad that might be possible when inhibitors bind to the enzyme cofactor complex, were constructed. The binding of the two alcohol competitive inhibitors were studied using automatic docking by the Internal Coordinate Mechanics program, molecular dynamic (MD) simulations with the AMBER program package, calculation of the free energy of ligand binding by the linear interaction energy method, and the hydropathic interactions force field. The calculations indicated that deprotonated Tyr acts as a strong base in the binary enzyme-NAD⁺ complex. Molecular dynamic simulations for 5 ns confirmed that deprotonated Tyr is essential for anchoring and orientating the inhibitors at the active site, which might be a general trend for the family of SDRs. The findings here have implications for the development of therapeutically important SDR inhibitors.