Paenibacillus polymyxa β-glucosidase B (BglB), belongs to a GH family 1, is a monomeric enzyme that acts as an exo-β-glucosidase hydrolysing cellobiose and cellodextrins of higher degree of polymerization using retaining mechanism. A molecular dynamics (MD) simulation was performed at 300 K under periodic boundary condition for 5 ns using the complexes structure obtained from previous docking study, namely BglB-Beta-d-glucose and BglB-Cellobiose. From the root-mean-square deviation analysis, both enzyme complexes were reported to deviate from the initial structure in the early part of the simulation but it was stable afterwards. The root-mean-square fluctuation analysis revealed that the most flexible regions comprised of the residues from 26 to 29, 43 to 53, 272 to 276, 306 to 325 and 364 to 367. The radius of gyration analysis had shown the structure of BglB without substrate became more compact towards the end of the simulation compare to other two complexes. The residues His122 and Trp410 were observed to form stable hydrogen bond with occupancy higher than 10%. In conclusion, the behaviour of BglB enzyme towards the substrate binding was successfully explored via MD simulation approaches.
Indolizines are heteroaromatic compounds, and their synthetic analogues have reportedly showed promising pharmacological properties. In this study, a series of synthetic 7-methoxy-indolizine derivatives were synthesised, characterised and evaluated for in vitro whole-cell anti-tuberculosis (TB) screening against susceptible (H37Rv) and multi-drug-resistant (MDR) strains of Mycobacterium tuberculosis (MTB) using the resazurin microplate assay method. The cytotoxicity was evaluated using the MTT assay. In silico molecular-docking study was conducted for compounds 5a-j against enoyl-[acyl-carrier] protein reductase, a key enzyme of the type II fatty acid synthesis that has attracted much interest for the development of novel anti-TB compounds. Thereafter, molecular dynamic (MD) simulation was undertaken for the most active inhibitors. Compounds 5i and 5j with the methoxy functional group at the meta position of the benzoyl group, which was at the third position of the indolizine nucleus, demonstrated encouraging anti-TB activity against MDR strains of MTB at 16 μg/mL. In silico studies showed binding affinity within the range of 7.07-8.57 kcal/mol, with 5i showing the highest binding affinity. Hydrogen bonding, π-π- interactions, and electrostatic interactions were common with the active site. Most of these interactions occurred with the catalytic amino acids (Pro193, Tyr158, Phe149, and Lys165). MD simulation showed that 5j possessed the highest binding affinity toward the enzyme, according to the two calculation methods (MM/PBSA and MM/GBSA). The single-crystal X-ray studies of compounds 5c and 5d revealed that the molecular arrangements in these two structures were mostly guided by C-H···O hydrogen-bonded dimeric motifs and C-H···N hydrogen bonds, while various secondary interactions (such as π···π and C-H···F) also contributed to crystal formation. Compounds 5a, 5c, 5i, and 5j exhibited no toxicity up to 500 μg/mL. In conclusion, 5i and 5j are promising anti-TB compounds that have shown high affinity based on docking and MD simulation results.
The utilization of organic solvents as reaction media for enzymatic reactions provides numerous industrially attractive advantages. However, an adaptation of enzyme towards organic solvent is unpredictable and not fully understood because of limited information on the organic solvent tolerant enzymes. To understand how the enzyme can adapt to the organic solvent environment, structural and computational approaches were employed. A recombinant elastase from Pseudomonas aeruginosa strain K was an organic solvent tolerant zinc metalloprotease was successfully crystallized and diffracted up to 1.39 Å. Crystal structure of elastase from strain K showed the typical, canonical alpha-beta hydrolase fold consisting of 10-helices (118 residues), 10- β-strands (38 residues) and 142 residues were formed other secondary structure such as loop and coil to whole structure. The elastase from Pseusomonas aeruginosa strain K possess His-140, His-144 and Glu-164 served as a ligand for zinc ion. The conserved catalytic triad was composed of Glu-141, Tyr-155 and His-223. Three-dimensional structure features such as calcium-binding and presence of disulphide-bridge contribute to the stabilizing the elastase structure. Molecular dynamic (MD) simulation of elastase revealed that, amino acid residues located at the surface area and disulphide bridge in Cys-30 to Cys-58 were responsible for enzyme stability in organic solvents.
DNA polymerase β is a 39 kDa enzyme that is a major component of Base Excision Repair in human cells. The enzyme comprises two major domains, a 31 kDa domain responsible for the polymerase activity and an 8 kDa domain, which bind ssDNA and has a deoxyribose phosphate (dRP) lyase activity. DNA polymerase β was shown to be phosphorylated in vitro with protein kinase C (PKC) at serines 44 and 55 (S44 and S55), resulting in loss of its polymerase enzymic activity, but not its ability to bind ssDNA. In this study, we investigate the potential phosphorylation-induced structural changes for DNA polymerase β using molecular dynamics simulations. The simulations show drastic conformational changes of the polymerase structure as a result of S44 phosphorylation. Phosphorylation-induced conformational changes transform the closed (active) enzyme structure into an open one. Further analysis of the results points to a key hydrogen bond and newly formed salt bridges as potential drivers of these structural fluctuations. The changes observed with S55/44 and S55 phosphorylation were less dramatic and the integrity of the H-bond was not compromised. Thus the phosphorylation of S44 is the major contributor to structural fluctuations that lead to loss of enzymatic activity.
Both the inactive- and active-state CB1 receptor crystal structures have now been solved, allowing their application in various structure-based drug design methods. One potential method utilizing these crystal structures is the Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) method of predicting relative binding free.
Adamantyl-based compounds are clinically important for the treatments of type 2 diabetes and for their antiviral abilities, while many more are under development for other pharmaceutical uses. This study focused on the acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) inhibitory activities of adamantyl-based ester derivatives with various substituents on the phenyl ring using Ellman's colorimetric method. Compound 2e with a 2,4-dichloro electron-withdrawing substituent on the phenyl ring exhibited the strongest inhibition effect against AChE, with an IC50 value of 77.15 µM. Overall, the adamantyl-based ester with the mono-substituent at position 3 of the phenyl ring exhibited good AChE inhibition effects with an ascending order for the substituents: Cl < NO₂ < CH₃ < OCH₃. Furthermore, compounds with electron-withdrawing groups (Cl and NO₂) substituted at position 3 on their phenyl rings demonstrated stronger AChE inhibition effects, in comparison to their respective positional isomers. On the other hand, compound 2j with a 3-methoxyphenyl ring showed the highest inhibition effect against BChE, with an IC50 value of 223.30 µM. Molecular docking analyses were conducted for potential AChE and BChE inhibitors, and the results demonstrated that the peripheral anionic sites of target proteins were predominant binding sites for these compounds through hydrogen bonds and halogen interactions instead of hydrophobic interactions in the catalytic active site.
The crystallins are a family of monomeric proteins present in the mammalian lens and mutations in these proteins cause various forms of cataracts. The aim of our current study is to emphasize the structural characterization of aggregation propensity of mutation R58H on γD crystallin using molecular dynamics (MD) approach. MD result revealed that difference in the sequence level display a wide variation in the backbone atomic position, and thus exhibits rigid conformational dynamics. Changes in the flexibility of residues favoured to increase the number of intra-molecular hydrogen bonds in mutant R58H. Moreover, notable changes in the hydrogen bonding interaction resulted to cause the misfolding of mutant R58H by introducing α-helix. Principal component analysis (PCA) result suggested that mutant R58H showed unusual conformational dynamics along the two principal components when compared to the wild-type (WT)-γD crystallin. In a nutshell, the increased surface hydrophobicity could be the cause of self-aggregation of mutant R58H leading to aculeiform cataract.
The l-2-haloacid dehalogenases (EC 3.8.1.2) specifically cleave carbon-halogen bonds in the L-isomers of halogenated organic acids. These enzymes have potential applications for the bioremediation and synthesis of various industrial products. One such enzyme is DehL, the l-2-haloacid dehalogenase from Rhizobium sp. RC1, which converts the L-isomers of 2-halocarboxylic acids into the corresponding D-hydroxycarboxylic acids. However, its catalytic mechanism has not been delineated, and to enhance its efficiency and utility for environmental and industrial applications, knowledge of its catalytic mechanism, which includes identification of its catalytic residues, is required. Using ab initio fragment molecular orbital calculations, molecular mechanics Poisson-Boltzmann surface area calculations, and classical molecular dynamic simulation of a three-dimensional model of DehL-l-2-chloropropionic acid complex, we predicted the catalytic residues of DehL and propose its catalytic mechanism. We found that when Asp13, Thr17, Met48, Arg51, and His184 were individually replaced with an alanine in silico, a significant decrease in the free energy of binding for the DehL-l-2-chloropropionic acid model complex was seen, indicating the involvement of these residues in catalysis and/or structural integrity of the active site. Furthermore, strong inter-fragment interaction energies calculated for Asp13 and L-2-chloropropionic acid, and for a water molecule and His184, and maintenance of the distances between atoms in the aforementioned pairs during the molecular dynamics run suggest that Asp13 acts as the nucleophile and His184 activates the water involved in DehL catalysis. The results of this study should be important for the rational design of a DehL mutant with improved catalytic efficiency.
Tuberculosis (TB) still remains a global threat due to the emergence of a drug-resistant strain. Instead of focusing on the drug target of active stage TB, we are highlighting the isocitrate lyase (ICL) at the dormant stage TB. ICL is one of the persistent factors for Mycobacterium tuberculosis (MTB) to survive during the dormant phase. In addition, the absence of ICL in human has made ICL a potential drug target for TB therapy. However, the dynamic details of ICL which could give insights to the ICL-ligand interaction have yet to be solved. Therefore, a series of ICL dimer dynamics studies through molecular dynamics simulation were performed in this work. The ICL active site entrance gate closure is contributed to by hydrogen bonding and electrostatic interactions with the C-terminal. Analysis suggested that the open-closed behavior of the ICL active site entrance depends on the type of ligand present in the active site. We also observed four residues (Ser91, Asp108, Asp153, and Cys191) which could possibly be the nucleophiles for nucleophilic attack on the cleavage of isocitrate at the C2-C3bond. We hope that the elucidation of ICL dynamics can benefit future works such as lead identification or antibody design against ICL for TB therapeutics.
Palm oil-based esters (POEs) are unsaturated and non-ionic esters with a great potential to act as chemical penetration enhancers and drug carriers for transdermal drug nano-delivery. A ratio of palmitate ester and nonionic Tween80 with and without diclofenac acid was chosen from an experimentally determined phase diagram. Molecular dynamics simulations were performed for selected compositions over a period of 15 ns. Both micelles showed a prolate-like shape, while adding the drug produced a more compact micellar structure. Our results proposed that the drug could behave as a co-surfactant in our simulated model.
The distinctive morphology characteristics of microfold cells (M cells) allow the vaccine antigen not only to interact with immune cells directly, but also to effectively stimulate mucosal immune responses via receptors on its apical surface. Human prion protein, a transmembrane receptor for Brucella abortus Hsp60, is highly expressed on the M cell surface. Nonetheless, this protein tends to express in inclusion body in prokaryotic hosts. In this study, the shorter interacting regions of human prion protein were identified via computational methods such as docking and molecular dynamics simulations to minimize its aggregation tendency. The computational calculations revealed three novel human prion protein-interacting regions, namely PrP125, PrP174, and PrP180. In accordance with in silico prediction, the biologically synthesized peptides fusing with GST tag demonstrated their specific binding to Hsp60 protein via pull-down assay. Hence, this finding laid the groundwork for M-cell targeting candidate validation through these newly identified interacting regions.
Genome mining revealed a 1011 nucleotide-long fragment encoding a type I L-asparaginase (J15 asparaginase) from the halo-tolerant Photobacterium sp. strain J15. The gene was overexpressed in pET-32b (+) vector in E. coli strain Rosetta-gami B (DE3) pLysS and purified using two-step chromatographic methods: Ni(2+)-Sepharose affinity chromatography and Q-Sepharose anion exchange chromatography. The final specific activity and yield of the enzyme achieved from these steps were 20 U/mg and 49.2%, respectively. The functional dimeric form of J15-asparaginase was characterised with a molecular weight of ~70 kDa. The optimum temperature and pH were 25°C and pH 7.0, respectively. This protein was stable in the presence of 1 mM Ni(2+) and Mg(2+), but it was inhibited by Mn(2+), Fe(3+) and Zn(2+) at the same concentration. J15 asparaginase actively hydrolysed its native substrate, l-asparagine, but had low activity towards l-glutamine. The melting temperature of J15 asparaginase was ~51°C, which was determined using denatured protein analysis of CD spectra. The Km, Kcat, Kcat/Km of J15 asparaginase were 0.76 mM, 3.2 s(-1), and 4.21 s(-1) mM(-1), respectively. Conformational changes of the J15 asparaginase 3D structure at different temperatures (25°C, 45°C, and 65°C) were analysed using Molecular Dynamic simulations. From the analysis, residues Tyr₂₄ , His₂₂, Gly₂₃, Val₂₅ and Pro₂₆ may be directly involved in the 'open' and 'closed' lid-loop conformation, facilitating the conversion of substrates during enzymatic reactions. The properties of J15 asparaginase, which can work at physiological pH and has low glutaminase activity, suggest that this could be a good candidate for reducing toxic effects during cancer treatment.
Halogenated compounds are recalcitrant environmental pollutants prevalent in agricultural fields, waste waters and industrial by-products, but they can be degraded by dehalogenase-containing microbes. Notably, 2-haloalkanoic acid dehalogenases are employed to resolve optically active chloropropionates, as exemplified by the d-specific dehalogenase from Rhizobium sp. RCI (DehD), which acts on d-2-chloropropionate but not on its l-enantiomer. The catalytic residues of this dehalogenase responsible for its affinity toward d-2-chloropropionate have not been experimentally determined, although its three-dimensional crystal structure has been solved. For this study, we performed in silico docking and molecular dynamic simulations of complexes formed by this dehalogenase and d- or l-2-chloropropionate. Arg134 of the enzyme plays the key role in the stereospecific binding and Arg16 is in a position that would allow it to activate a water molecule for hydrolytic attack on the d-2-chloropropionate chiral carbon for release of the halide ion to yield l-2-hydroxypropionate. We propose that within the DehD active site, the NH group of Arg134 can form a hydrogen bond with the carboxylate of d-2-chloropropionate with a strength of ∼4 kcal/mol that may act as an acid-base catalyst, whereas, when l-2-chloropropionate is present, this bond cannot be formed. The significance of the present work is vital for rational design of this dehalogenase in order to confirm the involvement of Arg16 and Arg134 residues implicated in hydrolysis and binding of d-2-chloropropionate in the active site of d-specific dehalogenase from Rhizobium sp. RC1.
Glycolipid, found commonly in membranes, is also a liquid crystal material which can self-assemble without the presence of a solvent. Here, the dielectric and conductivity properties of three synthetic glycolipid thin films in different thermotropic liquid crystal phases were investigated over a frequency and temperature range of (10(-2)-10(6) Hz) and (303-463 K), respectively. The observed relaxation processes distinguish between the different phases (smectic A, columnar/hexagonal, and bicontinuous cubic Q) and the glycolipid molecular structures. Large dielectric responses were observed in the columnar and bicontinuous cubic phases of the longer branched alkyl chain glycolipids. Glycolipids with the shortest branched alkyl chain experience the most restricted self-assembly dynamic process over the broad temperature range studied compared to the longer ones. A high frequency dielectric absorption (Process I) was observed in all samples. This is related to the dynamics of the hydrogen bond network from the sugar group. An additional low-frequency mechanism (Process II) with a large dielectric strength was observed due to the internal dynamics of the self-assembly organization. Phase sensitive domain heterogeneity in the bicontinuous cubic phase was related to the diffusion of charge carriers. The microscopic features of charge hopping were modelled using the random walk scheme, and two charge carrier hopping lengths were estimated for two glycolipid systems. For Process I, the hopping length is comparable to the hydrogen bond and is related to the dynamics of the hydrogen bond network. Additionally, that for Process II is comparable to the bilayer spacing, hence confirming that this low-frequency mechanism is associated with the internal dynamics within the phase.
We describe a server that allows the interrogation of the Protein Data Bank for hypothetical 3D side chain patterns that are not limited to known patterns from existing 3D structures. A minimal side chain description allows a variety of side chain orientations to exist within the pattern, and generic side chain types such as acid, base and hydroxyl-containing can be additionally deployed in the search query. Moreover, only a subset of distances between the side chains need be specified. We illustrate these capabilities in case studies involving arginine stacks, serine-acid group arrangements and multiple catalytic triad-like configurations. The IMAAAGINE server can be accessed at http://mfrlab.org/grafss/imaaagine/.
DAPK3 belongs to family of DAPK (death-associated protein kinases) and is involved in the regulation of progression of the cell cycle, cell proliferation, apoptosis and autophagy. It is considered as a tumor suppressor kinase, suggesting the loss of its function in case of certain specific mutations. The T112M, D161N and P216S mutations in DAPK3 have been observed in cancer patients. These DAPK3 mutants have been associated with very low kinase activity, which results in the cellular progression towards cancer. However, a clear understanding of the structural and biophysical variations that occur in DAPK3 with these mutations, resulting in the decreased kinase activity has yet not been deciphered. We performed a molecular dynamic simulation study to investigate such structural variations. Our results revealed that mutations caused a significant structural variation in DAPK3, majorly concentrated in the flexible loops that form part of the ATP binding pocket. Interestingly, D161N and P216S mutations collapsed the ATP binding pocket through flexible loops invasion, hindering ATP binding which resulted in very low kinase activity. On the contrary, T112M mutant DAPK3 reduces ATP binding potential through outward distortion of flexible loops. In addition, the mutant lacked characteristic features of the active protein kinase including proper interaction between HR/FD and DFG motifs, well structured hydrophobic spine and Lys42-Glu64 salt bridge interaction. These observations could possibly explain the underlying mechanism associated with the loss of kinase activity with T112M, D161N and P216S mutation in DAPK3.
The limited sequence similarity of protein sequences with known structures has led to an indispensable need for computational technology to predict their structures. Structural bioinformatics (SB) has become integral in elucidating the sequence-structure-function relationship of a protein. This report focuses on the applications of SB within the context of protein engineering including its limitation and future challenges.
Immobilisation of aptameric ligands on solid stationary supports for effective binding of target molecules requires understanding of the relationship between aptamer-polymer interactions and the conditions governing the mass transfer of the binding process. Herein, key process parameters affecting the molecular anchoring of a thrombin-binding aptamer (TBA) onto polymethacrylate monolith pore surface, and the binding characteristics of the resulting macroporous aptasensor were investigated. Molecular dynamics (MD) simulations of the TBA-thrombin binding indicated enhanced Guanine 4 (G4) structural stability of TBA upon interaction with thrombin in an ionic environment. Fourier-transform infrared spectroscopy and thermogravimetric analyses were used to characterise the available functional groups and thermo-molecular stability of the immobilised polymer generated with Schiff-base activation and immobilisation scheme. The initial degradation temperature of the polymethacrylate stationary support increased with each step of the Schiff-base process: poly(Ethylene glycol Dimethacrylate-co-Glycidyl methacrylate) or poly(EDMA-co-GMA) [196.0 °C (±1.8)]; poly(EDMA-co-GMA)-Ethylenediamine [235.9 °C (±6.1)]; poly(EDMA-co-GMA)-Ethylenediamine-Glutaraldehyde [255.4 °C (±2.7)]; and aptamer-modified monolith [273.7 °C (±2.5)]. These initial temperature increments reflected in the associated endothermic energies were determined with differential scanning calorimetry. The aptameric ligand density obtained after immobilisation was 480 pmol/μL. Increase in pH and ionic concentration affected the surface charge distribution and the binding characteristics of the aptamer-modified disk-monoliths, resulting in the optimum binding pH and ionic concentration of 8.0 and 5 mM Mg2+, respectively. These results are critical in understanding and setting parametric constraints indispensable to develop and enhance the performance of aptasensors.
Membrane-bound fatty acid desaturases perform oxygenated desaturation reactions to insert double bonds within fatty acyl chains in regioselective and stereoselective manners. The Δ9-fatty acid desaturase strictly creates the first double bond between C9 and 10 positions of most saturated substrates. As the three-dimensional structures of the bacterial membrane fatty acid desaturases are not available, relevant information about the enzymes are derived from their amino acid sequences, site-directed mutagenesis and domain swapping in similar membrane-bound desaturases. The cold-tolerantPseudomonassp. AMS8 was found to produce high amount of monounsaturated fatty acids at low temperature. Subsequently, an active Δ9-fatty acid desaturase was isolated and functionally expressed inEscherichia coli. In this paper we report homology modeling and docking studies of a Δ9-fatty acid desaturase from a Cold-tolerantPseudomonassp. AMS8 for the first time to the best of our knowledge. Three dimensional structure of the enzyme was built using MODELLER version 9.18 using a suitable template. The protein model contained the three conserved-histidine residues typical for all membrane-bound desaturase catalytic activity. The structure was subjected to energy minimization and checked for correctness using Ramachandran plots and ERRAT, which showed a good quality model of 91.6 and 65.0%, respectively. The protein model was used to preform MD simulation and docking of palmitic acid using CHARMM36 force field in GROMACS Version 5 and Autodock tool Version 4.2, respectively. The docking simulation with the lowest binding energy, -6.8 kcal/mol had a number of residues in close contact with the docked palmitic acid namely, Ile26, Tyr95, Val179, Gly180, Pro64, Glu203, His34, His206, His71, Arg182, Thr85, Lys98 and His177. Interestingly, among the binding residues are His34, His71 and His206 from the first, second, and third conserved histidine motif, respectively, which constitute the active site of the enzyme. The results obtained are in compliance with thein vivoactivity of the Δ9-fatty acid desaturase on the membrane phospholipids.
Aptamer has been long studied as a substitute of antibodies for many purposes. However, due to the exceeded length of the aptamers obtained in vitro, difficulties arise in its manipulation during its molecular conjugation on the matrix surfaces. Current study focuses on computational improvement for aptamers screening of hepatitis B surface antigen (HBsAg) through optimization of the length sequences obtained from SELEX. Three original aptamers with affinity against HBsAg were truncated into five short hairpin structured aptamers and their affinity against HBsAg was thoroughly studied by molecular docking, molecular dynamics (MD) simulation, and Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) method. The result shows that truncated aptamers binding on HBsAg "a" determinant region are stabilized by the dynamic H-bond formation between the active binding residues and nucleotides. Amino acids residues with the highest hydrogen bonds hydrogen bond interactions with all five aptamers were determined as the active binding residues and further characterized. The computational prediction of complexes binding will include validations through experimental assays in future studies. Current study will improve the current in vitro aptamers by minimizing the aptamer length for its easy manipulation.