A comprehensive and non-redundant database of protein domain movements

Collective dynamics of periplasmic glutamine binding protein upon domain closure

The glutamine binding protein is a vital component of the associated ATP binding cassette transport systems responsible for the uptake of glutamine into the cell. We have investigated the global movements of this protein by molecular dynamics simulations and principal component analysis (PCA). We confirm that the most dominant mode corresponds to the biological function of the protein, i.e., a hinge-type motion upon ligand binding. The closure itself was directly observed from two independent trajectories whereby PCA was used to elucidate the nature of this closing reaction. Two intermediary states are identified and described in detail. The ligand binding induces the structural change of the hinge regions from a discontinuous beta-sheet to a continuous one, which also enhances softness of the hinge and modifies the direction of hinge motion to enable closing. We also investigated the convergence behavior of PCA modes, which were found to converge rather quickly when the associated magnitudes of the eigenvalues are well separated.

A method for the analysis of domain movements in large biomolecular complexes

A new method for the analysis of domain movements in large, multichain, biomolecular complexes is presented. The method is applicable to any molecule for which two atomic structures are available that represent a conformational change indicating a possible domain movement. The method is blind to atomic bonding and atom type and can, therefore, be applied to biomolecular complexes containing different constituent molecules such as protein, RNA, or DNA. At the heart of the method is the use of blocks located at grid points spanning the whole molecule. The rotation vector for the rotation of atoms from each block between the two conformations is calculated. Treating components of these vectors as coordinates means that each block is associated with a point in a "rotation space" and that blocks with atoms that rotate together, perhaps as part of the same rigid domain, will have colocated points. Thus a domain can be identified from the clustering of points from blocks that span it. Domain pairs are accepted for analysis of their relative movements in terms of screw axes based upon a set of reasonable criteria. Here, we report on the application of the method to biomolecules covering a considerable size range: hemoglobin, liver alcohol dehydrogenase, S-Adenosylhomocysteine hydrolase, aspartate transcarbamylase, and the 70S ribosome. The results provide a depiction of the conformational change within each molecule that is easily understood, giving a perspective that is expected to lead to new insights. Of particular interest is the allosteric mechanism in some of these molecules. Results indicate that common boundaries between subunits and domains are good regions to focus on as movement in one subunit can be transmitted to another subunit through such interfaces.

Database of ligand-induced domain movements in enzymes

BACKGROUND

Conformational change induced by the binding of a substrate or coenzyme is a poorly understood stage in the process of enzyme catalysed reactions. For enzymes that exhibit a domain movement, the conformational change can be clearly characterized and therefore the opportunity exists to gain an understanding of the mechanisms involved. The development of the non-redundant database of protein domain movements contains examples of ligand-induced domain movements in enzymes, but this valuable data has remained unexploited.

DESCRIPTION

The domain movements in the non-redundant database of protein domain movements are those found by applying the DynDom program to pairs of crystallographic structures contained in Protein Data Bank files. For each pair of structures cross-checking ligands in their Protein Data Bank files with the KEGG-LIGAND database and using methods that search for ligands that contact the enzyme in one conformation but not the other, the non-redundant database of protein domain movements was refined down to a set of 203 enzymes where a domain movement is apparently triggered by the binding of a functional ligand. For these cases, ligand binding information, including hydrogen bonds and salt-bridges between the ligand and specific residues on the enzyme is presented in the context of dynamical information such as the regions that form the dynamic domains, the hinge bending residues, and the hinge axes.

CONCLUSION

The presentation at a single website of data on interactions between a ligand and specific residues on the enzyme alongside data on the movement that these interactions induce, should lead to new insights into the mechanisms of these enzymes in particular, and help in trying to understand the general process of ligand-induced domain closure in enzymes. The website can be found at: http://www.cmp.uea.ac.uk/dyndom/enzymeList.do.

The crystal structures of the open and catalytically competent closed conformation of Escherichia coli glycogen synthase

Escherichia coli glycogen synthase (EcGS, EC 2.4.1.21) is a retaining glycosyltransferase (GT) that transfers glucose from adenosine diphosphate glucose to a glucan chain acceptor with retention of configuration at the anomeric carbon. EcGS belongs to the GT-B structural superfamily. Here we report several EcGS x-ray structures that together shed considerable light on the structure and function of these enzymes. The structure of the wild-type enzyme bound to ADP and glucose revealed a 15.2 degrees overall domain-domain closure and provided for the first time the structure of the catalytically active, closed conformation of a glycogen synthase. The main chain carbonyl group of His-161, Arg-300, and Lys-305 are suggested by the structure to act as critical catalytic residues in the transglycosylation. Glu-377, previously thought to be catalytic is found on the alpha-face of the glucose and plays an electrostatic role in the active site and as a glucose ring locator. This is also consistent with the structure of the EcGS(E377A)-ADP-HEPPSO complex where the glucose moiety is either absent or disordered in the active site.

Dihydroorotase from the hyperthermophile Aquifex aeolicus is activated by stoichiometric association with aspartate transcarbamoylase and forms a one-pot reactor for pyrimidine biosynthesis

In prokaryotes, the first three enzymes in pyrimidine biosynthesis, carbamoyl phosphate synthetase (CPS), aspartate transcarbamoylase (ATC), and dihydroorotase (DHO), are commonly expressed separately and either function independently (Escherichia coli) or associate into multifunctional complexes (Aquifex aeolicus). In mammals the enzymes are expressed as a single polypeptide chain (CAD) in the order CPS-DHO-ATC and associate into a hexamer. This study presents the three-dimensional structure of the noncovalent hexamer of DHO and ATC from the hyperthermophile A. aeolicus at 2.3 A resolution. It is the first structure of any multienzyme complex in pyrimidine biosynthesis and is a possible model for the core of mammalian CAD. The structure has citrate, a near isosteric analogue of carbamoyl aspartate, bound to the active sites of both enzymes. Three active site loops that are intrinsically disordered in the free, inactive DHO are ordered in the complex. The reorganization also changes the peptide bond between Asp153, a ligand of the single zinc atom in DHO, and Gly154, to the rare cis conformation. In the crystal structure, six DHO and six ATC chains form a hollow dodecamer, in which the 12 active sites face an internal reaction chamber that is approximately 60 A in diameter and connected to the cytosol by narrow tunnels. The entrances and the interior of the chamber are both electropositive, which suggests that the architecture of this nanoreactor modifies the kinetics of the bisynthase, not only by steric channeling but also by preferential escape of the product, dihydroorotase, which is less negatively charged than its precursors, carbamoyl phosphate, aspartate, or carbamoyl aspartate.

Principal component analysis of native ensembles of biomolecular structures (PCA_NEST): insights into functional dynamics

MOTIVATION

To efficiently analyze the 'native ensemble of conformations' accessible to proteins near their folded state and to extract essential information from observed distributions of conformations, reliable mathematical methods and computational tools are needed.

RESULT

Examination of 24 pairs of structures determined by both NMR and X-ray reveals that the differences in the dynamics of the same protein resolved by the two techniques can be tracked to the most robust low frequency modes elucidated by principal component analysis (PCA) of NMR models. The active sites of enzymes are found to be highly constrained in these PCA modes. Furthermore, the residues predicted to be highly immobile are shown to be evolutionarily conserved, lending support to a PCA-based identification of potential functional sites. An online tool, PCA_NEST, is designed to derive the principal modes of conformational changes from structural ensembles resolved by experiments or generated by computations.

AVAILABILITY

http://ignm.ccbb.pitt.edu/oPCA_Online.htm

Dynamic energy landscape view of coupled binding and protein conformational change: induced-fit versus population-shift mechanisms

Allostery, the coupling between ligand binding and protein conformational change, is the heart of biological network and it has often been explained by two representative models, the induced-fit and the population-shift models. Here, we clarified for what systems one model fits better than the other by performing molecular simulations of coupled binding and conformational change. Based on the dynamic energy landscape view, we developed an implicit ligand-binding model combined with the double-basin Hamiltonian that describes conformational change. From model simulations performed for a broad range of parameters, we uncovered that each of the two models has its own range of applicability, stronger and longer-ranged interaction between ligand and protein favors the induced-fit model, and weaker and shorter-ranged interaction leads to the population-shift model. We further postulate that the protein binding to small ligand tends to proceed via the population-shift model, whereas the protein docking to macromolecules such as DNA tends to fit the induced-fit model.

pFlexAna: detecting conformational changes in remotely related proteins

The pFlexAna (protein flexibility analyzer) web server detects and displays conformational changes in remotely related proteins, without relying on sequence homology. To do so, it first applies a reliable statistical test to align core protein fragments that are structurally similar and then clusters these aligned fragment pairs into 'super-alignments', according to the similarity of geometric transformations that align them. The result is that the dominant conformational changes occur between the clusters, while the smaller conformational changes occur within a cluster. pFlexAna is available at http://bigbird.comp.nus.edu.sg/pfa2/.

Leishmania trypanothione synthetase-amidase structure reveals a basis for regulation of conflicting synthetic and hydrolytic activities

The bifunctional trypanothione synthetase-amidase catalyzes biosynthesis and hydrolysis of the glutathione-spermidine adduct trypanothione, the principal intracellular thiol-redox metabolite in parasitic trypanosomatids. These parasites are unique with regard to their reliance on trypanothione to determine intracellular thiol-redox balance in defense against oxidative and chemical stress and to regulate polyamine levels. Enzymes involved in trypanothione biosynthesis provide essential biological activities, and those absent from humans or for which orthologues are sufficiently distinct are attractive targets to underpin anti-parasitic drug discovery. The structure of Leishmania major trypanothione synthetase-amidase, determined in three crystal forms, reveals two catalytic domains. The N-terminal domain, a cysteine, histidine-dependent amidohydrolase/peptidase amidase, is a papain-like cysteine protease, and the C-terminal synthetase domain displays an ATP-grasp family fold common to C:N ligases. Modeling of substrates into each active site provides insight into the specificity and reactivity of this unusual enzyme, which is able to catalyze four reactions. The domain orientation is distinct from that observed in a related bacterial glutathionylspermidine synthetase. In trypanothione synthetase-amidase, the interactions formed by the C terminus, binding in and restricting access to the amidase active site, suggest that the balance of ligation and hydrolytic activity is directly influenced by the alignment of the domains with respect to each other and implicate conformational changes with amidase activity. The potential inhibitory role of the C terminus provides a mechanism to control relative levels of the critical metabolites, trypanothione, glutathionylspermidine, and spermidine in Leishmania.

HingeProt: automated prediction of hinges in protein structures

Proteins are highly flexible molecules. Prediction of molecular flexibility aids in the comprehension and prediction of protein function and in providing details of functional mechanisms. The ability to predict the locations, directions, and extent of molecular movements can assist in fitting atomic resolution structures to low-resolution EM density maps and in predicting the complex structures of interacting molecules (docking). There are several types of molecular movements. In this work, we focus on the prediction of hinge movements. Given a single protein structure, the method automatically divides it into the rigid parts and the hinge regions connecting them. The method employs the Elastic Network Model, which is very efficient and was validated against a large data set of proteins. The output can be used in applications such as flexible protein-protein and protein-ligand docking, flexible docking of protein structures into cryo-EM maps, and refinement of low-resolution EM structures. The web server of HingeProt provides convenient visualization of the results and is available with two mirror sites at http://www.prc.boun.edu.tr/appserv/prc/HingeProt3 and http://bioinfo3d.cs.tau.ac.il/HingeProt/.

Three-dimensional structures of apo- and holo-L-alanine dehydrogenase from Mycobacterium tuberculosis reveal conformational changes upon coenzyme binding

L-alanine dehydrogenase from Mycobacterium tuberculosis catalyzes the NADH-dependent reversible conversion of pyruvate and ammonia to L-alanine. Expression of the gene coding for this enzyme is up-regulated in the persistent phase of the organism, and alanine dehydrogenase is therefore a potential target for pathogen control by antibacterial compounds. We have determined the crystal structures of the apo- and holo-forms of the enzyme to 2.3 and 2.0 A resolution, respectively. The enzyme forms a hexamer of identical subunits, with the NAD-binding domains building up the core of the molecule and the substrate-binding domains located at the apical positions of the hexamer. Coenzyme binding stabilizes a closed conformation where the substrate-binding domains are rotated by about 16 degrees toward the dinucleotide-binding domains, compared to the open structure of the apo-enzyme. In the structure of the abortive ternary complex with NAD+ and pyruvate, the substrates are suitably positioned for hydride transfer between the nicotinamide ring and the C2 carbon atom of the substrate. The approach of the nucleophiles water and ammonia to pyruvate or the reaction intermediate iminopyruvate, respectively, is, however, only possible through conformational changes that make the substrate binding site more accessible. The crystal structures identified the conserved active-site residues His96 and Asp270 as potential acid/base catalysts in the reaction. Amino acid replacements of these residues by site-directed mutagenesis led to inactive mutants, further emphasizing their essential roles in the enzymatic reaction mechanism.

What is the relationship between the global structures of apo and holo proteins?

It is well known that ligand binding and release may induce a wide range of structural changes in a receptor protein, varying from small movements of loops or side chains in the binding pocket to large-scale domain hinge-bending and shear motions or even partial unfolding that facilitates the capture and release of a ligand. An interesting question is what in general are the conformational changes triggered by ligand binding? The aim of this work is analyze the magnitude of structural changes in a protein resulting from ligand binding to assess if the state of ligand binding needs to be included in template-based protein structure prediction algorithms. To address this issue, a nonredundant dataset of 521 paired protein structures in the ligand-free and ligand-bound form was created and used to estimate the degree of both local and global structure similarity between the apo and holo forms. In most cases, the proteins undergo relatively small conformational rearrangements of their tertiary structure upon ligand binding/release (most root-mean-square-deviations from native, RMSD, are <1 A). However, a clear difference was observed between single- and multiple-domain proteins. For the latter, RMSD changes greater than 1 A and sometimes larger were found for almost 1/3 of the cases; these are mainly associated with large-scale hinge-bending movements of entire domains. The changes in the mutual orientation of individual domains in multiple-domain proteins upon ligand binding were investigated using a mechanistic model based on mass-weighted principal axes as well as interface buried surface calculations. Some preferences toward the anticipated mechanism of protein domain movements are predictable based on the examination of just the ligand-free structural form. These results have applications to protein structure prediction, particularly in the context of protein domain assembly, if additional information concerning ligand binding is exploited.

Essential dynamics sampling study of adenylate kinase: Comparison to citrate synthase and implication for the hinge and shear mechanisms of domain motions

Essential dynamics sampling simulations of the domain conformations of unliganded Escherichia coli adenylate kinase have been performed to determine whether the ligand-induced closed-domain conformation is accessible to the open unliganded enzyme. Adenylate kinase is a three- domain protein with a central CORE domain and twoflanking domains, the LID and the NMPbind domains. The sampling simulations were applied to the CORE and NMPbind domain pair and the CORE and LID domain pair separately. One aim is to compare the results to those of a similar study on the enzyme citrate synthase to determine whether a similar domain-locking mechanism operates in adenylate kinase. Although for adenylate kinase the simulations suggest that the closed-domain conformation of the unliganded enzyme is at a slightly higher free energy than the open for both domain pairs, the results are radically different to those found for citrate synthase. In adenylate kinase the targeted domain conformations could always be achieved, whereas this was not the case in citrate synthase due to an apparent free-energy barrier between the open and closed conformations. Adenylate kinase has been classified as a protein that undergoes closure through a hinge mechanism, whereas citrate synthase has been assigned to the shear mechanism. This was quantified here in terms of the change in the number of interdomain contacting atoms upon closure which showed a considerable increase in adenylate kinase. For citrate synthase this number remained largely the same, suggesting that the domain faces slide over each other during closure. This suggests that shear and hinge mechanisms of domain closure may relate to the existence or absence of an appreciable barrier to closure for the unliganded protein, as the latter can hinge comparatively freely, whereas the former must follow a more constrained path. In general though it appears a bias toward keeping the unliganded enzyme in the open-domain conformation may be a common feature of domain enzymes.

Dichotomic Phylogenetic Tree of the Pyruvate Kinase Family

K+ dependence was assumed to be a feature of all pyruvate kinases until it was discovered that some enzymes express K+ -independent activity. Almost all the K+ -independent pyruvate kinases have Lys at position 117, instead of the Glu present in the K+ -dependent muscle enzyme. Mutagenesis studies show that the internal positive charge substitutes for the K+ requirement (Laughlin, L. T. & Reed, G. H. (1997) Arch. Biochem. Biophys. 348, 262-267). In this work a phylogenetic analysis of pyruvate kinase was performed to ascertain the abundance of K+ -independent activities and to explore whether the K+ activating effect is related to the evolutionary history of the enzyme. Of the 230 studied sequences, 46% have Lys at position 117, and the rest have Glu. Pyruvate kinases with Lys117 and Glu117 are separated in two clusters. All of the enzymes of the Glu117 cluster that have been characterized are K+ -dependent, whereas those of the Lys117 cluster are K+ -independent. Thus, there is a strict correlation between the dichotomy of the tree and the dependence of activity on K+. 77% of the pyruvate kinases that possess Lys117 have Lys113/Gln114; they also have Ile, Val, or Leu at position 120. These residues are replaced by Glu117 and Thr113/Lys114/Thr120 in 80% of K+ -dependent pyruvate kinases. Structural analysis indicates that these residues are in a hinge region involved in the acquisition of the catalytic conformation of the enzyme. The route of conversion from K+ -independent to K+ -dependent pyruvate kinases is described. A plausible explanation of how enzymes developed K+ dependence is put forth.

Amyloid Formation May Involve α- to β Sheet Interconversion via Peptide Plane Flipping

The toxic component of amyloid is not the mature fiber but a soluble prefibrillar intermediate. It has been proposed, from molecular dynamics simulations, that the precursor is composed of alpha sheet, which converts into the beta sheet of mature amyloid via peptide plane flipping. alpha sheet, not seen in proteins, occurs as isolated stretches of polypeptide. We show that the alpha- to beta sheet transition can occur by the flipping of alternate peptide planes. The flip can be described as alphaRalphaL<-->betabeta. A search conducted within sets of closely related protein crystal structures revealed that these flips are common, occurring in 8.5% of protein families. The average "alphaL" conformation found is in an adjacent and less populated region of the Ramachandran plot, as expected if the flanking peptide planes, being hydrogen bonded, are restricted in their movements. This work provides evidence for flips allowing direct alpha- to beta sheet interconversion.

Molecular dynamics simulations of NAD+-induced domain closure in horse liver alcohol dehydrogenase

Horse liver alcohol dehydrogenase is a homodimer, the protomer having a coenzyme-binding domain and a catalytic domain. Using all available x-ray structures and 50 ns of molecular dynamics simulations, we investigated the mechanism of NAD+-induced domain closure. When the well-known loop at the domain interface was modeled to its conformation in the closed structure, the NAD+-induced domain closure from the open structure could be simulated with remarkable accuracy. Native interactions in the closed structure between Arg369, Arg47, His51, Ala317, Phe319, and NAD+ were seen to form at different stages during domain closure. Removal of the Arg369 side-chain charge resulted in the loss of the tendency to close, verifying that specific interactions do help drive the domains closed. Further simulations and a careful analysis of x-ray structures suggest that the loop prevents domain closure in the absence of NAD+, and a cooperative mechanism operates between the subunits for domain closure. This cooperative mechanism explains the role of the loop as a block to closure because in the absence of NAD+ it would prevent the occurrence of an unliganded closed subunit when the other subunit closes on NAD+. Simulations that started with one subunit open and one closed supported this.