The effect of double (S238F/W159H) mutations on the structure and dynamics of PET degrading enzyme.

Polyethylene terephthalate (PET) is one of the highly produced synthetic polymers worldwide and had acquired attention due to its impact resistance, high clarity, and light weight. PET has become the first choice in making disposable bottles, leading to massive scales of production resulting in very high utilization across various facets of our daily life. Unfortunately, PET accumulates as waste and is highly resistant to biodegradation, thus presenting a serious threat to the ecosystem. Degradation of PET by enzymatic hydrolysis is a promising strategy to depolymerize the PET into its monomers. In recent studies, a plastic-degrading enzyme known as PETase (IsPETase) from the Ideonella sakaiensis has been identified to hydrolyze PET. The wild-type enzyme from Ideonella sp., has been engineered to improve the catalytic activity. While the IsPETase and its variants have been the subject of extensive structural and biochemical studies, the corresponding computational studies to support the improved activity of the mutant enzyme is not fully understood. In this work, we employed all-atom classical molecular dynamics simulations of the wild-type and double mutant IsPETase enzymes to investigate the underlying reason for the improved catalytic activity in the double mutant by means of structure-dynamics-function relationship. Our results show that the engineered mutations reshape the active site structure, volume, and dynamics of the protein loops which is crucial for substrate binding. We also demonstrate that addition of aromatic and hydrogen bond-forming residues near catalytic site improves binding affinity. This work will enable the rational design of mutants for enhanced PET degrading activity.Communicated by Ramaswamy H. Sarma.

Exploring the Conformational Impact of Glycine Receptor TM1-2 Mutations Through Coarse-Grained Analysis and Atomistic Simulations.

Pentameric ligand-gated ion channels (PLGICs) are a family of proteins that convert chemical signals into ion fluxes through cellular membranes. Their structures are highly conserved across all kingdoms from bacteria to eukaryotes. Beyond their classical roles in neurotransmission and neurological disorders, PLGICs have been recently related to cell proliferation and cancer. Here, we focus on the best characterized eukaryotic channel, the glycine receptor (GlyR), to investigate its mutational patterns in genomic-wide tumor screens and compare them with mutations linked to hyperekplexia (HPX), a Mendelian neuromotor disease that disrupts glycinergic currents. Our analysis highlights that cancer mutations significantly accumulate across TM1 and TM2, partially overlapping with HPX changes. Based on 3D-clustering, conservation, and phenotypic data, we select three mutations near the pore, expected to impact GlyR conformation, for further study by molecular dynamics (MD). Using principal components from experimental GlyR ensembles as framework, we explore the motions involved in transitions from the human closed and desensitized structures and how they are perturbed by mutations. Our MD simulations show that WT GlyR spontaneously explores opening and re-sensitization transitions that are significantly impaired by mutations, resulting in receptors with altered permeability and desensitization properties in agreement with HPX functional data.

Temperature-Dependent Kinetic Isotope Effects in R67 Dihydrofolate Reductase from Path-Integral Simulations.

Calculation of temperature-dependent kinetic isotope effects (KIE) in enzymes presents a significant theoretical challenge. Additionally, it is not trivial to identify enzymes with available experimental accurate intrinsic KIEs in a range of temperatures. In the current work, we present a theoretical study of KIEs in the primitive R67 dihydrofolate reductase (DHFR) enzyme and compare with experimental work. The advantage of R67 DHFR is its significantly lower kinetic complexity compared to more evolved DHFR isoforms. We employ mass-perturbation-based path-integral simulations in conjunction with umbrella sampling and a hybrid quantum mechanics-molecular mechanics Hamiltonian. We obtain temperature-dependent KIEs in good agreement with experiments and ascribe the temperature-dependent KIEs primarily to zero-point energy effects. The active site in the primitive enzyme is found to be poorly preorganized, which allows excessive water access to the active site and results in loosely bound reacting ligands.

Modeling the Role of a Flexible Loop and Active Site Side Chains in Hydride Transfer Catalyzed by Glycerol-3-phosphate Dehydrogenase.

Glycerol-3-phosphate dehydrogenase is a biomedically important enzyme that plays a crucial role in lipid biosynthesis. It is activated by a ligand-gated conformational change that is necessary for the enzyme to reach a catalytically competent conformation capable of efficient transition-state stabilization. While the human form (hlGPDH) has been the subject of extensive structural and biochemical studies, corresponding computational studies to support and extend experimental observations have been lacking. We perform here detailed empirical valence bond and Hamiltonian replica exchange molecular dynamics simulations of wild-type hlGPDH and its variants, as well as providing a crystal structure of the binary hlGPDH·NAD R269A variant where the enzyme is present in the open conformation. We estimated the activation free energies for the hydride transfer reaction in wild-type and substituted hlGPDH and investigated the effect of mutations on catalysis from a detailed structural study. In particular, the K120A and R269A variants increase both the volume and solvent exposure of the active site, with concomitant loss of catalytic activity. In addition, the R269 side chain interacts with both the Q295 side chain on the catalytic loop, and the substrate phosphodianion. Our structural data and simulations illustrate the critical role of this side chain in facilitating the closure of hlGPDH into a catalytically competent conformation, through modulating the flexibility of a key catalytic loop (292-LNGQKL-297). This, in turn, rationalizes a tremendous 41,000 fold decrease experimentally in the turnover number, k cat, upon truncating this residue, as loop closure is essential for both correct positioning of key catalytic residues in the active site, as well as sequestering the active site from the solvent. Taken together, our data highlight the importance of this ligand-gated conformational change in catalysis, a feature that can be exploited both for protein engineering and for the design of allosteric inhibitors targeting this biomedically important enzyme.

EnzyDock: Protein-Ligand Docking of Multiple Reactive States along a Reaction Coordinate in Enzymes.

Enzymes play a pivotal role in all biological systems. These biomachines are the most effective catalysts known, dramatically enhancing the rate of reactions by more than 10 orders of magnitude relative to the uncatalyzed reactions in solution. Predicting the correct, mechanistically appropriate binding modes for substrate and product, as well as all reaction intermediates and transition states, along a reaction pathway is immensely challenging and remains an unsolved problem. In the present work, we developed an effective methodology for identifying probable binding modes of multiple ligand states along a reaction coordinate in an enzyme active site. The program is called EnzyDock and is a CHARMM-based multistate consensus docking program that includes a series of protocols to predict the chemically relevant orientation of substrate, reaction intermediates, transition states, product, and inhibitors. EnzyDock is based on simulated annealing molecular dynamics and Monte Carlo sampling and allows ligand, as well as enzyme side-chain and backbone flexibility. The program can employ many user-defined constraints and restraints and classical force field potentials, as well as a range of hybrid quantum mechanics-molecular mechanics potentials. Herein, we apply EnzyDock to several different kinds of problems. First, we study two terpene synthase reactions, namely bornyl diphosphate synthase and the bacterial diterpene synthase CotB2. Second, we use EnzyDock to predict reaction coordinate states in a pair of Diels-Alder reactions in the enzymes spirotetronate AbyU and LepI. Third, we study a couple of racemases: the cofactor-dependent serine racemase and the cofactor independent proline racemase. Finally, we study several cases of covalent docking involving the Michael addition reaction. For all systems we predict binding modes that are consistent with available experimental observations, as well as with theoretical modeling studies from the literature. EnzyDock provides a platform for generating mechanistic insight into enzyme reactions, useful and reliable starting points for in-depth multiscale modeling projects, and rational design of noncovalent and covalent enzyme inhibitors.

Chemical and Biochemical Approaches for the Synthesis of Substituted Dihydroxybutanones and Di- and Tri-Hydroxypentanones.

Polyhydroxylated compounds are building blocks for the synthesis of carbohydrates and other natural products. Their synthesis is mainly achieved by different synthetic versions of aldol-coupling reactions, catalyzed either by organocatalysts, enzymes, or metal-organic catalysts. We have investigated the formation of 1,4-substituted 2,3-dihydroxybutan-1-one derivatives from para- and meta-substituted phenylacetaldehydes by three distinctly different strategies. The first involved a direct aldol reaction with hydroxyacetone, dihydroxyacetone, or 2-hydroxyacetophenone, catalyzed by the cinchona derivative cinchonine. The second was reductive cross-coupling with methyl- or phenylglyoxal promoted by SmI2, resulting in either 5-substituted 3,4-dihydroxypentan-2-ones or 1,4 bis-phenyl-substituted butanones, respectively. Finally, in the third case, aldolase catalysis was employed for synthesis of the corresponding 1,3,4-trihydroxylated pentan-2-one derivatives. The organocatalytic route with cinchonine generated distereomerically enriched syn-products (de = 60-99%), with moderate enantiomeric excesses (ee = 43-56%) but did not produce aldols with either hydroxyacetone or dihydroxyacetone as donor ketones. The SmI2-promoted reductive cross-coupling generated product mixtures with diastereomeric and enantiomeric ratios close to unity. This route allowed for the production of both 1-methyl- and 1-phenyl-substituted 2,3-dihydroxybutanones at yields between 40-60%. Finally, the biocatalytic approach resulted in enantiopure syn-(3 R,4 S) 1,3,4-trihydroxypentan-2-ones.

Evolutionary Effects on Bound Substrate p Ka in Dihydrofolate Reductase.

In the present study, we address the effect of active site structure and dynamics of different dihydrofolate reductase (DHFR) isoforms on the p Ka of the bound substrate 7,8-dihydrofolate, in an attempt to understand possible evolutionary trends. We apply a hybrid QM/MM free energy perturbation method to estimate the p Ka of the N5 position of the bound substrate. We observe a gradual increase in N5 basicity as we move from primitive to more evolved DHFR isoforms. Structural analysis of these isoforms reveals a gradual sequestering of water molecules from the active site in the more evolved enzymes, thereby modulating the local dielectric environment near the substrate. Furthermore, the present study reveals a clear correlation between active site hydration and the N5 p Ka of the substrate. We emphasize the role of the M20 loop in controlling the active site hydration level, via a preorganized active site with a more hydrophobic environment and reduced loop flexibility as evolution progresses from bacterial to the human enzyme.

Towards a comprehensive understanding of the structural dynamics of a bacterial diterpene synthase during catalysis.

Terpenes constitute the largest and structurally most diverse natural product family. Most terpenoids exhibit a stereochemically complex macrocyclic core, which is generated by C-C bond forming of aliphatic oligo-prenyl precursors. This reaction is catalysed by terpene synthases (TPSs), which are capable of chaperoning highly reactive carbocation intermediates through an enzyme-specific reaction. Due to the instability of carbocation intermediates, the proteins' structural dynamics and enzyme:substrate interactions during TPS catalysis remain elusive. Here, we present the structure of the diterpene synthase CotB2, in complex with an in crystallo cyclised abrupt reaction product and a substrate-derived diphosphate. We captured additional snapshots of the reaction to gain an overview of CotB2's catalytic mechanism. To enhance insights into catalysis, structural information is augmented with multiscale molecular dynamic simulations. Our data represent fundamental TPS structure dynamics during catalysis, which ultimately enable rational engineering towards tailored terpene macrocycles that are inaccessible by conventional chemical synthesis.

Effect of Asp122 Mutation on the Hydride Transfer in E. coli DHFR Demonstrates the Goldilocks of Enzyme Flexibility.

Dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate (THF) in the presence of NADPH. The key hydride transfer step in the reaction is facilitated by a combination of enzyme active site preorganization and correlated protein motions in the Michaelis-Menten (E:NADPH:DHF) complex. The present theoretical study employs mutagenesis to examine the relation between structural and functional properties of the enzyme. We mutate Asp122 in Escherichia coli DHFR, which is a conserved amino acid in the DHFR family. The consequent effect of the mutation on enzyme catalysis is examined from an energetic, structural and short-time dynamic perspective. Our investigations suggest that the structural and short-time dynamic perturbations caused by Asp122X mutations (X = Asn, Ser, Ala) are along the reaction coordinate and lower the rate of hydride transfer. Importantly, analysis of the correlated and principle component motions in the enzyme suggest that the mutation alters the coupled motions that are present in the wild-type enzyme. In the case of D122N and D122S, the mutations inhibit coupled motion, whereas in the case of D122A, the mutation enhances coupled motion, although all mutations result in similar rate reduction. These results emphasize a Goldilocks principle of enzyme flexibility, that is, enzymes should neither be too rigid nor too flexible.

The role of the Met20 loop in the hydride transfer in Escherichia coli dihydrofolate reductase.

A key question concerning the catalytic cycle of Escherichia coli dihydrofolate reductase (ecDHFR) is whether the Met20 loop is dynamically coupled to the chemical step during catalysis. A more basic, yet unanswered question is whether the Met20 loop adopts a closed conformation during the chemical hydride transfer step. To examine the most likely conformation of the Met20 loop during the chemical step, we studied the hydride transfer in wild type (WT) ecDHFR using hybrid quantum mechanics-molecular mechanics free energy simulations with the Met20 loop in a closed and disordered conformation. Additionally, we investigated three mutant forms (I14X; X = Val, Ala, Gly) of the enzyme that have increased active site flexibility and donor-acceptor distance dynamics in closed and disordered Met20 loop states. We found that the conformation of the Met20 loop has a dramatic effect on the ordering of active site hydration, although the Met20 loop conformation only has a moderate effect on the hydride transfer rate and donor-acceptor distance dynamics. Finally, we evaluated the pKa of the substrate N5 position in closed and disordered Met20 loop states and found a strong correlation between N5 basicity and the conformation of the Met20 loop.

Functional stability and structural transitions of Kallikrein: spectroscopic and molecular dynamics studies.

Kallikrein, a physiologically vital serine protease, was investigated for its functional and conformational transitions during chemical (organic solvents, Gdn-HCl), thermal, and pH induced denaturation using biochemical and biophysical techniques and molecular dynamics (MD) simulations approach. The enzyme was exceptionally stable in isopropanol and ethanol showing 110% and 75% activity, respectively, after 96 h, showed moderate tolerance in acetonitrile (45% activity after 72 h) and much lower stability in methanol (40% activity after 24 h) (all the solvents [90% v/v]). Far UV CD and fluorescence spectra indicated apparent reduction in compactness of KLKp structure in isopropanol system. MD simulation studies of the enzyme in isopropanol revealed (1) minimal deviation of the structure from native state (2) marginal increase in radius of gyration and solvent accessible surface area (SASA) of the protein and the active site, and (3) loss of density barrier at the active site possibly leading to increased accessibility of substrate to catalytic triad as compared to methanol and acetonitrile. Although kallikrein was structurally stable up to 90 °C as indicated by secondary structure monitoring, it was functionally stable only up to 45 °C, implicating thermolabile active site geometry. In GdnHCl [1.0 M], 75% of the activity of KLKp was retained after incubation for 4 h, indicating its denaturant tolerance. A molten globule-like structure of KLKp formed at pH 1.0 was more thermostable and exhibited interesting structural transitions in organic solvents. The above results provide deeper understanding of functional and structural stability of the serine proteases at molecular level.

Free Energy of Bare and Capped Gold Nanoparticles Permeating through a Lipid Bilayer.

Herein, we study the permeation free energy of bare and octane-thiol-capped gold nanoparticles (AuNPs) translocating through a lipid membrane. To investigate this, we have pulled the bare and capped AuNPs from bulk water to the membrane interior and estimated the free energy cost. The adsorption of the bare AuNP on the bilayer surface is energetically favorable but further loading inside it requires energy. However, the estimated free-energy barrier for loading the capped AuNP into the lipid membrane is much higher compared to bare AuNP. We also demonstrate the details of the permeation process of bare and capped AuNPs. Bare AuNP induces the curvature in the lipid membrane whereas capped AuNP creates an opening in the interacting monolayer and get inserted into the membrane. The insertion of capped AuNP induces a partial unzipping of the lipid bilayer, which results in the ordering of the local lipids interacting with the nanoparticle. However, bare AuNP disrupts the lipid membrane by pushing the lipid molecules inside the membrane. We also analyze pore formation due to the insertion of capped AuNP into the membrane, which results in water molecules penetrating the hydrophobic region.

Nucleoside-2',3'/3',5'-bis(thio)phosphate antioxidants are also capable of disassembly of amyloid beta42-Zn(ii)/Cu(ii) aggregates via Zn(ii)/Cu(ii)-chelation.

Currently, there is an urgent need for biocompatible metal-ion chelators capable of antioxidant activity and disassembly of amyloid beta (Aß)-aggregates as potential therapeutics for Alzheimer's disease (AD). We recently demonstrated the promising antioxidant activity of adenine/guanine 2',3' or 3',5'-bis(thio)phosphate analogues, 2'-dA/G3'5'PO/S and A2'3'PO/S, and their affinity to Zn(ii)-ions. These findings encouraged us to evaluate them as agents for the dissolution of Aß42-Zn(ii)/Cu(ii) aggregates. Specifically, we explored their ability to bind Cu(ii)/Zn(ii)-ions, the geometry and stoichiometry of these complexes, Cu(ii)/Zn(ii)-binding-sites and binding mode, and the ability of these analogues to dissolve Aß42-Zn(ii)/Cu(ii) aggregates, as well as their effect on the secondary structure of those aggregates. Finally, we identified the most promising agents for dissolution of Aß42-Zn(ii)/Cu(ii) aggregates. Specifically, we observed the formation of a 1 : 1 complex between 2'-dG3'5'PO and Cu(ii), involving O4 ligands. Zn(ii) was coordinated by both thiophosphate groups of 2'-dA3'5'PS and A2'3'PS involving O2S2 ligands in a 1 : 1 stoichiometry. A2'3'PS dissolves Aß42-Zn(ii) and Aß42-Cu(ii) aggregates as effectively as, and 2.5-fold more effectively than EDTA, respectively. Furthermore, 2'-dG3'5'PS and A2'3'PS reverted the Aß42-M(ii) structure, back to that of the free Aß42. Finally, cryo-TEM and TEM images confirmed the disassembly of Aß42 and Aß42-M(ii) aggregates by A2'3'PS. Hence, 2'-dG3'5'PS and A2'3'PS may serve as promising scaffolds for new AD therapeutics, acting as both effective antioxidants and agents for solubilization of Aß42-Cu(ii)/Zn(ii) aggregates.

Probing the ATP-induced conformational flexibility of the PcrA helicase protein using molecular dynamics simulation.

Helicases are enzymes that unwind double-stranded DNA (dsDNA) into its single-stranded components. It is important to understand the binding and unbinding of ATP from the active sites of helicases, as this knowledge can be used to elucidate the functionality of helicases during the unwinding of dsDNA. In this work, we investigated the unbinding of ATP and its effect on the active-site residues of the helicase PcrA using molecular dynamic simulations. To mimic the unbinding process of ATP from the active site of the helicase, we simulated the application of an external force that pulls ATP from the active site and computed the free-energy change during this process. We estimated an energy cost of ~85 kJ/mol for the transformation of the helicase from the ATP-bound state (1QHH) to the ATP-free state (1PJR). Unbinding led to conformational changes in the residues of the protein at the active site. Some of the residues at the ATP-binding site were significantly reoriented when the ATP was pulled. We observed a clear competition between reorientation of the residues and energy stabilization by hydrogen bonds between the ATP and active-site residues. We also checked the flexibility of the PcrA protein using a principal component analysis of domain motion. We found that the ATP-free state of the helicase is more flexible than the ATP-bound state.

Self-assembly of phospholipids on flat supports.

The current study deals with the self-assembly of phospholipids on flat supports using the Martini coarse grain model. We reported here the effect of the hydrophilic and hydrophobic nature of the solid supports on the lipid self-assembly. The hydrophilic and hydrophobic supports were modeled on the basis of water droplet simulations. The present work addresses the self-assembly mechanism of lipids on eight different supports with different strengths of hydrophilicity and hydrophobicity. We demonstrated how interplay between the interactions of lipid and water with the support can guide the lipid self-assembly process. Thereafter, we calculated the energetics of the components of the system to quantify the competitions between water and a lipid head-group with hydrophilic supports. Finally, the properties of the self-assembled bilayers were also analyzed and reported here.

Effect of gold nanoparticle on structure and fluidity of lipid membrane.

This paper deals with the effect of different size gold nanoparticles on the fluidity of lipid membrane at different regions of the bilayer. To investigate this, we have considered significantly large bilayer leaflets and incorporated only one nanoparticle each time, which was subjected to all atomistic molecular dynamics simulations. We have observed that, lipid molecules located near to the gold nanoparticle interact directly with it, which results in deformation of lipid structure and slower dynamics of lipid molecules. However, lipid molecules far away from the interaction site of the nanoparticle get perturbed, which gives rise to increase in local ordering of the lipid domains and decrease in fluidity. The bilayer thickness and area per head group in this region also get altered. Similar trend, but with different magnitude is also observed when different size nanoparticle interact with the bilayer.