Triggering loops and enzyme function: identification of loops that trigger and modulate movements

Remote loop evolution reveals a complex biological function for chitinase enzymes beyond the active site

Loops are small secondary structural elements that play a crucial role in the emergence of new enzyme functions. However, the evolutionary molecular mechanisms how proteins acquire these loop elements and obtain new function is poorly understood. To address this question, we study glycoside hydrolase family 19 (GH19) chitinase-an essential enzyme family for pathogen degradation in plants. By revealing the evolutionary history and loops appearance of GH19 chitinase, we discover that one loop which is remote from the catalytic site, is necessary to acquire the new antifungal activity. We demonstrate that this remote loop directly accesses the fungal cell wall, and surprisingly, it needs to adopt a defined structure supported by long-range intramolecular interactions to perform its function. Our findings prove that nature applies this strategy at the molecular level to achieve a complex biological function while maintaining the original activity in the catalytic pocket, suggesting an alternative way to design new enzyme function.

Loop dynamics and the evolution of enzyme activity

In the early 2000s, Tawfik presented his 'New View' on enzyme evolution, highlighting the role of conformational plasticity in expanding the functional diversity of limited repertoires of sequences. This view is gaining increasing traction with increasing evidence of the importance of conformational dynamics in both natural and laboratory evolution of enzymes. The past years have seen several elegant examples of harnessing conformational (particularly loop) dynamics to successfully manipulate protein function. This Review revisits flexible loops as critical participants in regulating enzyme activity. We showcase several systems of particular interest: triosephosphate isomerase barrel proteins, protein tyrosine phosphatases and β-lactamases, while briefly discussing other systems in which loop dynamics are important for selectivity and turnover. We then discuss the implications for engineering, presenting examples of successful loop manipulation in either improving catalytic efficiency, or changing selectivity completely. Overall, it is becoming clearer that mimicking nature by manipulating the conformational dynamics of key protein loops is a powerful method of tailoring enzyme activity, without needing to target active-site residues.

Increasing the O2 Resistance of the [FeFe]-Hydrogenase CbA5H through Enhanced Protein Flexibility

The high turnover rates of [FeFe]-hydrogenases under mild conditions and at low overpotentials provide a natural blueprint for the design of hydrogen catalysts. However, the unique active site (H-cluster) degrades upon contact with oxygen. The [FeFe]-hydrogenase fromClostridium beijerinckii (CbA5H) is characterized by the flexibility of its protein structure, which allows a conserved cysteine to coordinate to the active site under oxidative conditions. Thereby, intrinsic cofactor degradation induced by dioxygen is minimized. However, the protection from O₂ is only partial, and the activity of the enzyme decreases upon each exposure to O₂. By using site-directed mutagenesis in combination with electrochemistry, ATR-FTIR spectroscopy, and molecular dynamics simulations, we show that the kinetics of the conversion between the oxygen-protected inactive state (cysteine-bound) and the oxygen-sensitive active state can be accelerated by replacing a surface residue that is very distant from the active site. This sole exchange of methionine for a glutamate residue leads to an increased resistance of the hydrogenase to dioxygen. With our study, we aim to understand how local modifications of the protein structure can have a crucial impact on protein dynamics and how they can control the reactivity of inorganic active sites through outer sphere effects.

Structural and biochemical insight into mode of action and subsite specificity of a chitosan degrading enzyme from Bacillus spec. MN

Chitosans, partially de-N-acetylated derivatives of chitin, are multifunctional biopolymers. In nature, biological activities of partially acetylated chitosan polymers are mediated in part by their oligomeric breakdown products, which are generated in situ by the action of chitosanolytic enzymes. Understanding chitosanolytic enzymes, therefore, can lead to the production of chitosan oligomers with fully defined structures that may confer specific bioactivities. To address whether defined oligomer products can be produced via chitosanolytic enzymes, we here characterized a GH8 family chitosanase from Bacillus spec. MN, determining its mode of action and product profiles. We found that the enzyme has higher activity towards polymers with lower degree of acetylation. Oligomeric products were dominated by GlcN₃, GlcN₃GlcNAc₁, and GlcN₄GlcNAc₁. The product distribution from oligomers were GlcN₃ > GlcN₂. Modeling and simulations show that the binding site comprises subsites ranging from (-3) to (+3), and a putative (+4) subsite, with defined preferences for GlcN or GlcNAc at each subsite. Flexible loops at the binding site facilitate enzyme-substrate interactions and form a cleft at the active site which can open and close. The detailed insight gained here will help to engineer enzyme variants to produce tailored chitosan oligomers with defined structures that can then be used to probe their specific biological activities.

pH-Induced Rotation of Lidless Lipase LipA from Bacillus subtilis at Lipase-Detergent Interface

Lipases exhibit a unique process during the catalysis of the hydrolysis of triglyceride substrates called interfacial activation. Surfactants are used as cosolvents with water not only to offer a less polar environment to the lipases needed for their interfacial activation but also to solvate the substrate which are poorly soluble in water. However, the presence of detergent in the medium can affect both the lipase and the substrate, making the construction of a microkinetic model for lipase activity in the presence of the detergent difficult. Herein, we study the interfacial activation of a lidless lipase LipA from Bacillus subtilis using extensive atomistic molecular dynamics simulations at different concentrations of the surfactant, Thesit (C12E8), at two pH values. Residues which bind to the monomeric detergent are found to be the same as the ones which have been reported earlier to bind to the substrate. Very importantly, a pH-induced rotation of the enzyme with respect to surfactant aggregate has been observed which not only explains the experimentally observed pH-dependent enzymatic activity of this lidless lipase, but also suggests its reorientation at an aqueous-lipodophilic interface.

A consensus-guided approach yields a heat-stable alkane-producing enzyme and identifies residues promoting thermostability

Aldehyde-deformylating oxygenase (ADO) is an essential enzyme for production of long-chain alkanes as drop-in biofuels, which are compatible with existing fuel systems. The most active ADOs are present in mesophilic cyanobacteria, especially Nostoc punctiforme Given the potential applications of thermostable enzymes in biorefineries, here we generated a thermostable (Cts)-ADO based on a consensus of ADO sequences from several thermophilic cyanobacterial strains. Using an in silico design pipeline and a metagenome library containing 41 hot-spring microbial communities, we created Cts-ADO. Cts-ADO displayed a 3.8-fold increase in pentadecane production on raising the temperature from 30 to 42 °C, whereas ADO from N. punctiforme (Np-ADO) exhibited a 1.7-fold decline. 3D structure modeling and molecular dynamics simulations of Cts- and Np-ADO at different temperatures revealed differences between the two enzymes in residues clustered on exposed loops of these variants, which affected the conformation of helices involved in forming the ADO catalytic core. In Cts-ADO, this conformational change promoted ligand binding to its preferred iron, Fe2, in the di-iron cluster at higher temperature, but the reverse was observed in Np-ADO. Detailed mapping of residues conferring Cts-ADO thermostability identified four amino acids, which we substituted individually and together in Np-ADO. Among these substitution variants, A161E was remarkably similar to Cts-ADO in terms of activity optima, kinetic parameters, and structure at higher temperature. A161E was located in loop L6, which connects helices H5 and H6, and supported ligand binding to Fe2 at higher temperatures, thereby promoting optimal activity at these temperatures and explaining the increased thermostability of Cts-ADO.

Engineered control of enzyme structural dynamics and function

Enzymes undergo a range of internal motions from local, active site fluctuations to large-scale, global conformational changes. These motions are often important for enzyme function, including in ligand binding and dissociation and even preparing the active site for chemical catalysis. Protein engineering efforts have been directed towards manipulating enzyme structural dynamics and conformational changes, including targeting specific amino acid interactions and creation of chimeric enzymes with new regulatory functions. Post-translational covalent modification can provide an additional level of enzyme control. These studies have not only provided insights into the functional role of protein motions, but they offer opportunities to create stimulus-responsive enzymes. These enzymes can be engineered to respond to a number of external stimuli, including light, pH, and the presence of novel allosteric modulators. Altogether, the ability to engineer and control enzyme structural dynamics can provide new tools for biotechnology and medicine.

Dynamics of the active site loops in catalyzing aminoacylation reaction in seryl and histidyl tRNA synthetases

Aminoacylation reaction is the first step of protein biosynthesis. The catalytic reorganization at the active site of aminoacyl tRNA synthetases (aaRSs) is driven by the loop motions. There remain lacunae of understanding concerning the catalytic loop dynamics in aaRSs. We analyzed the functional loop dynamics in seryl tRNA synthetase from Methanopyrus kandleri (^mkSerRS) and histidyl tRNA synthetases from Thermus thermophilus (^ttHisRS), respectively, using molecular dynamics. Results confirm that the motif 2 loop and other active site loops are flexible spots within the catalytic domain. Catalytic residues of the loops form a network of interaction with the substrates to form a reactive state. The loops undergo transitions between closed state and open state and the relaxation of the constituent residues occurs in femtosecond to nanosecond time scale. Order parameters are higher for constituent catalytic residues which form a specific network of interaction with the substrates to form a reactive state compared to the Gly residues within the loop. The development of interaction is supported from mutation studies where the catalytic domain with mutated loop exhibits unfavorable binding energy with the substrates. During the open-close motion of the loops, the catalytic residues make relaxation by ultrafast librational motion as well as fast diffusive motion and subsequently relax rather slowly via slower diffusive motion. The Gly residues act as a hinge to facilitate the loop closing and opening by their faster relaxation behavior. The role of bound water is analyzed by comparing implicit solvent-based and explicit solvent-based simulations. Loops fail to form catalytically competent geometry in absence of water. The present result, for the first time reveals the nature of the active site loop dynamics in aaRS and their influence on catalysis.

Two strategies to engineer flexible loops for improved enzyme thermostability

Flexible sites are potential targets for engineering the stability of enzymes. Nevertheless, the success rate of the rigidifying flexible sites (RFS) strategy is still low due to a limited understanding of how to determine the best mutation candidates. In this study, two parallel strategies were applied to identify mutation candidates within the flexible loops of Escherichia coli transketolase (TK). The first was a “back to consensus mutations” approach, and the second was computational design based on ΔΔG calculations in Rosetta. Forty-nine single variants were generated and characterised experimentally. From these, three single-variants I189H, A282P, D143K were found to be more thermostable than wild-type TK. The combination of A282P with H192P, a variant constructed previously, resulted in the best all-round variant with a 3-fold improved half-life at 60 °C, 5-fold increased specific activity at 65 °C, 1.3-fold improved k_cat and a T_m increased by 5 °C above that of wild type. Based on a statistical analysis of the stability changes for all variants, the qualitative prediction accuracy of the Rosetta program reached 65.3%. Both of the two strategies investigated were useful in guiding mutation candidates to flexible loops, and had the potential to be used for other enzymes.

Network-based modelling and percolation analysis of conformational dynamics and activation in the CDK2 and CDK4 proteins: dynamic and energetic polarization of the kinase lobes may determine divergence of the regulatory mechanisms

Network modeling and percolation analysis of conformational dynamics and energetics of regulatory mechanisms in cyclin-dependent kinases.

Protein structure based prediction of catalytic residues

BACKGROUND

Worldwide structural genomics projects continue to release new protein structures at an unprecedented pace, so far nearly 6000, but only about 60% of these proteins have any sort of functional annotation.

RESULTS

We explored a range of features that can be used for the prediction of functional residues given a known three-dimensional structure. These features include various centrality measures of nodes in graphs of interacting residues: closeness, betweenness and page-rank centrality. We also analyzed the distance of functional amino acids to the general center of mass (GCM) of the structure, relative solvent accessibility (RSA), and the use of relative entropy as a measure of sequence conservation. From the selected features, neural networks were trained to identify catalytic residues. We found that using distance to the GCM together with amino acid type provide a good discriminant function, when combined independently with sequence conservation. Using an independent test set of 29 annotated protein structures, the method returned 411 of the initial 9262 residues as the most likely to be involved in function. The output 411 residues contain 70 of the annotated 111 catalytic residues. This represents an approximately 14-fold enrichment of catalytic residues on the entire input set (corresponding to a sensitivity of 63% and a precision of 17%), a performance competitive with that of other state-of-the-art methods.

CONCLUSIONS

We found that several of the graph based measures utilize the same underlying feature of protein structures, which can be simply and more effectively captured with the distance to GCM definition. This also has the added the advantage of simplicity and easy implementation. Meanwhile sequence conservation remains by far the most influential feature in identifying functional residues. We also found that due the rapid changes in size and composition of sequence databases, conservation calculations must be recalibrated for specific reference databases.

Coupling between catalytic loop motions and enzyme global dynamics

Catalytic loop motions facilitate substrate recognition and binding in many enzymes. While these motions appear to be highly flexible, their functional significance suggests that structure-encoded preferences may play a role in selecting particular mechanisms of motions. We performed an extensive study on a set of enzymes to assess whether the collective/global dynamics, as predicted by elastic network models (ENMs), facilitates or even defines the local motions undergone by functional loops. Our dataset includes a total of 117 crystal structures for ten enzymes of different sizes and oligomerization states. Each enzyme contains a specific functional/catalytic loop (10–21 residues long) that closes over the active site during catalysis. Principal component analysis (PCA) of the available crystal structures (including apo and ligand-bound forms) for each enzyme revealed the dominant conformational changes taking place in these loops upon substrate binding. These experimentally observed loop reconfigurations are shown to be predominantly driven by energetically favored modes of motion intrinsically accessible to the enzyme in the absence of its substrate. The analysis suggests that robust global modes cooperatively defined by the overall enzyme architecture also entail local components that assist in suitable opening/closure of the catalytic loop over the active site.

Protein loops have critical roles in ligand binding and catalysis. An unresolved issue in this context is the extent to which the intrinsic dynamics of proteins predispose loops to perform their molecular function. In this work, we (i) critically examine the structural changes undergone by functional/catalytic loops based on a set of enzyme crystal structures in the presence/absence of a ligand, and (ii) examine to what extent those motions are correlated with, or driven by, the global modes that are predictable using simplified, physics-based models. Using a dataset of 117 structures for ten enzymes of different sizes and oligomerization states, we show that the collective modes defined by the protein topology favor loop rearrangements in reasonable agreement with those experimentally observed upon activation. These results suggest that simple but robust motions encoded by the entire architecture, not the local binding site only, assist in binding of the ligand, positioning of the catalytic loop, and/or sequestration of the catalytic site, which in turn, enable efficient catalysis.

Cooperativity in monomeric enzymes with single ligand-binding sites

Cooperativity is widespread in biology. It empowers a variety of regulatory mechanisms and impacts both the kinetic and thermodynamic properties of macromolecular systems. Traditionally, cooperativity is viewed as requiring the participation of multiple, spatially distinct binding sites that communicate via ligand-induced structural rearrangements; however, cooperativity requires neither multiple ligand binding events nor multimeric assemblies. An underappreciated manifestation of cooperativity has been observed in the non-Michaelis-Menten kinetic response of certain monomeric enzymes that possess only a single ligand-binding site. In this review, we present an overview of kinetic cooperativity in monomeric enzymes. We discuss the primary mechanisms postulated to give rise to monomeric cooperativity and highlight modern experimental methods that could offer new insights into the nature of this phenomenon. We conclude with an updated list of single subunit enzymes that are suspected of displaying cooperativity, and a discussion of the biological significance of this unique kinetic response.

A comparison of the dynamics of pantothenate synthetase from M. tuberculosis and E. coli: computational studies

Pantothenate synthetase (PS) catalyzes the final step of the pantothenate pathway, in which pantothenate is formed from pantoate and β-alanine in an ATP-dependent reaction. Mycobacterium tuberculosis PS (MTB PS) is functionally a dimer and a potential target for novel antitubercular drugs. Molecular dynamics simulations show that the functional dynamics of the enzyme are dominated by motions of a flexible gate loop in the N-terminal domain and of the C-terminal domain. The gate loop motions dominate in MTB PS while the C-terminal domain motion dominates in Escherichia coli PS. Simulations also show that the correlated motions of the domains are severely compromised in the monomeric forms. Mutations that reduce the mobility of the gate loop in MTB PS and increased it in E. coli PS were designed and validated through simulations.

Structural basis for evolution of product diversity in soybean glutathione biosynthesis

The redox active peptide glutathione is ubiquitous in nature, but some plants also synthesize glutathione analogs in response to environmental stresses. To understand the evolution of chemical diversity in the closely related enzymes homoglutathione synthetase (hGS) and glutathione synthetase (GS), we determined the structures of soybean (Glycine max) hGS in three states: apoenzyme, bound to gamma-glutamylcysteine (gammaEC), and with hGSH, ADP, and a sulfate ion bound in the active site. Domain movements and rearrangement of active site loops change the structure from an open active site form (apoenzyme and gammaEC complex) to a closed active site form (hGSH*ADP*SO(4)(2-) complex). The structure of hGS shows that two amino acid differences in an active site loop provide extra space to accommodate the longer beta-Ala moiety of hGSH in comparison to the glycinyl group of glutathione. Mutation of either Leu-487 or Pro-488 to an Ala improves catalytic efficiency using Gly, but a double mutation (L487A/P488A) is required to convert the substrate preference of hGS from beta-Ala to Gly. These structures, combined with site-directed mutagenesis, reveal the molecular changes that define the substrate preference of hGS, explain the product diversity within evolutionarily related GS-like enzymes, and reinforce the critical role of active site loops in the adaptation and diversification of enzyme function.

Modeling of solvent-dependent conformational transitions in Burkholderia cepacia lipase

BACKGROUND

The characteristic of most lipases is the interfacial activation at a lipid interface or in non-polar solvents. Interfacial activation is linked to a large conformational change of a lid, from a closed to an open conformation which makes the active site accessible for substrates. While for many lipases crystal structures of the closed and open conformation have been determined, the pathway of the conformational transition and possible bottlenecks are unknown. Therefore, molecular dynamics simulations of a closed homology model and an open crystal structure of Burkholderia cepacia lipase in water and toluene were performed to investigate the influence of solvents on structure, dynamics, and the conformational transition of the lid.

RESULTS

The conformational transition of B. cepacia lipase was dependent on the solvent. In simulations of closed B. cepacia lipase in water no conformational transition was observed, while in three independent simulations of the closed lipase in toluene the lid gradually opened during the first 10-15 ns. The pathway of conformational transition was accessible and a barrier was identified, where a helix prevented the lid from opening to the completely open conformation. The open structure in toluene was stabilized by the formation of hydrogen bonds.In simulations of open lipase in water, the lid closed slowly during 30 ns nearly reaching its position in the closed crystal structure, while a further lid opening compared to the crystal structure was observed in toluene. While the helical structure of the lid was intact during opening in toluene, it partially unfolded upon closing in water. The closing of the lid in water was also observed, when with eight intermediate structures between the closed and the open conformation as derived from the simulations in toluene were taken as starting structures. A hydrophobic beta-hairpin was moving away from the lid in all simulations in water, which was not observed in simulations in toluene. The conformational transition of the lid was not correlated to the motions of the beta-hairpin structure.

CONCLUSION

Conformational transitions between the experimentally observed closed and open conformation of the lid were observed by multiple molecular dynamics simulations of B. cepacia lipase. Transitions in both directions occurred without applying restraints or external forces. The opening and closing were driven by the solvent and independent of a bound substrate molecule.

Including Functional Annotations and Extending the Collection of Structural Classifications of Protein Loops (ArchDB)

Loops represent an important part of protein structures. The study of loop is critical for two main reasons: First, loops are often involved in protein function, stability and folding. Second, despite improvements in experimental and computational structure prediction methods, modeling the conformation of loops remains problematic. Here, we present a structural classification of loops, ArchDB, a mine of information with application in both mentioned fields: loop structure prediction and function prediction. ArchDB ( http://sbi.imim.es/archdb ) is a database of classified protein loop motifs. The current database provides four different classification sets tailored for different purposes. ArchDB-40, a loop classification derived from SCOP40, well suited for modeling common loop motifs. Since features relevant to loop structure or function can be more easily determined on well-populated clusters, we have developed ArchDB-95, a loop classification derived from SCOP95. This new classification set shows a ~40% increase in the number of subclasses, and a large 7-fold increase in the number of putative structure/function-related subclasses. We also present ArchDB-EC, a classification of loop motifs from enzymes, and ArchDB-KI, a manually annotated classification of loop motifs from kinases. Information about ligand contacts and PDB sites has been included in all classification sets. Improvements in our classification scheme are described, as well as several new database features, such as the ability to query by conserved annotations, sequence similarity, or uploading 3D coordinates of a protein. The lengths of classified loops range between 0 and 36 residues long. ArchDB offers an exhaustive sampling of loop structures. Functional information about loops and links with related biological databases are also provided. All this information and the possibility to browse/query the database through a web-server outline an useful tool with application in the comparative study of loops, the analysis of loops involved in protein function and to obtain templates for loop modeling.

Comparative analysis of structural and dynamic properties of the loaded and unloaded hemophore HasA: functional implications

A heme-acquisition system present in several Gram-negative bacteria requires the secretion of hemophores. These extracellular carrier proteins capture heme and deliver it to specific outer membrane receptors. The Serratia marcescens HasA hemophore is a monodomain protein that binds heme with a very high affinity. Its alpha/beta structure, as that of its binding pocket, has no common features with other iron- or heme-binding proteins. Heme is held by two loops L1 and L2 and coordinated to iron by an unusual ligand pair, H32/Y75. Two independent regions of the hemophore beta-sheet are involved in HasA-HasR receptor interaction. Here, we report the 3-D NMR structure of apoHasA and the backbone dynamics of both loaded and unloaded hemophore. While the overall structure of HasA is very similar in the apo and holo forms, the hemophore presents a transition from an open to a closed form upon ligand binding, through a large movement, of up to 30 A, of loop L1 bearing H32. Comparison of loaded and unloaded HasA dynamics on different time scales reveals striking flexibility changes in the binding pocket. We propose a mechanism by which these structural and dynamic features provide the dual function of heme binding and release to the HasR receptor.

Substrate-induced conformational changes and dynamics of UDP-N-acetylgalactosamine:polypeptide N-acetylgalactosaminyltransferase-2

O-Glycan biosynthesis is initiated by the transfer of N-acetylgalactosamine (GalNAc) from a nucleotide sugar donor (UDP-GalNAc) to Ser/Thr residues of an acceptor substrate. The detailed transfer mechanism, catalyzed by the UDP-GalNAc polypeptide:N-acetyl-alpha-galactosaminyltransferases (ppGalNAcTs), remains unclear despite structural information available for several isoforms in complex with substrates at various stages along the catalytic pathway. We used all-atom molecular dynamics simulations with explicit solvent and counterions to study the conformational dynamics of ppGalNAcT-2 in several enzymatic states along the catalytic pathway. ppGalNAcT-2 is simulated both in the presence and in the absence of substrates and reaction products to examine the role of conformational changes in ligand binding. In multiple 40-ns-long simulations of more than 600 ns total run time, we studied systems ranging from 45,000 to 95,000 atoms. Our simulations accurately identified dynamically active regions of the protein, as previously revealed by the X-ray structures, and permitted a detailed, atomistic description of the conformational changes of loops near the active site and the characterization of the ensemble of structures adopted by the transferase complex on the transition pathway between the ligand-bound and ligand-free states. In particular, the conformational transition of a functional loop adjacent to the active site from closed (active) to open (inactive) is correlated with the rotameric state of the conserved residue W331. Analysis of water dynamics in the active site revealed that internal water molecules have an important role in enhancing the enzyme flexibility. We also found evidence that charged side chains in the active site rearrange during site opening to facilitate ligand binding. Our results are consistent with the single-displacement transfer mechanism previously proposed for ppGalNAcTs based on X-ray structures and mutagenesis data and provide new evidence for possible functional roles of certain amino acids conserved across several isoforms.

Identification of carotenoid cleavage dioxygenases from Nostoc sp. PCC 7120 with different cleavage activities

Carotenoid cleavage dioxygenases (CCDs) are a class of enzymes that oxidatively cleave carotenoids into apocarotenoids. Dioxygenases have been identified in plants and animals and produce a wide variety of cleavage products. Despite what is known about apocarotenoids in higher organisms, very little is known about apocarotenoids and CCDs in microorganisms. This study surveyed cleavage activities of ten putative carotenoid cleavage dioxygenases from five different cyanobacteria in recombinant Escherichia coli cells producing different carotenoid substrates. Three CCD homologs identified in Nostoc sp. PCC 7120 were purified, and their cleavage activities were investigated. Two of the three enzymes showed cleavage of beta,beta-carotene at the 9,10 and 15,15' positions, respectively. The third enzyme did not cleave full-length carotenoids but cleaved the apocarotenoid beta-apo-8'-carotenal at the 9,10 position. 9,10-Apocarotenoid cleavage specificity has previously not been described. The diversity of carotenoid cleavage activities identified in one cyanobacteria suggests that CCDs not only facilitate the degradation of photosynthetic pigments but generate apocarotenals with yet to be determined biological roles in microorganisms.

Is allostery an intrinsic property of all dynamic proteins?

Allostery involves coupling of conformational changes between two widely separated binding sites. The common view holds that allosteric proteins are symmetric oligomers, with each subunit existing in "at least" two conformational states with a different affinity for ligands. Recent observations such as the allosteric behavior of myoglobin, a classical example of a nonallosteric protein, call into question the existing allosteric dogma. Here we argue that all (nonfibrous) proteins are potentially allosteric. Allostery is a consequence of re-distributions of protein conformational ensembles. In a nonallosteric protein, the binding site shape may not show a concerted second-site change and enzyme kinetics may not reflect an allosteric transition. Nevertheless, appropriate ligands, point mutations, or external conditions may facilitate a population shift, leading a presumably nonallosteric protein to behave allosterically. In principle, practically any potential drug binding to the protein surface can alter the conformational redistribution. The question is its effectiveness in the redistribution of the ensemble, affecting the protein binding sites and its function. Here, we review experimental observations validating this view of protein allostery.

Evidence of a double-lid movement in Pseudomonas aeruginosa lipase: insights from molecular dynamics simulations

Pseudomonas aeruginosa lipase is a 29-kDa protein that, following the determination of its crystal structure, was postulated to have a lid that stretched between residues 125 and 148. In this paper, using molecular dynamics simulations, we propose that there exists, in addition to the above-mentioned lid, a novel second lid in this lipase. We further show that the second lid, covering residues 210–222, acts as a triggering lid for the movement of the first. We also investigate the role of hydrophobicity in the movement of the lids and show that two residues, Phe214 and Ala217, play important roles in lid movement. To our knowledge, this is the first time that a double-lid movement of the type described in our manuscript has been presented to the scientific community. This work also elucidates the interplay of hydrophobic interactions in the dynamics, and hence the function, of an enzyme.

Lipases hydrolyse long-chain fatty acid esters at water-oil interfaces through the mechanism of interfacial activation mediated by the movement of a lid subdomain that covers the active site. Studying lid movement is an area of active research in the field of protein dynamics. The lipase from Pseudomonas aeruginosa is a 29-kDa protein that was previously crystallized in the open conformation, and as expected, an approximately 20-residue lid subdomain was identified. In the present study, the authors report extensive molecular dynamics simulations of the P. aeruginosa lipase. They show that this protein has two lids covering the substrate-binding pocket. The first lid is the one proposed from the known crystal structure. The second lid, a much shorter one, lies over the binding pocket facing the first lid. Furthermore, using position-restrained simulations, these authors show that movement of the second lid may actually be a trigger for the movement of the first, and that this triggering action is driven by hydrophobic contacts between the two lids. This computational study paves a way for experimentalists to study the structure and dynamics of this protein in greater detail in order to understand coupled subdomain movements in a comprehensive fashion.

Changing the metal ion selectivity of rabbit muscle enolase by mutagenesis: effects of the G37A and G41A mutations

During the reaction catalyzed by enolase, a mobile loop, residues 36-45, closes over the active site. In order to probe the role of this loop movement in catalysis, the glycines at positions 37 and 41 of rabbit muscle enolase (beta beta) have been mutated to alanines. The mutant forms-G37A and G41A-of enolase are both active, but have altered selectivity for divalent cations. G37A, when assayed with Mg(2+), has 12% the activity of the wild type. However, it is twice as active as wild type when assayed with Mn(2+), Zn(2+), or Co(2+). G41A has 4% the activity of the wild type with Mg(2+), is more active than wild type with Co(2+), and slightly less active than wild type with Mn(2+) and Zn(2+). The kinetic isotope effect for both mutants is greater than that of the wild type with all 4 divalent cations. These results indicate that the flexibility of this loop has subtle effects on catalytic activity.

Dioxygenases without requirement for cofactors and their chemical model reaction: compulsory order ternary complex mechanism of 1H-3-hydroxy-4-oxoquinaldine 2,4-dioxygenase involving general base catalysis by histidine 251 and single-electron oxidation of the substrate dianion

1H-3-Hydroxy-4-oxoquinaldine 2,4-dioxygenase (Hod) is a cofactor-less dioxygenase belonging to the alpha/beta hydrolase fold family, catalyzing the cleavage of 1H-3-hydroxy-4-oxoquinaldine (I) and 1H-3-hydroxy-4-oxoquinoline (II) to N-acetyl- and N-formylanthranilate, respectively, and carbon monoxide. Bisubstrate steady-state kinetics and product inhibition patterns of HodC, the C69A protein variant of Hod, suggested a compulsory-order ternary-complex mechanism, in which binding of the organic substrate precedes dioxygen binding, and carbon monoxide is released first. The specificity constants, k(cat)/K(m,A) and k(cat)/K(m,O)()2, were 1.4 x 10(8) and 3.0 x 10(5) M(-1) s(-1) with I and 1.2 x 10(5) and 0.41 x 10(5) M(-1) s(-1) with II, respectively. Whereas HodC catalyzes formation of the dianion of its organic substrate prior to dioxygen binding, HodC-H251A does not, suggesting that H251, which aligns with the histidine of the catalytic triad of the alpha/beta hydrolases, acts as general base in catalysis. Investigation of base-catalyzed dioxygenolysis of I by electron paramagnetic resonance (EPR) spectroscopy revealed formation of a resonance-stabilized radical upon exposure to dioxygen. Since in D(2)O spectral properties are not affected, exchangeable protons are not involved, confirming that the dianion is the reactive intermediate that undergoes single-electron oxidation. We suggest that in the ternary complex of the enzyme, direct single-electron transfer from the substrate dianion to dioxygen may occur, resulting in a radical pair. Based on the estimated spin distribution within the radical anion (observed in the model reaction of I), radical recombination may produce a C4- or C2-hydroperoxy(di)anion. Subsequent intramolecular attack would result in the 2,4-endoperoxy (di)anion that may collapse to the reaction products.

Structure and substrate-binding mechanism of human Ap4A hydrolase

Asymmetric diadenosine 5',5'''-P(1),P(4)-tetraphosphate (Ap(4)A) hydrolases play a major role in maintaining homeostasis by cleaving the metabolite diadenosine tetraphosphate (Ap(4)A) back into ATP and AMP. The NMR solution structures of the 17-kDa human asymmetric Ap(4)A hydrolase have been solved in both the presence and absence of the product ATP. The adenine moiety of the nucleotide predominantly binds in a ring stacking arrangement equivalent to that observed in the x-ray structure of the homologue from Caenorhabditis elegans. The binding site is, however, markedly divergent to that observed in the plant/pathogenic bacteria class of enzymes, opening avenues for the exploration of specific therapeutics. Binding of ATP induces substantial conformational and dynamic changes that were not observed in the C. elegans structure. In contrast to the C. elegans homologue, important side chains that play a major role in substrate binding do not have to reorient to accommodate the ligand. This may have important implications in the mechanism of substrate recognition in this class of enzymes.

GDSL family of serine esterases/lipases

GDSL esterases and lipases are hydrolytic enzymes with multifunctional properties such as broad substrate specificity and regiospecificity. They have potential for use in the hydrolysis and synthesis of important ester compounds of pharmaceutical, food, biochemical, and biological interests. This new subclass of lipolytic enzymes possesses a distinct GDSL sequence motif different from the GxSxG motif found in many lipases. Unlike the common lipases, GDSL enzymes do not have the so called nucleophile elbow. Studies show that GDSL hydrolases have a flexible active site that appears to change conformation with the presence and binding of the different substrates, much like the induced fit mechanism proposed by Koshland. Some of the GDSL enzymes have thioesterase, protease, arylesterase, and lysophospholipase activity, yet they appear to be the same protein with similar molecular weight ( approximately 22-60 kDa for most esterases), although some have multiple glycosylation sites with higher apparent molecular weight. GDSL enzymes have five consensus sequence (I-V) and four invariant important catalytic residues Ser, Gly, Asn, and His in blocks I, II, III, and V, respectively. The oxyanion structure led to a new designation of these enzymes as SGNH-hydrolase superfamily or subfamily. Phylogenetic analysis revealed that block IIA which belonged to the SGNH-hydrolases was found only in clade I. Therefore, this family of hydrolases represents a new example of convergent evolution of lipolytic enzymes. These enzymes have little sequence homology to true lipases. Another important differentiating feature of GDSL subfamily of lipolytic enzymes is that the serine-containing motif is closer to the N-terminus unlike other lipases where the GxSxG motif is near the center. Since the first classification of these subclass or subfamily of lipases as GDSL(S) hydrolase, progress has been made in determining the consensus sequence, crystal structure, active site and oxyanion residues, secondary structure, mechanism of catalysis, and understanding the conformational changes. Nevertheless, much still needs to be done to gain better understanding of in vivo biological function, 3-D structure, how this group of enzymes evolved to utilize many different substrates, and the mechanism of reactions. Protein engineering is needed to improve the substrate specificity, enantioselectivity, specific activity, thermostability, and heterologous expression in other hosts (especially food grade microorganisms) leading to eventual large scale production and applications. We hope that this review will rekindle interest among researchers and the industry to study and find uses for these unique enzymes.

Modulating Functional Loop Movements: The Role of Highly Conserved Residues in the Correlated Loop Motions

Loop flexibility in enzymes plays a vital role in correctly positioning catalytically important residues. This strong relationship between enzyme flexibility and function provides an opportunity to engineer new substrates and inhibitors. It further allows the design of site-directed mutagenesis experiments to explore enzymatic activity through the control of flexibility of a functional loop. Earlier, we described a novel mechanism in which a small loop triggers the motions of a functional loop in three enzymes (beta-1,4-galactosyltransferase, lipase, and enolase) unrelated in sequence, structure, or function. Here, we further address the question of how the interactions between various flexible loops modulate the movements of the functional loop. We examine beta-1,4-galactosyltransferase as a model system in which a Long loop undergoes a large conformational change (moves in space up to 20 A) upon substrate binding in addition to a small loop (Trp loop) that shows a considerably smaller conformational change. Our molecular-dynamics simulations carried out in implicit and explicit solvent show that, in addition to these two loops, two other neighboring loops are also highly flexible. These loops are in contact with either the Long loop or the Trp loop. Analysis of the covariance of the spatial displacement of the residues reveals that coupled motions occur only in one of these two loops. Sequence analysis indicates that loops correlated in their motions also have highly conserved residues involved in the loop-loop interactions. Further, analysis of crystal structures and simulations in explicit water open the possibility that the Trp loop that triggers the movement of the Long loop in the unbound conformation may also play the same role in the substrate-bound conformation through its contact with the conserved and correlated third loop. Our proposition is supported by the observation that four of the five conserved positions in the third loop are at the interface with the Trp loop. Evolution appears to select residues that drive the functional Long loop to a large conformational change. These observations suggest that altering selected loop-loop interactions might modulate the movements of the functional loop.