scholarly journals Molecular Dynamics Analysis of Glucose Oxidase Stability against Temperature

2021 ◽  
Vol 12 (3) ◽  
pp. 4062-4073

Glucose oxidase (GOD) from local isolated Aspergillus niger IPBCC.08.610 shows a widespread application, specifically as a bioanode in glucose-based biofuel cells. Enzymes with adequate thermal stability are necessary for enhancing product efficiency. Also, evaluating the structural dynamics to improve temperature helps to determine the residue. The molecular dynamics simulation of GOD_IPBCC_1CF3 at temperatures of 300, 400, and 500 K was carried out to analyze important amino acid residues for thermal stability. The results showed that the amino acid residues responsible for thermal stability were dispersed into several essential regions, including D576 at the C terminal, E266-R250, and E38-R237 in the FAD-binding domain E485-R470 in the substrate-binding antiparallel beta system. However, the FAD molecular flexibility against temperature depends on conserve E48 by stabilizing the ribose sugar moiety.

2021 ◽  
Vol 12 ◽  
Author(s):  
Zheng Zequn ◽  
Lian Jiangfang

The rapidly activating delayed rectifier K+ current generated by the cardiac hERG potassium channel encoded by KCNH2 is the most important reserve current for cardiac repolarization. The unique inward rectification characteristics of the hERG channel depend on the gating regulation, which involves crucial structural domains and key single amino acid residues in the full-length hERG channel. Identifying critical molecules involved in the regulation of gating kinetics for the hERG channel requires high-resolution structures and molecular dynamics simulation models. Based on the latest progress in hERG structure and molecular dynamics simulation research, summarizing the molecules involved in the changes in the channel state helps to elucidate the unique gating characteristics of the channel and the reason for its high affinity to cardiotoxic drugs. In this review, we aim to summarize the significant advances in understanding the voltage gating regulation of the hERG channel based on its structure obtained from cryo-electron microscopy and computer simulations, which reveal the critical roles of several specific structural domains and amino acid residues.


Antibiotics ◽  
2021 ◽  
Vol 10 (5) ◽  
pp. 474
Author(s):  
Shazi Shakil ◽  
Syed M. Danish Rizvi ◽  
Nigel H. Greig

Background: Multidrug resistant bacteria are a major therapeutic challenge. CTX-M-type enzymes are an important group of class A extended-spectrum β-lactamases (ESBLs). ESBLs are the enzymes that arm bacterial pathogens with drug resistance to an array of antibiotics, notably the advanced-generation cephalosporins. The current need for an effective CTX-M-inhibitor is high. Objective: The aim of the current study was to identify a promising anti-CTX-M-15 ligand whose chemical skeleton could be used as a ‘seed-molecule’ for future drug design against resistant bacteria. Methods: Virtual screening of 5,000,000 test molecules was performed by ‘MCULE Drug Discovery Platform’. ‘ADME analyses’ was performed by ‘SWISS ADME’. TOXICITY CHECKER of MCULE was employed to predict the safety profile of the test molecules. The complex of the ‘Top inhibitor’ with the ‘bacterial CTX-M-15 enzyme’ was subjected to 102.25 ns molecular dynamics simulation. This simulation was run for 3 days on a HP ZR30w workstation. Trajectory analyses were performed by employing the macro ‘md_analyze.mcr’ of YASARA STRUCTURE version 20.12.24.W.64 using AMBER14 force field. YANACONDA macro language was used for complex tasks. Figures, including RMSD and RMSF plots, were generated. Snapshots were acquired after every 250 ps. Finally, two short videos of ‘41 s’ and ‘1 min and 22 s’ duration were recorded. Results: 5-Amino-1-(2H-[1,2,4]triazino[5,6-b]indol-3-yl)-1H-pyrazole-4-carbonitrile, denoted by the MCULE-1352214421-0-56, displayed the most efficient binding with bacterial CTX-M-15 enzyme. This screened molecule significantly interacted with CTX-M-15 via 13 amino acid residues. Notably, nine amino acid residues were found common to avibactam binding (the reference ligand). Trajectory analysis yielded 410 snapshots. The RMSD plot revealed that around 26 ns, equilibrium was achieved and, thereafter, the complex remained reasonably stable. After a duration of 26 ns and onwards until 102.25 ns, the backbone RMSD fluctuations were found to be confined within a range of 0.8–1.4 Å. Conclusion: 5-Amino-1-(2H-[1,2,4]triazino[5,6-b]indol-3-yl)-1H-pyrazole-4-carbonitrile could emerge as a promising seed molecule for CTX-M-15-inhibitor design. It satisfied ADMET features and displayed encouraging ‘simulation results’. Advanced plots obtained by trajectory analyses predicted the stability of the proposed protein-ligand complex. ‘Hands on’ wet laboratory validation is warranted.


2011 ◽  
Vol 10 (02) ◽  
pp. 245-259 ◽  
Author(s):  
SHANG-ZHI PU ◽  
WEN-HUA ZHANG ◽  
BI SHI

Molecular dynamics simulations were carried out to investigate the effect of pH on structure and stability of collagen-like peptide. All simulations were performed using the consistent valence force field (CVFF) molecular mechanical force field and isothermal-isobaric ensemble (NPT). The initial geometries of the collagen-like peptide were from an X-ray crystallographic structure. Some analyses from the molecular dynamics trajectories have been completed. The results show that the diameter of collagen-like peptide increases and the volume swells obviously in basic environment; however, the size of peptide changes slightly in acidic environment. The stability of collagen-like peptide decreases in acid and basic environment comparing to neutral environment based on root mean square deviation (RMSD). The number of hydrogen bond formed by peptide has a tendency to decrease in both acidic and basic environment. The average of intra-molecular H-bond is minimal under basic condition, and the average of inter-molecular H-bond between amino acid residues and water molecules is minimal under acid condition. The radial distribution function (RDF) shows that side-chain oxygen atoms are easier to form hydrogen bonds with water than side-chain nitrogen atoms. The interaction of various amino acid residues with water is position dependent. Distance between two triple helices increases markedly under highly basic condition, but changes slightly under highly acidic condition.


Author(s):  
Muntahi Mourin ◽  
Arittra Bhattacharjee ◽  
Alvan Wai ◽  
George Hausner ◽  
Joe O'Neil ◽  
...  

The genome of Vibrio cholerae contains three structural genes for the NhaP-type cation-proton antiporter paralogues, Vc-NhaP1, 2 and 3 mediating exchange of K+ and or Na+ for protons across the membrane. Based on phenotype analysis of chromosomal Vc-NhaP1, 2 and 3 triple deletion mutants we suggested that Vc-NhaP paralogues might play a role in the Acid Tolerance Response (ATR) of V. cholerae as it passes through the gastric acid barrier of the stomach. Comparison of the biochemical properties of Vc-NhaP isoforms revealed that Vc-NhaP2 is the most active among all three paralogues. Therefore, Vc-NhaP2 antiporter is a plausible therapeutic target for developing novel inhibitors targeting these ion exchangers. Our structural and mutational analysis of Vc-NhaP2 identified a putative cation binding pocket formed by antiparallel extended regions of two transmembrane segments (TMSs V/XII) along with TMS VI. Molecular Dynamics (MD) simulations suggested that the flexibility of TMS-V/XII is crucial for the intra-molecular conformational events in Vc-NhaP2. In this study, we developed some putative Vc-NhaP2 inhibitors from Amiloride analogs (AAs). Amiloride is a potent inhibitor of human Na+/H+ exchanger-1 (NHE1). Based on the pharmacokinetic properties and potential binding affinity scores we chose six AAs showing high binding affinity scores to Vc-NhaP2. In silico, the six AAs interacted with the functionally important amino acid residues located in TMSs III, IV, V, VI, VIIII and IX either from the cytoplasmic side (three AAs) or the periplasmic side (three AAs) of Vc-NhaP2. Four AAs were modified to reduce their toxicity profile compared to the original AAs. Molecular docking of the modified AAs revealed promising binding. The four selected drugs interacted with functionally important amino acid residues located on the cytoplasmic side of TMS VI, the extended chain region of TMS V and TMS XII and the loop region between TMSs VIIII and IX. Molecular dynamics simulations revealed that binding of the selected drugs destabilized the Vc-NhaP2 and altered the flexibility of functionally important TMS VI.


2020 ◽  
Vol 16 (4) ◽  
pp. 451-459 ◽  
Author(s):  
Fortunatus C. Ezebuo ◽  
Ikemefuna C. Uzochukwu

Background: Sulfotransferase family comprises key enzymes involved in drug metabolism. Oxamniquine is a pro-drug converted into its active form by schistosomal sulfotransferase. The conformational dynamics of side-chain amino acid residues at the binding site of schistosomal sulfotransferase towards activation of oxamniquine has not received attention. Objective: The study investigated the conformational dynamics of binding site residues in free and oxamniquine bound schistosomal sulfotransferase systems and their contribution to the mechanism of oxamniquine activation by schistosomal sulfotransferase using molecular dynamics simulations and binding energy calculations. Methods: Schistosomal sulfotransferase was obtained from Protein Data Bank and both the free and oxamniquine bound forms were subjected to molecular dynamics simulations using GROMACS-4.5.5 after modeling it’s missing amino acid residues with SWISS-MODEL. Amino acid residues at its binding site for oxamniquine was determined and used for Principal Component Analysis and calculations of side-chain dihedrals. In addition, binding energy of the oxamniquine bound system was calculated using g_MMPBSA. Results: The results showed that binding site amino acid residues in free and oxamniquine bound sulfotransferase sampled different conformational space involving several rotameric states. Importantly, Phe45, Ile145 and Leu241 generated newly induced conformations, whereas Phe41 exhibited shift in equilibrium of its conformational distribution. In addition, the result showed binding energy of -130.091 ± 8.800 KJ/mol and Phe45 contributed -9.8576 KJ/mol. Conclusion: The results showed that schistosomal sulfotransferase binds oxamniquine by relying on hybrid mechanism of induced fit and conformational selection models. The findings offer new insight into sulfotransferase engineering and design of new drugs that target sulfotransferase.


RSC Advances ◽  
2014 ◽  
Vol 4 (89) ◽  
pp. 48621-48631 ◽  
Author(s):  
Eleanor R. Turpin ◽  
Sam Mulholland ◽  
Andrew M. Teale ◽  
Boyan B. Bonev ◽  
Jonathan D. Hirst

Sign in / Sign up

Export Citation Format

Share Document