Mapping the Dynamics Landscape of Conformational Transitions in Enzyme: The Adenylate Kinase Case

Conformational transition describes the essential dynamics and mechanism of enzymes in pursuing their various functions. The fundamental and practical challenge to researchers is to quantitatively describe the roles of large-scale dynamic transitions for regulating the catalytic processes. In this study, we tackled this challenge by exploring the pathways and free energy landscape of conformational changes in adenylate kinase (AdK), a key ubiquitous enzyme for cellular energy homeostasis. Using explicit long-timescale (up to microseconds) molecular dynamics and bias-exchange metadynamics simulations, we determined at the atomistic level the intermediate conformational states and mapped the transition pathways of AdK in the presence and absence of ligands. There is clearly chronological operation of the functional domains of AdK. Specifically in the ligand-free AdK, there is no significant energy barrier in the free energy landscape separating the open and closed states. Instead there are multiple intermediate conformational states, which facilitate the rapid transitions of AdK. In the ligand-bound AdK, the closed conformation is energetically most favored with a large energy barrier to open it up, and the conformational population prefers to shift to the closed form coupled with transitions. The results suggest a perspective for a hybrid of conformational selection and induced fit operations of ligand binding to AdK. These observations, depicted in the most comprehensive and quantitative way to date, to our knowledge, emphasize the underlying intrinsic dynamics of AdK and reveal the sophisticated conformational transitions of AdK in fulfilling its enzymatic functions. The developed methodology can also apply to other proteins and biomolecular systems.     

Zipping and Unzipping of Adenylate Kinase: Atomistic Insights into the Ensemble of Open ↔ Closed Transitions

Adenylate kinase (AdK), a phosphotransferase enzyme, plays an important role in cellular energy homeostasis. It undergoes a large conformational change between an open and a closed state, even in the absence of substrate. We investigate the apo-AdK transition at the atomic level both with free-energy calculations and with our new dynamic importance sampling (DIMS) molecular dynamics method. DIMS is shown to sample biologically relevant conformations as verified by comparing an ensemble of hundreds of DIMS transitions to AdK crystal structure intermediates. The simulations reveal in atomic detail how hinge regions partially and intermittently unfold during the transition. Conserved salt bridges are seen to have important structural and dynamic roles; in particular, four ionic bonds that open in a sequential, zipper-like fashion and, thus, dominate the free-energy landscape of the transition are identified. Transitions between the closed and open conformations only have to overcome moderate free-energy barriers. Unexpectedly, the closed state and the open state encompass broad free-energy basins that contain conformations differing in domain hinge motions by up to 40°. The significance of these extended states is discussed in relation to recent experimental Förster resonance energy transfer measurements. Taken together, these results demonstrate how a small number of cooperative key interactions can shape the overall dynamics of an enzyme and suggest an “all-or-nothing” mechanism for the opening and closing of AdK. Our efficient DIMS molecular dynamics computer simulation approach can provide a detailed picture of a functionally important macromolecular transition and thus help to interpret and suggest experiments to probe the conformational landscape of dynamic proteins such as AdK.     

Plasticity of cytochrome P450 2B4 as investigated by hydrogen-deuterium exchange mass spectrometry and X-ray crystallography

Crystal structures of the xenobiotic metabolizing cytochrome P450 2B4 have demonstrated markedly different conformations in the presence of imidazole inhibitors or in the absence of ligand. However, knowledge of the plasticity of the enzyme in solution has remained scant. Thus, hydrogen-deuterium exchange mass spectrometry (DXMS) was utilized to probe the conformations of ligand-free P450 2B4 and the complex with 4-(4-chlorophenyl)imidazole (4-CPI) or 1-biphenyl-4-methyl-1H-imidazole (1-PBI). The results of DXMS indicate that the binding of 4-CPI slowed the hydrogen-deuterium exchange rate over the B'- and C-helices and portions of the F-G-helix cassette compared with P450 2B4 in the absence of ligands. In contrast, there was little difference between the ligand-free and 1-PBI-bound exchange sets. In addition, DXMS suggests that the ligand-free P450 2B4 is predominantly open in solution. Interestingly, a new high resolution structure of ligand-free P450 2B4 was obtained in a closed conformation very similar to the 4-CPI complex. Molecular dynamics simulations performed with the closed ligand-free structure as the starting point were used to probe the energetically accessible conformations of P450 2B4. The simulations were found to equilibrate to a conformation resembling the 1-PBI-bound P450 2B4 crystal structure. The results indicate that conformational changes observed in available crystal structures of the promiscuous xenobiotic metabolizing cytochrome P450 2B4 are consistent with its solution structural behavior.     

Investigation of the mechanism of domain closure in citrate synthase by molecular dynamics simulation

Six, 2 ns molecular dynamics simulations have been performed on the homodimeric enzyme citrate synthase. In three, both monomers were started from the open, unliganded X-ray conformation. In the remaining three, both monomers started from a closed, liganded X-ray conformation, with the ligands removed. Projecting the motion from the simulations onto the experimental domain motion revealed that the free-energy profile is rather flat around the open conformation, with steep sides. The most closed conformations correspond to hinge-bending angles of 12-14 compared to the 20 degrees that occurs upon the binding of oxaloacetate. It is also found that the open, unliganded X-ray conformation is situated at the edge of the steep rise in free energy, although conformations that are about 5 degrees more open were sampled. A rigid-body essential dynamics analysis of the combined open trajectories has shown that domain motions in the direction of the closed X-ray conformation are compatible with the natural domain motion of the unliganded protein, which has just two main degrees of freedom. The simulations starting from the closed conformation suggest a free-energy profile with a small barrier in going from the closed to open conformation. A combined essential dynamics and hinge-bending analysis of a trajectory that spontaneously converts from the closed to open state shows an almost exact correspondence to the experimental transition that occurs upon ligand binding. The simulations support the conclusion from an earlier analysis of the experimental transition that the beta-hairpin acts as a mechanical hinge by attaching the small domain to the large domain through a conserved main-chain hydrogen bond and salt-bridges, and allowing rotation to occur via its two flexible termini. The results point to a mechanism of domain closure in citrate synthase that has analogy to the process of closing a door.     

Reaction path and free energy calculations of the transition between alternate conformations of HIV-1 protease

Two different structures of ligand-free HIV protease have been determined by X-ray crystallography. These structures differ in the position of two 12 residue, β-hairpin regions (or “flaps”) which cap the active site. The movements of the flaps must be involved in the binding of substrates since, in either conformation, the flaps block the binding site. One of these structures is similar to structures of the ligand-bound enzyme; however, the importance of both structures to enzyme function is unclear. This transformation takes place on a time scale too long for conventional molecular dynamics simulations, so the process was studied by first identifying a reaction path between the two structures and then calculating the free energy along this path using umbrella sampling. For the ligand-free enzyme, it is found that the two structures are nearly equally stable, with the ligand-bound-type structure being less stable, consistent with X-ray crystallography data. The more stable open structure does not have a lower potential energy, but is stabilized by entropy. The transition occurs through a collapse and reformation of the β-sheet structure of the conformationally flexible, glycine-rich flap ends. Additionally, some problems in studying conformational changes in proteins through the use of a single reaction path are addressed. Proteins 32:7–16, 1998. © 1998 Wiley-Liss, Inc.     

Molecular Dynamics Study of the Opening Mechanism for DNA Polymerase I

During DNA replication, DNA polymerases follow an induced fit mechanism in order to rapidly distinguish between correct and incorrect dNTP substrates. The dynamics of this process are crucial to the overall effectiveness of catalysis. Although X-ray crystal structures of DNA polymerase I with substrate dNTPs have revealed key structural states along the catalytic pathway, solution fluorescence studies indicate that those key states are populated in the absence of substrate. Herein, we report the first atomistic simulations showing the conformational changes between the closed, open, and ajar conformations of DNA polymerase I in the binary (enzyme∶DNA) state to better understand its dynamics. We have applied long time-scale, unbiased molecular dynamics to investigate the opening process of the fingers domain in the absence of substrate for B. stearothermophilis DNA polymerase in silico. These simulations are biologically and/or physiologically relevant as they shed light on the transitions between states in this important enzyme. All closed and ajar simulations successfully transitioned into the fully open conformation, which is known to be the dominant binary enzyme-DNA conformation from solution and crystallographic studies. Furthermore, we have detailed the key stages in the opening process starting from the open and ajar crystal structures, including the observation of a previously unknown key intermediate structure. Four backbone dihedrals were identified as important during the opening process, and their movements provide insight into the recognition of dNTP substrate molecules by the polymerase binary state. In addition to revealing the opening mechanism, this study also demonstrates our ability to study biological events of DNA polymerase using current computational methods without biasing the dynamics.     

Investigating the accessibility of the closed domain conformation of citrate synthase using essential dynamics sampling

A molecular dynamics study of pig heart citrate synthase is presented that aims to directly address the question of whether, for this enzyme, the ligand-induced closed domain conformation is accessible to the open unliganded enzyme. The approach utilises the technique of essential dynamics sampling, which is used in two modes. In exploring mode, the enzyme is encouraged to explore domain conformations it might not normally sample in free molecular dynamics simulation. In targeting mode, the enzyme is encouraged to adopt the domain conformation of a target structure. Using both modes extensively, it has been found that when the enzyme is prepared from a crystallographic open-domain structure and is in the unliganded state, it is unable to adopt the crystallographic closed-domain conformation of the liganded enzyme. Likewise, when the enzyme is prepared from the crystallographic closed liganded conformation with the ligands removed, it is unable to adopt the crystallographic open domain conformation. Structural investigations point to a common structural difference that is the source of this energy barrier; namely, the shift of alpha-helix 328-341 along its own axis relative to the large domain. Without this shift, the domains are unable to close or open fully. The charged substrate, oxaloacetate, binds near the base of this helix in the large domain and the interaction of Arg329 at the base of the helix with oxaloacetate is one that is consistent with the shift of this helix in going from the crystallographic open to closed structure. Therefore, the results suggest that without the substrate the enzyme remains in a partially open conformation ready to receive the substrate. In this way, the efficiency of the enzyme should be increased over one that is closed part of the time, with its binding site inaccessible to the substrate.     

Structural basis for the altered activity of Gly794 variants of Escherichia coli beta-galactosidase

The open-closed conformational switch in the active site of Escherichia coli beta-galactosidase was studied by X-ray crystallography and enzyme kinetics. Replacement of Gly794 by alanine causes the apoenzyme to adopt the closed rather than the open conformation. Binding of the competitive inhibitor isopropyl thio-beta-D-galactoside (IPTG) requires the mutant enzyme to adopt its less favored open conformation, weakening affinity relative to wild type. In contrast, transition-state inhibitors bind to the enzyme in the closed conformation, which is favored for the mutant, and display increased affinity relative to wild type. Changes in affinity suggest that the free energy difference between the closed and open forms is 1-2 kcal/mol. By favoring the closed conformation, the substitution moves the resting state of the enzyme along the reaction coordinate relative to the native enzyme and destabilizes the ground state relative to the first transition state. The result is that the rate constant for galactosylation is increased but degalactosylation is slower. The covalent intermediate may be better stabilized than the second transition state. The substitution also results in better binding of glucose to both the free and the galactosylated enzyme. However, transgalactosylation with glucose to produce allolactose (the inducer of the lac operon) is slower with the mutant than with the native enzyme. This suggests either that the glucose is misaligned for the reaction or that the galactosylated enzyme with glucose bound is stabilized relative to the transition state for transgalactosylation.     

A highly strained nuclear conformation of the exportin Cse1p revealed by molecular dynamics simulations

To investigate the stability of the open nuclear state of the exportin Cse1p and its closing mechanism at the atomic level, we have performed multiple molecular dynamics simulations. The simulations revealed a strikingly fast transition of Cse1p from the open conformation to the closed cytoplasmic form, consistent with the proposal that Cse1p represents a "spring-loaded molecule." The structure of the ring-shaped state obtained in the simulations is remarkably close to the crystal structure of the cytoplasmic state, though the open nuclear structure was used as the only input. The conformational change is initially driven by release of strain due to RanGTP/importin-alpha binding. Subsequently, a stable closed state is formed, driven by attraction of electrostatically complementary interfaces. These results are consistent with and extend previous proposals. Reverse-charge and neutral mutants remained in an open state. The simulations predict a detailed reaction pathway and resolve the role of suggested hinge regions.     

Molecular dynamics study of the energetic, mechanistic, and structural implications of a closed phosphate tube in ncd

The switch 1 region of myosin forms a lid over the nucleotide phosphates as part of a structure known as the phosphate-tube. The homologous region in kinesin-family motors is more open, not interacting with the nucleotide. We used molecular dynamics (MD) simulations to examine a possible displacement of switch 1 of the microtubule motor, ncd, from the open conformation to the closed conformation seen in myosin. MD simulations were done of both the open and the closed conformations, with either MgADP or MgATP at the active site. All MD structures were stable at 300 K for 500 ps, implying that the open and closed conformers all represented local minima on a global free energy surface. Free energy calculations indicated that the open structure was energetically favored with MgADP at the active site, suggesting why only the open structure has been captured in crystallographic work. With MgATP, the closed and open structures had roughly equal energies. Simulated annealing MD showed the transformation from the closed phosphate-tube ncd structure to an open configuration. The MD simulations also showed that the coordination of switch 1 to the nucleotide dramatically affected the position of both the bound nucleotide and switch 2 and that a closed phosphate-tube may be necessary for catalysis.     

Exploring the role of large conformational changes in the fidelity of DNA polymerase beta

The relationships between the conformational landscape, nucleotide insertion catalysis and fidelity of DNA polymerase beta are explored by means of computational simulations. The simulations indicate that the transition states for incorporation of right (R) and wrong (W) nucleotides reside in substantially different protein conformations. The protein conformational changes that reproduce the experimentally observed fidelity are significantly larger than the small rearrangements that usually accompany motions from the reactant state to the transition state in common enzymatic reactions. Once substrate binding has occurred, different constraints imposed on the transition states for insertion of R and W nucleotides render it highly unlikely that both transition states can occur in the same closed structure, because the predicted fidelity would then be many orders of magnitude too large. Since the conformational changes reduce the transition state energy of W incorporation drastically they decrease fidelity rather than increase it. Overall, a better agreement with experimental data is attained when the R is incorporated through a transition state in a closed conformation and W is incorporated through a transition state in one or perhaps several partially open conformations. The generation of free energy surfaces for R and W also allow us to analyze proposals about the relationship between induced fit and fidelity.     

Crystal structures of a ligand-free and malonate-bound human caspase-1: implications for the mechanism of substrate binding

Caspase-1, a mediator of the posttranslational processing of IL-1beta and IL-18, requires an aspartic acid in the P1 position of its substrates. The mechanisms of caspase-1 activation remain poorly understood despite numerous structures of the enzyme complexed with aspartate-based inhibitors. Here we report a crystal structure of ligand-free caspase-1 that displays dramatic rearrangements of loops defining the active site to generate a closed conformation that is incompatible with substrate binding. A structure of the enzyme complexed with malonate shows the protein in its open (active-site ligand-bound) conformation in which malonate reproduces the hydrogen bonding network observed in structures with covalent inhibitors. These results illustrate the essential function of the obligatory aspartate recognition element that opens the active site of caspase-1 to substrates and may be the determinant responsible for the conformational changes between ligand-free and -bound forms of the enzyme, and suggest a new approach for identifying novel aspartic acid mimetics.     

Conformational flexibility of the leucine binding protein examined by protein domain coarse-grained molecular dynamics

Periplasmic binding proteins are the initial receptors for the transport of various substrates over the inner membrane of gram-negative bacteria. The binding proteins are composed of two domains, and the substrate is entrapped between these domains. For several of the binding proteins it has been established that a closed-up conformation exists even without substrate present, suggesting a highly flexible apo-structure which would compete with the ligand-bound protein for the transporter interaction. For the leucine binding protein (LBP), structures of both open and closed conformations are known, but no closed-up structure without substrate has been reported. Here we present molecular dynamics simulations exploring the conformational flexibility of LBP. Coarse grained models based on the MARTINI force field are used to access the microsecond timescale. We show that a standard MARTINI model cannot maintain the structural stability of the protein whereas the ELNEDIN extension to MARTINI enables simulations showing a stable protein structure and nanosecond dynamics comparable to atomistic simulations, but does not allow the simulation of conformational flexibility. A modification to the MARTINI-ELNEDIN setup, referred to as domELNEDIN, is therefore presented. The domELNEDIN setup allows the protein domains to move independently and thus allows for the simulation of conformational changes. Microsecond domELNEDIN simulations starting from either the open or the closed conformations consistently show that also for LBP, the apo-structure is flexible and can exist in a closed form.

The dynamics of the MgATP-driven closure of MalK, the energy-transducing subunit of the maltose ABC transporter

The nucleotide binding domains (NBDs) are the energy supplying subunits of ATP-binding cassette (ABC) proteins. They power transport by binding and hydrolyzing ATP. Tracing the pathway between different conformational states of the NBDs during ATP binding, hydrolysis, and release has, however, proven difficult. We have used molecular dynamics simulations to study the ATP-driven association of the NBDs of the maltose ABC transporter, MalK, based on the crystal structures of its open and semiopen dimers. When MgATP was introduced into the binding pockets, the semiopen dimer transitioned to a closed conformation, whereas the open dimer evolved to a semiopen state. In the absence of docked MgATP, however, the twin NBDs of both the open and semiopen starting configurations drifted further apart. Both the presence of MgATP and direct cross-interface protein-protein hydrogen bonds, primarily involving the D-loop, quite likely play a key role in initiating closure. The simulations of the MgATP-docked semiopen form indicate that completion of closure is driven mainly by cross-interface contacts between the gamma-phosphate of ATP and residues in the signature motif. Our simulations also give insight into possible interactions of MalK with the regulatory proteins MalT and enzyme IIA(glc).     

High count rates with total internal reflection fluorescence correlation spectroscopy

We have modeled the structure of KirBac1.1 in an open state using as a starting point the structure of KirBac1.1 in its closed conformation (Protein Data Bank 1P7B). To test the validity of the open-state model, molecular dynamics simulations in octane, a lipid bilayer mimetic, were carried out. Simulations of the closed conformer were used for comparison purposes. The total simulation time was approximately 138 ns. The initial open model was refined by using projection maps obtained from electron microscopy experiments on two-dimensional crystals of the inwardly rectifying K+ channel KirBac3.1 from Magentospirillum magnetotacticum captured in its open state (C. Vénien-Bryan, unpublished data). Significant movements of the outer helices take place in going from the closed to the open model in agreement with structural and biochemical data in potassium channels, which suggests that gating is accomplished by a conformational change that takes place in the transmembrane domain upon an external stimulus. The motion of the inner helices is mainly achieved by bending at conserved glycine residues that have been previously reported to act as molecular hinges. Overall, these simulations suggest that the open conformer is stable, providing a plausible all-atom model that will enable the study of potential gating mechanisms in more detail.     

Conformational plasticity and dynamics in the generic protein folding catalyst SlyD unraveled by single-molecule FRET

The relation between conformational dynamics and chemistry in enzyme catalysis recently has received increasing attention. While, in the past, the mechanochemical coupling was mainly attributed to molecular motors, nowadays, it seems that this linkage is far more general. Single-molecule fluorescence methods are perfectly suited to directly evidence conformational flexibility and dynamics. By labeling the enzyme SlyD, a member of peptidyl-prolyl cis-trans isomerases of the FK506 binding protein type with an inserted chaperone domain, with donor and acceptor fluorophores for single-molecule fluorescence resonance energy transfer, we directly monitor conformational flexibility and conformational dynamics between the chaperone domain and the FK506 binding protein domain. We find a broad distribution of distances between the labels with two main maxima, which we attribute to an open conformation and to a closed conformation of the enzyme. Correlation analysis demonstrates that the conformations exchange on a rate in the 100 Hz range. With the aid from Monte Carlo simulations, we show that there must be conformational flexibility beyond the two main conformational states. Interestingly, neither the conformational distribution nor the dynamics is significantly altered upon binding of substrates or other known binding partners. Based on these experimental findings, we propose a model where the conformational dynamics is used to search the conformation enabling the chemical step, which also explains the remarkable substrate promiscuity connected with a high efficiency of this class of peptidyl-prolyl cis-trans isomerases.     

How maltose influences structural changes to bind to maltose‐binding protein: Results from umbrella sampling simulation

A well‐studied periplasmic‐binding protein involved in the abstraction of maltose is maltose‐binding protein (MBP), which undergoes a ligand‐induced conformational transition from an open (ligand‐free) to a closed (ligand‐bound) state. Umbrella sampling simulations have been us to estimate the free energy of binding of maltose to MBP and to trace the potential of mean force of the unbinding event using the center‐of‐mass distance between the protein and ligand as the reaction coordinate. The free energy thus obtained compares nicely with the experimentally measured value justifying our theoretical basis. Measurement of the domain angle (N‐terminal‐domain – hinge – C‐terminal‐domain) along the unbinding pathway established the existence of three different states. Starting from a closed state, the protein shifts to an open conformation during the initial unbinding event of the ligand then resides in a semi‐open conformation and later resides predominantly in an open‐state. These transitions along the ligand unbinding pathway have been captured in greater depth using principal component analysis. It is proposed that in mixed‐model, both conformational selection and an induced‐fit mechanism combine to the ligand recognition process in MBP. Proteins 2013. © 2012 Wiley Periodicals, Inc.