Unusual Role of the 3-OH Group of Oligosaccharide Substrates in the Mechanism of Bacillus 1,3-1,4-β-glucanase

Abstract

In the mechanism of retaining β-glycosidases, the 2-hydroxyl group of the substrate in the monosaccharyl unit involved in catalysis (subsite -1) is beleived to play an important role through hydrogen bonding interactions with protein residues that are optimized at the transition state. Commonly, removal of the 2-OH group of the substrate results in a 10–12 kcal·mol^-1 transition state destabilization. However, this effect seems not to be general as reported here for Bacillus 1,3-1,4-β-glucanase, a family 16 retaining endo-glycosidase. A p-nitrophenol 2-deosxy tetrasaccharide substrate was synthesized to probe the involvement of the 2-OH group in catalysis. Comparative kinetics with wild-type and subsite +1 mutants show that the 2-deoxy analog is a better substrate than the corresponding 2-hydroxy substrate. It is tentatively proposed that the 2-deoxy analog adopts a different conformation upon binding that compensates for the lack of the 2-OH substituent.     

Comparison of Full-atomic and Coarse-grained Models to Examine the Molecular Fluctuations of c-AMP Dependent Protein Kinase

Abstract

Molecular fluctuations of the native conformation of c-AMP dependent protein kinase (cAPK) have been investigated with three different approaches. The first approach is the full atomic normal mode analysis (NMA) with empirical force fields. The second and third approaches are based on a coarse-grained model with a single single-parameter- harmonic potential between close residues in the crystal structure of the molecule without any residue specificity. The second method calculates only the magnitude of fluctuations whereas the third method is developed to find the directionality of the fluctuations which are essential to understand the functional importance of biological molecules. The aim, in this study, is to determine whether using such coarse-grained models are appropriate for elucidating the global dynamic characteristics of large proteins which reduces the size of the system at least by a factor of ten. The mean-square fluctuations of C^α atoms and the residue cross-correlations are obtained by three approaches. These results are then compared to test the results of coarse grained models on the overall collective motions. AH three of the approaches show that highly flexible regions correspond to the activation and solvent exposed loops, whereas the conserved residues (especially in substrate binding regions) exhibit almost no flexibility, adding stability to the structure. The anti-correlated motions of the two lobes of the catalytic core provide flexibility to the molecule. High similarities among the results of these methods indicate that the slowest modes governing the most global motions are preserved in the coarse grained models for proteins. This finding may suggest that the general shapes of the structures are representative of their dynamic characteristics and the dominant motions of protein structures are robust at coarse-grained levels.     

Identification and characterization of Helicobacter pylori O‐acetylserine‐dependent cystathionine β‐synthase,a distinct member of the PLP‐II family

O‐acetylserine sulfhydrylase (OASS) and cystathionine β‐synthase (CBS) are members of the PLP‐II family, and involved in L‐cysteine production. OASS produces L‐cysteine via a de novo pathway while CBS participates in the reverse transsulfuration pathway. O‐acetylserine‐dependent CBS (OCBS) was previously identified as a new member of the PLP‐II family, which are predominantly seen in bacteria. The bacterium Helicobacter pylori possess only one OASS (hp0107) gene and we showed that the protein coded by this gene actually functions as an OCBS and utilizes L‐homocysteine and O‐acetylserine (OAS) to produce cystathionine. HpOCBS did not show CBS activity with the substrate L‐serine and required OAS exclusively. The HpOCBS structure in complex with methionine showed a closed cleft state, explaining the initial mode of substrate binding. Sequence and structural analyses showed differences between the active sites of OCBS and CBS, and explain their different substrate preferences. We identified three hydrophobic residues near the active site of OCBS, corresponding to one serine and two tyrosine residues in CBSs. Mutational studies were performed on HpOCBS and Saccharomyces cerevisiae CBS. A ScCBS double mutant (Y158F/Y226V) did not display activity with L‐serine, indicating indispensability of these polar residues for selecting substrate L‐serine, however, did show activity with OAS.     

Kinetics of protein folding: Nucleation mechanism,time scales,and pathways

The kinetics and thermodynamics of protein folding is investigated using low friction Langevin simulation of minimal continuum mode of proteins. We show that the model protein has two characteristic temperatures: (a) T_θ, at which the chain undergoes a collapse transition from an extended conformation; (b) T_f(< T_θ), at which a finite size first-order transition to the folded state takes place. The kinetics of approach to the native state from initially denatured conformations is probed by several novel correlation functions. We find that the overall kinetics of approach to the native conformation occurs via a three-stage multiple pathway mechanism. The initial stage, characterized by a series of local dihedral angle transitions, eventually results in the compaction of the protein. Subsequently, the molecule acquires native-like structures during the second stage of folding. The final stage of folding involves activated transitions from one of the native-like structures to the native conformation. The first two stages are characterized by a multiplicity of pathways while relatively few paths are involved in the final stage. A detailed analysis of the dynamics of individual trajectories reveals a novel picture of protein folding. We find that afraction of the initial population reaches the native conformation without the formation of any detectable intermediates. This pathway is associated with a nucleation mechanism, i.e., once a critical number of tertiary contacts are established then the native state is reached rapidly. The remaining fraction of molecules become trapped in misfolded structures (stabilized by incorrect tertiary contacts). The slow folding involves transitions over barriers from these structures to the native conformation. The theoretical predictions are compared with recent experiments that probe protein folding kinetics by hydrogen exchange labeling technique. © 1995 John Wiley & Sons, Inc.     

Kinetics and thermodynamics of the rate-limiting conformational change in the actomyosin V mechanochemical cycle

We used FRET to examine the kinetics and thermodynamics of structural changes associated with ADP release in myosin V, which is thought to be a strain-sensitive step in many muscle and non-muscle myosins. We also explored essential dynamics using FIRST/FRODA starting with three different myosin V X-ray crystal structures to examine intrinsic flexibility and correlated motions. Our steady-state and time-resolved FRET analysis demonstrates a temperature-dependent reversible conformational change in the nucleotide-binding pocket (NBP). Our kinetic results demonstrate that the NBP goes from a closed to an open conformation prior to the release of ADP, while the actin-binding cleft remains closed. Interestingly, we find that the temperature dependence of the maximum actin-activated myosin V ATPase rate is similar to the pocket opening step, demonstrating that this is the rate-limiting structural transition in the ATPase cycle. Thermodynamic analysis demonstrates that the transition from the open to closed NBP conformation is unfavorable because of a decrease in entropy. The intrinsic flexibility analysis is consistent with conformational entropy playing a role in this transition as the MV.ADP structure is highly flexible compared to the MV.APO structure. Our experimental and modeling studies support the conclusion of a novel post-power-stroke actomyosin.ADP state in which the NBP and actin-binding cleft are closed. The novel state may be important for strain sensitivity as the transition from the closed to open NBP conformation may be altered by lever arm position.     

Investigation of Domain Motions in Bacteriophage T4 Lysozyme

Abstract

Hinge-bending in T4 lysozyme has been inferred from single amino acid mutant crystalline allomorphs by Matthews and coworkers. This raises an important question: are the different conformers in the unit cell artifacts of crystal packing forces, or do they represent different solution state structures? The objective of this theoretical study is to determine whether domain motions and hinge-bending could be simulated in T4 lysozyme using molecular dynamics. An analysis of a 400 ps molecular dynamics simulation of the 164 amino acid enzyme T4 lysozyme is presented. Molecular dynamics calculations were computed using the Discover software package (Biosym Technologies). All hydrogen atoms were modeled explicitly with the inclusion of all 152 crystallographic waters at a temperature of 300 K. The native T4 lysozyme molecular dynamics simulation demonstrated hinge-bending in the protein. Relative domain motions between the N-terminal and C-terminal domains were evident. The enzyme hinge bending sites resulted from small changes in backbone atom conformations over several residues rather than rotation about a single bound. Two hinge loci were found in the simulation. One locus comprises residues 8–14 near the C-terminal of the A helix; the other site, residues 77–83 near the C-terminal of the C helix. Comparison of several snapshot structures from the dynamics trajectory clearly illustrates domain motions between the two lysozyme lobes. Time correlated atomic motions in the protein were analyzed using a dynamical cross-correlation map. We found a high degree of correlated atomic motions in each of the domains and, to a lesser extent, anticorrelated motions between the two domains. We also found that the hairpin loop in the N-terminal lobe (residues 19–24) acted as a mobile ‘flap’ and exhibited highly correlated dynamic motions across the cleft of the active site, especially with residue 142.     

Smyd3 open & closed lock mechanism for substrate recruitment: The hinge motion of C-terminal domain inferred from μ-second molecular dynamics simulations

BackgroundThe human lysine methyltransferase Smyd3, a member of the SET and MYND domain containing protein family, harbors methylation activity on both histone and non-histone targets in a tightly regulated manner. The mechanism of how Smyd3 dynamically regulates substrate recognition is still not fully unveiled.MethodsHere, we employed molecular dynamics simulations on full length human Smyd3, performed to a total of 1.2 μ-second, in the presence (holo) and absence (apo) of the S-Adenosyl methionine (AdoMet) cofactor. The dynamical features of Smyd3 in apo and holo states have been examined and compared via examining geometrical and electrostatic properties.ResultsThe results show a distinct dynamics of the C-terminal domain (CTD) in the two states. In the apo state, the CTD undergoes a large hinge like motion and samples more opened configurations, thus acting like a loosened clamp and resulting in expanded substrate binding crevice. In the holo state, the CTD exhibits a restricted motion while the overall structure remains compact, mimicking a closed clamp. This leads to a localized increase in the negative potential at the substrate binding cleft. Further, solvent accessibility of critical residues at the target lysine access channel, important for methylation activity, is increased.ConclusionsWe postulate that AdoMet cofactor acts like a key and locks Smyd3 in a closed conformation. In effect, the cofactor binding restricts the elasticity of the CTD, presenting a compact substrate binding cleft with high negative potential, which may have implications on substrate recruitment via long range electrostatics.General significanceThe deletion of the CTD from Smyd3 has been shown to abolish the basal histone methylation activity. Our study highlights the importance of the CTD elasticity in shaping the substrate binding site for recognition and supports the previously proposed role of the CTD in stabilizing the active site for methylation activity.     

Dissecting the Conformational Dynamics of the Bile Acid Transporter Homologue ASBTNM

Apical sodium-dependent bile acid transporter (ASBT) catalyses uphill transport of bile acids using the electrochemical gradient of Na⁺ as the driving force. The crystal structures of two bacterial homologues ASBT_NM and ASBT_Yf have previously been determined, with the former showing an inward-facing conformation, and the latter adopting an outward-facing conformation accomplished by the substitution of the critical Na⁺-binding residue glutamate-254 with an alanine residue. While the two crystal structures suggested an elevator-like movement to afford alternating access to the substrate binding site, the mechanistic role of Na⁺ and substrate in the conformational isomerization remains unclear. In this study, we utilized site-directed alkylation monitored by in-gel fluorescence (SDAF) to probe the solvent accessibility of the residues lining the substrate permeation pathway of ASBT_NM under different Na⁺ and substrate conditions, and interpreted the conformational states inferred from the crystal structures. Unexpectedly, the crosslinking experiments demonstrated that ASBT_NM is a monomer protein, unlike the other elevator-type transporters, usually forming a homodimer or a homotrimer. The conformational dynamics observed by the biochemical experiments were further validated using DEER measuring the distance between the spin-labelled pairs. Our results revealed that Na⁺ ions shift the conformational equilibrium of ASBT_NM toward the inward-facing state thereby facilitating cytoplasmic uptake of substrate. The current findings provide a novel perspective on the conformational equilibrium of secondary active transporters.     

Docking Simulation and Competitive Experiments Validate the Interaction Between the 2,5-Xylidine Inhibitor and Rigidoporus lignosus Laccase

Abstract

Laccases are polyphenol oxidases which oxidize a broad range of reducing substrates, preferably phenolic compounds, and their use in biotechnological applications is increasing.

Recently, the first X-ray structure of active laccase from white rot fungus Rigidoporus lignosus has been reported containing a full complement of copper ions. Comparison among selected fungal laccases of known 3D structure has shown that the Rigidoporus lignosus laccase has a very high similarity with the Trametes versicolor laccase that, being co-crystallized with 2,5-xylidine, shows a well defined binding pocket for the substrate. Global sequence alignment between Rigidoporus lignosus and Trametes versicolor laccases shows 73% of identity but, surprisingly, there is no identity and neither conservative substitutions between the residues composing the loops directly contacting the 2,5-xylidine. Moreover the structural alignment of these two enzymes identifies in these loops a striking structural similarity proposing the question if 2,5-xylidine may bind in same enzyme pocket.

Here we report the protein-ligand docking simulation of 3D structure of Rigidoporus lignosus laccase and 2,5-xylidine. Docking simulation analyses show that spatial conformation of the two 2,5-xylidine binding pockets, despite differences in the residues directly contacting the ligand, may arrange a similar pocket that allows a comparable accommodation of the inhibitor. To validate these results the binding of 2,5-xylidine in the substrate cavity has been confirmed by kinetic competitive experiments.     

Poisson-Boltzmann model studies of molecular electrostatic properties of the cAMP-dependent protein kinase

Protonation equilibria of residues important in the catalytic mechanism of a protein kinase were analyzed on the basis of the Poisson-Boltzmann electrostatic model along with a cluster-based treatment of the multiple titration state problem. Calculations were based upon crystallographic structures of the mammalian cAMP-dependent protein kinase, one representing the so called closed form of the enzyme and the other representing an open conformation. It was predicted that at pH 7 the preferred form of the phosphate group at the catalytically essential threonine 197 (P-Thr197) in the closed form is dianionic, whereas in the open form a monoanionic ionization state is preferred. This dianionic state of P-Thr197, in the closed form, is stabilized by interactions with ionizable residues His87, Arg165, and Lys189. Our calculations predict that the hydroxyl of the Ser residue in the peptide substrate is very difficult to ionize, both in the closed and open structures of the complex. Also, the supposed catalytic base, Asp166, does not seem to have a pK _a appropriate to remove the hydroxyl group proton of the peptide substrate. However, when Ser of the peptide substrate is forced to remain ionized, the predicted pK _a of Asp166 increases strongly, which suggests that the Asp residue is a likely candidate to attract the proton if the Ser residue becomes deprotonated, possibly during some structural change preceding formation of the transition state. Finally, in accord with suggestions made on the basis of the pH-dependence of kinase kinetics, our calculations predict that Glu230 and His87 are the residues responsible for the molecular pK _a values of 6.2 and 8.5, observed in the experiment. Received: 19 October 1998 / Revised version: 12 April 1999 / Accepted: 15 April 1999     

Protonation of Glu135 Facilitates the Outward-to-Inward Structural Transition of Fucose Transporter

Major facilitator superfamily (MFS) transporters typically need to alternatingly sample the outward-facing and inward-facing conformations, in order to transport the substrate across membrane. To understand the mechanism, in this work, we focused on one MFS member, the L-fucose/H⁺ symporter (FucP), whose crystal structure exhibits an outward-open conformation. Previous experiments imply several residues critical to the substrate/proton binding and structural transition of FucP, among which Glu¹³⁵, located in the periplasm-accessible vestibule, is supposed as being involved in both proton translocation and conformational change of the protein. Here, the structural transition of FucP in presence of substrate was investigated using molecular-dynamics simulations. By combining the equilibrium and accelerated simulations as well as thermodynamic calculations, not only was the large-scale conformational change from the outward-facing to inward-facing state directly observed, but also the free energy change during the structural transition was calculated. The simulations confirm the critical role of Glu¹³⁵, whose protonation facilitates the outward-to-inward structural transition both by energetically favoring the inward-facing conformation in thermodynamics and by reducing the free energy barrier along the reaction pathway in kinetics. Our results may help the mechanistic studies of both FucP and other MFS transporters.     

Lid closure dynamics of porcine pancreatic lipase in aqueous solution

BackgroundPancreatic lipases hydrolyze fatty acids in dietary pathway. The activity of porcine pancreatic lipase (PPL) is controlled by lid domain along with a coenzyme, colipase. The active open-state conformation of the protein could be induced by detergents or bile salts which would be further stabilized by binding of colipase. In the absence of these interactions, the lid preferably attains a closed conformation in water.MethodsMolecular dynamic simulation was used to monitor the lid movement of PPL in open and closed conformations in water. Free energy surface was constructed from the simulation. Energy barriers and major structural changes during lid opening were evaluated.ResultsThe lid closure of PPL in water from its open conformation might be initiated by columbic interactions which initially move the lid away from domain 1. This is followed by major dihedral changes on the lid residues which alter the trajectory of motion. The lid then swirls back towards domain 1 to attain closed conformation. This is accompanied with conformational changes around β5- and β9-loops as well. However, PPL in closed conformation shows only the domain movements and the lid remains in its closed conformation.ConclusionsPPL in closed conformation is stable in water and the open conformation is driven towards closed state. The lid follows a swirling trajectory during the closure.General significanceConformational state of the lid regulates the activity and substrate specificity of PPL. Hence, it is essential to understand the lid dynamics and the role of specific amino acid residues involved.     

Influence of an extrinsic cross-link on the folding pathway of ribonuclease A. Kinetics of folding-unfolding

The kinetics of folding/unfolding of cross-linked Lys7-dinitrophenylene-Lys41-ribonuclease A were studied and compared to those of unmodified ribonuclease A (RNase A) at various concentrations of guanidine hydrochloride. The folding of the denatured cross-linked protein involved one fast-folding species (22 +/- 4%) and two slow-folding species, as observed in unmodified ribonuclease A. Also, a nativelike intermediate, analogous to that reported previously for unmodified ribonuclease A [Cook, K. H., Schmid, F. X., & Baldwin, R. L. (1979) Proc. Natl. Acad. Sci. U.S.A. 76, 6157], has been detected on the folding pathway of cross-linked ribonuclease A. The extrinsic cross-link between Lys7 and Lys41 did not affect the rate constants for the folding kinetics of these three species. The cross-link did, however, significantly affect the rate constant for unfolding of the native protein. The conformation of the protein in the transition state of the unfolding pathway was deduced from an analysis of the kinetic data. It appears that the 41 N-terminal residues are unfolded in the transition state of the unfolding pathway. Thus, the unfolding pathway of RNase A is sequential in that further unfolding (after the transition state) follows the unfolding of the 41 N-terminal residues. Also, the conformation of the 41 N-terminal residues does not play a role in the folding pathway. Presumably, if the cross-link were introduced instead between two other residues that are in the segment(s) involved in the rate-limiting step(s), it could increase the refolding rate constants and possibly the concentration of fast-folding species.     

Structural Insights into the Substrate Specificity of (S)-Ureidoglycolate Amidohydrolase and Its Comparison with Allantoate Amidohydrolase

In plants, the ureide pathway is a metabolic route that converts the ring nitrogen atoms of purine into ammonia via sequential enzymatic reactions, playing an important role in nitrogen recovery. In the final step of the pathway, (S)-ureidoglycolate amidohydrolase (UAH) catalyzes the conversion of (S)-ureidoglycolate into glyoxylate and releases two molecules of ammonia as by-products. UAH is homologous in structure and sequence with allantoate amidohydrolase (AAH), an upstream enzyme in the pathway with a similar function as that of an amidase but with a different substrate. Both enzymes exhibit strict substrate specificity and catalyze reactions in a concerted manner, resulting in purine degradation. Here, we report three crystal structures of Arabidopsis thaliana UAH (bound with substrate, reaction intermediate, and product) and a structure of Escherichia coli AAH complexed with allantoate. Structural analyses of UAH revealed a distinct binding mode for each ligand in a bimetal reaction center with the active site in a closed conformation. The ligand directly participates in the coordination shell of two metal ions and is stabilized by the surrounding residues. In contrast, AAH, which exhibits a substrate-binding site similar to that of UAH, requires a larger active site due to the additional ureido group in allantoate. Structural analyses and mutagenesis revealed that both enzymes undergo an open-to-closed conformational transition in response to ligand binding and that the active-site size and the interaction environment in UAH and AAH are determinants of the substrate specificities of these two structurally homologous enzymes.     

Holo‐BtuF stabilizes the open conformation of the vitamin B12 ABC transporter BtuCD

While there is evidence that other ABC transporters can tell between empty and loaded substrate binding protein, reconstitution experiments suggest otherwise for the Escherichia coli vitamin B12 importer BtuCD‐F. Here, we address the question of BtuCD‐F substrate sensitivity in a combined protein–protein docking and molecular dynamics simulation approach. Starting from the BtuCD and holo‐BtuF crystal structures, we model two holo‐BtuCD‐F docking complexes differing by a 180° orientation of BtuF. One of these is similar to the apo‐BtuCD‐F crystal structure. Both docking complexes were embedded in a lipid/water environment to investigate their dynamics and BtuCD's conformational response to the presence and absence of BtuF, vitamin B12, and Mg‐ATP in a series of 28 independent MD simulations. We find holo‐BtuF stabilizing the open conformation of BtuCD, whereas the transporter begins to close again when BtuF or vitamin B12 is removed—suggesting BtuCD‐F is capable of substrate sensitivity. We identified BtuC transmembrane helices 3 and 5, the L‐loops and the adjacent helices comprised of BtuC residues 170–180 as hotspots of conformational change. We propose the latter to act as substrate sensors. BtuF‐Trp44 appears to act as a lid on the vitamin B12 binding cleft in BtuF X‐ray structures and protrudes into the BtuCD transport channel in one of our simulations, which might represent an initial step in vitamin B12 uptake. On an average, we observe subunit motions where the nucleotide binding domains approach each other while the transmembrane domains display an opening trend toward the periplasm. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Statistical Mechanical Approach for Predicting the Transition to Non-B DNA Structures in Supercoiled DNA

Abstract

Supercoiling causes global twist of DNA structure and the supercoiled state has wide influence on conformational transition. A statistical mechanical approach was made for prediction of the transition probability to non-B DNA structures under torsional stress. A conditional partition function was defined as the sum over all possible states of the DNA sequence with basepair 1 and basepair n being in B-form helix and a recurrence formula was developed which expressed the partition function for basepair n with those for less number of pairs. This new definition permits a quick enumeration of every configuration of secondary structures. Energetic parameters of all conformations concerned, involving B-form, interior loop, cruciform and Z-form, were included in the equation. The probability of transition to each non-B conformation could be derived from these conditional partition functions. For treatment of effects of superhelicity, supercoiling energy was considered, and a twist of each conformation was determined to minimize the supercoiling energy. As the twist itself affects the transition probability, the whole scheme of equations was solved by renormalization technique. The present method permits a simultaneous treatment of serveral types of conformations under a common torsional stress.

A set of energetic parameters of DNA secondary structures has been chosen for calculation. Some DNA sequences were submitted to the calculation, and all the sequences that we submitted gave stable convergence. Some of them have been investigated the critical supercoil density for the transition to non-B DNA structures. Even though the reliability of the set of parameters was not enough, the prediction of secondary structure transition showed good agreement with reported observation. Hence, the present algorithm can estimate the probability of local conformational change of DNA under a given supercoil density, and also be employed to predict some specific sequences in which conformational change is sensitive to superhelicity.     

Influence of long-range contacts and surrounding residues on the transition state structures of proteins

Understanding the parameters influencing the formation of transition state structures in proteins is an important problem in protein folding and kinetics. In this work, we have analyzed the structure-based parameters, surrounding hydrophobicity, secondary structure, solvent accessibility, number of medium- and long-range contacts, and surrounding residues for understanding the transition state structures of 15 proteins. The analysis of Φ-values shows that 29% of the studied 378 mutants have a Φ-value of more than 0.5. The combination of different structure-based parameters could discriminate the residues that have a Φ-value cutoff of more than 0.5 with a 5-fold cross-validation accuracy of 68%, which indicates that the surrounding residues and contacts play important roles in the formation of transition state structures. Systematic analysis on different proteins reveals that the proteins azurin, cold shock protein, and C-terminal domain of ribosomal protein L9 are influenced by the number of medium- and long-range proteins, whereas barnase, FK506 binding protein, and IM9 are influenced by surrounding residues. The discrimination accuracy lies in the ranges of 81–95% and 74–85% for these respective classes of protein. Furthermore, the combination of surrounding residues and contacts improved the accuracy up to 24% in other considered proteins. We suggest that the structure-based parameters along with noncovalent interactions and conservation of residues may aid in identifying the potential residues in the formation of transition state structures in proteins.     

Insight Derived from Molecular Dynamics Simulation into Substrate-Induced Changes in Protein Motions of Proteinase K

Abstract

Because of the significant industrial, agricultural and biotechnological importance of serine protease proteinase K, it has been extensively investigated using experimental approaches such as X-ray crystallography, site-directed mutagenesis and kinetic measurement. However, detailed aspects of enzymatic mechanism such as substrate binding, release and relevant regulation remain unstudied. Molecular dynamics (MD) simulations of the proteinase K alone and in complex with the peptide substrate AAPA were performed to investigate the effect of substrate binding on the dynamics/molecular motions of proteinase K. The results indicate that during simulations the substrate-complexed proteinase K adopt a more compact and stable conformation than the substrate-free form. Further essential dynamics (ED) analysis reveals that the major internal motions are confined within a subspace of very small dimension. Upon substrate binding, the overall flexibility of the protease is reduced; and the noticeable displacements are observed not only in substrate-binding regions but also in regions opposite the substrate-binding groove/pockets. The dynamic pockets caused by the large concerted motions are proposed to be linked to the substrate recognition, binding, orientation and product release; and the significant displacements in regions opposite the binding groove/pockets are considered to play a role in modulating the dynamics of enzyme-substrate interaction. Our simulation results complement the biochemical and structural studies, highlighting the dynamic mechanism of the functional properties of proteinase K.     

Structural enzymological studies of 2-enoyl thioester reductase of the human mitochondrial FAS II pathway: new insights into its substrate recognition properties

Structural and kinetic properties of the human 2-enoyl thioester reductase [mitochondrial enoyl-coenzyme A reductase (MECR)/ETR1] of the mitochondrial fatty acid synthesis (FAS) II pathway have been determined. The crystal structure of this dimeric enzyme (at 2.4 Å resolution) suggests that the binding site for the recognition helix of the acyl carrier protein is in a groove between the two adjacent monomers. This groove is connected via the pantetheine binding cleft to the active site. The modeled mode of NADPH binding, using molecular dynamics calculations, suggests that Tyr94 and Trp311 are critical for catalysis, which is supported by enzyme kinetic data. A deep, water-filled pocket, shaped by hydrophobic and polar residues and extending away from the catalytic site, was recognized. This pocket can accommodate a fatty acyl tail of up to 16 carbons. Mutagenesis of the residues near the end of this pocket confirms the importance of this region for the binding of substrate molecules with long fatty acyl tails. Furthermore, the kinetic analysis of the wild-type MECR/ETR1 shows a bimodal distribution of catalytic efficiencies, in agreement with the notion that two major products are generated by the mitochondrial FAS II pathway.     

Identification of Residues Surrounding the Active Site of Type A Botulinum Neurotoxin Important for Substrate Recognition and Catalytic Activity

Type A botulinum neurotoxin is one of the most lethal of the seven serotypes and is increasingly used as a therapeutic agent in neuromuscular dysfunctions. Its toxic function is related to zinc-endopeptidase activity of the N-terminal light chain (LC) on synaptosome-associated protein-25 kDa (SNAP-25) of the SNARE complex. To understand the determinants of substrate specificity and assist the development of strategies for effective inhibitors, we used site-directed mutagenesis to investigate the effects of 13 polar residues of the LC on substrate binding and catalysis. Selection of the residues for mutation was based on a computational analysis of the three-dimensional structure of the LC modeled with a 17-residue substrate fragment of SNAP-25. Steady-state kinetic parameters for proteolysis of the substrate fragment were determined for a set of 16 single mutants. Of the mutated residues non-conserved among the serotypes, replacement of Arg-230 and Asp-369 by polar or apolar residues resulted in drastic lowering of the catalytic rate constant (k _cat), but had less effect on substrate affinity (K _m). Substitution of Arg-230 with Lys decreased the catalytic efficiency (k _cat/K _m) by 50-fold, whereas replacement by Leu yielded an inactive protein. Removal of the electrostatic charge at Asp-369 by mutation to Asn resulted in 140-fold decrease in k _cat/K _m. Replacement of other variable residues surrounding the catalytic cleft (Glu-54, Glu-63, Asn-66, Asp-130, Asn-161, Glu-163, Glu-170, Glu-256), had only marginal effect on decreasing the catalytic efficiency, but unexpectedly the substitution of Lys-165 with Leu resulted in fourfold increase in k _cat/K _m. For comparison purposes, two conserved residues Arg-362 and Tyr-365 were investigated with substitutions of Leu and Phe, respectively, and their catalytic efficiency decreased 140- and 10-fold, respectively, whereas substitution of the tyrosine ring with Asn abolished activity. The altered catalytic efficiencies of the mutants were not due to any significant changes in secondary or tertiary structures, or in zinc content and thermal stability. We suggest that, despite the large minimal substrate size for catalysis, only a few non-conserved residues surrounding the active site are important to render the LC competent for catalysis or provide conformational selection of the substrate.