Domain motion and interdomain hot spots in a multidomain enzyme

The aim of this article is to analyze conformational changes by comparing 10 different structures of Pseudomonas aeruginosa phosphomannomutase/phosphoglucomutase (PMM/PGM), a four‐domain enzyme in which both substrate binding and catalysis require substantial movement of the C‐terminal domain. We focus on changes in interdomain and active site crevices using a method called computational solvent mapping rather than superimposing the structures. The method places molecular probes (i.e., small organic molecules containing various functional groups) around the protein to find hot spots. One of the most important hot spots is in the active site, consistent with the ability of the enzyme to bind both glucose and mannose phosphosugar substrates. The protein has eight additional hot spots at domain‐domain interfaces and hinge regions. The locations and nature of six of these hot spots vary between the open, half‐open, and closed conformers of the enzyme, in good agreement with the ligand‐induced conformational changes. In the closed structures the number of probe clusters at the hinge region significantly depends on the position of the phosphorylated oxygen in the substrate (e.g., glucose 1‐phosphate versus glucose 6‐phosphate), but the protein remains almost unchanged in terms of the overall RMSD, indicating that computational solvent mapping is a more sensitive approach to detect changes in binding sites and interdomain crevices. Focusing on multidomain proteins we show that the subresolution conformational differences revealed by the mapping are in fact significant, and present a general statistical method of analysis to determine the significance of rigid body domain movements in X‐ray structures.     

Domain motions and the open-to-closed conformational transition of an enzyme: a normal mode analysis of S-adenosyl-L-homocysteine hydrolase

The structure and fluctuations of the enzyme S-adenosyl-L-homocysteine hydrolase (SAHH) are analyzed in an effort to explain its biological function. Besides the previously identified open structure, characteristic of the substrate-free enzyme, we find two distinct structures in enzyme-inhibitor complexes, the closed and closed-twisted conformers. Both closed conformers differ from the open form by a hinge bending motion of two large domains within each subunit, which isolate the inhibitor bound in the active site from the bulk solvent. The closed-twisted form further differs from the closed form by a rigid body twist of the two-subunit dimers. The local structural fluctuations of SAHH are analyzed by performing block normal mode analysis of the tetrameric enzyme in its three forms. For the open form, we find that the four lowest-frequency normal modes, corresponding to the collective motions of the protein with the largest amplitudes, are essentially combinations of the hinge bending deformations of the individual subunits. Thus, the mechanical properties of the open structure of SAHH lead to the presence of structural fluctuations in the direction of the open-to-closed conformational transition. A candidate for such a motion has been observed in previous fluorescence depolarization studies of the enzyme. Both structural and normal mode analyses indicate that residues 180-190 and 350-356 form hinge regions, connecting large domains which tend to move as rigid bodies in response to interactions with substrate, intermediates, and the product of the enzymatic reactions. We propose that these hinge regions play a crucial role in the enzymatic mechanism of SAHH. In contrast to the open form, normal mode calculations for the closed conformations show strong coupling of the hinge bending motions of the individual subunits to each other and to other low-frequency vibrations. Thus, information about structural changes related to reaction progress in one active site may be mechanically transmitted to other subunits of the protein, explaining the cooperativity found in the enzyme kinetics.     

An interaction-based analysis of calcium-induced conformational changes in Ca2+ sensor proteins.

Calcium sensor proteins translate transient increases in intracellular calcium levels into metabolic or mechanical responses, by undergoing dramatic conformational changes upon Ca2+ binding. A detailed analysis of the calcium binding-induced conformational changes in the representative calcium sensors calmodulin (CaM) and troponin C was performed to obtain insights into the underlying molecular basis for their response to the binding of calcium. Distance difference matrices, analysis of interresidue contacts, comparisons of interhelical angles, and inspection of structures using molecular graphics were used to make unbiased comparisons of the various structures. The calcium-induced conformational changes in these proteins are dominated by reorganization of the packing of the four helices within each domain. Comparison of the closed and open conformations confirms that calcium binding causes opening within each of the EF-hands. A secondary analysis of the conformation of the C-terminal domain of CaM (CaM-C) clearly shows that CaM-C occupies a closed conformation in the absence of calcium that is distinct from the semi-open conformation observed in the C-terminal EF-hand domains of myosin light chains. These studies provide insight into the structural basis for these changes and into the differential response to calcium binding of various members of the EF-hand calcium-binding protein family. Factors contributing to the stability of the Ca2+-loaded open conformation are discussed, including a new hypothesis that critical hydrophobic interactions stabilize the open conformation in Ca2+ sensors, but are absent in "non-sensor" proteins that remain closed upon Ca2+ binding. A role for methionine residues in stabilizing the open conformation is also proposed.     

Interdomain dynamics and ligand binding: molecular dynamics simulations of glutamine binding protein

Periplasmic binding proteins from Gram-negative bacteria possess a common architecture, comprised of two domains linked by a hinge region, a fold which they share with the neurotransmitter-binding domains of ionotropic glutamate receptors (GluRs). Glutamine-binding protein (GlnBP) is one such protein, whose crystal structure has been solved in both open and closed forms. Multi-nanosecond molecular dynamics simulations have been used to explore motions about the hinge region and how they are altered by ligand binding. Glutamine binding is seen to significantly reduce inter-domain motions about the hinge region. Essential dynamics analysis of inter-domain motion revealed the presence of both hinge-bending and twisting motions, as has been reported for a related sugar-binding protein. Significantly, the influence of the ligand on GlnBP dynamics is similar to that previously observed in simulations of rat glutamate receptor (GluR2) ligand-binding domain. The essential dynamics analysis of GlnBP also revealed a third class of motion which suggests a mechanism for signal transmission in GluRs.     

Conformational Changes and Ligand Recognition of Escherichia colid-Xylose Binding Protein Revealed

ATP binding cassette transport systems account for most import of necessary nutrients in bacteria. The periplasmic binding component (or an equivalent membrane-anchored protein) is critical to recognizing cognate ligand and directing it to the appropriate membrane permease. Here we report the X-ray structures of d-xylose binding protein from Escherichia coli in ligand-free open form, ligand-bound open form, and ligand-bound closed form at 2.15 Å, 2.2 Å, and 2.2 Å resolutions, respectively. The ligand-bound open form is the first such structure to be reported at high resolution; the combination of the three different forms from the same protein furthermore gives unprecedented details concerning the conformational changes involved in binding protein function. As is typical of the structural family, the protein has two similar globular domains, which are connected by a three-stranded hinge region. The open liganded structure shows that xylose binds first to the C-terminal domain, with only very small conformational changes resulting. After a 34° closing motion, additional interactions are formed with the N-terminal domain; changes in this domain are larger and serve to make the structure more ordered near the ligand. An analysis of the interactions suggests why xylose is the preferred ligand. Furthermore, a comparison with the most closely related proteins in the structural family shows that the conformational changes are distinct in each type of binding protein, which may have implications for how the individual proteins act in concert with their respective membrane permeases.     

Mechanisms of domain closure in proteins

Certain enzymes respond to the binding of substrates and coenzymes by the closure of an active site that lies in a cleft between two domains. We have examined the mechanism of the domain closure in citrate synthase, for which atomic co-ordinates are available for "open" and "closed" forms. We show that the mechanism of domain closure involves small shifts and rotations of packed helices within the two domains and at their interface. Large motions of distant segments of the structure are the cumulative effect of the small relative shifts in intervening pairs of packed segments. These shifts are accommodated not by changes in packing but rather by small conformational changes in side-chains. We call this the helix interface shear mechanism of domain closure. The relative movements of packed helices follow the principles suggested by our recent study of insulin. This mechanism of domain closure is quite different from the hinge mechanisms that allow the rigid body movements of domains in immunoglobulins. The large interface between the domains of citrate synthase precludes a simple hinge mechanism for its conformational change. The helix interface shear mechanism of conformational change occurs in other enzymes that contain extensive domain-domain interfaces.     

Trapping a 96 degrees domain rotation in two distinct conformations by engineered disulfide bridges

Engineering disulfide bridges is a common technique to lock a protein movement in a defined conformational state. We have designed two double mutants of Escherichia coli 5'-nucleotidase to trap the enzyme in both an open (S228C, P513C) and a closed (P90C, L424C) conformation by the formation of disulfide bridges. The mutant proteins have been expressed, purified, and crystallized, to structurally characterize the designed variants. The S228C, P513C is a double mutant crystallized in two different crystal forms with three independent conformers, which differ from each other by a rotation of up to 12 degrees of the C-terminal domain with respect to the N-terminal domain. This finding, as well as an analysis of the domain motion in the crystal, indicates that the enzyme still exhibits considerable residual domain flexibility. In the double mutant that was designed to trap the enzyme in the closed conformation, the structure analysis reveals an unexpected intermediate conformation along the 96 degrees rotation trajectory between the open and closed enzyme forms. A comparison of the five independent conformers analyzed in this study shows that the domain movement of the variant enzymes is characterized by a sliding movement of the residues of the domain interface along the interface, which is in contrast to a classical closure motion where the residues of the domain interface move perpendicular to the interface.     

Domain closure mechanism in transferrins: new viewpoints about the hinge structure and motion as deduced from high resolution crystal structures of ovotransferrin N-lobe

The crystal structure of holo hen ovotransferrin N-lobe refined at 1.65 A resolution has been obtained. The final model gave an R-factor of 0.173 in the resolution range between 10.0 and 1.65 A. The comparison of the structure with previous high-resolution apo and Fe(3+)-loaded, domain-opened intermediate structures provides new viewpoints on the domain closure mechanism upon Fe(3+) uptake in ovotransferrin N-lobe. Overall, conformational transition follows the common mechanism that has been first demonstrated for lactoferrin N-lobe; the domains 1 and 2 rotate 49.7 degrees as rigid bodies with a translation of 2.1 A around a screw-axis that passes through the two interdomain beta-strands (89-94 and 244-249). It is generally believed that the two strands display a hinge-like motion. Here, the latter strand indeed displays an ideal hinge nature: the segments 244-246 and 248-249 behave as a part of the rigid body of domain 2 and that of domain 1, respectively, and a sharp bend upon the domain closure is largely accounted for by the changes in the torsion angles phi and psi of Val247. We find, however, that the mode of the conformational change in the first beta-strand is much more complex. Two of the five inter beta-strand hydrogen bonds undergo crucial exchanges: from Ser91-N...Val247-O and Thr89-O...Ala249-N in the open apo and intermediate structures into Tyr92-N...Val247-O and Thr90-O...Ala249-N in the closed holo structure. These exchanges, which may be triggered in the intermediate state by modulation in the topological relation between the Fe(3+)-ligated hinge residue Tyr92-OH and the anion anchor residues of helix 5, are accompanied by a large conformational change and extensive hydrogen bonding rearrangements in a long stretch of segment of Glu82 to Tyr92. Such structural transition would work as a driving force for the domain closure, which highlights a "door closer"-like role, in addition to the canonical-hinge role, for the interdomain polypeptide segment pair. As an alternative hinge that secures the correct domain motion by being placed on a significant distance from the beta-strand hinge, we point out the participation of the van der Waals contacts formed between domain 1 residue of Met331 and domain 2 residues of Trp125, Ile129 and Trp140.     

Molecular dynamics simulations of the conformational changes of the glutamate receptor ligand-binding core in the presence of glutamate and kainate

Excitatory synaptic transmission is mediated by ionotropic glutamate receptors (iGluRs) through the induced transient opening of transmembrane ion channels. The three-dimensional structure of the extracellular ligand-binding core of iGluRs shares the overall features of bacterial periplasmic binding proteins (PBPs). In both families of proteins, the ligand-binding site is arranged in two domains separated by a cleft and connected by two peptide stretches. PBPs undergo a typical hinge motion of the two domains associated with ligand binding that leads to a conformational change from an open to a closed form. The common architecture suggests a similar closing mechanism in the ligand-binding core of iGluRs induced by the binding of specific agonists. Starting from the experimentally determined kainate-bound closed form of the S1S2 GluR2 construct, we have studied by means of molecular dynamics simulations the opening motion of the ligand-binding core in the presence and in the absence of both glutamate and kainate. Our results suggest that the opening/closing interdomain hinge motions are coupled to conformational changes in the insertion region of the transmembrane segments. These changes are triggered by the interaction of the agonists with the essential Glu 209 residue. A plausible mechanism for the coupling of agonist binding to channel gating is discussed.     

KcsA closed and open: modelling and simulation studies

Bacterial homologues of mammalian potassium channels provide structures of two states of a gated K channel. Thus, the crystal structure of KcsA represents a closed state whilst that of MthK represents an open state. Using homology modelling and molecular dynamics simulations we have built a model of the transmembrane domain of KcsA in an open state and have compared its conformational stability with that of the same domain of KcsA in a closed state. Approximate Born energy calculations of monovalent cations within the two KcsA channel states suggest that the intracellular hydrophobic gate in the closed state provides a barrier of height ~5 kT to ion permeation, whilst in the open state the barrier is absent. Simulations (10 ns duration) in an octane slab (a simple membrane mimetic) suggest that closed- and open-state models are of comparable conformational stability, both exhibiting conformational drifts of ~3.3 Å C RMSD relative to the respective starting models. Substantial conformational fluctuations are observed in the intracellular gate region during both simulations (closed state and open state). In the simulation of open-state KcsA, rapid (<5 ns) exit of all three K⁺ ions occurs through the intracellular mouth of the channel. Helix kink and swivel motion is observed at the molecular hinge formed by residue G99 of the M2 helix. This motion is more substantial for the open- than for the closed-state model of the channel.     

Molecular dynamics simulations reveal that apo‐HisJ can sample a closed conformation

The Escherichia coli histidine binding protein HisJ is a type II periplasmic binding protein (PBP) that preferentially binds histidine and interacts with its cytoplasmic membrane ABC transporter, HisQMP₂, to initiate histidine transport. HisJ is a bilobal protein where the larger Domain 1 is connected to the smaller Domain 2 via two linking strands. Type II PBPs are thought to undergo “Venus flytrap” movements where the protein is able to reversibly transition from an open to a closed conformation. To explore the accessibility of the closed conformation to the apo state of the protein, we performed a set of all‐atom molecular dynamics simulations of HisJ starting from four different conformations: apo‐open, apo‐closed, apo‐semiopen, and holo‐closed. The simulations reveal that the closed conformation is less dynamic than the open one. HisJ experienced closing motions and explored semiopen conformations that reverted to closed forms resembling those found in the holo‐closed state. Essential dynamics analysis of the simulations identified domain closing/opening and twisting as main motions. The formation of specific inter‐hinge strand and interdomain polar interactions contributed to the adoption of the closed apo‐conformations although they are up to 2.5‐fold less prevalent compared with the holo‐closed simulations. The overall sampling of the closed form by apo‐HisJ provides a rationale for the binding of unliganded PBPs with their cytoplasmic membrane ABC transporters. Proteins 2014; 82:386–398. © 2013 Wiley Periodicals, Inc.     

The crystal structure of iron-free human serum transferrin provides insight into inter-lobe communication and receptor binding

Serum transferrin reversibly binds iron in each of two lobes and delivers it to cells by a receptor-mediated, pH-dependent process. The binding and release of iron result in a large conformational change in which two subdomains in each lobe close or open with a rigid twisting motion around a hinge. We report the structure of human serum transferrin (hTF) lacking iron (apo-hTF), which was independently determined by two methods: 1) the crystal structure of recombinant non-glycosylated apo-hTF was solved at 2.7-A resolution using a multiple wavelength anomalous dispersion phasing strategy, by substituting the nine methionines in hTF with selenomethionine and 2) the structure of glycosylated apo-hTF (isolated from serum) was determined to a resolution of 2.7A by molecular replacement using the human apo-N-lobe and the rabbit holo-C1-subdomain as search models. These two crystal structures are essentially identical. They represent the first published model for full-length human transferrin and reveal that, in contrast to family members (human lactoferrin and hen ovotransferrin), both lobes are almost equally open: 59.4 degrees and 49.5 degrees rotations are required to open the N- and C-lobes, respectively (compared with closed pig TF). Availability of this structure is critical to a complete understanding of the metal binding properties of each lobe of hTF; the apo-hTF structure suggests that differences in the hinge regions of the N- and C-lobes may influence the rates of iron binding and release. In addition, we evaluate potential interactions between apo-hTF and the human transferrin receptor.     

Hinge-bending motion in citrate synthase arising from normal mode calculations

A normal mode analysis of the closed form of dimeric citrate synthase has been performed. The largest-amplitude collective motion predicted by this method compares well with the crystallographically observed hinge-bending motion. Such a result supports those obtained previously in the case of hinge-bending motions of smaller systems, such as lysozyme or hexokinase. Taken together, all these results suggest that low-frequency normal modes may become useful for determining a first approximation of the conformational path between the closed and open forms of these proteins. © 1995 Wiley-Liss, Inc.     

Pivot residue: an analysis of domain motion in proteins

In this study, we present an approach to identify some residues that represent the pivot points to experience conformational changes between open (unligand) and closed (ligand) forms of a protein. First, an angle, , formed by 4 consecutive Ca atoms in polypeptide backbones was introduced. The difference of this angle, , from the equivalent residues between the open and the closed form was used to represent the local torsion changes in the protein structure, and the residue with the maximum among was identified to be a pivot residue. We demonstrate the ability of our method by identifying the pivot residues from five proteins, Lysozyme mutates, Lactoferrin, Lay/Arg/Orn-binding protein, Calmodulin and Catabolit gene activator protein. These pivot residues are located at the hinges in the proteins, they are hinge points for the domain motion. These examples also show that the pivot residues are useful to distinguish the mechanism between shear motion and hinge motion in a protein     

Protein conformer selection by sequence-dependent packing contacts in crystals of 3-phosphoglycerate kinase

In several crystal structures of 3-phosphoglycerate kinase (PGK), the two domains occupy different relative positions. It is intriguing that the two extreme (open and closed) conformations have never been observed for the enzyme from the same species. Furthermore, in certain cases, these different crystalline conformations represent the enzyme-ligand complex of the same composition, such as the ternary complex containing either the substrate 3-phosphoglycerate (3-PG) and beta,gamma-imido-adenosine-5'-triphosphate (AMP-PNP), an analogue of the substrate MgATP, or 3-PG and the product MgADP. Thus, the protein conformation in the crystal is apparently determined by the origin of the isolated enzyme: PGK from pig muscle has only been crystallized in open conformation, whereas PGK from either Thermotoga maritima or Trypanosoma brucei has only been reported in closed conformations. A systematic analysis of the underlying sequence differences at the crucial hinge regions of the molecule and in the protein-protein contact surfaces in the crystal, in two independent pairs of open and closed states, have revealed that 1) sequential differences around the molecular hinges do not explain the appearance of fundamentally different conformations and 2) the species-specific intermolecular contacts between the nonconserved residues are responsible for stabilizing one conformation over the other in the crystalline state. A direct relationship between the steric position of the contacts in the three-dimensional structure and the conformational state of the protein has been demonstrated.     

Induced fit in guanidino kinases--comparison of substrate-free and transition state analog structures of arginine kinase

Arginine kinase (AK) is a member of the guanidino kinase family that plays an important role in buffering ATP concentration in cells with high and fluctuating energy demands. The AK specifically catalyzes the reversible phosphoryl transfer between ATP and arginine. We have determined the crystal structure of AK from the horseshoe crab (Limulus polyphemus) in its open (substrate-free) form. The final model has been refined at 2.35 A with a final R of 22.3% (R(free) = 23.7%). The structure of the open form is compared to the previously determined structure of the transition state analog complex in the closed form. Classically, the protein would be considered two domain, but dynamic domain (DynDom) analysis shows that most of the differences between the two structures can be considered as the motion between four rigid groups of nonsequential residues. ATP binds near a cluster of positively charged residues of a fixed dynamic domain. The other three dynamic domains close the active site with separate hinge rotations relative to the fixed domain. Several residues of key importance for the induced motion are conserved within the phosphagen kinase family, including creatine kinase. Substantial conformational changes are induced in different parts of the enzyme as intimate interactions are formed with both substrates. Thus, although induced fit occurs in a number of phosphoryl transfer enzymes, the conformational changes in phosphagen kinases appear to be more complicated than in prior examples.     

Conformational plasticity of the catalytic subunit of protein kinase CK2 and its consequences for regulation and drug design

At the first glance CK2α, the catalytic subunit of protein kinase CK2, is a rigid molecule: in contrast to many eukaryotic protein kinases in CK2α the canonical regulatory key elements like the activation segment occur exclusively in their typical active conformations. This observation fits well to the constitutive activity of the enzyme, meaning, its independence from phosphorylation or other characteristic control factors. Most CK2α structures are based on the enzyme from Zea mays, supplemented by an increasing number of human CK2α structures. In the latter a surprising plasticity of important ATP-binding elements – the interdomain hinge region and the glycine-rich loop – was discovered. In fully active CK2α the hinge region is open and does not anchor the ATP ribose, but alternatively it can adopt a closed conformation, form hydrogen bonds to the ribose moiety and thus retract the γ-phospho group from its functional position. In addition to this partially inactive state human CK2α was recently found in a fully inactive conformation. It is incompatible with ATP-binding due to a combination of a closed hinge and a collapse of the glycine-rich loop into the ATP cavity. These conformational transitions are apparently correlated with the occupation state of a remote docking site located at the interface to the non-catalytic subunit CK2β: if CK2β blocks this site, the fully active conformation of CK2α is stabilized, while the binding of certain small molecule seems to favour the partially and fully inactive states. This observation may be exploited to design effective and selective CK2 inhibitors.     

A mass weighted chemical elastic network model elucidates closed form domain motions in proteins

An elastic network model (ENM), usually Cα coarse‐grained one, has been widely used to study protein dynamics as an alternative to classical molecular dynamics simulation. This simple approach dramatically saves the computational cost, but sometimes fails to describe a feasible conformational change due to unrealistically excessive spring connections. To overcome this limitation, we propose a mass‐weighted chemical elastic network model (MWCENM) in which the total mass of each residue is assumed to be concentrated on the representative alpha carbon atom and various stiffness values are precisely assigned according to the types of chemical interactions. We test MWCENM on several well‐known proteins of which both closed and open conformations are available as well as three α‐helix rich proteins. Their normal mode analysis reveals that MWCENM not only generates more plausible conformational changes, especially for closed forms of proteins, but also preserves protein secondary structures thus distinguishing MWCENM from traditional ENMs. In addition, MWCENM also reduces computational burden by using a more sparse stiffness matrix.