Flavin Adenine Dinucleotide Binding Site Research Articles

Magnetic orientation in juvenile Atlantic herring (Clupea harengus) could involve cryptochrome 4 as a potential magnetoreceptor.

The Earth's magnetic field can provide reliable directional information, allowing migrating animals to orient themselves using a magnetic compass or estimate their position relative to a target using map-based orientation. Here we show for the first time that young, inexperienced herring (Clupea harengus, Ch) have a magnetic compass when they migrate hundreds of kilometres to their feeding grounds. In birds, such as the European robin (Erithacus rubecula), radical pair-based magnetoreception involving cryptochrome 4 (ErCRY4) was demonstrated; the molecular basis of magnetoreception in fish is still elusive. We show that cry4 expression in the eye of herring is upregulated during the migratory season, but not before, indicating a possible use for migration. The amino acid structure of herring ChCRY4 shows four tryptophans and a flavin adenine dinucleotide-binding site, a prerequisite for a magnetic receptor. Using homology modelling, we successfully reconstructed ChCRY4 of herring, DrCRY4 of zebrafish (Danio rerio) and StCRY4 of brown trout (Salmo trutta) and showed that ChCRY4, DrCRY4 and ErCRY4a, but not StCRY4, exhibit very comparable dynamic behaviour. The electron transfer could take place in ChCRY4 in a similar way to ErCRY4a. The combined behavioural, transcriptomic and simulation experiments provide evidence that CRY4 could act as a magnetoreceptor in Atlantic herring.

The protein family of pyruvate:quinone oxidoreductases: Amino acid sequence conservation and taxonomic distribution

Pyruvate:quinone oxidoreductases (PQOs) catalyse the oxidative decarboxylation of pyruvate to acetate and concomitant reduction of quinone to quinol with the release of CO2. They are thiamine pyrophosphate (TPP) and flavin-adenine dinucleotide (FAD) containing enzymes, which interact with the membrane in a monotopic way. PQOs are considered as part of alternatives to most recognized pyruvate catabolizing pathways, and little is known about their taxonomic distribution and structural/functional relationship.In this bioinformatics work we tackled these gaps in PQO knowledge. We used the KEGG database to identify PQO coding genes, performed a multiple sequence analysis which allowed us to study the amino acid conservation on these enzymes, and looked at their possible cellular function. We observed that PQOS are enzymes exclusively present in prokaryotes with most of the sequences identified in bacteria. Regarding the amino acid sequence conservation, we found that 75 amino acid residues (out of 570, on average) have a conservation over 90 %, and that the most conserved regions in the protein are observed around the TPP and FAD binding sites. We systematized the presence of conserved features involved in Mg2+, TPP and FAD binding, as well as residues directly linked to the catalytic mechanism. We also established the presence of a new motif named “HEH lock”, possibly involved in the dimerization process. The results here obtained for the PQO protein family contribute to a better understanding of the biochemistry of these respiratory enzymes.

FAD-BERT: Improved prediction of FAD binding sites using pre-training of deep bidirectional transformers

The electron transport chain is a series of protein complexes embedded in the process of cellular respiration, which is an important process to transfer electrons and other macromolecules throughout the cell. Identifying Flavin Adenine Dinucleotide (FAD) binding sites in the electron transport chain is vital since it helps biological researchers precisely understand how electrons are produced and are transported in cells. This study distills and analyzes the contextualized word embedding from pre-trained BERT models to explore similarities in natural language and protein sequences. Thereby, we propose a new approach based on Pre-training of Bidirectional Encoder Representations from Transformers (BERT), Position-specific Scoring Matrix profiles (PSSM), Amino Acid Index database (AAIndex) to predict FAD-binding sites from the transport proteins which are found in nature recently. Our proposed approach archives 85.14% accuracy and improves accuracy by 11%, with Matthew's correlation coefficient of 0.39 compared to the previous method on the same independent set. We also deploy a web server that identifies FAD-binding sites in electron transporters available for academics at http://140.138.155.216/fadbert/.

A topologically distinct class of photolyases specific for UV lesions within single-stranded DNA.

Photolyases are ubiquitously occurring flavoproteins for catalyzing photo repair of UV-induced DNA damages. All photolyases described so far have a bilobal architecture with a C-terminal domain comprising flavin adenine dinucleotide (FAD) as catalytic cofactor and an N-terminal domain capable of harboring an additional antenna chromophore. Using sequence-similarity network analysis we discovered a novel subgroup of the photolyase/cryptochrome superfamily (PCSf), the NewPHLs. NewPHL occur in bacteria and have an inverted topology with an N-terminal catalytic domain and a C-terminal domain for sealing the FAD binding site from solvent access. By characterizing two NewPHL we show a photochemistry characteristic of other PCSf members as well as light-dependent repair of CPD lesions. Given their common specificity towards single-stranded DNA many bacterial species use NewPHL as a substitute for DASH-type photolyases. Given their simplified architecture and function we suggest that NewPHL are close to the evolutionary origin of the PCSf.

Exploration of 5-cyano-6-phenylpyrimidin derivatives containing an 1,2,3-triazole moiety as potent FAD-based LSD1 inhibitors

Histone lysine specific demethylase 1 (LSD1) has become a potential therapeutic target for the treatment of cancer. Discovery and develop novel and potent LSD1 inhibitors is a challenge, although several of them have already entered into clinical trials. Herein, for the first time, we reported the discovery of a series of 5-cyano-6-phenylpyrimidine derivatives as LSD1 inhibitors using flavin adenine dinucleotide (FAD) similarity-based designing strategy, of which compound 14q was finally identified to repress LSD1 with IC50 = 183 nmol/L. Docking analysis suggested that compound 14q fitted well into the FAD-binding pocket. Further mechanism studies showed that compound 14q may inhibit LSD1 activity competitively by occupying the FAD binding sites of LSD1 and inhibit cell migration and invasion by reversing epithelial to mesenchymal transition (EMT). Overall, these findings showed that compound 14q is a suitable candidate for further development of novel FAD similarity-based LSD1 inhibitors.

Crystal Structure of the Active Site Mutant Form of Soluble Fumarate Reductase, Osm1

Soluble fumarate reductase is essential for survival under anaerobic conditions. This enzyme can maintain the redox balance in the cell by catalyzing the reduction of fumarate to succinate. Although the overall reaction mechanism of soluble fumarate reductase in yeast, Osm1, has been proposed by a previous structural study, the details of the underlying mechanism are not completely elucidated. The present study provides the structural information regarding the active site mutant form of Osm1 (R326A), thus, revealing that R326A mutation does not affect the substrate binding. Structural alterations of the residues surrounding the active site, and the missing 2nd flavin adenine dinucleotide (FAD) in the previously defined 2nd FAD binding site, were observed as characteristic features of the Osm1 R326A crystal structure. Based on these findings, we provided a clue that can explain the loss of activity of Osm1 R326A.

Prediction of FAD binding sites in electron transport proteins according to efficient radial basis function networks and significant amino acid pairs.

BackgroundCellular respiration is a catabolic pathway for producing adenosine triphosphate (ATP) and is the most efficient process through which cells harvest energy from consumed food. When cells undergo cellular respiration, they require a pathway to keep and transfer electrons (i.e., the electron transport chain). Due to oxidation-reduction reactions, the electron transport chain produces a transmembrane proton electrochemical gradient. In case protons flow back through this membrane, this mechanical energy is converted into chemical energy by ATP synthase. The convert process is involved in producing ATP which provides energy in a lot of cellular processes. In the electron transport chain process, flavin adenine dinucleotide (FAD) is one of the most vital molecules for carrying and transferring electrons. Therefore, predicting FAD binding sites in the electron transport chain is vital for helping biologists understand the electron transport chain process and energy production in cells.ResultsWe used an independent data set to evaluate the performance of the proposed method, which had an accuracy of 69.84 %. We compared the performance of the proposed method in analyzing two newly discovered electron transport protein sequences with that of the general FAD binding predictor presented by Mishra and Raghava and determined that the accuracy of the proposed method improved by 9–45 % and its Matthew’s correlation coefficient was 0.14–0.5. Furthermore, the proposed method enabled reducing the number of false positives significantly and can provide useful information for biologists.ConclusionsWe developed a method that is based on PSSM profiles and SAAPs for identifying FAD binding sites in newly discovered electron transport protein sequences. This approach achieved a significant improvement after we added SAAPs to PSSM features to analyze FAD binding proteins in the electron transport chain. The proposed method can serve as an effective tool for predicting FAD binding sites in electron transport proteins and can help biologists understand the functions of the electron transport chain, particularly those of FAD binding sites. We also developed a web server which identifies FAD binding sites in electron transporters available for academics.

Association of C677T transition of the human methylenetetrahydrofolate reductase (MTHFR) gene with male infertility.

The human methylenetetrahydrofolate reductase (MTHFR) gene encodes one of the key enzymes in folate metabolism. This gene is located on chromosome 1 (1p36.3), which has 12 exons. The aim of the present study was to investigate the possible association of the two (C677T and A1298C) polymorphisms of this gene with male infertility. In a case-control study, 250 blood samples were collected from IVF centres in Sari and Babol (Iran): 118 samples were from oligospermic men and 132 were from controls. Two single nucleotide polymorphisms of the MTHFR genotype were detected using polymerase chain reaction-restriction fragment length polymorphism. There was no association found between the A1298C variant and male infertility. However, carriers of the 677T allele (CT and TT genotypes) were at a higher risk of infertility than individuals with other genotypes (odds ratio 1.84; 95% confidence interval 1.11-3.04; P=0.0174). Structural analysis of human MTHFR flavoprotein showed that C677T transition played an important role in the change in affinity of the MTHFR-Flavin adenine dinucleotide binding site. Based on our results, we suggest that C677T transition in MTHFR may increase the risk of male infertility, and detection of the C677T polymorphism biomarker may be helpful in the screening of idiopathic male infertility.

Is the recovery of (photo) respiratory CO2 and intermediates minimal?

Is the recovery of (photo) respiratory <scp>CO</scp>₂ and intermediates minimal?

Computational analysis of a novel mutation in ETFDH gene highlights its long-range effects on the FAD-binding motif

BackgroundMultiple acyl-coenzyme A dehydrogenase deficiency (MADD) is an autosomal recessive disease caused by the defects in the mitochondrial electron transfer system and the metabolism of fatty acids. Recently, mutations in electron transfer flavoprotein dehydrogenase (ETFDH) gene, encoding electron transfer flavoprotein:ubiquinone oxidoreductase (ETF:QO) have been reported to be the major causes of riboflavin-responsive MADD. To date, no studies have been performed to explore the functional impact of these mutations or their mechanism of disrupting enzyme activity.ResultsHigh resolution melting (HRM) analysis and sequencing of the entire ETFDH gene revealed a novel mutation (p.Phe128Ser) and the hotspot mutation (p.Ala84Thr) from a patient with MADD. According to the predicted 3D structure of ETF:QO, the two mutations are located within the flavin adenine dinucleotide (FAD) binding domain; however, the two residues do not have direct interactions with the FAD ligand. Using molecular dynamics (MD) simulations and normal mode analysis (NMA), we found that the p.Ala84Thr and p.Phe128Ser mutations are most likely to alter the protein structure near the FAD binding site as well as disrupt the stability of the FAD binding required for the activation of ETF:QO. Intriguingly, NMA revealed that several reported disease-causing mutations in the ETF:QO protein show highly correlated motions with the FAD-binding site.ConclusionsBased on the present findings, we conclude that the changes made to the amino acids in ETF:QO are likely to influence the FAD-binding stability.

Interaction between the reductase Tah18 and highly conserved Fe‐S containing Dre2 C‐terminus is essential for yeast viability

Tah18-Dre2 is a recently identified yeast protein complex, which is highly conserved in human and has been implicated in the regulation of oxidative stress induced cell death and in cytosolic Fe-S proteins synthesis. Tah18 is a diflavin oxido-reductase with binding sites for flavin mononucleotide, flavin adenine dinucleotide and nicotinamide adenine dinucleotide phosphate, which is able to transfer electrons to Dre2 Fe-S clusters. In this work we characterized in details the interaction between Tah18 and Dre2, and analysed how it conditions yeast viability. We show that Dre2 C-terminus interacts in vivo and in vitro with the flavin mononucleotide- and flavin adenine dinucleotide-binding sites of Tah18. Neither the absence of the electron donor nicotinamide adenine dinucleotide phosphate-binding domain in purified Tah18 nor the absence of Fe-S in aerobically purified Dre2 prevents the binding in vitro. In vivo, when this interaction is affected in a dre2 mutant, yeast viability is reduced. Conversely, enhancing artificially the interaction between mutated Dre2 and Tah18 restores cellular viability despite still reduced cytosolic Fe-S cluster biosynthesis. We conclude that Tah18-Dre2 interaction in vivo is essential for yeast viability. Our study may provide new insight into the survival/death switch involving this complex in yeast and in human cells.

In vitro evolution of styrene monooxygenase from Pseudomonas putida CA-3 for improved epoxide synthesis

The styAB genes from Pseudomonas putida CA-3, which encode styrene monooxygenase, were subjected to three rounds of in vitro evolution using error-prone polymerase chain reaction with a view to improving the rate of styrene oxide and indene oxide formation. Improvements in styrene monooxygenase activity were monitored using an indole to indigo conversion assay. Each round of random mutagenesis generated variants improved in indigo formation with third round variants improved nine- to 12-fold over the wild type enzyme. Each round of in vitro evolution resulted in two to three amino acid substitutions in styrene monooxygenase. While the majority of mutations occurred in styA (oxygenase), mutations were also observed in styB (reductase). A mutation resulting in the substitution of valine with isoleucine at amino acid residue 303 occurred near the styrene and flavin adenine dinucleotide binding site of styrene monooxygenase. One mutation caused a shift in the reading frame in styA and resulted in a StyA variant that is 19 amino acids longer than the wild-type protein. Whole cells expressing the best styrene monooxygenase variants (round 3) exhibited eight- and 12-fold improvements in styrene and indene oxidation rates compared to the wild-type enzyme. In all cases, a single enantiomer, (S)-styrene oxide, was formed from styrene while (1S,2R)-indene oxide was the predominant enantiomer (e.e. 97%) formed from indene. The average yield of styrene oxide and indene oxide from their respective alkene substrates was 65% and 90%, respectively.

Functional and antigenic properties of GlpO from Mycoplasma mycoides subsp. mycoides SC: characterization of a flavin adenine dinucleotide-binding site deletion mutant.

L-α-glycerophosphate oxidase (GlpO) plays a central role in virulence of Mycoplasma mycoides subsp. mycoides SC, a severe bacterial pathogen causing contagious bovine pleuropneumonia (CBPP). It is involved in production and translocation of toxic H2O2 into the host cell, causing inflammation and cell death. The binding site on GlpO for the cofactor flavin adenine dinucleotide (FAD) has been identified as Gly 12−Gly13−Gly 14−Ile15−Ile16−Gly 17. Recombinant GlpO lacking these six amino acids (GlpOΔFAD) was unable to bind FAD and was also devoid of glycerophosphate oxidase activity, in contrast to non-modified recombinant GlpO that binds FAD and is enzymatically active. Polyclonal monospecific antibodies directed against GlpOΔFAD, similarly to anti-GlpO antibodies, neutralised H2O2 production of M. mycoides subsp. mycoides SC grown in the presence of glycerol, as well as cytotoxicity towards embryonic calf nasal epithelial (ECaNEp) cells. The FAD-binding site of GlpO is therefore suggested as a valuable target site for the future construction of deletion mutants to yield attenuated live vaccines of M. mycoides subsp. mycoides SC necessary to efficiently combat CBPP.

Structure and DNA binding of alkylation response protein AidB

Exposure of Escherichia coli to alkylating agents activates expression of AidB in addition to DNA repair proteins Ada, AlkA, and AlkB. AidB was recently shown to possess a flavin adenine dinucleotide (FAD) cofactor and to bind to dsDNA, implicating it as a flavin-dependent DNA repair enzyme. However, the molecular mechanism by which AidB acts to reduce the mutagenic effects of specific DNA alkylators is unknown. We present a 1.7-A crystal structure of AidB, which bears superficial resemblance to the acyl-CoA dehydrogenase superfamily of flavoproteins. The structure reveals a unique quaternary organization and a distinctive FAD active site that provides a rationale for AidB's limited dehydrogenase activity. A highly electropositive C-terminal domain not present in structural homologs was identified by mutational analysis as the DNA binding site. Structural analysis of the DNA and FAD binding sites provides evidence against AidB-catalyzed DNA repair and supports a model in which AidB acts to prevent alkylation damage by protecting DNA and destroying alkylating agents that have yet to reach their DNA target.

Cloning and Expression of Three Formate Oxidase Genes fromDebaryomyces vanrijiaeMH201

Debaryomyces vanrijiae MH201 produces formate oxidase (FOD) at estimated pI values by native isoelectric focusing of 5.1, 5.4, and 5.9. We cloned and expressed three formate oxidase cDNAs, FOD1, FOD2, and FDO3, of the yeast using Escherichia coli. The open reading frames of FOD1, FOD2, and FDO3 were 1,731 bp long, and encoded 576-amino acid polypeptides with molecular masses of 64,142, 63,794, and 63,836 Da respectively. Expression of FOD1, FOD2, and FOD3 resulted in the production of three isozymes, with pI values of 5.1, 5.9, and 5.9 respectively. Co-expression of FOD1 and FOD2 and of FOD1 and FOD3 resulted in the production of additional isozymes with pI values, of 5.4. The three amino acid sequences of FOD1, FOD2, and FOD3 contained a consensus motif of a flavin adenine dinucleotide binding site in their N-terminal parts and a glucose-methanol-choline oxidoreductase signature pattern, suggesting that formate oxidase ought to be classified in the glucose-methanol-choline oxidoreductase family.

Conformational Switching in the Fungal Light Sensor Vivid

The Neurospora crassa photoreceptor Vivid tunes blue-light responses and modulates gating of the circadian clock. Crystal structures of dark-state and light-state Vivid reveal a light, oxygen, or voltage Per-Arnt-Sim domain with an unusual N-terminal cap region and a loop insertion that accommodates the flavin cofactor. Photoinduced formation of a cystein-flavin adduct drives flavin protonation to induce an N-terminal conformational change. A cysteine-to-serine substitution remote from the flavin adenine dinucleotide binding site decouples conformational switching from the flavin photocycle and prevents Vivid from sending signals in Neurospora. Key elements of this activation mechanism are conserved by other photosensors such as White Collar-1, ZEITLUPE, ENVOY, and flavin-binding, kelch repeat, F-BOX 1 (FKF1).

Environmental stressors (salinity, heavy metals, H 2O 2) modulate expression of glutathione reductase (GR) gene from the intertidal copepod Tigriopus japonicus

Glutathione reductase (GR) plays an essential role in cell defense against reactive oxygen metabolites by sustaining the reduced status of an important antioxidant, glutathione. To address the effect of oxidative stresses on the intertidal copepod Tigriopus japonicus, we exposed specimens to hydrogen peroxide, heavy metals and different salinity levels, cloned and sequenced the oxidative stress-related GR gene. T. japonicus GR gene ( Tigriopus GR) cDNA contained 1526 bp including an open reading frame (ORF) encoding 458 amino acids with a theoretical p I of 6.58 and a calculated molecular weight of 49.6 kDa. Tigriopus GR showed a high similarity to frog Xenopus laevis GR (identity 57%) and the filarial parasite, Onchocerca volvulus GR (identity 57%). Specific motifs such as flavin adenine dinucleotide-binding site (LVLGGGSGGIASARRAAEF), pyridine nucleotide-disulphide oxidoreductases class-I active site (GGTCVNVGCVP), and NADPH binding motif (GxGYIAx18Rx5R) were highly conserved in the deduced amino acid sequence of Tigriopus GR. Interestingly, its expression and enzyme characteristics were different from GR homologue of filarial parasite O. volvulus. To investigate the biochemical and enzymatic characteristics of Tigriopus GR protein, we constructed the expression vector, pCRT7/TOPO NT containing Tigriopus GR. Tigriopus pCRT7/TOPO NT/GR was expressed in Escherichia coli, and the soluble protein was purified by 6× His-tag chromatography. The recombinant Tigriopus GR enzyme was found to make homodimer complexes of approximately 108 kDa on 12% native gel electrophoresis and showed enzymatic activity with NADPH and oxidized glutathione (GSSG) as substrates. To analyze the gene expression of Tigriopus GR against different environmental stresses (hydrogen peroxide, salinity, and heavy metals), we performed real-time reverse transcriptase-polymerase chain reaction (real-time RT-PCR). Slight down-regulation in the expression of Tigriopus GR at the initial stage was observed upon exposure to hydrogen peroxide. The expression recovered in 2 h, while there were significant changes upon heavy metal (Cu and Mn) exposures in a time-dependent manner. Also, Tigriopus GR expression was significantly increased with moderately high salt stress (24 and 40 ppt). In the case of low salt stress (0 and 12 ppt) the expression was found to be down-regulated. These findings provide a better understanding of cellular protection mechanisms in the intertidal copepod T. japonicus against the environmental stressors caused by non-optimal salt levels.

Uncommon missense and splice mutations and resulting biochemical phenotypes in German patients with X-linked chronic granulomatous disease

Chronic granulomatous disease is an inherited disease characterized by the inability of phagocytes to generate normal amounts of superoxide, leaving patients susceptible to opportunistic, life-threatening infections. In the majority of cases, cytochrome b 558 is absent in the X-chromosomal form of CGD. However, the neutrophils from six of nine X-linked CGD patients, reported here, expressed normal or decreased amounts of this cytochrome and are referred to as “variant” forms. In three of these six variant patients, a roughly proportional decrease in cytochrome b 558 expression and production of H 2O 2 were found. In two cases this phenotype could be well explained by special splice mutations, whereas in the third case it was caused by a missense mutation, predicting Ser 193 → Phe. In the other three variant patients, cytochrome b 558 expression and H 2O 2 production were clearly disproportionate as the generation of H 2O 2 was much more decreased than cytochrome expression. Missense mutations also were found in these cases. One of these mutations, predicting Leu 546 → Pro and affecting the putative nicotinamide adenine dinucleotide phosphate binding site, led to normal levels of cytochrome b 558 expression and reduced H 2O 2 production. In the other two mutations, predicting Pro 339 → His and His 338 → Tyr, the putative flavin adenine dinucleotide binding site was affected. This could explain the corresponding uncommon phenotypes, characterized by zero or trace amounts of H 2O 2 production and the expression of relatively high amounts of nonfunctional or low functional cytochrome b 558, respectively. The only missense mutation found that prevented the expression of any cytochrome b 558 was caused by a predicted His 222 → Arg exchange in one of the three classic cases. The two other classic phenotypes were caused by splice mutations.

Comamonas testosteroni 3-ketosteroid-delta 4(5 alpha)-dehydrogenase: gene and protein characterization

Comamonas testosteroni delta 4(5 alpha)- and delta1-dehydrogenases [delta4(5alpha)- and delta1DH] are key enzymes in the degradation of steroids having an A:B ring fusion in a trans configuration. We previously reported the isolation of the delta1dh gene (P. Plesiat, M. Grandguillot, S. Harayama, S. Vragar, and Y. Michel Briand, J. Bacteriol. 173:7219-7227, 1991). In this study, the gene encoding delta 4(5 alpha)DH was cloned in Escherichia coli on a 16-kbp BamHI fragment by screening a genomic bank of C. testosteroni ATCC 17410 with a probe derived from delta1dh. Subcloning experiments in plasmid pUC19 mapped delta 4(5 alpha)dh immediately downstream of delta1dh. The enzyme was overexpressed 18-fold in cells of E. coli JM109 carrying a 2.5-kbp cloned fragment (plasmid pXE25). However, much higher levels of enzymatic activity (264-fold) were obtained in Pseudomonas putida KT2440, using pMMB208 as an expression vector. Studies with crude lysates of KT2440 showed that delta4(5alpha)DH exhibits higher specificity and higher activity toward delta l-androstene-3,17-dione than toward the saturated derivative 5 alpha-androstane-3,17-dione. The reaction was found to be irreversible and to use efficiently typical flavoprotein electron acceptors; optimal conditions for the enzyme activity were pH 8 and 40 degrees C. Analysis of the nucleotide sequence of the insert of pXE25 revealed an open reading frame of 1,593 bp preceded by a putative ribosome-binding site and followed by a potential transcription terminator. The amino acid sequence of the deduced peptide showed a typical flavin adenine dinucleotide-binding site in its N-terminal region, confirming the flavoproteinic structure of delta 4(5 alpha)DH. The predicted molecular mass was consistent with that of the enzyme expressed in a T7 polymerase system (60 kDa). Alignment between delta 4(5 alpha)dh and delta1dh indicated that both genes, though coding for functionally related enzymes, do not derive from a common ancestor.

A molecular dynamics study of glutathione reductase

A 200 ps molecular dynamics trajectory for the 8634 atoms corresponding to the enzyme glutathione reductase were calculated with the GROMOS program. This enzyme is formed by two identical subunits, with 478 residues and a flavinadenine dinucleotide (FAD) molecule each. Four regions are well distinguished in its crystallographic structure, the dimer interface, the FAD binding site, the active site and the crystal contact region. The analysis shows an equilibrated trajectory after 40 ps. The r.m.s. difference between the X-ray structure and the conformation along the simulation is 2.4 Å when considering the r.m.s. of the whole molecule. The structure is locally well conserved. Values of r.m.s. lower than 2.0 Å are found for all regions, with the exception of the crystal contact region. The results give evidence of a fairly stable fold, especially for the catalytically active regions even in the absence of explicit solvent collisions.