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.

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.

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.

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.