An in silico analysis on the photoproteins Mnemiopsin 1 and Mnemiopsin 2 to explain the experimental results

Mnemiopsin 1 (Mn1) and Mnemiopsin 2 (Mn2) are photoproteins found in Mnemiopsis leidyi. We have tried to answer the question of whether the structural features of photoproteins can explain the observed activity data. According to the activity measurements data, they have the same characteristic wavelength. However, the initial intensity of Mn2 is significantly higher than that of Mn1, and decay time of Mn1 (0.92 s-1 ) is lower than that of Mn2 (1.46 s-1 ). The phylogenetic analysis demonstrates that, compared with Obelin and Aequorin from Obelia longissima and Aequorea victoria, respectively, a gene modification event may have caused the expansion of the N-terminal side of all photoproteins from M. leidyi. An in silico study has shown that the stability of the photoprotein-substrate complex of Mn2 is higher than that of Mn1, indicating a higher affinity of the substrate for Mn2 compared with Mn1. It was revealed that the active EF-hand loops 1 and III in Mn2 is locally more rigid compared with those in Mn1. We concluded that different stability of the photoprotein complexes leads to different initial intensity. While different patterns of the local dynamics of loops I and III may influence the decay rate.

Dynamic Behavior and Ligand Binding Properties of the Wild Type Deoxycytidine Kinase and its Characterized Gemcitabine-Resistant Variant: A Bioinformatics and Computational Study

The main challenge in using gemcitabine, a nucleoside analogue anti-cancer, includes resistance of some patients to this compound due to the genetic mutations on deoxycytidine kinase (dCK) that are found in a fraction of the human population. Here, we investigated the dynamics behavior and ligand binding properties of the wild type (WT) dCK and its characterized double mutant using a combination of bioinformatics tools and molecular dynamics (MD) simulation studies. Root-mean-square deviation (RMSD) values of the WT and mutant enzyme in complex with both ligands (deoxycytidine/gemcitabine) demonstrated that the WT enzyme forms a more stable complex with gemcitabine, as compared with its natural ligand (deoxycytidine). However, the stability of the double mutant-deoxycytidine complex is greater, when compared with gemcitabine. It was also found that Arg131 as a critical residue can affect the binding pattern of the enzyme to ligands in a manner that the WT enzyme can interact with both ligands, while, the mutant enzyme cannot establish efficient interaction with gemcitabine, and shows high affinity to deoxycytidine. These data together indicate that gemcitabine can interact with the WT enzyme in a competitive manner, while double mutant can be considered as a resistant variant against gemcitabine.

Effects of representative point mutations on dynamic behavior of the DISC1-Ndel1 complex: a molecular dynamics study

It has been found that the development of schizophrenia and some other psychiatric disorders is related to defects in the normal functioning of Disrupted-In-Schizophrenia 1 (DISC1). It is a large-sized protein containing 855 residues and acts as an active hub at the core of many interactions with various proteins. On the other hand, NudE Neurodevelopment Protein 1 Like 1 (Ndel1) plays a role in nervous system development via interaction with the DISC1. It was shown that some point mutations on DISC1 have clinical implications. In line with these reports, here we have used the NMR structure of the wild-type (WT) C-terminal tail of DISC1 in complex with the N-terminal fragment of Ndel1, and have constructed the three-dimensional structures of L62Q and L29Q mutants, as the pathologic variants of the complex. The time-dependent interaction of DISC1 with Ndel1 in the WT complex and mutants was simulated by performing molecular dynamics (MD) simulation using programs in the GROMACS package. It was found that the flexibility of residues in some regions of the protein chains increases, and secondary structural changes from ordered toward unordered one leads to destabilizing of the complex in mutants. Destabilization of the complex upon substitution of Leu by Gln was also confirmed by analysis of the contact map plot.Communicated by Ramaswamy H. Sarma.

Bioinformatics and Molecular Dynamics Studies on the Human DISC1 in Complex with the Ndel1

In this study we analyzed the sequence and structure of the human DISC1-Ndel1 complex using bioinformatics tools and molecular dynamics simulation studies. Multiple sequence alignment between the homologue protein sequences in primates shows that corresponding positions of residues in Ndel1 are highly conserved, while the DISC1 has variable conservation lines demonstrating its tolerability against various mutations during evolutionary time scale. In comparison with the mouse variant, structural analysis has shown that the evolutionary inserted charged residues in the human DISC1 (E[Formula: see text]-R[Formula: see text]) can establish intra-chain electrostatic interactions with the K[Formula: see text]-E[Formula: see text] dipeptide that may result in more stability of the DISC1 chain. According to MD simulation studies, the compactness for the human variant of the DISC1-Ndel1 is considerably lower than that of the mouse variant. Analysis of structural fluctuation shows that a fragment at the N-terminus side of the human DISC1 has more residual fluctuation. However, the Ndel1 chain of the human variant has globally more flexibility compared with the mouse variant. Considering the identical amino acid sequence of the Ndel1 chains of human and mouse, it concluded that there is a competition between the inter-chain and intra-chain electrostatic interaction in the human DISC1 that directs the complex to weaker inter-chain interactions with the expense of strengthening the intra-chain stabilizing interaction in complex.

PARS: a bioinformatic tool for Protein Assignment Regarding Secondary structure.. [DOI: 10.21203/rs.3.rs-1660859/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Based on the statistical observation, sequence motifs have been introduced as regular patterns of amino acids that convey a specific structural or functional feature to a group of proteins. But, there is no report about the frequency distribution of the short lengths regularities (such as dipeptides) regarding its structural feature. To find a systematic view about the issue as mentioned above, we developed a bioinformatics tool, namely, PARS (Protein Assignment Regarding Secondary Structure), to enumerate all possible di- and tri-peptides in a given protein structure accompanied by assigning secondary structural elements for the first and second residues. This program can extract structural information of a PDB file, and determines the frequency distribution of 16 and 64 possible secondary structural of possible dipeptides and tripeptides respectively. It is implemented in C# programming language using the dynamic programming method. Using the structural files of a group of proteins, one can construct a detailed structural dataset for further analysis. Other applications include systems biology studies and logical designing of mutations. As an example, a brief description of the published article will be introduced at the end of this paper. This program is freely available for download at: http://www.znu.ac.ir/members/newpage/702/en

Dynamic behavior and ligand binding properties of the wild type dCK and its characterized gemcitabine-resistant variant: a bioinformatics and computational study.. [DOI: 10.21203/rs.3.rs-1661054/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

In order to control the cell proliferation in various types of cancers, gemcitabine as the analogue of deoxycytidine is being used to perturb the DNA synthesis. The main challenge in using gemcitabine includes resistance of some patients to this compound due to the genetic mutations on deoxycytidine kinase (dCK) that are found in a fraction of human population. Here, we investigated the dynamics behavior and ligand binding properties of the wild type (WT) dCK and its characterized double mutant using a combination of the bioinformatics tools and molecular dynamics (MD) simulation studies. The crystal structure of the WT enzyme was trimmed by the MODELLER program and that of the double mutant was made by this program, and the structures were used as input for MD simulation. Root-mean-square deviation (RMSD) values of the WT and mutant enzyme in complex with both ligands, demonstrated that WT enzyme forms more stable complex with gemcitabine, as compared with its natural ligand (deoxycytidine). However, the stability of the double mutant-deoxycytidine complex is greater, when compared with gemcitabine. It was also found that Arg131 as a critical residue can affect the binding pattern of the enzyme to ligands in a manner that the WT enzyme can interact with both ligands, while, mutant enzyme cannot establish efficient interaction with gemcitabine, and shows high affinity to deoxycytidine. These data together indicate that gemcitabine can interact with the WT enzyme in a competitive manner, while double mutant can be considered as a resistant variant against gemcitabine.

Evolutionary adaptation of DHFR via expression of enzyme isoforms with various binding properties and dynamics behavior: a bioinformatics and computational study

We compared the binding properties and dynamics of three experimentally reviewed isoforms of human dihydrofolate reductase (DHFR). The cytoplasmic variants including isoforms1 and 2 (iso1 and iso2) are produced by alternative splicing; while the mitochondrial form is located in the mitochondria. The iso1 as the canonical sequence contains 187 residues, and iso2 differs from the iso1, where it has 1-52 residues missing at the N-terminus of canonical sequence. Here, the structural models of the iso2 and mitochondrial forms were constructed by the MODELLER program using the crystal structure of the iso1 as the template. Bioinformatics analysis on ligand-bearing structures demonstrates that mitochondrial variant forms more stable complex with ligands compared with iso1 and 2, indicating their different binding properties. The root mean square fluctuation (RMSF) data suggest that C-terminus of iso1 contains two representative highly flexible fragments, while iso2 contains a highly flexible fragment at N-terminus end. Interestingly, both ends of mitochondrial variant have a degree of rigidity. Finally, the observation of differences in structural dynamics and binding properties predicts that the simultaneous existence of enzyme isoforms is a way to increase the speed of the enzyme maneuver in response to various environmental conditions. This prediction needs to be tested experimentally.

Thermostability of Ctenophore and Coelenterate Ca2+-Regulated Apo-photoproteins: A Comparative Study

The stabilities of Ca²⁺-regulated ctenophore and coelenterate apo-photoproteins, apo-mnemiopsin (apo-Mne) and apo-aequorin (apo-Aeq), respectively, were compared biochemically, biophysically, and structurally. Despite high degrees of structural and functional conservation, drastic variations in stability and structural dynamics were found between the two proteins. Irreversible thermoinactivation experiments were performed upon incubation of apo-photoproteins at representative temperatures. The inactivation rate constants (k_inact) at 50 °C were determined to be 0.001 and 0.004 min^-1 for apo-Mne and apo-Aeq, respectively. Detailed analysis of the inactivation process suggests that the higher thermostability of apo-Mne is due to the higher activation energy (E_a) and subsequently higher values of ΔH* and ΔG* at a given temperature. According to molecular dynamics simulation studies, the higher hydrogen bond, electrostatic, and van der Waals energies in apo-Mne can validate the relationship between the thermal adaptation of apo-Mne and the energy barrier for the inactivation process. Our results show that favorable residues for protein thermostability such as hydrophobic, charged, and adopted α-helical structure residues are more frequent in the apo-Mne structure. Although the effect of acrylamide on fluorescence quenching suggests that the local flexibility in regions around Trp and Tyr residues of apo-Aeq is higher than that of apo-Mne, which results in it having a better ability to penetrate acrylamide molecules, the root-mean-square fluctuation of helix A in apo-Mne is higher than that in apo-Aeq. It seems that the greater flexibility of apo-Mne in these regions may be considered as a determining factor, affecting the thermal stability of apo-Mne through a balance between structural rigidity and flexibility.

RCOVID19: Recurrence-based SARS-CoV-2 features using chaos game representation

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the COVID-19 pandemic. It was first detected in China and was rapidly spread to other countries. Several thousands of whole genome sequences of SARS-CoV-2 have been reported and it is important to compare them and identify distinctive evolutionary/mutant markers. Utilizing chaos game representation (CGR) as well as recurrence quantification analysis (RQA) as a powerful nonlinear analysis technique, we proposed an effective process to extract several valuable features from genomic sequences of SARS-CoV-2. The represented features enable us to compare genomic sequences with different lengths. The provided dataset involves totally 18 RQA-based features for 4496 instances of SARS-CoV-2.

An evolution-based designing and characterization of mutants of cyclomaltodextrinase: Molecular modeling and spectroscopic studies

Cyclomaltodextrinase (CDase) is a member of the alpha-amylase family GH13, the subfamily GH13_20. In addition to CDase and neopullulanase, this subfamily also contains maltogenic amylase. They have common structural features, but different substrate specificity. In current work, a combination of bioinformatics and experimental tools were used for designing and constructions of single and double mutants of a new variant of CDase from Anoxybacillus flavithermus. Considering the evolutionary variable positions 123 and 127 at the dimer interface of subunits in the alpha-amylase family, these positions in CDase were modified and three mutants, including A123V, C127Q and A123V/C127Q were constructed. The tertiary structure of WT and mutants were made with the MODELLER program, and the phylogenetic tree of homologous protein sequences was built with selected programs in Phylip package. Enzyme kinetic studies revealed that the catalytic efficiency of mutants, especially double one, is lower than the WT enzyme. Heat-induced denaturation experiments were monitored by measuring the UV/Vis signal at 280 nm, and it was found that WT protein is structurally more stable at 25 °C. However, it is more susceptible to changes in temperature compared to the double mutant. It was concluded that the positions 123 and 127 at the dimeric interface of CDase, not only could affect the conformational stability; but also; the catalytic properties of the enzyme by setting up the active site configuration in the dimeric state.

Structural and functional consequences of replacement of His403 with Arg near the catalytic site of Anoxybacillus flavithermus cyclomaltodextrinase

The hydrolytic activity of a thermophilic cyclomaltodextrinase (CMD) from Anoxybacillus flavithermus ZNU-NGA and a representative single mutant were investigated against soluble substrates including α-, β- and γ-cyclomaltodestrines (CDs). Based on the occurrence of arginine (Arg) at position 403 in some homologue proteins, His⁴⁰³ in Wild-type (WT) CMD was replaced with Arg (H403R variant) with site-directed mutagenesis procedures. According to bioinformatics data, Arg⁴⁰³ in mutant protein is located near Glu³⁵⁷ as one of the catalytic residues in a manner that they are able to create a medium-range attractive electrostatistic interaction. Structural studies by Far UV-CD showed that this mutation is accompanied by conversion of a small fraction of α-helix to β-form structure. Fluorescence data reveals that, the hydrophobic regions at the surface of protein, as the binding sites for ANS (8-Anilinonaphthalene-1-sulfonic acid) increase in mutant protein, demonstrating relative inflation of H403R variant compared with WT protein. However, the polarity of microenvironment around chromophores did not change upon mutation. Activity measurement in different ranges of pH and temperatures showed that the optimum values of pH and temperature in mutant enzyme is the same as WT enzyme, however; the activity at optimum points increased in H403R variant. It was also revealed that the H403R variant had slightly improved catalytic efficiency for γ-CD. The same value of activation parameters for both protein variants indicates that mutation does not alter the mechanism of catalysis during enzyme-substrate formation.

Application of Chaotic Laws to Improve Haplotype Assembly Using Chaos Game Representation

Sequence data are deposited in the form of unphased genotypes and it is not possible to directly identify the location of a particular allele on a specific parental chromosome or haplotype. This study employed nonlinear time series modeling approaches to analyze the haplotype sequences obtained from the NGS sequencing method. To evaluate the chaotic behavior of haplotypes, we analyzed their whole sequences, as well as several subsequences from distinct haplotypes, in terms of the SNP distribution on their chromosomes. This analysis utilized chaos game representation (CGR) followed by the application of two different scaling methods. It was found that chaotic behavior clearly exists in most haplotype subsequences. For testing the applicability of the proposed model, the present research determined the alleles in gap positions and positions with low coverage by using chromosome subsequences in which 10% of each subsequence's alleles are replaced by gaps. After conversion of the subsequences' CGR into the coordinate series, a Local Projection (LP) method predicted the measure of ambiguous positions in the coordinate series. It was discovered that the average reconstruction rate for all input data is more than 97%, demonstrating that applying this knowledge can effectively improve the reconstruction rate of given haplotypes.

Designing and construction of novel variants of Chondroitinase ABC I to reduce aggregation rate

Chondroitinase ABC I (cABC I) can degrade inhibitory molecules for axon regrowth at the site of damage after spinal cord injury (SCI). One of the main problems in the practical application is the possibility of structural changes that lead to the inactivation of the enzyme. In current work, three variants of cABC I was designed and constructed by manipulation of a short helix conformation (Gln⁶⁷⁸-Leu⁶⁷⁹-Ser⁶⁸⁰-Gln⁶⁸¹); where Gln residues were converted to Glu. According to the enzyme kinetics studies, the catalytic efficiency of the Q681E and double mutant (Q678E/Q681E) increases in comparison with WT enzyme; while that of Q678E decreases. It was also shown that the rate of the inactivation of the enzyme variants over time is greater in WT and Q678E variants than that of the Q681E and double mutant. Negative values of entropy change of thermal inactivation measurements; demonstrate that inactivation of the WT and Q678E variants are mainly originated from aggregation. These observations can be explained by considering the repulsive electrostatic interaction between enzyme molecules that prevents protein aggregation over time. It is concluded that increasing the solubility of the Q681E and double mutant via favorable interactions of surface-exposed charged residues with dipole momentum of water molecules accompanied by the presence of intermolecular repulsive electrostatic interaction leads to decreasing the rate of aggregation in both long-term storage and heat-induced structural changes.

Molecular mechanisms governing the evolutionary conservation of Glycine in the 6th position of loops ΙΙΙ and ΙV in photoprotein mnemiopsin 2

Photoproteins in their functional form are complexed noncovalently with 2-hydroperoxycoelenterazine. A conformational change upon coordination of Ca⁺² ions with their EF-hand loops leads to oxidation of substrate and emission of light. In all photoproteins, EF-hand loops Ι, ΙΙΙ and ΙV have standard sequence for binding to Ca⁺² ion, however the second one is not able for Ca⁺² coordination. Sequence analysis of Mnemiopsin 2 and other known photoproteins shows that Glutamate (Glu) is occurred in the 6th position of its first EF-hand loop, but this position in other loops of mnemiopsin 2 and all functional loops of other photoproteins is occupied by Glycine (Gly). Here we designed and made single and double mutants where Gly residue at the 6th positions of loops ΙΙΙ and ΙV of mnemiopsin 2 was replaced with Glu. According to the activity measurements, wild-type (WT) and G142E variants have more initial luminescence intensity than G176E and double mutants; while WT and G176E have higher values of half decay time when compared with G142E and double mutants. According to the isothermal denaturation experiments, all protein variants are structurally more stable than WT mnemiopsin 2 and that the stabilizing effects of single mutants are paired resulting in more stability of double mutant against urea denaturation. We concluded that simultaneous occurrence of Gly in the 6th position of loops ΙΙΙ and ΙV is essential for evolutionary adjustment of initial intensity and decay rate of luminescence emission via affecting the interaction of the core structure of photoprotein with coelenteramide and binding affinity of Ca⁺² to the corresponding loops, respectively.

Temperature dependent dynamics in highly homologous adenylate kinases

To correlate the structural features of enzymes to temperature adaptation, we studied psychrophile, mesophile, and thermophile adenylate kinases as model enzymes using bioinformatics and computational tools. Phylogenetic analysis revealed that mesophile and thermophile variants are clustered in one stem of phylogenetic tree and are close to contemporary time, while psychrophile enzyme is more close to their common ancestor. This finding is in good agreement with the process of environmental changes from ice age toward current warm conditions on the earth. We also performed Molecular Dynamics simulation at corresponding temperatures of all enzyme variants including 308, 318, and 328 K. It was found that mesophile enzyme has no distinct deviation of Root Mean Square Deviation (RMSD) and Radius of Gyration (Rg) values from equilibrium states at operating temperature of thermophile enzyme as well as its own optimum temperature. However, psychrophile enzyme undergoes more fluctuations with higher amplitude of change; particularly at 328 K. It was also found that initial increasing of RMSD and Rg for Psychrophile enzyme at all temperatures is occurred gradually; while, the increment of this structural parameters for thermophile enzyme at 328 K is occurred in a highly cooperative and switching manner demonstrating snap structural change of thermophile enzyme in its own temperature. By analysis of Root Mean Square Fluctuation values at different temperatures, we identified two flexible fragments in adenylate kinases so that different dynamic behavior of these regions in mesophile enzyme against operating temperatures of psychrophile and thermophile variants is critical in compensation of flexibility challenges at respective temperatures.

Negative net charge of EF-hand loop I can affect both calcium sensitivity and substrate binding pattern in mnemiopsin 2

Mnemiopsin 2 from Mnemiopsis leidy has three Ca²⁺-binding motifs and has luminescence properties in the presence of calcium and coelenterazine.

Distribution of dipeptides in different protein structural classes: an effort to find new similarities

Finding any regularity in the sequences of proteins and determining their correlation with structural features are of great interest for an understanding of molecular biology. We statistically analyzed the relative frequencies of all 400 possible dipeptides in a data set containing randomly selected proteins of different defined structural classes including all-alpha, all-beta, alpha + beta and alpha/beta families. We found that the distribution of dipeptides is not the same for different structural classes, and some of them are significantly far from a random distribution. A tendency of a given amino acid to localize in the first or second position of a dipeptide depending on the structural class of protein was also found. Interestingly, some amino acids may be substituted for each other in the first or second positions of specific dipeptides in each structural class. This finding apparently contrasts with the routine expectation from the viewpoint of amino acid properties, as classically understood.

Evolutionary conservation of EF-hand ΙΙ loop in aequorin: Priority of intensity to decay rate in bioluminescence emission

As a Ca²⁺-regulated photoprotein, aequorin (Aeq) contains four EF-hand motifs, the second one lacks the standard sequence for Ca²⁺ coordination and doesn't bind to Ca²⁺. Here, we replaced this loop with a functional loop. According to structural studies, although the global stability of modified aequorin (4EFAeq) is higher than that of Aeq; increasing the local flexibility accompanied by internal structural rearrangements in 4EFAeq result in its penetrability to urea and acrylamide. A fast decay rate was observed for 4EFAeq. Assuming the presence of intermediate states in the luminescent reaction, this observation indicate that the loop replacement leads to the lowering of the half-life of intermediate states which results in increasing the rate of conformational switching of 4EFAeq to light emitting form. However, considerable reduction in initial luminescence intensity of 4EFAeq suggests that the number of functional complexes is reduced. Our findings demonstrate that the conformational effects of the second loop in Aeq elicit a delicate balance between local flexibility and global stability which may be considered as an important functional parameter in photoproteins. It was also concluded that evolutionary conservation of EF-hand ΙΙ in the current form is a consequence of priority of intensity to decay rate in bioluminescent organisms.

OCT4 spliced variants are highly expressed in brain cancer tissues and inhibition of OCT4B1 causes G2/M arrest in brain cancer cells

The new claim about the origin of cancer known as Cancer Stem Cell theory states that a somatic differentiated cell can dedifferentiated or reprogrammed for regaining the cancer cell features. It has been recently shown that expression of stemness factors such as Oct4, Sox2, Nanog and Klf4, in a variety of somatic cancers can leads to development of tumorogenesis. Here, the expression of Oct4 variants were evaluated in brain tumor tissues by quantitative RT-PCR and immunohistochemical (IHC) analysis. In next phase of our study, the expression of Oct4B1 was knock-down in brain cancer cell lines and its effect on cell cycle was assessed. Finally, in order to get insights into sequence-structure-function relationships of Oct4 isofroms, their sequences were analysed using bioinformatic tools. Our data revealed that all three variants of Oct4 are expressed in different types of brain cancer. The expression level of Oct4B1, in contast to Oct4B, was much higher in high-grade brain tumors compared with low-grade ones. In line with qPCR, the expression of Oct4A and B isofroms was confirmed with IHC in different types of brain tumors. Moreover, as a result of the suppression of Oct4B1 expression, the brain cancer cells were arrested in G2/M phase of cell cycle. Bioinfromatics data indicated that the predicted Oct4B1 protein have DNA binding properties. All together, our findings suggest that Oct4B1 has a potential role in tumorigenesis of brain cancer and can be considered as a new tumor marker with potential value in diagnosis and treatment of brain cancer.

Determination of structural elements on the folding reaction of mnemiopsin by spectroscopic techniques

Mnemiopsin 1 is a member of photoprotein family, made up of 206 amino acid residues. These Ca(2+)-regulated photoproteins are responsible for light emission in a variety of marine cnidarians and ctenophores. They composed of an apoprotein, a single polypeptide chain of 25 kDa, molecular oxygen and the non-covalently bound chromophore. In this study, we examined whether three mutations, namely R39K, S128G and V183T affect the thermodynamic stability as well as refolding and unfolding kinetics of mnemiopsin 1. Conformational stability measurements using fluorescence and far-UV CD spectroscopies revealed that all variants unfold in multi-step manner in which the secondary and tertiary structures are lost in different steps. However kinetic studies showed that point mutation S128G destabilizes both kinetic intermediate and native conformation; while, these structural elements are stabilized in V183T. We also found that the stability of folded and intermediate states increases in R39K. We concluded that the initial packing of helical segments within the protein structure is more facilitated when Lys with smaller side chain is present in the protein chain.

The effect of charge alteration and flexibility on the function and structural stability of sweet-tasting brazzein

To identify the structural and functional roles of Glu⁹ located at the first loop of sweet-tasting brazzein, two different mutants (E9K and E9G) of the minor form of brazzein were designed and constructed.

Critical Role of Glu175 on Stability and Folding of Bacterial Luciferase: Stopped-flow Fluorescence Study

Bacterial luciferase is a heterodimeric enzyme, which catalyzes the light emission reaction, utilizing reduced FMN (FMNH2), a long chain aliphatic aldehyde and O(2), to produce green-blue light. This enzyme can be readily classed as slow or fast decay based on their rate of luminescence decay in a single turnover. Mutation of Glu175 in alpha subunit to Gly converted slow decay Xenorhabdus Luminescence luciferase to fast decay one. The following studies revealed that changing the luciferase flexibility and lake of Glu-flavin interactions are responsible for the unusual kinetic properties of mutant enzyme. Optical and thermodynamics studies have caused a decrease in free energy and anisotropy of mutant enzyme. Moreover, the role of Glu175 in transition state of folding pathway by use of stopped-flow fluorescence technique has been studied which suggesting that Glu175 is not involved in transition state of folding and appears as surface residue of the nucleus or as a member of one of a few alternative folding nuclei. These results suggest that mutation of Glu175 to Gly extended the structure of Xenorhabdus Luminescence luciferase, locally.