Molecular dynamics study of in silico mutations in the auto-inhibitory loop of human endothelial nitric oxide synthase FMN sub-domain.

Structural flexibility of the peptide linker connecting two domains is essential for the functioning of multi-domain complex. Nitric oxide synthase (NOS) isoforms contain the oxygenase and the reductase domains connected by calmodulin binding linker (CBL) region. Additionally, the endothelial NOS (eNOS) isoform contain an auto-inhibitory loop (AI) in the FMN reductase sub-domain which represses the inter-domain electron transfer process. Binding of Ca2+-Calmodulin complex on the CBL region relieves the AI loop repression and facilitates electron transfer from FMN in the reductase domain to the heme in the oxygenase domain. Few experimental studies have reported that in vitro mutation of Serine-615 (S615D) and Serine-633 (S633D) in the FMN reductase sub-domain to aspartic acid increased NO production and increased Ca2+ sensitivity. To understand the role of AI loop in eNOS repression and activation in serine mutants (S615D and S633D), we modelled the FMN reductase sub-domain of human eNOS protein with and without the CBL region. Molecular dynamics simulations performed indicated that the mutant protein AI loop structure was stabilized by salt bridge formed between D615 and R602. It was also found that mutation increased the flexibility of C-terminal residues of eNOS CBL region. The hinge-like movement of the AI loop allowed rotation of the FMN sub-domain clockwise which may favour electron-transfer in the mutant protein. This study provides insight on mutation (S615D and S633D) induced changes in AI loop and increased flexibility of CBL region which may lead to the protein activation and may also facilitate Calmodulin binding at physiological Ca2+ concentration. Graphical Abstract Mutation of amino acid residues contribute to structural changes at molecular level leading to alteration in protein dynamics and its function. Serine-615 and Serine-633 in the auto-inhibitory loop of human eNOS reductase model was mutated to aspartic acid in silico and molecular dynamics simulations of the protein showed that steric hindrance due to mutation altered the auto-inhibitory loop rearrangement and the FMN sub-domain movement favouring electron transfer.

Molecular dynamics simulation study of Plasmodium falciparum and Escherichia coli SufA: Exploration of conformational changes possibly involved in iron-sulfur cluster transfer.

Iron-sulfur (Fe-S) clusters are one of the earliest known metal complexes in biological molecules. Suf system is one of the Fe-S biogenesis pathways. SufA belongs to the Suf pathway. It is an A-type carrier protein that transfers Fe-S clusters from the scaffold to target proteins. Structural studies were performed for the Suf pathway protein, SufA, in order to explore the conformational changes that probably aid in the transfer of Fe-S clusters to target proteins. Three-dimensional (3D) structure of Plasmodium falciparum (Pf) SufA homodimer was obtained by homology modeling using 3D structure of Escherichia coli (Ec) SufA as template. Molecular dynamics (MD) simulation of Pf SufA and Ec SufA homodimers followed by trajectory and pocket analyses were carried out. A co-ordinated displacement of the homodimeric chains in the interfacial region, resembling a swinging trapeze-like movement was observed. Potential involvement of this swinging trapeze-like movement of the residues belonging to the interfacial region has been proposed as a probable mechanism that assists in the transfer of Fe-S cluster from SufA to apo proteins. This was substantiated by protein-protein interaction studies in Pf SufA by performing molecular docking of 3D conformations of Pf SufA obtained from MD trajectory at every 1 ns interval with Pf ferredoxin.Communicated by Ramaswamy H. Sarma.

Molecular dynamics of the membrane interaction and localisation of prodigiosin.

The tripyrrolic antibiotic prodigiosin causes diverse reactions on its targets like energy spilling, membrane leakage, loss of motility and phototoxicity. It has bacteriostatic, bactericidal, anti-fungal, anti-cancer and immunosuppressive properties. Most of the functions suggest the role of prodigiosin in membrane disruption but the exact mechanism remains unknown. A molecular dynamics study was performed to understand the interactions of prodigiosin with the membrane. It was seen that prodigiosin from the solvent enters the membrane immediately either individually or as small clusters. Prodigiosin clusters with more than eight molecules do not appear to enter the membrane. Upon entry, the molecules orient themselves along the membrane-water interface with the pyrrole rings interacting with lipid head groups and with water. This orientation is stabilised by hydrogen bonding and hydrophobic interactions. The presence of prodigiosin molecules in the membrane changes the local lipid architecture and reduces the solvent accessibility of the membrane. The membrane fluidity, thickness or area per lipid head are largely unaffected. This suggests that prodigiosin could cause most damage in the vicinity of a membrane protein and thus could also explain the reason for varied effects on the targets.

Free enzyme dynamics of CmaA3 and CmaA2 cyclopropane mycolic acid synthases from Mycobacterium tuberculosis: Insights into residues with potential significance in cyclopropanation.

Mycolic acids are long chain alpha-alkyl beta-hydroxy fatty acids that are major constituents of the cell wall of Mycobacterium tuberculosis. M. tuberculosis produces three main types of mycolic acids, alpha mycolic acids and keto and methoxy mycolic acids. Cycloproponated mycolic acids make the cell wall less permeable, contribute to antibiotic resistance and host immunomodulation and protect from injury. Cyclopropanation is catalyzed by enzymes of the Cyclopropane Mycolic Acid Synthase (CMAS) family. In the current study, we addressed two CMAS enzymes, proximal alpha cyclopropane mycolic acid synthase (PcaA/CmaA3) and keto cyclopropane Mycolic acid synthase (CmaA2). All-atom Molecular Dynamics (MD) simulations were performed for these enzymes for a timeframe of 100ns each (in triplicate), using GROMACS. Based on the PDB structures of apo and holo states of related CMAS enzymes, we generated a framework which helped us correlate active or inactive states of the enzymes to different conformations sampled by the enzymes during MD simulations. Dynamics suggested that the free or unbound enzymes have intrinsic memory and sample different states of catalysis even in the absence of the substrate/cofactor. Additionally, we find that F200, P201 and W204 may have functional significance. MD simulation of CmaA2 was performed with the objective of gaining insights into the putative role of a loop insert. Analysis showed that acidic residues of this loop possibly play an important role during the active state by forming salt bridges. The insights gained in this study can potentially be utilized for design of effective inhibitors against CMAS enzymes.

Insights into the pH-dependent catalytic mechanism of Sulfolobus solfataricus ß-glycosidase: A molecular dynamics study.

Sulfolobus solfataricus ß-glycosidase (SS-ßGly) belongs to Glycosyl Hydrolase family1 (GH1) with broad substrate specificity. SS-ßGly catalyzes both hydrolysis and transglycosylation reactions. SS-ßGly is commonly used to synthesize variety of galacto-oligosaccharides. A comparison of SS-ßGly with bacterial and eukaryotic homologs, using DALI search, revealed unique inserts. Free enzyme molecular dynamics (MD) simulation was performed under two different pH conditions (pH 6.5 and 2.5) at a constant temperature of 65 °C using GROMACS. A probable active-site loop (residues 331-364) in SS-ßGly was identified. Dynamics of substrate binding cavity revealed that it was buried and inaccessible during most timeframes at pH 6.5 whereas open and accessible at pH 2.5. New cavities identified during both simulations may act as probable water channel or product egress path. Analyses of docked complexes of 3D structures obtained at every 1ns interval with compounds, involved in hydrolysis and tranglycosylation reactions, demonstrated that conformational states sampled by SS-ßGly during free enzyme dynamics mimic the stages in enzyme catalysis thereby providing a mechanistic perspective. Current study revealed that conformational changes were conducive for hydrolysis at pH 6.5 and multiple cycles of transglycosylation at pH 2.5. Probable role of salt-bridge interactions in determining the type of reaction mechanism was also explored.

Thalidomide and Its Analogs Differentially Target Fibroblast Growth Factor Receptors: Thalidomide Suppresses FGFR Gene Expression while Pomalidomide Dampens FGFR2 Activity.

Thalidomide is an infamous teratogen and it is continuously being explored for its anticancer properties. Fibroblast growth factor receptors (FGFRs) are implicated in embryo development and cancer pathophysiology. With striking similarities observed between FGFR implicated conditions and thalidomide embryopathy, we hypothesized thalidomide targets FGFRs. We utilized three different cell lines and chicken embryo model to investigate the effects of thalidomide and analogs on FGFR expression. We performed molecular docking, KINOMEscan analysis, and kinase activity assays to study the drug-protein interactions. The expression of FGFR1 and FGFR2 was differentially regulated by all the three drugs in cells as well as in developing organs. Transcriptome analysis of thalidomide-treated chick embryo strongly suggests the modulation of FGFR signaling and key transcription factors. Corroboration with previous studies suggests that thalidomide might affect FGFR expression through the transcription factor, E2F1. At the protein level, molecular docking predicted all three analogs to interact with lysine residue at 517th and 508th positions of FGFR2 and FGFR3, respectively. This lysine coordinates the ATP binding site of FGFR, thus hinting at the possible perturbation of FGFR activity by thalidomide. Kinome analysis revealed that kinase activities of FGFR2 and FGFR3 (G697C) reduced by 31% and 65%, respectively, in the presence of 10 µM thalidomide. Further, we checked and confirmed that the analogs inhibited the FGFR2 kinase activity in a dose-dependent manner. This study suggests that FGFRs could be potential targets of thalidomide and the two analogs, and also endorses the link between the teratogenicity and antitumor activities of the drugs.

Potential role of salt-bridges in the hinge-like movement of apicomplexa specific ß-hairpin of Plasmodium and Toxoplasma profilins: A molecular dynamics simulation study.

Profilin is one of the actin-binding proteins that regulate dynamics of actin polymerization. It plays a key role in cell motility and invasion. It also interacts with several other proteins notably through its poly-L-proline (PLP) binding site. Profilin in apicomplexa is characterized by a unique mini-domain consisting of a large ß-hairpin extension and an acidic loop which is relatively longer in Plasmodium species. Profilin is essential for the invasive blood stages of Plasmodium falciparum. In the current study, unbound profilins from Plasmodium falciparum (Pf), Toxoplasma gondii (Tg), and Homo sapiens (Hs) were subjected to molecular dynamics (MD) simulations for a timeframe of 100 ns each to understand the conformational dynamics of these proteins. It was found that the ß-hairpin of profilins from Pf and Tg shows a hinge-like movement. This movement in Pf profilin may possibly be driven by the loss of a salt-bridge within profilin. The impact of this conformational change on actin binding was assessed by docking three dimensional (3D) structures of profilin from Pf and Tg with their corresponding actins using ClusPro2.0. The stability of docked Pf profilin-actin complex was assessed through a 50 ns MD simulation. As Hs profilin I does not have the apicomplexa specific mini-domain, MD simulation was performed for this protein and its dynamics was compared to that of profilins from Pf and Tg. Using an immunoinformatics approach, potential epitope regions were predicted for Pf profilin. This has a potential application in the design of vaccines as they mapped to its unique mini-domain.

A pipeline for proteome-scale identification and studies on hormone sensitive lipases in Mycobacterium tuberculosis.

Hormone sensitive lipases (HSLs) play an important role in the survival of M. tuberculosis during dormancy. They help in the utilization of fatty acids from stored lipids. The objective of the current study was to identify all HSLs from the proteome of M. tuberculosis H37Rv. We have developed a novel HSL identification pipeline, based on amino acid sequence homology, presence of conserved motifs and other sequence features deciphered from known HSL dataset. Through this pipeline, we identified 10 proteins as putative HSLs in M. tuberculosis. We have annotated a lipase LipT, as putative p-nitrobenzyl esterase and also identified a new motif "PGG" which is a possible characteristic motif of a subfamily of HSLs.

Comparative pan genome analysis of oral Prevotella species implicated in periodontitis.

Prevotella is part of the oral bacterial community implicated in periodontitis. Pan genome analyses of eight oral Prevotella species, P. dentalis, P. enoeca, P. fusca, P. melaninogenica, P. denticola, P. intermedia 17, P. intermedia 17-2 and P. sp. oral taxon 299 are presented in this study. Analysis of the Prevotella pan genome revealed features such as secretion systems, resistance to oxidative stress and clustered regularly interspaced short palindromic repeat (CRISPR)-Cas systems that enable the bacteria to adapt to the oral environment. We identified the presence of type VI secretion system (T6SS) in P. fusca and P. intermedia strains. For some VgrG and Hcp proteins which were not part of the core T6SS loci, we used gene neighborhood analysis and identified putative effector proteins and putative polyimmunity loci in P. fusca and polymorphic toxin systems in P. intermedia strains. Earlier studies have identified the presence of Por secretion system (PorSS) in P. gingivalis, P. melaninogenica and P. intermedia. We noted the presence of their homologs in six other oral Prevotella studied here. We suggest that in Prevotella, PorSS is used to secrete cysteine proteases such as interpain and C-terminal domain containing proteins with a "Por_secre_tail" domain. We identified subtype I-B CRISPR-Cas system in P. enoeca. Putative CRISPR-Cas system subtypes for 37 oral Prevotella and 30 non-oral Prevotella species were also predicted. Further, we performed a BLASTp search of the Prevotella proteins which are also conserved in the red-complex pathogens, against the human proteome to identify potential broad-spectrum drug targets. In summary, the use of a pan genome approach enabled identification of secretion systems and defense mechanisms in Prevotella that confer adaptation to the oral cavity.

Functional insights from a comparative study on the dynamics of Antigen85 proteins and MPT51 from Mycobacterium tuberculosis.

Antigen85 (Ag85) proteins of Mycobacterium tuberculosis are mycolyl transferases that aid in cell wall biosynthesis. MPT51 (Ag85D) is closely related to Ag85 proteins. We have performed a comparative molecular dynamics (MD) simulation study of Ag85 proteins (Ag85A, Ag85B, and Ag85C) and MPT51. We observe that helix α5, ß7-α9 loop, and N-terminal region of helix α9 of Ag85 proteins are mobile, suggestive of lid like movement over the active site. Further, in Ag85B, we observe the proposed scooting mode of the hydrophobic gating residue Phe232. Our simulations also show a similar scooting mode for Phe232 of Ag85A and Trp158 of Ag85C. We also found aromatic residue clusters at the ends of the hydrophobic channel of Ag85 proteins, which may have functional significance. Although MPT51 lacks the tunnel, it has the aromatic clusters. The aromatic cluster region has the ability to bind trehalose. From an immunoinformatics study, a promiscuous linear epitope was identified in MPT51 which could be useful in subunit vaccine studies. Recent studies have shown that a mycobacterial protein HupB, interacts with Ag85 proteins and has a regulatory role in cell wall biogenesis, with implications in growth rate and latency. We performed molecular docking studies of HupB protein with Ag85 proteins and predicted potential sites of interaction in Ag85 proteins. The insights gained through the current study can potentially pave way for newer therapeutic interventions. Graphical Abstract Dynamics of antigen85 proteins and MPT51 from Mycobacterium tuberculosis.

Cavities create a potential back door in epoxide hydrolase Rv1938 from Mycobacterium tuberculosis-A molecular dynamics simulation study.

Mycobacterium tuberculosis (Mtb) is the causative organism of tuberculosis. Extensively drug resistant strains and latency have posed formidable challenges in the treatment of tuberculosis. The current study addresses an alpha/beta hydrolase fold bearing enzyme, epoxide hydrolase Rv1938 from Mtb. Epoxide hydrolases are involved in detoxification processes, catabolism and regulation of signaling molecules. Using GROMACS, a 100ns Molecular Dynamics (MD) simulation was performed for Rv1938. Cavities were identified within the protein at various time frames of the simulation and their volumes were computed. During MD simulation, in addition to the substrate binding cavity, opening of two new cavities located behind the active site was observed. These cavities may be similar to the backdoor proposed for acetylcholinesterase. Structural superimposition of epoxide hydrolase from Mtb with the epoxide hydrolase of Agrobacterium radiobacter1 AD1 (Ephy) indicates that cavity1 in Mtb lies at an identical position to that of the water tunnel in Ephy. Further, docking of the substrate and an inhibitor with protein structures obtained from MD simulation at various time frames was also performed. The potential role of these cavities is discussed.

Mechanistic insights from molecular dynamic simulation of Rv0045c esterase in Mycobacterium tuberculosis.

Rv0045c is an esterase involved in lipid metabolism of Mycobacterium tuberculosis. It belongs to the α/ß hydrolase family. In the current study, we performed sequence- and structure-based analysis of Rv0045c followed by molecular dynamics (MD) simulation for 100 ns to investigate conformational changes in the enzyme. Sequence analysis revealed that this enzyme is possibly a hormone-sensitive lipase. Further, through structural analysis, a putative catalytic tetrad containing "Ser-Asp-Ser-His" and residues involved in the formation of an oxyanion hole were identified. MD simulation of Rv0045c revealed a conformational transition from an open to a closed state. The active site pocket was found to be gated by four loops. The potential role of the cap domain and the mobile histidine is discussed. From the simulation, we see that the conformational changes mimic the different stages in the reaction mechanism of Rv0045c. These results support the hypothesis that free enzyme simulation encompasses all the conformations necessary for the different stages of catalysis. Our findings add to the growing knowledge of an important family of esterases in Mycobacterium tuberculosis.

Molecular dynamics simulations of lectin domain of FimH and immunoinformatics for the design of potential vaccine candidates.

Adhesion of uropathogenic E. coli (UPEC) to uroepithelial cell receptors is facilitated through the lectin domain of FimH adhesin. In the current study, Molecular Dynamics (MD) simulations were performed for the lectin domain of FimH from UPEC J96. The high affinity state lectin domain was found to be stable and rigid during the simulations. Further, based on conserved subsequences around one of the disulfide forming cysteines, two sequence motifs were designed. An immunoinformatics approach was utilized to identify linear and discontinuous epitopes for the lectin domain of FimH. We propose that the accessibility of predicted epitopes should also be assessed in a dynamic aqueous environment to evaluate the potential of vaccine candidates. Since MD simulation data enables assessing the accessibility in a dynamic environment, we evaluated the accessibility of the top ranked discontinuous and linear epitopes using structures obtained at every nanosecond (ns) in the 1-20 ns MD simulation timeframe. Knowledge gained in this study has a potential utility in the design of vaccine candidates for Urinary Tract Infection (UTI).

Molecular Dynamics simulations of Inhibitor of Apoptosis Proteins and identification of potential small molecule inhibitors.

Chemotherapeutic resistance due to over expression of Inhibitor of Apoptosis Proteins (IAPs) XIAP, survivin and livin has been observed in various cancers. In the current study, Molecular Dynamics (MD) simulations were carried out for all three IAPs and a common ligand binding scaffold was identified. Further, a novel sequence based motif specific to these IAPs was designed. SMAC is an endogenous inhibitor of IAPs. Screening of ChemBank for compounds similar to lead SMAC-non-peptidomimetics yielded a cemadotin related compound NCIMech_000654. Cemadotin is a derivative of natural anti-tumor peptide dolastatin-15; hence these compounds were docked against all three IAPs. Based on our analysis, we propose that NCIMech_000654/dolastatin-15/cemadotin derivatives may be investigated for their potential in inhibiting XIAP, survivin and livin.

Understanding the lid movements of LolA in Escherichia coli using molecular dynamics simulation and in silico point mutation.

The Lol system in Escherichia coli is involved in localization of lipoproteins and hence is essential for growth of the organism. LolA is a periplasmic chaperone that binds to outer-membrane specific lipoproteins and transports them from inner membrane to outer membrane through LolB. The hydrophobic lipid-binding cavity of LolA consists of α-helices which act as a lid in regulating the transfer of lipoproteins from LolA to LolB. The current study aims to investigate the structural changes observed in LolA during the transition from open to closed conformation in the absence of lipoprotein. Molecular dynamics (MD) simulations were carried out for two LolA crystal structures; LolA(R43L), and in silico mutated MsL43R for a simulation time of 50 ns in water environment. We have performed an in silico point mutation of leucine to arginine in MsL43R to evaluate the importance of arginine to induce structural changes and impact the stability of protein structure. A complete dynamic analysis of open to closed conformation reveals the existence of two distinct levels; closing of lid and closing of entrance of hydrophobic cavity. Our analysis reveals that the structural flexibility of LolA is an important factor for its role as a periplasmic chaperone.

A matrix based algorithm for Protein-Protein Interaction prediction using Domain-Domain Associations.

Protein-Protein Interactions (PPI) are vital to many cellular processes. The availability of high-throughput protein interaction data has provided us with an opportunity to assess domain associations in interacting proteins using computational approaches. High throughput PPI data, wherein the interaction status of every protein in the dataset has been experimentally tested against all the other proteins in the dataset contains information not only on protein interactions but also on proteins which do not interact with each other. We call such datasets "all against all" datasets. In the current study, using these datasets and the Pfam domain composition of the proteins in the sets, we have developed a matrix based method for predicting PPI. We infer positive and negative Domain-Domain Associations (DDA) by our method. We have generated more than a million domain association values which can be utilized for predicting new PPI. The performance of the algorithm was evaluated against a test set and the sensitivity and specificity was found to be 68.1% and 65.3%, respectively. The overall prediction accuracy of the algorithm with individual test sets from IntAct, DIP, 3did, iPfam databases and a literature curated set from Saccharomyces cerevisiae was found to be around 70%. The insights gained in the study have a potential application in providing leads for experimental interaction studies and understanding host pathogen interactions amongst others.

A meta-metabolome network of carbohydrate metabolism: interactions between gut microbiota and host.

With the current knowledge of the multitude of microbes that inhabit the human body, it is increasingly clear that they constitute an integral component of the host. The gut microbiota community is principally involved in the metabolism of dietary constituents such as carbohydrates which account for majority of the energy intake from diet. Diet has gained an important role in shaping the composition of gut inhabitants. The quantity and type of food consumed is recognized as a causal factor for metabolic disorders such as obesity and diabetes. Analysis of host-microbe interactions can thus contribute to the understanding of such metabolic disorders. In this study, data from Kyoto Encyclopedia of Genes and Genomes and Carbohydrate Active EnZYmes Database was utilized as a starting point. Enzyme information from the host Homo sapiens coupled with details of the three predominant phyla of gut bacteria, namely Firmicutes, Bacteroidetes and Actinobacteria were used in the creation of a comprehensive metabolic network, which we refer to as 'meta-metabolome'. This 'meta-metabolome' provides a perspective of the degree to which microbes influence carbohydrate metabolism, in conjunction with host specific enzymes. Analysis of reactions in the network reveals the amplification of monosaccharide content brought about by microbial enzyme activity. The framework outlined in this study provides a holistic approach to assess host-microbe symbiosis. It also provides us with a means of analyzing how diet can be modulated to provide beneficial effects to the host or how probiotics can potentially be used to relieve certain metabolic disorders.

In silico analysis of DosR regulon proteins of Mycobacterium tuberculosis.

One of the challenges faced by Mycobacterium tuberculosis (M. tuberculosis) in dormancy is hypoxia. DosR/DevR of M. tuberculosis is a two component dormancy survival response regulator which induces the expression of 48 genes. In this study, we have used DosR regulon proteins of M. tuberculosis H37Rv as the query set and performed a comprehensive homology search against the non-redundant database. Homologs were found in environmental mycobacteria, environmental bacteria and archaebacteria. Analysis of genomic context of DosR regulon revealed that they are distributed as nine blocks in the genome of M. tuberculosis with many transposases and integrases in their vicinity. Further, we classified DosR regulon proteins into eight functional categories. One of the hypothetical proteins Rv1998c could probably be a methylisocitrate lyase or a phosphonomutase. Another hypothetical protein, Rv0572 was found only in mycobacteria. Insights gained in this study can potentially aid in the development of novel therapeutic interventions.

Molecular dynamics simulations of human and dog gastric lipases: insights into domain movements.

Mammalian gastric lipases are stable and active under acidic conditions and also in the duodenal lumen. There has been considerable interest in acid stable lipases owing to their potential application in the treatment of pancreatic exocrine insufficiency. In order to gain insights into the domain movements of these enzymes, molecular dynamics simulations of human gastric lipase was performed at an acidic pH and under neutral conditions. For comparative studies, simulation of dog gastric lipase was also performed at an acidic pH. Analyses show, that in addition to the lid region, there is another region of high mobility in these lipases. The potential role of this novel region is discussed.

Investigations on domain movements of N-acetyltransferase in nano scale.

Molecular Dynamics simulations were performed on aminoglycoside N-acetyltransferase ((AAC) (2')-Ic) of Mycobacterium tuberculosis, to gain insight into enzyme flexibility and conformation around the ligand and acetyl-CoA binding sites in nano scale. Simulations were performed in water and methanol to further study the effect of solvation on the enzyme. AAC(2')-Ic of M. tuberculosis consists of negatively charged residues of aspartic and glutamic acids, which are believed to form a docking platform for the positively charged aminoglycosides. The acetyl-CoA binding site involves a wide groove surrounded by helices alphal, alpha3, alpha4 and strand beta4. Simulation results revealed the enzyme to be stable in water and the enzyme was found to be flexible around the docking platform (ligand binding site) with Asp 35, Asp 40 and Asp 179 exhibiting maximum movement. An interesting observation was the flexibility and conformational change at beta4, hydrophobic pocket of the acetyl-CoA binding site a prerequisite for catalysis to occur with Leu 95 exhibiting maximum displacement with a root mean square fluctuation of 0.18 nm. The enzyme was found to be very unstable in methanol with root mean square deviation not stabilizing even at 10 ns.