Peptidomics of Three Bothrops Snake Venoms: Insights Into the Molecular Diversification of Proteomes and Peptidomes

Snake venom proteomes/peptidomes are highly complex and maintenance of their integrity within the gland lumen is crucial for the expression of toxin activities. There has been considerable progress in the field of venom proteomics, however, peptidomics does not progress as fast, because of the lack of comprehensive venom sequence databases for analysis of MS data. Therefore, in many cases venom peptides have to be sequenced manually by MS/MS analysis or Edman degradation. This is critical for rare snake species, as is the case of Bothrops cotiara (BC) and B. fonsecai (BF), which are regarded as near threatened with extinction. In this study we conducted a comprehensive analysis of the venom peptidomes of BC, BF, and B. jararaca (BJ) using a combination of solid-phase extraction and reversed-phase HPLC to fractionate the peptides, followed by nano-liquid chromatography-tandem MS (LC-MS/MS) or direct infusion electrospray ionization-(ESI)-MS/MS or MALDI-MS/MS analyses. We detected marked differences in the venom peptidomes and identified peptides ranging from 7 to 39 residues in length by de novo sequencing. Forty-four unique sequences were manually identified, out of which 30 are new peptides, including 17 bradykinin-potentiating peptides, three poly-histidine-poly-glycine peptides and interestingly, 10 L-amino acid oxidase fragments. Some of the new bradykinin-potentiating peptides display significant bradykinin potentiating activity. Automated database search revealed fragments from several toxins in the peptidomes, mainly from L-amino acid oxidase, and allowed the determination of the peptide bond specificity of proteinases and amino acid occurrences for the P4-P4' sites. We also demonstrate that the venom lyophilization/resolubilization process greatly increases the complexity of the peptidome because of the imbalance caused to the venom proteome and the consequent activity of proteinases on venom components. The use of proteinase inhibitors clearly showed different outcomes in the peptidome characterization and suggested that degradomic-peptidomic analysis of snake venoms is highly sensitive to the conditions of sampling procedures. Molecular & Cellular Proteomics 11: 10.1074/mcp.M112.019331, 1245-1262, 2012.

Combined snake venomics and venom gland transcriptomic analysis of Bothropoides pauloensis

Unraveling the repertoire of venom toxins of Bothropoides pauloensis was assessed by snake venomics and venom gland transcriptomic surveys. Both approaches yielded converging overall figures, pointing to metalloproteinases (similar to 37%), PLA(2)s (26-32%), and vasoactive (bradykinin-potentiating) peptides (12-17%) as the major toxin classes. The high occurrence of SVMPs, PLA(2) molecules, vasoactive peptides, along with serine proteinases, explains the local and systemic effects observed in envenomations by B. pauloensis. Minor (<3%) C-type lectin, serine proteinase, L-amino acid oxidase, nerve growth factor, and CRISP molecules were also identified in the transcriptome and the proteome. Low abundance (0.3%) EST singletons coding for vascular endothelial growth factor (svVEGF), ohanin, hyaluronidase, and 5' nucleotidase were found only in the venom gland cDNA library. At the molecular level, the transcriptomic and proteomic datasets display low compositional concordance. In particular, although there is good agreement between transcriptome and proteome in the identity of BPPs, PLA(2) molecules and L-amino acid oxidase, both datasets strongly depart in their C-type lectin and SVMP complements. These data support the view that venom composition is influenced by transcriptional and translational mechanisms and emphasize the value of combining proteomic and transcriptomic approaches to acquire a more complete understanding of the toxinological profile and natural history of the snake venom. (C) 2012 Elsevier B.V. All rights reserved.

Interaction between effects of genes coding for dopamine and glutamate transmission on striatal and parahippocampal function

The genes for the dopamine transporter (DAT) and the D-Amino acid oxidase activator (DAOA or G72) have been independently implicated in the risk for schizophrenia and in bipolar disorder and/or their related intermediate phenotypes. DAT and G72 respectively modulate central dopamine and glutamate transmission, the two systems most robustly implicated in these disorders. Contemporary studies have demonstrated that elevated dopamine function is associated with glutamatergic dysfunction in psychotic disorders. Using functional magnetic resonance imaging we examined whether there was an interaction between the effects of genes that influence dopamine and glutamate transmission (DAT and G72) on regional brain activation during verbal fluency, which is known to be abnormal in psychosis, in 80 healthy volunteers. Significant interactions between the effects of G72 and DAT polymorphisms on activation were evident in the striatum, parahippocampal gyrus, and supramarginal/angular gyri bilaterally, the right insula, in the right pre-/postcentral and the left posterior cingulate/retrosplenial gyri (P < 0.05, FDR-corrected across the whole brain). This provides evidence that interactions between the dopamine and the glutamate system, thought to be altered in psychosis, have an impact in executive processing which can be modulated by common genetic variation.

D-serine, an endogenous synaptic modulator: localization to astrocytes and glutamate-stimulated release.

Using an antibody highly specific for D-serine conjugated to glutaraldehyde, we have localized endogenous D-serine in rat brain. Highest levels of D-serine immunoreactivity occur in the gray matter of the cerebral cortex, hippocampus, anterior olfactory nucleus, olfactory tubercle, and amygdala. Localizations of D-serine immunoreactivity correlate closely with those of D-serine binding to the glycine modulatory site of the N-methyl-D-aspartate (NMDA) receptor as visualized by autoradiography and are inversely correlated to the presence of D-amino acid oxidase. D-Serine is enriched in process-bearing glial cells in neuropil with the morphology of protoplasmic astrocytes. In glial cultures of rat cerebral cortex, D-serine is enriched in type 2 astrocytes. The release of D-serine from these cultures is stimulated by agonists of non-NMDA glutamate receptors, suggesting a mechanism by which astrocyte-derived D-serine could modulate neurotransmission. D-Serine appears to be the endogenous ligand for the glycine site of NMDA receptors.

Polymorphisms in the trace amine receptor 4 (TRAR4) gene on chromosome 6q23.2 are associated with susceptibility to schizophrenia

Several linkage studies across multiple population groups provide convergent support for a susceptibility locus for schizophrenia - and, more recently, for bipolar disorder - on chromosome 6q13-q26. We genotyped 192 European-ancestry and African American (AA) pedigrees with schizophrenia from samples that previously showed linkage evidence to 6q13-q26, focusing on the MOXD1-STX7-TRARs gene cluster at 6q23.2, which contains a number of prime candidate genes for schizophrenia. Thirty-one screening single-nucleotide polymorphisms (SNPs) were selected, providing a minimum coverage of at least 1 SNP/20 kb. The association observed with rs4305745 (P = .0014) within the TRAR4 (trace amine receptor 4) gene remained significant after correction for multiple testing. Evidence for association was proportionally stronger in the smaller AA sample. We performed database searches and sequenced genomic DNA in a 30-proband subsample to obtain a high-density map of 23 SNPs spanning 21.6 kb of this gene. Single-SNP analyses and also haplotype analyses revealed that rs4305745 and/or two other polymorphisms in perfect linkage disequilibrium (LD) with rs4305745 appear to be the most likely variants underlying the association of the TRAR4 region with schizophrenia. Comparative genomic analyses further revealed that rs4305745 and/or the associated polymorphisms in complete LD with rs4305745 could potentially affect gene expression. Moreover, RT-PCR studies of various human tissues, including brain, confirm that TRAR4 is preferentially expressed in those brain regions that have been implicated in the pathophysiology of schizophrenia. These data provide strong preliminary evidence that TRAR4 is a candidate gene for schizophrenia; replication is currently being attempted in additional clinical samples.

Identification and analysis of venom gland-specific genes from the coastal taipan (Oxyuranus scutellatus) and related species

Australian terrestrial elapid snakes contain amongst the most potently toxic venoms known. However, despite the well-documented clinical effects of snake bite, little research has focussed on individual venom components at the molecular level. To further characterise the components of Australian elapid venoms, a complementary (cDNA) microarray was produced from the venom gland of the coastal taipan (Oxyuranus scutellatus) and subsequently screened for venom gland-specific transcripts. A number of putative toxin genes were identified, including neurotoxins, phospholipases, a pseudechetoxin-like gene, a venom natriuretic peptide and a nerve growth factor together with other genes involved in cellular maintenance. Venom gland-specific components also included a calglandulin-like protein implicated in the secretion of toxins from the gland into the venom. These toxin transcripts were subsequently identified in seven other related snake species, producing a detailed comparative analysis at the cDNA and protein levels. This study represents the most detailed description to date of the cloning and characterisation of different genes associated with envenomation from Australian snakes.

Astroglial d-serine is the endogenous co-agonist at the presynaptic NMDA receptor in rat entorhinal cortex

Presynaptic NMDA receptors facilitate the release of glutamate at excitatory cortical synapses and are involved in regulation of synaptic dynamics and plasticity. At synapses in the entorhinal cortex these receptors are tonically activated and provide a positive feedback modulation of the level of background excitation. NMDA receptor activation requires obligatory occupation of a co-agonist binding site, and in the present investigation we have examined whether this site on the presynaptic receptor is activated by endogenous glycine or d-serine. We used whole-cell patch clamp recordings of spontaneous AMPA receptor-mediated synaptic currents from rat entorhinal cortex neurones in vitro as a monitor of presynaptic glutamate release. Addition of exogenous glycine or d-serine had minimal effects on spontaneous release, suggesting that the co-agonist site was endogenously activated and likely to be saturated in our slices. This was supported by the observation that a co-agonist site antagonist reduced the frequency of spontaneous currents. Depletion of endogenous glycine by enzymatic breakdown with a bacterial glycine oxidase had little effect on glutamate release, whereas d-serine depletion with a yeast d-amino acid oxidase significantly reduced glutamate release, suggesting that d-serine is the endogenous agonist. Finally, the effects of d-serine depletion were mimicked by compromising astroglial cell function, and this was rescued by exogenous d-serine, indicating that astroglial cells are the provider of the d-serine that tonically activates the presynaptic NMDA receptor. We discuss the significance of these observations for the aetiology of epilepsy and possible targeting of the presynaptic NMDA receptor in anticonvulsant therapy. © 2014 Elsevier Ltd. All rights reserved.

Caracterização funcional da BaltMTx, uma miotoxina isolada da peçonha de Bothrops alternatus

CHAPTER II: Snake venoms are a complex mixture of organic and inorganic compounds, proteins and peptides such as aminotransferases, acetylcholinesterase, hyaluronidases, L-amino acid oxidase, phospholipase A2, metalloproteases, serine proteases, lectins, disintegrins, and others. Phospholipase A2 directly or indirectly influence the pathophysiological effect on envenomation, as well as their participation in the digestion of the prey. They have several other activities such as hemolytic indirect action, cardiotoxicity, aggregating of platelets, anticoagulant, edema, myotoxic and inflammatory activities. In this work, we describe the functional characterization of BaltMTx, a PLA2 from Bothrops alternatus that inhibits platelet aggregation and present bactericidal effect. The purification of BaltMTx was carried out through three chromatographic steps (ion-exchange on a DEAE-Sephacel column, followed by hydrophobic chromatography on Phenyl–Sepharose and affinity chromatography on HiTrap™ Heparin HP). The protein was purified to homogeneity as judged by its migration profile in SDS–PAGE stained with coomassie blue, and showed a molecular mass of about 15 kDa under reducing conditions and approximately 25 kDa in non-reducing conditions. BaltMTx showed a rather specific inhibitory effect on platelet aggregation induced by epinephrine in human platelet-rich plasma in a dose-dependent manner, whereas it had little or no effect on platelet aggregation induced by collagen or adenosine diphosphate. BaltMTx also showed antibacterial activity against Staphylococcus aureus and Escherichia coli. High concentrations of BatlMTx stimulated the proliferation of Leishmania (Leishmania) infantum and Leishmania (Viania) braziliensis. BaltMTx induced production of inflammatory mediators such as IL-10, IL-12, TNF-α and NO. BaltMTx could be of medical interest as a new tool for the development of novel therapeutic agents for the prevention and treatment of thrombotic disorders as well as bactericidal agent.

Effects of dissolved oxygen availability and culture biomass at induction upon the intracellular expression of monoamine oxidase by recombinant E. coli in fed batch bioprocesses

The effects of oxygen availability and induction culture biomass upon production of an industrially important monoamine oxidase (MAO) were investigated in fed-batch cultures of a recombinant E. coli. For each induction cell biomass 2 different oxygenation methods were used, aeration and oxygen enriched air. Induction at higher biomass levels increased the culture demand for oxygen, leading to fermentative metabolism and accumulation of high levels of acetate in the aerated cultures. Paradoxically, despite an almost eight fold increase in acetate accumulation to levels widely reported to be highly detrimental to protein production, when induction wet cell weight (WCW) rose from 100% to 137.5%, MAO specific activity in these aerated processes showed a 3 fold increase. By contrast, for oxygenated cultures induced at WCW's 100% and 137.5% specific activity levels were broadly similar, but fell rapidly after the maxima were reached. Induction at high biomass levels (WCW 175%) led to very low levels of specific MAO activity relative to induction at lower WCW's in both aerated and oxygenated cultures. Oxygen enrichment of these cultures was a useful strategy for boosting specific growth rates, but did not have positive effects upon specific enzyme activity. Based upon our findings, consideration of the amino acid composition of MAO and previous studies on related enzymes, we propose that this effect is due to oxidative damage to the MAO enzyme itself during these highly aerobic processes. Thus, the optimal process for MAO production is aerated, not oxygenated, and induced at moderate cell density, and clearly represents a compromise between oxygen supply effects on specific growth rate/induction cell density, acetate accumulation, and high specific MAO activity. This work shows that the negative effects of oxygen previously reported in free enzyme preparations, are not limited to these acellular environments but are also discernible in the sheltered environment of the cytosol of E. coli cells.

Preparative deracemization of unnatural amino acids

Unnatural amino acids are a growing class of intermediates required for pharmaceuticals, agrochemicals and other industrial products. However, no single method has proven sufficiently versatile to prepare these compounds broadly at scale. To address this need, we have developed a general chemoenzymatic process to prepare enantiomerically pure L- and D-amino acids in high yield by deracemization of racemic starting materials. This method involves the concerted action of an enantioselective oxidase biocatalyst and a non-selective chemical reducing agent to effect the stereoinversion of one enantiomer and can result in an enantiomeric excess of >99% from the starting racemate, and product yields of over 90%. This approach compares very favourably with resolution processes, which have a maximum single-pass yield of 50%. We have developed efficient methods to adapt the process towards new target compounds and to optimize key factors that influence process efficiency and offer competitive economics at scale.

Studies on crystalline complexes involving amino acids. V. The structure of L-serine-L-ascorbic acid

L-Serine-L-ascorbic acid, C3HTNOa. C6HsO6, a 1:1 complex between the amino acid serine and the vitamin ascorbic acid, crystallizes in the orthorhombic space group P2~2~2~ with four formula units in a cell of dimensions a = 5.335(3), b = 8.769(2), c = 25.782 (5) A. The structure was solved by direct methods and refined by full-matrix least squares to an R of 0.036 for 951 observed reflections. Both molecules are neutral in the structure. The conformation of the serine molecule is different from that observed in the crystal structures of L-serine, DL-serine and L-serine monohydrate. The enediol group in the ascorbic acid molecule is planar, whereas significant departures from planarity are observed in the lactone group. The conformation of this molecule is similar to that observed in arginine ascorbate. The unlike molecules aggregate into separate columns in the crystal structure. The columns are held together by hydrogen bonds. Among these, a pair of hydrogen bonds between the enediol group of ascorbic acid and the carboxylate group of serine provides a possible model for a specific interaction between ascorbic acid and a carboxylate ion.

Redox-linked proton transfer by cytochrome c oxidase

The respiratory chain is found in the inner mitochondrial membrane of higher organisms and in the plasma membrane of many bacteria. It consists of several membrane-spanning enzymes, which conserve the energy that is liberated from the degradation of food molecules as an electrochemical proton gradient across the membrane. The proton gradient can later be utilized by the cell for different energy requiring processes, e.g. ATP production, cellular motion or active transport of ions. The difference in proton concentration between the two sides of the membrane is a result of the translocation of protons by the enzymes of the respiratory chain, from the negatively charged (N-side) to the positively charged side (P-side) of the lipid bilayer, against the proton concentration gradient. The endergonic proton transfer is driven by the flow of electrons through the enzymes of the respiratory chain, from low redox-potential electron donors to acceptors of higher potential, and ultimately to oxygen. Cytochrome c oxidase is the last enzyme in the respiratory chain and catalyzes the reduction of dioxygen to water. The redox reaction is coupled to proton transport across the membrane by a yet unresolved mechanism. Cytochrome c oxidase has two proton-conducting pathways through which protons are taken up to the interior part of the enzyme from the N-side of the membrane. The K-pathway transfers merely substrate protons, which are consumed in the process of water formation at the catalytic site. The D-pathway transfers both substrate protons and protons that are pumped to the P-side of the membrane. This thesis focuses on the role of two conserved amino acids in proton translocation by cytochrome c oxidase, glutamate 278 and tryptophan 164. Glu278 is located at the end of the D-pathway and is thought to constitute the branching point for substrate and pumped protons. In this work, it was shown that although Glu278 has an important role in the proton transfer mechanism, its presence is not an obligatory requirement. Alternative structural solutions in the area around Glu278, much like the ones present in some distantly related heme-copper oxidases, could in the absence of Glu278 support the formation of a long hydrogen-bonded water chain through which proton transfer from the D-pathway to the catalytic site is possible. The other studied amino acid, Trp164, is hydrogen bonded to the ∆-propionate of heme a3 of the catalytic site. Mutation of this amino acid showed that it may be involved in regulation of proton access to a proton acceptor, a pump site, from which the proton later is expelled to the P-side of the membrane. The ion pair that is formed by the ∆-propionate of heme a3 and arginine 473 is likely to form a gate-like structure, which regulates proton mobility to the P-side of the membrane. The same gate may also be part of an exit path through which water molecules produced at the catalytically active site are removed towards the external side of the membrane. Time-resolved optical and electrometrical experiments with the Trp164 to phenylalanine mutant revealed a so far undetected step in the proton pumping mechanism. During the A to PR transition of the catalytic cycle, a proton is transferred from Glu278 to the pump site, located somewhere in the vicinity of the ∆-propionate of heme a3. A mechanism for proton pumping by cytochrome c oxidase is proposed on the basis of the presented results and the mechanism is discussed in relation to some relevant experimental data. A common proton pumping mechanism for all members of the heme-copper oxidase family is moreover considered.

Characterisation, cloning and production of industrially interesting enzymes: gluconolactone oxidase of Penicillium cyaneo-fulvum and gluconate 5-dehydrogenase of Gluconobacter suboxydans

The work covered in this thesis is focused on the development of technology for bioconversion of glucose into D-erythorbic acid (D-EA) and 5-ketogluconic acid (5-KGA). The task was to show on proof-of-concept level the functionality of the enzymatic conversion or one-step bioconversion of glucose to these acids. The feasibility of both studies to be further developed for production processes was also evaluated. The glucose - D-EA bioconversion study was based on the use of a cloned gene encoding a D-EA forming soluble flavoprotein, D-gluconolactone oxidase (GLO). GLO was purified from Penicillium cyaneo-fulvum and partially sequenced. The peptide sequences obtained were used to isolate a cDNA clone encoding the enzyme. The cloned gene (GenBank accession no. AY576053) is homologous to the other known eukaryotic lactone oxidases and also to some putative prokaryotic lactone oxidases. Analysis of the deduced protein sequence of GLO indicated the presence of a typical secretion signal sequence at the N-terminus of the enzyme. No other targeting/anchoring signals were found, suggesting that GLO is the first known lactone oxidase that is secreted rather than targeted to the membranes of the endoplasmic reticulum or mitochondria. Experimental evidence supports this analysis, as near complete secretion of GLO was observed in two different yeast expression systems. Highest expression levels of GLO were obtained using Pichia pastoris as an expression host. Recombinant GLO was characterised and the suitability of purified GLO for the production of D-EA was studied. Immobilised GLO was found to be rapidly inactivated during D-EA production. The feasibility of in vivo glucose - D-EA conversion using a P. pastoris strain co-expressing the genes of GLO and glucose oxidase (GOD, E.C. 1.1.3.4) of A. niger was demonstrated. The glucose - 5-KGA bioconversion study followed a similar strategy to that used in the D-EA production research. The rationale was based on the use of a cloned gene encoding a membrane-bound pyrroloquinoline quinone (PQQ)-dependent gluconate 5-dehydrogenase (GA 5-DH). GA 5-DH was purified to homogeneity from the only source of this enzyme known in literature, Gluconobacter suboxydans, and partially sequenced. Using the amino acid sequence information, the GA 5-DH gene was cloned from a genomic library of G. suboxydans. The cloned gene was sequenced (GenBank accession no. AJ577472) and found to be an operon of two adjacent genes encoding two subunits of GA 5-DH. It turned out that GA 5-DH is a rather close homologue of a sorbitol dehydrogenase from another G. suboxydans strain. It was also found that GA 5-DH has significant polyol dehydrogenase activity. The G. suboxydans GA 5-DH gene was poorly expressed in E. coli. Under optimised conditions maximum expression levels of GA 5-DH did not exceed the levels found in wild-type G. suboxydans. Attempts to increase expression levels resulted in repression of growth and extensive cell lysis. However, the expression levels were sufficient to demonstrate the possibility of bioconversion of glucose and gluconate into 5-KGA using recombinant strains of E. coli. An uncharacterised homologue of GA 5-DH was identified in Xanthomonas campestris using in silico screening. This enzyme encoded by chromosomal locus NP_636946 was found by a sequencing project of X. campestris and named as a hypothetical glucose dehydrogenase. The gene encoding this uncharacterised enzyme was cloned, expressed in E. coli and found to encode a gluconate/polyol dehydrogenase without glucose dehydrogenase activity. Moreover, the X. campestris GA 5-DH gene was expressed in E. coli at nearly 30 times higher levels than the G. suboxydans GA 5-DH gene. Good expressability of the X. campestris GA-5DH gene makes it a valuable tool not only for 5-KGA production in the tartaric acid (TA) bioprocess, but possibly also for other bioprocesses (e.g. oxidation of sorbitol into L-sorbose). In addition to glucose - 5-KGA bioconversion, a preliminary study of the feasibility of enzymatic conversion of 5-KGA into TA was carried out. Here, the efficacy of the first step of a prospective two-step conversion route including a transketolase and a dehydrogenase was confirmed. It was found that transketolase convert 5-KGA into TA semialdehyde. A candidate for the second step was suggested to be succinic dehydrogenase, but this was not tested. The analysis of the two subprojects indicated that bioconversion of glucose to TA using X. campestris GA 5-DH should be prioritised first and the process development efforts in future should be focused on development of more efficient GA 5-DH production strains by screening a more suitable production host and by protein engineering.

Oxygen reduction and proton translocation by cytochrome c oxidase

Energy conversion by living organisms is central dogma of bioenergetics. The effectiveness of the energy extraction by aerobic organisms is much greater than by anaerobic ones. In aerobic organisms the final stage of energy conversion occurs in respiratory chain that is located in the inner membrane of mitochondria or cell membrane of some aerobic bacteria. The terminal complex of the respiratory chain is cytochrome c oxidase (CcO) - the subject of this study. The primary function of CcO is to reduce oxygen to water. For this, CcO accepts electrons from a small soluble enzyme cytochrome c from one side of the membrane and protons from another side. Moreover, CcO translocates protons across the membrane. Both oxygen reduction and proton translocation contributes to generation of transmembrane electrochemical gradient that is used for ATP synthesis and different types of work in the cell. Although the structure of CcO is defined with a relatively high atomic resolution (1.8 Å), its function can hardly be elucidated from the structure. The electron transfer route within CcO and its steps are very well defined. Meanwhile, the proton transfer roots were predicted from the site-specific mutagenesis and later proved by X-ray crystallography, however, the more strong proof of the players of the proton translocation machine is still required. In this work we developed new methods to study CcO function based on FTIR (Fourier Transform Infrared) spectroscopy. Mainly with use of these methods we answered several questions that were controversial for many years: [i] the donor of H+ for dioxygen bond splitting was identified and [ii] the protolytic transitions of Glu-278 one of the key amino acid in proton translocation mechanism was shown for the first time.

Molecular effects of encapsulation of glucose oxidase dimer by graphene

Knowing the nature of the enzyme-graphene interface is critical for a design of graphene-based biosensors. Extensive contacts between graphene and enzyme could be obtained by employing a suitable encapsulation which does not impede its enzymatic reaction. We have performed molecular dynamics simulations to obtain an insight on many forms of contact between glucose oxidase dimer and the single-layer graphene nano-sheets. The unconnected graphene sheets tended to form a flat stack regardless of their initial positions around the enzyme, whereas the same graphene sheets linked together formed a flower-like shape engendering different forms of wrapping of the enzyme. During the encapsulation no core hydrophobic residues of the enzyme were exposed. Since the polar and charged amino acids populated the enzyme's surface we also estimated, using DFT calculations, the interaction energies of individual polar and charged amino acid residues with graphene. It was found that the negatively charged residues can bind to graphene unexpectedly strongly; however, the main effect of encapsulation comes from the overlap of adjacent edges of graphene sheets.