Probing the Effect of Glucose on the Activity and Stability of β-Glucosidase: An All-Atom Molecular Dynamics Simulation Investigation

β-Glucosidase (EC 3.2.1.21) plays an essential role in the removal of glycosyl residues from disaccharide cellobiose to produce glucose during the hydrolysis of lignocellulosic biomass. Although there exist a few β-glucosidase that are tolerant to large concentrations of glucose, these enzymes are typically prone to glucose inhibition. Understanding the basis of this inhibition is important for the production of cheaper biofuels from lignocellulose. In this study, all-atom molecular dynamics simulation at different temperatures and glucose concentrations was used to understand the molecular basis of glucose inhibition of GH1 β-glucosidase (B8CYA8) from Halothermothrix orenii. Our results show that glucose induces a broadening of the active site tunnel through residues lining the tunnel and facilitates the accumulation of glucose. In particular, we observed that glucose accumulates at the tunnel entrance and near the catalytic sites to block substrate accessibility and inhibit enzyme activity. The reduction of enzyme activity was also confirmed experimentally through specific activity measurements in the presence of 0-2.5 M glucose. We also show that the increase in glucose concentrations leads to a decrease in the number of water molecules inside the tunnel to affect substrate hydrolysis. Overall, the results help in understanding the role of residues along the active site tunnel for the engineering of glucose-tolerant β-glucosidase.

Effect of double mutations T790M/L858R on conformation and drug-resistant mechanism of epidermal growth factor receptor explored by molecular dynamics simulations

Epidermal growth factor receptor (EGFR) is one of the most promising targets for the treatment of cancers. Double mutations T790M/L858R lead to different degrees of drug resistance toward inhibitors. In this study, molecular dynamics (MD) simulations followed by principal component analysis are performed to study the conformational changes of EGFR induced by T790M/L858R. The results suggest that T790M/L858R cause obvious disturbance of the structural stability of EGFR. Molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) and residue-based free energy decomposition methods are integrated to explore the drug-resistant mechanism of T790M/L858R toward inhibitors. The results show that the decrease in van der Waals interactions of inhibitors with the mutated EFGR relative to the wild-type (WT) one is the main force inducing drug resistance of T790M/L858R toward inhibitors TAK-285, while drug resistance toward W2P and HKI-272 is dominated by the decrease in van der Waals interactions and the increase in polar interactions. We expect that the information obtained from this study can aid rational design of effective drugs to relieve drug resistance of EGFR induced by T790M/L858R.

Different Behaviors of a Substrate in P450 Decarboxylase and Hydroxylase Reveal Reactivity-Enabling Actors

Biological routes to the production of fuels from renewable feedstocks hold significant promise in our efforts towards a sustainable future. The fatty acid decarboxylase enzyme (OleT_JE) is a cytochrome P450 enzyme that converts long and medium chain fatty acids to terminal alkenes and shares significant similarities in terms of structure, substrate scope and mechanism with the hydroxylase cytochrome P450 (P450_BSβ). Recent reports have demonstrated that catalytic pathways in these enzymes bifurcate when the heme is in its iron-hydroxo (compound II) state. In spite of significant similarities, the fundamental underpinnings of their different characteristic wild-type reactivities remain ambiguous. Here, we develop point charges, modified parameters and report molecular simulations of this crucial intermediate step. Water occupancies and substrate mobility at the active site are observed to be vital differentiating aspects between the two enzymes in the compound II state and corroborate recent experimental hypotheses. Apart from increased substrate mobility in the hydroxylase, which could have implications for enabling the rebound mechanism for hydroxylation, OleT_JE is characterized by much stronger binding of the substrate carboxylate group to the active site arginine, implicating it as an important enabling actor for decarboxylation.

In silico insights into the solvation characteristics of the ionic liquid 1-methyltriethoxy-3-ethylimidazolium acetate for cellulosic biomass

Lignocellulosic biomass is a domestically grown, sustainable, and potentially carbon-neutral feedstock for the production of liquid fuels and other value added chemicals. This underutilized renewable feedstock has the potential to alleviate some of the current socio-economic dependence on foreign petroleum supplies while stimulating rural economies. Unfortunately, the potential of biomass has largely been underdeveloped due to the recalcitrant nature of lignocellulosic materials. Task-specific ionic liquids (ILs) have shown considerable promise as an alternative non-aqueous solvent for solvation and deconstruction of lignocellulose in the presence of metal chloride catalyst or enzymes. Recently it has been hypothesized that adding oxygen atoms to the tail of an imidazolium cation would alleviate some of the negative characteristics of the ILs by increasing mass transport properties, and decreasing IL deactivation of enzymes, while at the same time retaining favorable solvation characteristics for lignocellulose. Reported here are fully atomistic molecular dynamic simulations of 1-methyltriethoxy-3-ethylimidazolium acetate ([Me-(OEt)3-Et-IM(+)] [OAc(-)]) that elucidate promising molecular-level details pertaining to the solvation characteristics of model compounds of cellulose, and IL-induced side-chain and ring puckering conformations. It is found that the anion interactions with the saccharide induce alternate ring puckering conformations from those seen in aqueous environments (i.e.(1)C4), while the cation interactions are found to influence the conformation of the ω dihedral. These perturbations in saccharide structures are discussed in the context of their contribution to the disruption of hydrogen bonding in cellulosic architecture and their role in solvation.