Why are some Enzymes Dimers? Flexibility and Catalysis in Thermotoga Maritima Dihydrofolate Reductase

Cooperative Conformational Transitions Underpin the Activation Heat Capacity in the Temperature Dependence of Enzyme Catalysis

Many enzymes display non-Arrhenius behavior with curved Arrhenius plots in the absence of denaturation. There has been significant debate about the origin of this behavior and recently the role of the activation heat capacity (ΔCP⧧) has been widely discussed. If enzyme-catalyzed reactions occur with appreciable negative values of ΔCP⧧ (arising from narrowing of the conformational space along the reaction coordinate), then curved Arrhenius plots are a consequence. To investigate these phenomena in detail, we have collected high precision temperature-rate data over a wide temperature interval for a model glycosidase enzyme MalL, and a series of mutants that change the temperature-dependence of the enzyme-catalyzed rate. We use these data to test a range of models including macromolecular rate theory (MMRT) and an equilibrium model. In addition, we have performed extensive molecular dynamics (MD) simulations to characterize the conformational landscape traversed by MalL in the enzyme-substrate complex and an enzyme-transition state complex. We have crystallized the enzyme in a transition state-like conformation in the absence of a ligand and determined an X-ray crystal structure at very high resolution (1.10 Å). We show (using simulation) that this enzyme-transition state conformation has a more restricted conformational landscape than the wildtype enzyme. We coin the term "transition state-like conformation (TLC)" to apply to this state of the enzyme. Together, these results imply a cooperative conformational transition between an enzyme-substrate conformation (ES) and a transition-state-like conformation (TLC) that precedes the chemical step. We present a two-state model as an extension of MMRT (MMRT-2S) that describes the data along with a convenient approximation with linear temperature dependence of the activation heat capacity (MMRT-1L) that can be used where fewer data points are available. Our model rationalizes disparate behavior seen for MalL and previous results for a thermophilic alcohol dehydrogenase and is consistent with a raft of data for other enzymes. Our model can be used to characterize the conformational changes required for enzyme catalysis and provides insights into the role of cooperative conformational changes in transition state stabilization that are accompanied by changes in heat capacity for the system along the reaction coordinate. TLCs are likely to be of wide importance in understanding the temperature dependence of enzyme activity and other aspects of enzyme catalysis.

Electrostatics as a Guiding Principle in Understanding and Designing Enzymes

Enzyme design faces challenges related to the implementation of the basic principles that govern the catalytic activity in natural enzymes. In this work, we revisit basic electrostatic concepts that have been shown to explain the origin of enzymatic efficiency like preorganization and reorganization. Using magnitudes such as the electrostatic potential and the electric field generated by the protein, we explain how these concepts work in different enzymes and how they can be used to rationalize the consequences of point mutations. We also discuss examples of protein design in which electrostatic effects have been implemented. For the near future, molecular simulations, coupled with the use of machine learning methods, can be used to implement electrostatics as a guiding principle for enzyme designs.

Predicting DNA-binding protein and coronavirus protein flexibility using protein dihedral angle and sequence feature

The flexibility of protein structure is related to various biological processes, such as molecular recognition, allosteric regulation, catalytic activity, and protein stability. At the molecular level, protein dynamics and flexibility are important factors to understand protein function. DNA-binding proteins and Coronavirus proteins are of great concern and relatively unique proteins. However, exploring the flexibility of DNA-binding proteins and Coronavirus proteins through experiments or calculations is a difficult process. Since protein dihedral rotational motion can be used to predict protein structural changes, it provides key information about protein local conformation. Therefore, this paper introduces a method to improve the accuracy of protein flexibility prediction, DihProFle (Prediction of DNA-binding proteins and Coronavirus proteins flexibility introduces the calculated dihedral Angle information). Based on protein dihedral Angle information, protein evolution information, and amino acid physical and chemical properties, DihProFle realizes the prediction of protein flexibility in two cases on DNA-binding proteins and Coronavirus proteins, and assigns flexibility class to each protein sequence position. In this study, compared with the flexible prediction using sequence evolution information, and physicochemical properties of amino acids, the flexible prediction accuracy based on protein dihedral Angle information, sequence evolution information and physicochemical properties of amino acids improved by 2.2% and 3.1% in the nonstrict and strict conditions, respectively. And DihProFle achieves better performance than previous methods for protein flexibility analysis. In addition, we further analyzed the correlation of amino acid properties and protein dihedral angles with residues flexibility. The results show that the charged hydrophilic residues have higher proportion in the flexible region, and the rigid region tends to be in the angular range of the protein dihedral angle (such as the ψ angle of amino acid residues is more flexible than rigid in the range of 91°-120°). Therefore, the results indicate that hydrophilic residues and protein dihedral angle information play an important role in protein flexibility.

The Activation Parameters of a Cold-Adapted Short Chain Dehydrogenase Are Insensitive to Enzyme Oligomerization

The structural principles of enzyme cold adaptation are of fundamental interest both for understanding protein evolution and for biotechnological applications. It has become clear in recent years that structural flexibility plays a major role in tuning enzyme activity at low temperatures, which is reflected by characteristic changes in the thermodynamic activation parameters for psychrophilic enzymes, compared to those of mesophilic and thermophilic ones. Hence, increased flexibility of the enzyme surface has been shown to lead to a lower enthalpy and a more negative entropy of activation, which leads to higher activity in the cold. This immediately raises the question of how enzyme oligomerization affects the temperature dependence of catalysis. Here, we address this issue by computer simulations of the catalytic reaction of a cold-adapted bacterial short chain dehydrogenase in different oligomeric states. Reaction free energy profiles are calculated at different temperatures for the tetrameric, dimeric, and monomeric states of the enzyme, and activation parameters are obtained from the corresponding computational Arrhenius plots. The results show that the activation free energy, enthalpy, and entropy are remarkably insensitive to the oligomeric state, leading to the conclusion that assembly of the subunit interfaces does not compromise cold adaptation, even though the mobilities of interfacial residues are indeed affected.

Module function analysis of a full-length κ-carrageenase from Pseudoalteromonas sp. ZDY3

κ-Carrageenan oligosaccharides with many excellent biological properties could be produced by κ-carrageenases selectively. In this study, based on the encoding gene of full length κ-carrageenase obtained from Pseudoalteromonas sp. ZDY3 and the reported mature secreted κ-carrageenase composed of 275 amino acid residues (N26-T300), CgkPZ_GH16 was expressed in E. coli, but no soluble active protein could be detected. Fortunately, the signal peptide of wild-type κ-carrageenase was recognized, and cleaved in the soluble and folding form in E. coli, the K_m and k_cat values of CgkPZ_SP_GH16 was 1.007 mg/mL and 362.8 s^-1. By molecular dynamics simulations, it was showed that YjdB domain might affect the activity of κ-carrageenase. Due to the absence of mature processing modification system in E. coli, YjdB was remained in recombinant full length κ-carrageenase, and the lost catalytic efficiency of CgkPZ was compensated by expression level and thermal stability. Interestingly, CgkPZ_GH16_YjdB was expressed soluble without the signal peptide, which indicated that YjdB could contribute to the expression and folding of κ-carrageenase. These results provide new insight into the effects of different modules of κ-carrageenase on the expression and properties of enzyme.

Dimerization of PHGDH via the catalytic unit is essential for its enzymatic function

Human D-3-phosphoglycerate dehydrogenase (PHGDH), a key enzyme in de novo serine biosynthesis, is amplified in various cancers and serves as a potential target for anticancer drug development. To facilitate this process, more information is needed on the basic biochemistry of this enzyme. For example, PHGDH was found to form tetramers in solution and the structure of its catalytic unit (sPHGDH) was solved as a dimer. However, how the oligomeric states affect PHGDH enzyme activity remains elusive. We studied the dependence of PHGDH enzymatic activity on its oligomeric states. We found that sPHGDH forms a mixture of monomers and dimers in solution with a dimer dissociation constant of ∼0.58 μM, with the enzyme activity depending on the dimer content. We computationally identified hotspot residues at the sPHGDH dimer interface. Single-point mutants at these sites disrupt dimer formation and abolish enzyme activity. Molecular dynamics simulations showed that dimer formation facilitates substrate binding and maintains the correct conformation required for enzyme catalysis. We further showed that the full-length PHGDH exists as a dynamic mixture of monomers, dimers, and tetramers in solution with enzyme concentration-dependent activity. Mutations that can completely disrupt the sPHGDH dimer show different abilities to interrupt the full-length PHGDH tetramer. Among them, E108A and I121A can also disrupt the oligomeric structures of the full-length PHGDH and abolish its enzyme activity. Our study indicates that disrupting the oligomeric structure of PHGDH serves as a novel strategy for PHGDH drug design and the hotspot residues identified can guide the design process.

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.

Structure-based mechanisms: On the way to apply alcohol dehydrogenases/reductases to organic-aqueous systems

Alcohol dehydrogenases/reductases catalyze enantioselective syntheses of versatile chiral compounds relying on direct hydride transfer from cofactor to substrates, or to an intermediate and then to substrates. Since most of the substrates catalyzed by alcohol dehydrogenases/reductases are insoluble in aqueous solutions, increasing interest has been turning to organic-aqueous systems. However, alcohol dehydrogenases/reductases are normally instable in organic solvents, leading to the unsatisfied enantioselective synthesis efficiency. The behaviors of these enzymes in organic solvents at an atomic level are unclear, thus it is of great importance to understand its structure-based mechanisms in organic-aqueous systems to improve their relative stability. Here, we summarized the accessible structures of alcohol dehydrogenases/reductases in Protein Data Bank crystallized in organic-aqueous systems, and compared the structures of alcohol dehydrogenases/reductases which have different tolerance towards organic solvents. By understanding the catalytic behaviors and mechanisms of these enzymes in organic-aqueous systems, the efficient enantioselective syntheses mediated by alcohol dehydrogenases/reductases and further challenges are also discussed through solvent engineering and enzyme-immobilization in the last decade.

Differences in thermal structural changes and melting between mesophilic and thermophilic dihydrofolate reductase enzymes

The thermal resistance of two homolog enzymes is investigated, with an emphasis on their local stability and flexibility, and on the possible implications regarding their reactivity.