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.     

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.     

Conformational Dynamics of a Ligand-Free Adenylate Kinase

Adenylate kinase (AdK) is a phosphoryl-transfer enzyme with important physiological functions. Based on a ligand-free open structure and a ligand-bound closed structure solved by crystallography, here we use molecular dynamics simulations to examine the stability and dynamics of AdK conformations in the absence of ligands. We first perform multiple simulations starting from the open or the closed structure, and observe their free evolutions during a simulation time of 100 or 200 nanoseconds. In all seven simulations starting from the open structure, AdK remained stable near the initial conformation. The eight simulations initiated from the closed structure, in contrast, exhibited large variation in the subsequent evolutions, with most (seven) undergoing large-scale spontaneous conformational changes and approaching or reaching the open state. To characterize the thermodynamics of the transition, we propose and apply a new sampling method that employs a series of restrained simulations to calculate a one-dimensional free energy along a curved pathway in the high-dimensional conformational space. Our calculated free energy profile features a single minimum at the open conformation, and indicates that the closed state, with a high (∼13 kcal/mol) free energy, is not metastable, consistent with the observed behaviors of the unrestrained simulations. Collectively, our simulations suggest that it is energetically unfavorable for the ligand-free AdK to access the closed conformation, and imply that ligand binding may precede the closure of the enzyme.     

Insights into the Structure of Invisible Conformations of Large Methyl Group Labeled Molecular Machines from High Pressure NMR

Most proteins are highly flexible and can adopt conformations that deviate from the energetically most favorable ground state. Structural information on these lowly populated, alternative conformations is often lacking, despite the functional importance of these states. Here, we study the pathway by which the Dcp1:Dcp2 mRNA decapping complex exchanges between an autoinhibited closed and an open conformation. We make use of methyl Carr–Purcell–Meiboom–Gill (CPMG) NMR relaxation dispersion (RD) experiments that report on the population of the sparsely populated open conformation as well as on the exchange rate between the two conformations. To obtain volumetric information on the open conformation as well as on the transition state structure we made use of RD measurements at elevated pressures. We found that the open Dcp1:Dcp2 conformation has a lower molecular volume than the closed conformation and that the transition state is close in volume to the closed state. In the presence of ATP the volume change upon opening of the complex increases and the volume of the transition state lies in-between the volumes of the closed and open state. These findings show that ATP has an effect on the volume changes that are associated with the opening-closing pathway of the complex. Our results highlight the strength of pressure dependent NMR methods to obtain insights into structural features of protein conformations that are not directly observable. As our work makes use of methyl groups as NMR probes we conclude that the applied methodology is also applicable to high molecular weight complexes.     

The Effect of Tensile Stress on the Conformational Free Energy Landscape of Disulfide Bonds

Disulfide bridges are no longer considered to merely stabilize protein structure, but are increasingly recognized to play a functional role in many regulatory biomolecular processes. Recent studies have uncovered that the redox activity of native disulfides depends on their C–C–S–S dihedrals, and . Moreover, the interplay of chemical reactivity and mechanical stress of disulfide switches has been recently elucidated using force–clamp spectroscopy and computer simulation. The and angles have been found to change from conformations that are open to nucleophilic attack to sterically hindered, so–called closed states upon exerting tensile stress. In view of the growing evidence of the importance of C–C–S–S dihedrals in tuning the reactivity of disulfides, here we present a systematic study of the conformational diversity of disulfides as a function of tensile stress. With the help of force-clamp metadynamics simulations, we show that tensile stress brings about a large stabilization of the closed conformers, thereby giving rise to drastic changes in the conformational free energy landscape of disulfides. Statistical analysis shows that native TDi, DO and interchain Ig protein disulfides prefer open conformations, whereas the intrachain disulfide bridges in Ig proteins favor closed conformations. Correlating mechanical stress with the distance between the two –carbons of the disulfide moiety reveals that the strain of intrachain Ig protein disulfides corresponds to a mechanical activation of about 100 pN. Such mechanical activation leads to a severalfold increase of the rate of the elementary redox reaction step. All these findings constitute a step forward towards achieving a full understanding of functional disulfides.     

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.     

Protein folding and binding in confined spaces and in crowded solutions

Simple theoretical models are presented to illustrate the effects of spatial confinement and macromolecular crowding on the equilibria and rates of protein folding and binding. Confinement is expected to significantly stabilize the folded state, but for crowding only a marginal effect on protein stability is expected. In confinement the unfolded chain is restricted to a cage but in crowding the unfolded chain may explore different interstitial voids. Because confinement and crowding eliminate the more expanded conformations of the unfolded state, folding from the compact unfolded state is expected to speed up. Crowding will shift the binding equilibrium of proteins toward the bound state. The significant slowing down in protein diffusion by crowding, perhaps beneficial for chaperonin action, could result in a decrease in protein binding rates.     

Direct Observation of T4 Lysozyme Hinge-Bending Motion by Fluorescence Correlation Spectroscopy

Bacteriophage T4 Lysozyme (T4L) catalyzes the hydrolysis of the peptidoglycan layer of the bacterial cell wall late in the infection cycle. It has long been postulated that equilibrium dynamics enable substrate access to the active site located at the interface between the N- and C-terminal domains. Crystal structures of WT-T4L and point mutants captured a range of conformations that differ by the hinge-bending angle between the two domains. Evidence of equilibrium between open and closed conformations in solution was gleaned from distance measurements between the two domains but the nature of the equilibrium and the timescale of the underlying motion have not been investigated. Here, we used fluorescence fluctuation spectroscopy to directly detect T4L equilibrium conformational fluctuations in solution. For this purpose, Tetramethylrhodamine probes were introduced at pairs of cysteines in regions of the molecule that undergo relative displacement upon transition from open to closed conformations. Correlation analysis of Tetramethylrhodamine intensity fluctuations reveals hinge-bending motion that changes the relative distance and orientation of the N- and C-terminal domains with ≅15 μs relaxation time. That this motion involves interconversion between open and closed conformations was further confirmed by the dampening of its amplitude upon covalent substrate trapping. In contrast to the prevalent two-state model of T4L equilibrium, molecular brightness and number of particles obtained from cumulant analysis suggest that T4L populates multiple intermediate states, consistent with the wide range of hinge-bending angles trapped in the crystal structure of T4L mutants.     

Gating Transition of Pentameric Ligand-Gated Ion Channels

Pentameric ligand-gated ion channels are an important family of membrane proteins and play key roles in physiological processes, including signal transduction at chemical synapses. Here, we study the conformational changes associated with the opening and closing of the channel pore. Based on recent crystal structures of two prokaryotic members of the family in open and closed states, respectively, mixed elastic network models are constructed for the transmembrane domain. To explore the conformational changes in the gating transition, a coarse-grained transition path is computed that smoothly connects the closed and open conformations of the channel. We find that the conformational transition involves no major rotations of the transmembrane helices, and is instead characterized by a concerted tilting of helices M2 and M3. In addition, helix M2 changes its bending state, which results in an early closure of the pore during the open-to-closed transition.     

Modulation of functionally significant conformational equilibria in adenylate kinase by high concentrations of trimethylamine oxide attributed to volume exclusion

The effect of an inert small molecule osmolyte, trimethyl amine N-oxide (TMAO), upon the conformational equilibria of Escherichia coli adenylate kinase was studied using time-resolved FRET. The relative populations of open and closed clefts between the LID and the CORE domains were measured as functions of the concentrations of the substrate ATP over the concentration range 0–18 mM and TMAO over the concentration range 0–4 M. A model was constructed according to which the enzyme exists in equilibrium among four conformational states, corresponding to combinations of open and closed conformations of the LID-CORE and AMP-CORE clefts. ATP is assumed to bind only to those conformations with the closed LID-CORE cleft, and TMAO is assumed to be differentially excluded as a hard spherical particle from each of the four conformations in accordance with calculations based upon x-ray crystallographic structures. This model was found to describe quantitatively the dependence of the fraction of the closed LID-CORE cleft upon the concentrations of both ATP and TMAO over the entire range of concentrations with just five undetermined parameters.     

Use of helper enzymes for ADP removal in infrared spectroscopic experiments: application to Ca2+-ATPase

Adenylate kinase (AdK) and apyrase were employed as helper enzymes to remove ADP in infrared spectroscopic experiments that study the sarcoplasmic reticulum Ca(2+)-ATPase. The infrared absorbance changes of their enzymatic reactions were characterized and used to monitor enzyme activity. AdK transforms ADP to ATP and AMP, whereas apyrase consumes ATP and ADP to generate AMP and inorganic phosphate. The benefits of using them as helper enzymes are severalfold: i), both remove ADP generated after ATP hydrolysis by ATPase, which enables repeat of ATP-release experiments several times with the same sample without interference by ADP; ii), AdK helps maintain the presence of ATP for a longer time by regenerating 50% of the initial ATP; iii), apyrase generates free P(i), which can help stabilize the ADP-insensitive phosphoenzyme (E2P); and iv), apyrase can be used to monitor ADP dissociation from transient enzyme intermediates with relatively high affinity to ADP, as shown here for ADP dissociation from the ADP-sensitive phosphoenzyme intermediate (Ca(2)E1P). The respective infrared spectra indicate that ADP dissociation relaxes the closed conformation immediately after phosphorylation partially back toward the open conformation of Ca(2)E1 but does not trigger the transition to E2P. The helper enzyme approach can be extended to study other nucleotide-dependent proteins.     

Crowding Induces Complex Ergodic Diffusion and Dynamic Elongation of Large DNA Molecules

Despite the ubiquity of molecular crowding in living cells, the effects of crowding on the dynamics of genome-sized DNA are poorly understood. Here, we track single, fluorescent-labeled large DNA molecules (11, 115 kbp) diffusing in dextran solutions that mimic intracellular crowding conditions (0–40%), and determine the effects of crowding on both DNA mobility and conformation. Both DNAs exhibit ergodic Brownian motion and comparable mobility reduction in all conditions; however, crowder size (10 vs. 500 kDa) plays a critical role in the underlying diffusive mechanisms and dependence on crowder concentration. Surprisingly, in 10-kDa dextran, crowder influence saturates at ∼20% with an ∼5× drop in DNA diffusion, in stark contrast to exponentially retarded mobility, coupled to weak anomalous subdiffusion, with increasing concentration of 500-kDa dextran. Both DNAs elongate into lower-entropy states (compared to random coil conformations) when crowded, with elongation states that are gamma distributed and fluctuate in time. However, the broadness of the distribution of states and the time-dependence and length scale of elongation length fluctuations depend on both DNA and crowder size with concentration having surprisingly little impact. Results collectively show that mobility reduction and coil elongation of large crowded DNAs are due to a complex interplay between entropic effects and crowder mobility. Although elongation and initial mobility retardation are driven by depletion interactions, subdiffusive dynamics, and the drastic exponential slowing of DNA, up to ∼300×, arise from the reduced mobility of larger crowders. Our results elucidate the highly important and widely debated effects of cellular crowding on genome-sized DNA.     

Toxin binding reveals two open state structures for one acid-sensing ion channel

Of the three principal conformations of acid-sensing ion channels (ASICs)—closed, open and desensitized—only the atomic structure of the desensitized conformation had been known. Two recent papers report the crystal structure of chicken ASIC1 in complex with the spider toxin psalmotoxin 1, and one of these studies finds that, depending on the pH, channels are in two different open conformations. Compared with the desensitized conformation, toxin binding induces only subtle structural changes in the lower part of the large extracellular domain but a complete rearrangement of the two transmembrane domains (TMDs), suggesting that desensitization gating (the transition from open to desensitized) is mainly associated with conformational rearrangements of the TMDs. Moreover, the study reveals how two different arrangements of the TMDs in the open state give rise to ion pores with different selectivity for monovalent cations.     

Conformational equilibria and free energy profiles for the allosteric transition of the ribose-binding protein

The ribose-binding protein (RBP) is a sugar-binding bacterial periplasmic protein whose function is associated with a large allosteric conformational change from an open to a closed conformation upon binding to ribose. The crystal structures of RBP in open and closed conformations have been solved. It has been hypothesized that the open and closed conformations exist in a dynamic equilibrium in solution, and that sugar binding shifts the population from open conformations to closed conformations. Here, we study by computer simulations the thermodynamic changes that accompany this conformational change, and model the structural changes that accompany the allosteric transition, using umbrella sampling molecular dynamics and the weighted histogram analysis method. The open state is comprised of a diverse ensemble of conformations; the open ribose-free X-ray crystal conformations being representative of this ensemble. The unligated open form of RBP is stabilized by conformational entropy. The simulations predict detectable populations of closed ribose-free conformations in solution. Additional interdomain hydrogen bonds stabilize this state. The predicted shift in equilibrium from the open to the closed state on binding to ribose is in agreement with experiments. This is driven by the energetic stabilization of the closed conformation due to ribose-protein interactions. We also observe a significant population of a hitherto unobserved ribose-bound partially open state. We believe that this state is the one that has been suggested to play a role in the transfer of ribose to the membrane-bound permease complex.     

Effects of Molecular Crowding on the Dynamics of Intrinsically Disordered Proteins

Inside cells, the concentration of macromolecules can reach up to 400 g/L. In such crowded environments, proteins are expected to behave differently than in vitro. It has been shown that the stability and the folding rate of a globular protein can be altered by the excluded volume effect produced by a high density of macromolecules. However, macromolecular crowding effects on intrinsically disordered proteins (IDPs) are less explored. These proteins can be extremely dynamic and potentially sample a wide ensemble of conformations under non-denaturing conditions. The dynamic properties of IDPs are intimately related to the timescale of conformational exchange within the ensemble, which govern target recognition and how these proteins function. In this work, we investigated the macromolecular crowding effects on the dynamics of several IDPs by measuring the NMR spin relaxation parameters of three disordered proteins (ProTα, TC1, and α-synuclein) with different extents of residual structures. To aid the interpretation of experimental results, we also performed an MD simulation of ProTα. Based on the MD analysis, a simple model to correlate the observed changes in relaxation rates to the alteration in protein motions under crowding conditions was proposed. Our results show that 1) IDPs remain at least partially disordered despite the presence of high concentration of other macromolecules, 2) the crowded environment has differential effects on the conformational propensity of distinct regions of an IDP, which may lead to selective stabilization of certain target-binding motifs, and 3) the segmental motions of IDPs on the nanosecond timescale are retained under crowded conditions. These findings strongly suggest that IDPs function as dynamic structural ensembles in cellular environments.     

Phosphorylated Phospholamban Stabilizes a Compact Conformation of?the Cardiac Calcium-ATPase

The sarcoendoplasmic reticulum calcium ATPase (SERCA) plays a key role in cardiac calcium handling and is considered a high-value target for the treatment of heart failure. SERCA undergoes conformational changes as it harnesses the chemical energy of ATP for active transport. X-ray crystallography has provided insight into SERCA structural substates, but it is not known how well these static snapshots describe in vivo conformational dynamics. The goals of this work were to quantify the direction and magnitude of SERCA motions as the pump performs work in live cardiac myocytes, and to identify structural determinants of SERCA regulation by phospholamban. We measured intramolecular fluorescence resonance energy transfer (FRET) between fluorescent proteins fused to SERCA cytoplasmic domains. We detected four discrete structural substates for SERCA expressed in cardiac muscle cells. The relative populations of these discrete states oscillated with electrical pacing. Low FRET states were most populated in low Ca (diastole), and were indicative of an open, disordered structure for SERCA in the E2 (Ca-free) enzymatic substate. High FRET states increased with Ca (systole), suggesting rigidly closed conformations for the E1 (Ca-bound) enzymatic substates. Notably, a special compact E1 state was observed after treatment with β-adrenergic agonist or with coexpression of phosphomimetic mutants of phospholamban. The data suggest that SERCA calcium binding induces the pump to undergo a transition from an open, dynamic conformation to a closed, ordered structure. Phosphorylated phospholamban stabilizes a unique conformation of SERCA that is characterized by a compact architecture.     

Two distinct gating mechanisms in gap junction channels: CO2-sensitive and voltage-sensitive.

The chemical gating of single-gap junction channels was studied by the dual whole-cell voltage-clamp method in HeLa cells transfected with connexin43 (HeLa43) and in fibroblasts from sciatic nerves. Junctional current (Ij), single-channel conductance, and Ij kinetics were studied in cell pairs during CO2 uncoupling and recoupling at small transjunctional voltages (Vj < 35 mV: Vj gating absent) and at high Vj (Vj > 40 mV: Vj gating strongly activated). In the absence of Vj gating, CO2 exclusively caused Ij slow transitions from open to closed channel states (mean transition time: approximately 10 ms), corresponding to a single-channel conductance of approximately 120 pS. At Vj > 40 mV, Vj gating induced fast Ij flickering between open, gamma j(main state), and residual, gamma j(residual), states (transition time: approximately 2 ms). The ratio gamma j(main state)/gamma j(residual) was approximately 4-5. No obvious correlation between Ij fast flickering and CO2 treatment was noticed. At high Vj, in addition to slow Ij transitions between open and closed states, CO2 induced slow transitions between residual and closed states. During recoupling, each channel reopened by a slow transition (mean transition time: approximately 10 ms) from closed to open state (rarely from closed to residual state). Fast Ij flickering between open and residual states followed. The data are in agreement with the hypothesis that gap junction channels possess two gating mechanisms, and indicate that CO2 induces channel gating exclusively by the slow gating mechanism.     

The blind search for the closed states of hinge-bending proteins

The hinge-bending proteins provide the most pronounced example of the large-amplitude slow motions in a number of proteins, which are critical for their functioning. They are often used as the test ground for developing novel approaches to the simulation of slow protein dynamics. In the present study, we present the algorithm, which allows physically-consistent simulations of slow protein dynamics in globular proteins. Our algorithm is based on the hierarchical clustering of the correlation patterns (HCCP) technique of domain identification, which allows subdividing the protein into the hierarchy of the rigid-body-like clusters. The clusters are allowed to rotate relative to one another on the automatically identified hinges. The clusters are found in the course of automated, objective and well-tested procedure. In the present communication, our technique is applied to 10 hinge-bending proteins. For each of the proteins, we performed the blind search for the closed conformation, staring from the open one. Resulting closed conformations are compared with the closed states observed in crystallographic structures. It is shown that our technique produces realistic closed conformations for 8 out of 10 studied proteins. This demonstrates that HCCP technique can be used for finding alternative protein conformations and for sampling the slow motions in proteins.     

End-point targeted molecular dynamics: large-scale conformational changes in potassium channels

Large-scale conformational changes in proteins that happen often on biological time scales may be relatively rare events on the molecular dynamics time scale. We have implemented an approach to targeted molecular dynamics called end-point targeted molecular dynamics that transforms proteins between two specified conformational states through the use of nonharmonic “soft” restraints. A key feature of the method is that the protein is free to discover its own conformational pathway through the plethora of possible intermediate states. The method is applied to the Shaker K_v1.2 potassium channel in implicit solvent. The rate of cycling between the open and closed states was varied to explore how slow the cycling rate needed to be to ensure that microscopic reversibility along the transition pathways was well approximated. Results specific to the K⁺ channel include: 1), a variation in backbone torsion angles of residues near the Pro-Val-Pro motif in the inner helix during both opening and closing; 2), the identification of possible occlusion sites in the closed channel located among Pro-Val-Pro residues and downstream; 3), a difference in the opening and closing pathways of the channel; and 4), evidence of a transient intermediate structural substate. The results also show that likely intermediate conformations during the opening-closing process can be generated in computationally tractable simulation times.