Correlated conformational fluctuations during enzymatic catalysis: Implications for catalytic rate enhancement

Correlated enzymatic conformational fluctuations are shown to contribute to the rate of enhancement achieved during catalysis. Cytidine deaminase serves as a model system. Crystallographic temperature factor data for this enzyme complexed with substrate analog, transition-state analog, and product are available, thereby establishing a measure of atomic scale conformational fluctuations along the (approximate) reaction coordinate. First, a neural network-based algorithm is used to visualize the decreased conformational fluctuations at the transition state. Second, a dynamic diffusion equation along the reaction coordinate is solved and shows that the flux velocity through the associated enzymatic conformation space is greatest at the transition state. These results suggest (1) that there are both dynamic and energetic restrictions to conformational fluctuations at the transition state, (2) that enzymatic catalysis occurs on a fluctuating potential energy surface, and (3) a form for the potential energy. The Michaelis-Menten equations are modified to describe catalysis on this fluctuating potential energy profile, leading to enhanced catalytic rates when fluctuations along the reaction coordinate are appropriately correlated. This represents a dynamic tuning of the enzyme for maximally effective transformation of the ES complex into EP.     

On the role of the conformational flexibility of the active-site lid on the allosteric kinetics of glucosamine-6-phosphate deaminase

The active site of glucosamine-6-phosphate deaminase from Escherichia coli (GlcN6P deaminase, EC 3.5.99.6) has a complex lid formed by two antiparallel beta-strands connected by a helix-loop segment (158-187). This motif contains Arg172, which is a residue involved in binding the substrate in the active-site, and three residues that are part of the allosteric site, Arg158, Lys160 and Thr161. This dual binding role of the motif forming the lid suggests that it plays a key role in the functional coupling between active and allosteric sites. Previous crystallographic work showed that the temperature coefficients of the active-site lid are very large when the enzyme is in its T allosteric state. These coefficients decrease in the R state, thus suggesting that this motif changes its conformational flexibility as a consequence of the allosteric transition. In order to explore the possible connection between the conformational flexibility of the lid and the function of the deaminase, we constructed the site-directed mutant Phe174-Ala. Phe174 is located at the C-end of the lid helix and its side-chain establishes hydrophobic interactions with the remainder of the enzyme. The crystallographic structure of the T state of Phe174-Ala deaminase, determined at 2.02 A resolution, shows no density for the segment 162-181, which is part of the active-site lid (PDB 1JT9). This mutant form of the enzyme is essentially inactive in the absence of the allosteric activator, N-acetylglucosamine-6-P although it recovers its activity up to the wild-type level in the presence of this ligand. Spectrometric and binding studies show that inactivity is due to the inability of the active-site to bind ligands when the allosteric site is empty. These data indicate that the conformational flexibility of the active-site lid critically alters the binding properties of the active site, and that the occupation of the allosteric site restores the lid conformational flexibility to a functional state.     

Crystal structure determination and dynamic studies of Mycobacterium tuberculosis Cytidine deaminase in complex with products

Cytidine deaminase (CDA) is a key enzyme in the pyrimidine salvage pathway. It is involved in the hydrolytic deamination of cytidine or 2′-deoxycytidine to uridine or 2′-deoxyuridine, respectively. Here we report the crystal structures of Mycobacterium tuberculosis CDA (MtCDA) in complex with uridine (2.4 Å resolution) and deoxyuridine (1.9 Å resolution). Molecular dynamics (MD) simulation was performed to analyze the physically relevant motions involved in the protein–ligand recognition process, showing that structural flexibility of some protein regions are important to product binding. In addition, MD simulations allowed the analysis of the stability of tetrameric MtCDA structure. These findings open-up the possibility to use MtCDA as a target in future studies aiming to the rational design of new inhibitor of MtCDA-catalyzed chemical reaction with potential anti-proliferative activity on cell growth of M. tuberculosis, the major causative agent of tuberculosis.     

Interaction with the CCT chaperonin complex limits APOBEC3A cytidine deaminase cytotoxicity

The APOBEC3 cytidine deaminases are implicated as the cause of a prevalent somatic mutation pattern found in cancer genomes. The APOBEC3 enzymes act as viral restriction factors by mutating viral genomes. Mutation of the cellular genome is presumed to be an off‐target activity of the enzymes, although the regulatory measures for APOBEC3 expression and activity remain undefined. It is therefore difficult to predict circumstances that enable APOBEC3 interaction with cellular DNA that leads to mutagenesis. The APOBEC3A (A3A) enzyme is the most potent deaminase of the family. Using proteomics, we evaluate protein interactors of A3A to identify potential regulators. We find that A3A interacts with the chaperonin‐containing TCP‐1 (CCT) complex, a cellular machine that assists in protein folding and function. Importantly, depletion of CCT results in A3A‐induced DNA damage and cytotoxicity. Evaluation of cancer genomes demonstrates an enrichment of A3A mutational signatures in cancers with silencing mutations in CCT subunit genes. Together, these data suggest that the CCT complex interacts with A3A, and that disruption of CCT function results in increased A3A mutational activity.     

Modulation of human cytidine deaminase by specific aminoacids involved in the intersubunit interactions

An investigation was made of the role exerted by some residues supposed to be involved in the intersubunit interaction and also in the catalytic site of homotetrameric human cytidine deaminase (T-CDA). Attention was focused on Y33, Y60, R68, and F137 residues that are a part of a conserved region in most T-CDAs. Hence, a series of site-directed mutagenesis experiments was set up obtaining seven mutants: Y60G, Y33G, Y33F Y33S, F137A, R68G, and R68Q. Each active purified mutant protein was characterized kinetically, with a series of substrates and inhibitors, and the effect of temperature on enzyme activity and stability was also investigated. Circular dichroism (CD) experiments at different temperatures and in presence of small amounts of sodium dodecyl sulphate (SDS) were performed in all the soluble mutant CDAs. The results obtained by site-directed mutagenesis studies were compared to the crystallographic data of B. subtilis CDA and E. coli CDA and to molecular modeling studies previously performed on human CDA. The mutation of Y60 to glycine produced an enzyme with a more compact quaternary structure with respect to the wild-type; this mutation did not have a dramatic effect on cytidine deamination, but it slightly affected the binding with the substrate. None of the mutant CDAs in Y33 showed enzymatic activity; they existed only as monomers, indicating that this residue, located at the intersubunit interface, may be responsible for the correct folding of human CDA. The insertion of an alanine instead of phenylalanine at position 137 led to a soluble but completely inactive enzyme unable to form a tetramer, suggesting that F137 residue may be important for the assembling of the tetramer and also for the arrangement of the CDA active site. Finally, R68G and R68Q mutations revealed that the presence of the amino group seems to be important for the catalytic process but not for substrate binding, as already shown in B. subtilis CDA. The quaternary structure of R68Q was not affected by the mutation, as shown by the SDS-induced dissociation experiments and CD studies, whereas R68G dissociated very easily in presence of small amounts of SDS. These experiments indicated that in the human CDA, the side chain of arginine 68 involved in the catalytic process in one subunit active site might come from another subunit. The data obtained from these studies confirmed the presence of a complicated set of intersubunit interactions in the active site of human CDA, as shown in other T-CDAs.     

Binding and structural studies of the complexes of type 1 ribosome inactivating protein from Momordica balsamina with cytosine,cytidine, and cytidine diphosphate

The type 1 ribosome inactivating protein from Momordica balsamina (MbRIP1) has been shown to interact with purine bases, adenine and guanine of RNA/DNA. We report here the binding and structural studies of MbRIP1 with a pyrimidine base, cytosine; cytosine containing nucleoside, cytidine; and cytosine containing nucleotide, cytidine diphosphate. All three compounds bound to MbRIP1 at the active site with dissociation constants of 10^?4 M–10^?7 M. As reported earlier, in the structure of native MbRIP1, there are 10 water molecules in the substrate binding site. Upon binding of cytosine to MbRIP1, four water molecules were dislodged from the substrate binding site while five water molecules were dislodged when cytidine bound to MbRIP1. Seven water molecules were dislocated when cytidine diphosphate bound to MbRIP1. This showed that cytidine diphosphate occupied a larger space in the substrate binding site enhancing the buried surface area thus making it a relatively better inhibitor of MbRIP1 as compared to cytosine and cytidine. The key residues involved in the recognition of cytosine, cytidine and cytidine diphosphate were Ile71, Glu85, Tyr111 and Arg163. The orientation of cytosine in the cleft is different from that of adenine or guanine indicating a notable difference in the modes of binding of purine and pyrimidine bases. Since adenine containing nucleosides/nucleotides are suitable substrates, the cytosine containing nucleosides/nucleotides may act as inhibitors.     

Structural basis of hierarchical multiple substates of a protein. II: Monte Carlo simulation of native thermal fluctuations and energy minimization

Conformational fluctuations in a globular protein, bovine pancreatic trypsin inhibitor, in the time range between picoseconds and nanoseconds are studied by a Monte Carlo simulation method. Multiple energy minima are derived from sampled conformations by minimizing their energy. They are distributed in clusters in the conformational space. A hierarchical structure is observed in the simulated dynamics. In the time range between 10(-14) and 10(-10) seconds dynamics is well represented by a superposition of vibrational motions within an energy well with transitions among minima within each cluster. Transitions among clusters take place in the time range of nanoseconds or longer.     

Effect of cytidine on the modification of phospholipid metabolism induced by ischemia

[1-¹⁴C-]Arachidonic acid was injected into the lateral ventricle of the gerbils (meriones unguiculatus) two hours before producing brain ischemia by the bilateral ligation of the carotid arteries. Ten minutes before the carotid ligation a group of animals received an additional intraventricular injection of cold cytidine (2.5 mol/brain). Control animals with and without cytidine, together with the ischemic group, were decapitated directly into liquid nitrogen ten minutes after carotid ligation or sham surgery. Cytidine is able to both stimulate arachidonic acid incorporation into lipids and noticeably correct the release of this acid from polar lipids induced by ischemia. Based on these findings, it is possible to assume that cytidine exerts an effect on the biosynthesis of phosphoglycerides as well as on their catabolic activities.     

Heat shock protein 70 (Hsp70)-stimulated deoxycytidine deaminases from a human lymphoma cell but not the activation-induced cytidine deaminase (AID) from Ramos 6.4 human Burkitt's lymphoma cells

Deoxycytidine deaminase enzyme activity was reduced in lysates of human leukemic THP1 cells 24 h after transfection with siRNA designed to inhibit cell synthesis of heat shock protein 70 (Hsp70)1a and Hsp701b. The cytidine deaminase enzyme activity from the cell lysates was purified from an affinity column which contained bound single-stranded oligodeoxycytidylic acid. Deficient enzyme activity in certain elution fractions from the siRNA-transfected cells was restored by including recombinant HSP 70 in the assays. Enzyme activity in some other fractions was increased after siRNA transfection. Activation-induced cytidine deaminase (AID) is a central factor in the immune response. A more specific assay for AID was used to study the influence of Hsp70 on AID activity. Unlike Hsp70's ability to stimulate certain enzymes of DNA base excision repair and other cytidine deaminases, it had little effect on AID activity in vitro, or was weakly inhibitory.     

A pH-dependent conformational change, rather than the chemical step, appears to be rate-limiting in the hammerhead ribozyme cleavage reaction.

We have investigated the chemical basis for a previously observed 7.8 A conformational change in the hammerhead ribozyme that positions the substrate for in-line attack. We have found that the conformational change can only be observed at or above pH 8.5 (in the presence of Co(2+)) and requires the presence of an ionizable 2'-OH at the cleavage site, and note that this observed apparent pK(a) of 8.5 for the conformational change is within experimental error (+/-0.5) of the previously reported apparent kinetic pK(a) of 8.5 for the hammerhead ribozyme in the presence of Co(2+). We have solved two crystal structures of hammerhead ribozymes having 2'-OCH(3) or 2'-F substitutions at the cleavage site and have found that these will not undergo a conformational change equivalent to that observed for the hammerhead ribozyme having an unmodified attacking nucleophile under otherwise identical conditions. We have also characterized the kinetics of cleavage in the crystal. In addition to verifying that the particular sequence of RNA that we crystallized cleaves faster in the crystal than in solution, we also find that the extent of cleavage in the crystal is complete, unlike in solution where this and most other hammerhead ribozyme substrates are cleaved only to about 70 % completion. The initial cleavage rate in the crystal obeys the expected log-linear relation between cleavage-rate and pH with a slope of 0.7, as has been observed for other hammerhead ribozyme sequences in solution, indicating that in both the crystal and in solution the pH-dependent step is rate-limiting. However, the cleavage rate in the crystal is biphasic, with the most dramatic distinction between initial (slower) and final (faster) phases appearing at pH 6.0. The initial phase corresponds to the pH-dependent cleavage rate observed in solution, but the second, faster phase is roughly pH-independent and closely parallels the cleavage rate observed at pH 8 (0.4/minute). This result is particularly remarkable because it entails that the rapidly cleaving phase at pH 6 is comparable to the cleavage rate for the fastest cleaving hammerhead ribozymes at pH 6. Based upon these observations, we conclude that the pH-dependent conformational change is the rate-determining step under standard conditions for the hammerhead ribozyme self-cleavage reaction, and that an ionizable 2'-proton at cleavage site is required for this conformational change. We further hypothesize that deprotonation of the cleavage-site 2'-oxygen drives this conformational change.     

A helical twist-induced conformational switch activates cleavage in the hammerhead ribozyme

We have captured the structure of a pre-catalytic conformational intermediate of the hammerhead ribozyme using a phosphodiester tether formed between I and Stem II. This phosphodiester tether appears to mimic interactions in the wild-type hammerhead RNA that enable switching between nuclease and ligase activities, both of which are required in the replicative cycles of the satellite RNA viruses from which the hammerhead ribozyme is derived. The structure of this conformational intermediate reveals how the attacking nucleophile is positioned prior to cleavage, and demonstrates how restricting the ability of Stem I to rotate about its helical axis, via interactions with Stem II, can inhibit cleavage. Analogous covalent crosslinking experiments have demonstrated that imposing such restrictions on interhelical movement can change the hammerhead ribozyme from a nuclease to a ligase. Taken together, these results permit us to suggest that switching between ligase and nuclease activity is determined by the helical orientation of Stem I relative to Stem II.     

Universal features of fluctuations in globular proteins

Using data from 2000 non‐homologous protein crystal structures, we show that the distribution of residue B factors of proteins collapses onto a single master curve. We show by maximum entropy arguments that this curve is a Gamma function whose order and dispersion are obtained from experimental data. The distribution for any given specific protein can be generated from the master curve by a linear transformation. Any perturbation of the B factor distribution of a protein, imposed at constant energy, causes a decrease in the entropy of the protein relative to that of the reference state. Proteins 2016; 84:721–725. © 2016 Wiley Periodicals, Inc.     

Symmetry-free cryo-EM structures of the chaperonin TRiC along its ATPase-driven conformational cycle

The eukaryotic group II chaperonin TRiC/CCT is a 16-subunit complex with eight distinct but similar subunits arranged in two stacked rings. Substrate folding inside the central chamber is triggered by ATP hydrolysis. We present five cryo-EM structures of TRiC in apo and nucleotide-induced states without imposing symmetry during the 3D reconstruction. These structures reveal the intra- and inter-ring subunit interaction pattern changes during the ATPase cycle. In the apo state, the subunit arrangement in each ring is highly asymmetric, whereas all nucleotide-containing states tend to be more symmetrical. We identify and structurally characterize an one-ring closed intermediate induced by ATP hydrolysis wherein the closed TRiC ring exhibits an observable chamber expansion. This likely represents the physiological substrate folding state. Our structural results suggest mechanisms for inter-ring-negative cooperativity, intra-ring-positive cooperativity, and protein-folding chamber closure of TRiC. Intriguingly, these mechanisms are different from other group I and II chaperonins despite their similar architecture.     

Selection of the Best Blood Compartment to Measure Cytidine Deaminase Activity to Stratify for Optimal Gemcitabine or Cytarabine Treatment

Cytidine deaminase (CDA) plays a crucial role in the degradation of cytidine analogs, such as gemcitabine and cytarabine. Several studies showed that a low CDA activity is associated with more toxicity but a higher efficacy, while a high activity will lead to a lower efficacy but less toxicity. A stratified dosing strategy based on the relative CDA activity would increase efficiency. In order to predict these events, a reliable measurement of CDA with a validated method is crucial. We aimed to determine which phenotype assay would be most suitable; a spectrophotometric assay using cytidine as a substrate, or an HPLC assay using gemcitabine as a substrate. In serum and whole blood of 26 volunteers, both assays showed an excellent correlation (R > 0.999), but not in plasma nor in red blood cells. Moreover, there was no difference between males and females.

In conclusion, the spectrophotometric assay seems the most simple and cost-effective test. It should be performed in serum, while it should be normalized on protein content as measured by the Bicinchoninic Acid.     

Length of the hypermutation motif DGYW/WRCH in the focus of statistical limits. Implications for a double-motif or extended motif recognition models

The motif DGYW/WRCH (Mh) and its frequently discussed simplified derivative GYW/WRC (Mhs) are involved in immunoglobulin (Ig) hypermutation. Both these motifs appear to be markedly shorter than the corresponding conventionally predicted minima of valid sequence lengths (MVSL). The same conclusion concerning both Mh and Mhs can also be obtained in the combined case including a less strict semi-empirically defined w-value and one nucleotide length tolerance related to MVSL. Such disagreement indicates considerably low information content in Mh and Mhs when evaluating these motifs as alphabetical structures (words). This fact raises a question of actually recognized structures (presumably longer than Mh and Mhs). Interestingly, both Mh and Mhs dimers or pairs of closely located Mh or Mhs achieve confirmation of length validity in the case of w=0.05, suggesting thus double-motif recognition as one of statistically consistent explanations. This possibility is also in agreement with the results of our model sequence study of mRNA derived from variable Ig gene sequences (rIgV) with respect to the most frequently occurring structures formed by motif overlaps in all model sequence sets. On the other hand, additional superior occurrence of motif pairs at a structurally important distance of a single DNA thread was found in the conserved domain (cd00099) related sequences of Elasmobranchii origin and less markedly in the corresponding human rIgV, but not in a randomly selected human subset of rIgV. The data are discussed with respect to statistical evaluation and structural properties of hypermutation motifs or the competent enzyme, i.e. activation-induced cytidine deaminase.     

Chromosomal location, cloning and nucleotide sequence of the Bacillus subtilis cdd gene encoding cytidine/deoxycytidine deaminase

Summary The Bacillus subtilis cdd gene encoding cytidine/2-deoxycytidine deaminase has been located by transduction at approximately 225 degrees on the chromosome, and the gene order rpC-lys-cdd-aroD was established. The gene was isolated from a library of B. subtilis DNA cloned in D69 by complementation of an Escherichia coli cdd mutation. Minicell experiments revealed a molecular mass of 14000 dalton for the cytidine deaminase subunit encoded by the cloned DNA fragment. The molecular weight of the native enzyme was determined to be 58000, suggesting that it consists of four identical subunits. The nucleotide sequence of 1170 bp, including the cdd gene, was determined. An open reading frame encoding a polypeptide with a calculated molecular mass of 14800 dalton was deduced to be the coding region for cdd. The deduced amino acid composition of the 136-amino acid-long subunit shows that it contains six cysteine residues. A computer search in the GenBank DNA sequence library revealed that the 476 bp HindIII fragment containing the putative promoter region and the first ten codons of cdd is identical to the P43 promoter-containing fragment previously isolated by Wang and Doi (1984). They showed that the fragment contained overlapping promoters transcribed by B. subtilis 43 and 37 RNA polymerase holoenzymes during growth and stationary phase.Abbreviations SDS sodium dodecyl sulphate - Ap ampicillin resistance - Tet^r tetracycline resistance - Km^r kanamycin resistance     

Role of conformational fluctuations in the enzymatic reaction of HIV-1 protease

The emergence of compensatory drug-resistant mutations in HIV-1 protease challenges the common view of the reaction mechanism of this enzyme. Here, we address this issue by performing classical and ab initio molecular dynamics simulations (MD) on a complex between the enzyme and a peptide substrate. The classical MD calculation reveals large-scale protein motions involving the flaps and the cantilever. These motions modulate the conformational properties of the substrate at the cleavage site. The ab initio calculations show in turn that substrate motion modulates the activation free energy barrier of the enzymatic reaction dramatically. Thus, the catalytic power of the enzyme does not arise from the presence of a pre-organized active site but from the protein mechanical fluctuations. The implications of this finding for the emergence of drug-resistance are discussed.     

Exploring the role of large conformational changes in the fidelity of DNA polymerase beta

The relationships between the conformational landscape, nucleotide insertion catalysis and fidelity of DNA polymerase beta are explored by means of computational simulations. The simulations indicate that the transition states for incorporation of right (R) and wrong (W) nucleotides reside in substantially different protein conformations. The protein conformational changes that reproduce the experimentally observed fidelity are significantly larger than the small rearrangements that usually accompany motions from the reactant state to the transition state in common enzymatic reactions. Once substrate binding has occurred, different constraints imposed on the transition states for insertion of R and W nucleotides render it highly unlikely that both transition states can occur in the same closed structure, because the predicted fidelity would then be many orders of magnitude too large. Since the conformational changes reduce the transition state energy of W incorporation drastically they decrease fidelity rather than increase it. Overall, a better agreement with experimental data is attained when the R is incorporated through a transition state in a closed conformation and W is incorporated through a transition state in one or perhaps several partially open conformations. The generation of free energy surfaces for R and W also allow us to analyze proposals about the relationship between induced fit and fidelity.     

Structural snapshots of MTA/AdoHcy nucleosidase along the reaction coordinate provide insights into enzyme and nucleoside flexibility during catalysis

MTA/AdoHcy nucleosidase (MTAN) irreversibly hydrolyzes the N9-C1' bond in the nucleosides, 5'-methylthioadenosine (MTA) and S-adenosylhomocysteine (AdoHcy) to form adenine and the corresponding thioribose. MTAN plays a vital role in metabolic pathways involving methionine recycling, biological methylation, polyamine biosynthesis, and quorum sensing. Crystal structures of a wild-type (WT) MTAN complexed with glycerol, and mutant-enzyme and mutant-product complexes have been determined at 2.0A, 2.0A, and 2.1A resolution, respectively. The WT MTAN-glycerol structure provides a purine-free model and in combination with the previously solved thioribose-free MTAN-ADE structure, we now have separate apo structures for both MTAN binding subsites. The purine and thioribose-free states reveal an extensive enzyme-immobilized water network in their respective binding subsites. The Asp197Asn MTAN-MTA and Glu12Gln MTAN-MTR.ADE structures are the first enzyme-substrate and enzyme-product complexes reported for MTAN, respectively. These structures provide representative snapshots along the reaction coordinate and allow insight into the conformational changes of the enzyme and the nucleoside substrate. A "catalytic movie" detailing substrate binding, catalysis, and product release is presented.