Structural and mutagenetic analyses of a 1,3–1,4-β-glucanase from Paecilomyces thermophila

The thermostable 1,3–1,4-β-glucanase PtLic16A from the fungus Paecilomyces thermophila catalyzes stringent hydrolysis of barley β-glucan and lichenan with an outstanding efficiency and has great potential for broad industrial applications. Here, we report the crystal structures of PtLic16A and an inactive mutant E113A in ligand-free form and in complex with the ligands cellobiose, cellotetraose and glucotriose at 1.80 Å to 2.25 Å resolution. PtLic16A adopts a typical β-jellyroll fold with a curved surface and the concave face forms an extended ligand binding cleft. These structures suggest that PtLic16A might carry out the hydrolysis via retaining mechanism with E113 and E118 serving as the nucleophile and general acid/base, respectively. Interestingly, in the structure of E113A/1,3–1,4-β-glucotriose complex, the sugar bound to the − 1 subsite adopts an intermediate-like (α-anomeric) configuration. By combining all crystal structures solved here, a comprehensive binding mode for a substrate is proposed. These findings not only help understand the 1,3–1,4-β-glucanase catalytic mechanism but also provide a basis for further enzymatic engineering.     

Computational study of IAG-nucleoside hydrolase: determination of the preferred ground state conformation and the role of active site residues

The mechanism of action of inosine-adenosine-guanosine nucleoside hydrolase (IAG-NH) has been investigated by long-term molecular dynamics (MD) simulation in TIP3P water using stochastic boundary conditions. Special attention has been given to the role of leaving group pocket residues and conformation of the bound substrate at the active site of IAG-NH. We also describe the positioning of the residues of an important flexible loop at the active site, which was previously unobservable by X-ray crystallography due to high B-factors. Five MD simulations have been performed with the Enzyme x Substrate complexes: Enzyme x anti-Adenosine with Asp40-COOH [E(40H) x Ade(a)], Enzyme x anti-Adenosine with Asp40-COO- [E(40) x Ade(a)], Enzyme x syn-Adenosine with Asp40-COOH [E(40H) x Ade(s)], Enzyme x syn-Adenosine with Asp40-COO- [E(40) x Ade(s)], and Enzyme x anti-Inosine with Asp40-COO- [E(40) x Ino(a)]. Overall, the structures generated from the MD simulation of E(40H) x Ade(s) preserve the catalytically important hydrogen bonds as well as electrostatic and hydrophobic interactions to provide a plausible catalytic structure. When deprotonated Asp40 (Asp4-COO-) is present, the active site is open to water solvent which interferes with the base stacking between Trp83 and nucleobase. A calculation using Poisson-Boltzmann equation module supports that Asp40 indeed has an elevated pK(app). Solvent accessible surface area (SASA) calculations on all the five MD structures shows that systems with protonated Asp40, namely, E(40H) x Ade(a) and E(40H) x Ade(s), have zero SASA. It is found that a water molecule is hydrogen-bonded to the N7 of the nucleobase and is probably the essential general acid to protonate the departing nucleobase anion. The N7-bonded water is in turn hydrogen-bonded to waters in a channel, held in place by the residues of the flexible loop, Tyr257, His247, and Cys245. Using normal-mode analysis with elastic network model, we find that the flexible loop explores a conformational space much larger than in the MD trajectory, leading to a "gating"-like motion with respect to the active site.     

Structure and function of a novel purine specific nucleoside hydrolase from Trypanosoma vivax

The purine salvage pathway of parasitic protozoa is currently considered as a target for drug development because these organisms cannot synthesize purines de novo. Insight into the structure and mechanism of the involved enzymes can aid in the development of potent inhibitors, leading to new curative drugs. Nucleoside hydrolases are key enzymes in the purine salvage pathway of Trypanosomatidae, and they are especially attractive because they have no equivalent in mammalian cells. We cloned, expressed and purified a nucleoside hydrolase from Trypanosoma vivax. The substrate activity profile establishes the enzyme to be a member of the inosine-adenosine-guanosine-preferring nucleoside hydrolases (IAG-NH). We solved the crystal structure of the enzyme at 1.6 A resolution using MAD techniques. The complex of the enzyme with the substrate analogue 3-deaza-adenosine is presented. These are the first structures of an IAG-NH reported in the literature. The T. vivax IAG-NH is a homodimer, with each subunit consisting of ten beta-strands, 12 alpha-helices and three small 3(10)-helices. Six of the eight strands of the central beta-sheet form a motif resembling the Rossmann fold. Superposition of the active sites of this IAG-NH and the inosine-uridine-preferring nucleoside hydrolase (IU-NH) of Crithidia fasciculata shows the molecular basis of the different substrate specificity distinguishing these two classes of nucleoside hydrolases. An "aromatic stacking network" in the active site of the IAG-NH, absent from the IU-NH, imposes the purine specificity. Asp10 is the proposed general base in the reaction mechanism, abstracting a proton from a nucleophilic water molecule. Asp40 (replaced by Asn39 in the IU-NH) is positioned appropriately to act as a general acid and to protonate the purine leaving group. The second general acid, needed for full enzymatic activity, is probably part of a flexible loop located in the vicinity of the active site.     

Molecular mechanism of ADP-ribose hydrolysis by human NUDT5 from structural and kinetic studies

Human NUDT5 (hNUDT5) is an ADP-ribose (ADPR) pyrophosphatase (ADPRase) that plays important roles in controlling the intracellular levels of ADPR and preventing non-enzymatic ADP-ribosylation of proteins by hydrolyzing ADPR to AMP and ribose 5′-phosphate. We report the crystal structure of hNUDT5 in complex with a non-hydrolyzable ADPR analogue, α,β-methyleneadenosine diphosphoribose, and three Mg^2 + ions representing the transition state of the enzyme during catalysis. Analysis of this structure and comparison with previously reported hNUDT5 structures identify key residues involved in substrate binding and catalysis. In the transition-state structure, three metal ions are bound at the active site and are coordinated by surrounding residues and water molecules. A conserved water molecule is at an ideal position for nucleophilic attack on the α-phosphate of ADPR. The side chain of Glu166 on loop L9 changes its conformation to interact with the conserved water molecule compared with that in the substrate-bound structure and appears to function as a catalytic base. Mutagenesis and kinetic studies show that Trp28 and Trp46 are important for the substrate binding; Arg51 is involved in both the substrate binding and the catalysis; and Glu112 and Glu116 of the Nudix motif, Glu166 on loop L9, and Arg111 are critical for the catalysis. The structural and biochemical data together reveal the molecular basis of the catalytic mechanism of ADPR hydrolysis by hNUDT5. Specifically, Glu166 functions as a catalytic base to deprotonate a conserved water molecule that acts as a nucleophile to attack the α-phosphate of ADPR, and three Mg^2 + ions are involved in the activation of the nucleophile and the binding of the substrate. Structural comparison of different ADPRases also suggests that most dimeric ADPRases may share a similar catalytic mechanism of ADPR hydrolysis.     

Substrate binding to a GH131 β-glucanase catalytic domain from Podospora anserina

β-Glucanases have been utilized widely in industry to treat various carbohydrate-containing materials. Recently, the Podospora anserina β-glucanase 131A (PaGluc131A) was identified and classified to a new glycoside hydrolases GH131 family. It shows exo-β-1,3/exo-β-1,6 and endo-β-1,4 glucanase activities with a broad substrate specificity for laminarin, curdlan, pachyman, lichenan, pustulan, and cellulosic derivatives. Here we report the crystal structures of the PaGluc131A catalytic domain with or without ligand (cellotriose) at 1.8 Å resolution. The cellotriose was clearly observed to occupy the +1 to +3 subsites in substrate binding cleft. The broadened substrate binding groove may explain the diverse substrate specificity. Based on our crystal structures, the GH131 family enzyme is likely to carry out the hydrolysis through an inverting catalytic mechanism, in which E99 and E139 are supposed to serve as the general base and general acid.     

ATP-triggered ADP release from the asymmetric chaperonin GroEL/GroES/ADP7 is not the rate-limiting step of the GroEL/GroES reaction cycle

The GroEL/GroES protein folding chamber is formed and dissociated by ATP binding and hydrolysis. ATP hydrolysis in the GroES-bound (cis) ring gates entry of ATP into the opposite unoccupied trans ring, which allosterically ejects cis ligands. While earlier studies suggested that hydrolysis of cis ATP is the rate-limiting step of the cycle (t_½ ∼ 10 s), a recent study suggested that ADP release from the cis ring may be rate-limiting (t_½ ∼ 15-20 s). Here we have measured ADP release using a coupled enzyme assay and observed a t_½ for release of ?4-5 s, indicating that this is not the rate-limiting step of the reaction cycle.     

Predictive modelling of topology and loop variations in dimeric DNA quadruplex structures

We have used a combination of simulated annealing (SA), molecular dynamics (MD) and locally enhanced sampling (LES) methods in order to predict the favourable topologies and loop conformations of dimeric DNA quadruplexes with T2 or T3 loops. This follows on from our previous MD simulation studies on the influence of loop lengths on the topology of intramolecular quadruplex structures [P. Hazel et al. (2004) J. Am. Chem. Soc., 126, 16405–16415], which provided results consistent with biophysical data. The recent crystal structures of d(G4T3G4)2 and d(G4BrUT2G4) (P. Hazel et al. (2006) J. Am. Chem. Soc., in press) and the NMR-determined topology of d(TG4T2G4T)2 [A.T. Phan et al. (2004) J. Mol. Biol., 338, 93–102] have been used in the present study for comparison with simulation results. These together with MM-PBSA free-energy calculations indicate that lateral T3 loops are favoured over diagonal loops, in accordance with the experimental structures; however, distinct loop conformations have been predicted to be favoured compared to those found experimentally. Several lateral and diagonal loop conformations have been found to be similar in energy. The simulations suggest an explanation for the distinct patterns of observed dimer topology for sequences with T3 and T2 loops, which depend on the loop lengths, rather than only on G-quartet stability.     

Mutagenesis and subsite mapping underpin the importance for substrate specificity of the aglycon subsites of glycoside hydrolase family 11 xylanases

Glycoside hydrolase family (GH) 11 xylanase A from Bacillus subtilis (BsXynA) was subjected to site-directed mutagenesis to probe the role of aglycon active site residues with regard to activity, binding of decorated substrates and hydrolysis product profile. Targets were those amino acids identified to be important by 3D structure analysis of BsXynA in complex with substrate bound in the glycon subsites and the + 1 aglycon subsite. Several aromatic residues in the aglycon subsites that make strong substrate–protein interactions and that are indispensable for enzyme activity, were also important for the specificity of the xylanase. In the + 2 subsite of BsXynA, Tyr65 and Trp129 were identified as residues that are involved in the binding of decorated substrates. Most interestingly, replacement of Tyr88 by Ala in the + 3 subsite created an enzyme able to produce a wider variety of hydrolysis products than wild type BsXynA. The contribution of the + 3 subsite to the substrate specificity of BsXynA was established more in detail by mapping the enzyme binding site of the wild type xylanase and mutant Y88A with labelled xylo-oligosaccharides. Also, the length of the cord – a long loop flanking the aglycon subsites of GH11 xylanases – proved to impact the hydrolytic action of BsXynA. The aglycon side of the active site cleft of BsXynA, therefore, offers great potential for engineering and design of xylanases with a desired specificity.     

Crystal Structure of an Essential Enzyme in Seed Starch Degradation: Barley Limit Dextrinase in Complex with Cyclodextrins

Barley limit dextrinase [Hordeum vulgare limit dextrinase (HvLD)] catalyzes the hydrolysis of α-1,6 glucosidic linkages in limit dextrins. This activity plays a role in starch degradation during germination and presumably in starch biosynthesis during grain filling. The crystal structures of HvLD in complex with the competitive inhibitors α-cyclodextrin (CD) and β-CD are solved and refined to 2.5 Å and 2.1 Å, respectively, and are the first structures of a limit dextrinase. HvLD belongs to glycoside hydrolase 13 family and is composed of four domains: an immunoglobulin-like N-terminal eight-stranded β-sandwich domain, a six-stranded β-sandwich domain belonging to the carbohydrate binding module 48 family, a catalytic (β/α)₈-like barrel domain that lacks α-helix 5, and a C-terminal eight-stranded β-sandwich domain of unknown function. The CDs are bound at the active site occupying carbohydrate binding subsites + 1 and + 2. A glycerol and three water molecules mimic a glucose residue at subsite − 1, thereby identifying residues involved in catalysis. The bulky Met440, a unique residue at its position among α-1,6 acting enzymes, obstructs subsite − 4. The steric hindrance observed is proposed to affect substrate specificity and to cause a low activity of HvLD towards amylopectin. An extended loop (Asp513-Asn520) between β5 and β6 of the catalytic domain also seems to influence substrate specificity and to give HvLD a higher affinity for α-CD than pullulanases. The crystal structures additionally provide new insight into cation sites and the concerted action of the battery of hydrolytic enzymes in starch degradation.     

Molecular dynamics simulations of the d(CCAACGTTGG)(2) decamer: influence of the crystal environment

Molecular dynamics (MD) simulations of the DNA duplex d(CCAACGTTGG)(2) were used to study the relationship between DNA sequence and structure. Two crystal simulations were carried out; one consisted of one unit cell containing two duplexes, and the other of two unit cells containing four duplexes. Two solution simulations were also carried out, one starting from canonical B-DNA and the other starting from the crystal structure. For many helicoidal parameters, the results from the crystal and solution simulations were essentially identical. However, for other parameters, in particular, alpha, gamma, delta, (epsilon - zeta), phase, and helical twist, differences between crystal and solution simulations were apparent. Notably, during crystal simulations, values of helical twist remained comparable to those in the crystal structure, to include the sequence-dependent differences among base steps, in which values ranged from 20 degrees to 50 degrees per base step. However, in the solution simulations, not only did the average values of helical twist decrease to approximately 30 degrees per base step, but every base step was approximately 30 degrees, suggesting that the sequence-dependent information may be lost. This study reveals that MD simulations of the crystal environment complement solution simulations in validating the applicability of MD to the analysis of DNA structure.     

Recognition of nucleoside monophosphate substrates by Haemophilus influenzae class C acid phosphatase

The e (P4) phosphatase from Haemophilus influenzae functions in a vestigial NAD⁺ utilization pathway by dephosphorylating nicotinamide mononucleotide to nicotinamide riboside. P4 is also the prototype of class C acid phosphatases (CCAPs), which are nonspecific 5′,3′-nucleotidases localized to the bacterial outer membrane. To understand substrate recognition by P4 and other class C phosphatases, we have determined the crystal structures of a substrate-trapping mutant P4 enzyme complexed with nicotinamide mononucleotide, 5′-AMP, 3′-AMP, and 2′-AMP. The structures reveal an anchor-shaped substrate-binding cavity comprising a conserved hydrophobic box that clamps the nucleotide base, a buried phosphoryl binding site, and three solvent-filled pockets that contact the ribose and the hydrogen-bonding edge of the base. The span between the hydrophobic box and the phosphoryl site is optimal for recognizing nucleoside monophosphates, explaining the general preference for this class of substrate. The base makes no hydrogen bonds with the enzyme, consistent with an observed lack of base specificity. Two solvent-filled pockets flanking the ribose are key to the dual recognition of 5′-nucleotides and 3′-nucleotides. These pockets minimize the enzyme's direct interactions with the ribose and provide sufficient space to accommodate 5′ substrates in an anti conformation and 3′ substrates in a syn conformation. Finally, the structures suggest that class B acid phosphatases and CCAPs share a common strategy for nucleotide recognition.     

Crystal structure analysis,covalent docking,and molecular dynamics calculations reveal a conformational switch in PhaZ7 PHB depolymerase

An open and a closed conformation of a surface loop in PhaZ7 extracellular poly(3‐hydroxybutyrate) depolymerase were identified in two high‐resolution crystal structures of a PhaZ7 Y105E mutant. Molecular dynamics (MD) simulations revealed high root mean square fluctuations (RMSF) of the 281–295 loop, in particular at residue Asp289 (RMSF 7.62 Å). Covalent docking between a 3‐hydroxybutyric acid trimer and the catalytic residue Ser136 showed that the binding energy of the substrate is significantly more favorable in the open loop conformation compared to that in the closed loop conformation. MD simulations with the substrate covalently bound depicted 1 Å RMSF higher values for the residues 281–295 in comparison to the apo (substrate‐free) form. In addition, the presence of the substrate in the active site enhanced the ability of the loop to adopt a closed form. Taken together, the analysis suggests that the flexible loop 281–295 of PhaZ7 depolymerase can act as a lid domain to control substrate access to the active site of the enzyme. Proteins 2017; 85:1351–1361. © 2017 Wiley Periodicals, Inc.     

Crystal structures of a type-1 ribosome inactivating protein from Momordica balsamina in the bound and unbound states

The ribosome inactivating proteins (RIPs) of type 1 are plant toxins that eliminate adenine base selectively from the single stranded loop of rRNA. We report six crystal structures, type 1 RIP from Momordica balsamina (A), three in complexed states with ribose (B), guanine (C) and adenine (D) and two structures of MbRIP-1 when crystallized with adenosine triphosphate (ATP) (E) and 2′-deoxyadenosine triphosphate (2′-dATP) (F). These were determined at 1.67 Å, 1.60 Å, 2.20 Å, 1.70 Å, 2.07 Å and 1.90 Å resolutions respectively. The structures contained, (A) unbound protein molecule, (B) one protein molecule and one ribose sugar, (C) one protein molecule and one guanine base, (D) one protein molecule and one adenine base, (E) one protein molecule and one ATP-product adenine molecule and (F) one protein molecule and one 2′-dATP-product adenine molecule. Three distinct conformations of the side chain of Tyr70 were observed with (i) χ¹ = − 66°and χ² = 165° in structures (A) and (B); (ii) χ¹ = − 95° and χ² = 70° in structures (C), (D) and (E); and (iii) χ¹ = − 163° and χ² = 87° in structure (F). The conformation of Tyr70 in (F) corresponds to the structure of a conformational intermediate. This is the first structure which demonstrates that the slow conversion of DNA substrates by RIPs can be trapped during crystallization.     

Modulation of advanced glycation endproduct synthesis by kynurenines in human lens proteins

Human lens proteins (HLP) become chemically modified by kynurenines and advanced glycation end products (AGEs) during aging and cataractogenesis. We investigated the effects of kynurenines on AGE synthesis in HLP. We found that incubation with 5 mM ribose or 5 mM ascorbate produced significant quantities of pentosidine, and this was further enhanced in the presence of two different kynurenines (200–500 µM): N-formylkynurenine (Nfk) and kynurenine (Kyn). Another related compound, 3-hydroxykynurenine (3OH-Kyn), had disparate effects; low concentrations (10–200 µM) promoted pentosidine synthesis, but high concentrations (200–500 µM) inhibited it. 3OH-Kyn showed similar effects on pentosidine synthesis from Amadori-enriched HLP or ribated lysine. Chelex-100 treatment of phosphate buffer reduced pentosidine synthesis from Amadori-enriched HLP by ∼ 90%, but it did not inhibit the stimulating effect of 3OH-Kyn and EDTA. 3OH-Kyn (100–500 μM) spontaneously produced copious amounts of H₂O₂ (10–25 μM), but externally added H₂O₂ had only a mild stimulating effect on pentosidine but had no effect on N^ε-carboxymethyl lysine (CML) synthesis in HLP from ribose and ascorbate. Further, human lens epithelial cells incubated with ribose and 3OH-Kyn showed higher intracellular pentosidine than cells incubated with ribose alone. CML synthesis from glycating agents was inhibited 30 to 50% by 3OH-Kyn at concentrations of 100–500 μM. Argpyrimidine synthesis from 5 mM methylglyoxal was slightly inhibited by all kynurenines at concentrations of 100–500 μM. These results suggest that AGE synthesis in HLP is modulated by kynurenines, and such effects indicate a mode of interplay between kynurenines and carbohydrates important for AGE formation during lens aging and cataract formation.     

Insights into the fatty acid chain length specificity of the carboxylesterase PA3859 from Pseudomonas aeruginosa: A combined structural,biochemical and computational study

The open reading frame PA3859 of Pseudomonas aeruginosa encodes an intracellular carboxylesterase belonging to a group of microbial enzymes (EC 3.1.1.1) that catalyze the hydrolysis of aliphatic and aromatic esters with a broad substrate specificity. With few exceptions, for this class of enzymes, belonging to the α/β-hydrolase fold superfamily, very little information is available regarding their biochemical activity and in vivo function. The X-ray crystal structure of recombinant PA3859 has been determined for two crystal forms (space groups P2₁ and P2₁2₁2). The kinetic properties of the enzyme were studied using p-nitrophenyl esters as substrates and data fitted to a surface dilution mixed micelle kinetic model. Enzymatic assays and computational docking simulations, pinpointed the enzyme’s preference for esters of palmitic and/or stearic acids and provided insights into the enzyme–substrate favorable binding modes.     

Low-mass molecular dynamics simulation: A simple and generic technique to enhance configurational sampling

CLN025 is one of the smallest fast-folding proteins. Until now it has not been reported that CLN025 can autonomously fold to its native conformation in a classical, all-atom, and isothermal–isobaric molecular dynamics (MD) simulation. This article reports the autonomous and repeated folding of CLN025 from a fully extended backbone conformation to its native conformation in explicit solvent in multiple 500-ns MD simulations at 277 K and 1 atm with the first folding event occurring as early as 66.1 ns. These simulations were accomplished by using AMBER forcefield derivatives with atomic masses reduced by 10-fold on Apple Mac Pros. By contrast, no folding event was observed when the simulations were repeated using the original AMBER forcefields of FF12SB and FF14SB. The results demonstrate that low-mass MD simulation is a simple and generic technique to enhance configurational sampling. This technique may propel autonomous folding of a wide range of miniature proteins in classical, all-atom, and isothermal–isobaric MD simulations performed on commodity computers—an important step forward in quantitative biology.     

Simulation of the substrate cavity dynamics of quercetinase

Molecular dynamics (MD) simulations have been performed on quercetin 2,3 dioxygenase (2,3QD) to study the mobility and flexibility of the substrate cavity. 2,3QD is the only firmly established Cu-containing dioxygenase known so far. It catalyses the breakage of the O-heterocycle of flavonols. The substrates occupy a shallow and overall hydrophobic cavity proximal to the metal centre of the homo-dimeric enzyme. The linker connecting the C-terminal and N-terminal domains in the monomer is partly disordered in the crystal structure and part of it forms a flexible lid at the entrance of the substrate cavity. This loop has been tentatively assigned a role in the enzyme mechanism: it helps lock the substrate into place. The dynamics of this loop has been investigated by MD simulation. The initial coordinates were taken from the crystal structure of 2,3QD in the presence of the substrate kaempferol (KMP). After equilibration and simulation over 7.2ns the substrate was removed and another equilibration and simulation of 7.2ns was performed. The results show that the structures of the free enzyme as well as of the enzyme-substrate complex are stable in MD simulation. The linker shows strongly enhanced mobility in the loop region that is close to the entrance to the substrate cavity (residues 154-169). Movement of the loop takes place on a timescale of 5-10ns. To confirm the conclusions about the loop dynamics drawn from the 7.2ns simulation, the simulation was extended with another 8ns. When substrate binds into the cavity the loop orders remarkably, although mobility is retained by residues 155-158. Some regions of the loop (residues 154-160 and 164-176) move over a considerable distance and approach the substrate closely, reinforcing the idea that they lock the substrate in the substrate cavity. The enthalpic component of the interaction of the loop with the protein and the KMP appears to favour the locking of the substrate. Two water molecules were found immobilised in the cavity, one of which exhibited rotation on the picosecond timescale. When the substrate is removed, the empty cavity fills up with water within 200ps.     

Energy cost for the proper ionization of active site residues in 6-phosphogluconate dehydrogenase from T. brucei

The catalytic mechanism of 6-phosphogluconate dehydrogenase requires the inversion of a Lys/Glu couple from its natural ionization state. The pK_a of these residues in free and substrate bound enzymes has been determined measuring by ITC the proton release/uptake induced by substrate binding at different pH values. Wt 6-phosphogluconate dehydrogenase from Trypanosoma brucei and two active site enzyme mutants, K185H and E192Q were investigated. Substrate binding was accompanied by proton release and was dependent on the ionization of a group with pK_a 7.07 which was absent in the E192Q mutant. Kinetic data highlighted two pK_a, 7.17 and 9.64, in the enzyme–substrate complex, the latter being absent in the E192Q mutant, suggesting that the substrate binding shifts Glu192 pK_a from 7.07 to 9.64. A comparison of wt and E192Q mutant appears to show that the substrate binding shifts Lys185 pK_a from 9.9 to 7.17. By comparing differences in proton release and the binding enthalpy of wt and mutant enzymes, the enthalpic cost of the change in the protonation state of Lys185 and Glu192 was estimated at ≈ 6.1 kcal/mol. The change in protonation state of Lys185 and Glu192 has little effect on Gibbs free energy, 240–325 cal/mol. However proton balance evidences the dissociation of other group(s) that can be collectively described by a single pK_a shift from 9.1 to 7.54. This further change in ionization state of the enzyme causes an increase of free energy with a total cost of 1.2–2.3 kcal/mol to set the enzyme into a catalytically competent form.     

Crystal contact-free conformation of an intrinsically flexible loop in protein crystal: Tim21 as the case study

BackgroundIn protein crystals, flexible loops are frequently deformed by crystal contacts, whereas in solution, the large motions result in the poor convergence of such flexible loops in NMR structure determinations. We need an experimental technique to characterize the structural and dynamic properties of intrinsically flexible loops of protein molecules.MethodsWe designed an intended crystal contact-free space (CCFS) in protein crystals, and arranged the flexible loop of interest in the CCFS. The yeast Tim 21 protein was chosen as the model protein, because one of the loops (loop 2) is distorted by crystal contacts in the conventional crystal.ResultsYeast Tim21 was fused to the MBP protein by a rigid α-helical linker. The space created between the two proteins was used as the CCFS. The linker length provides adjustable freedom to arrange loop 2 in the CCFS. We re-determined the NMR structure of yeast Tim21, and conducted MD simulations for comparison. Multidimensional scaling was used to visualize the conformational similarity of loop 2. We found that the crystal contact-free conformation of loop 2 is located close to the center of the ensembles of the loop 2 conformations in the NMR and MD structures.ConclusionsLoop 2 of yeast Tim21 in the CCFS adopts a representative, dominant conformation in solution.General significanceNo single powerful technique is available for the characterization of flexible structures in protein molecules. NMR analyses and MD simulations provide useful, but incomplete information. CCFS crystallography offers a third route to this goal.