Effects of environment on the structure of Pyrococcus furiosus rubredoxin: a molecular dynamics study

Molecular dynamics simulations based on a 0.95-A resolution crystal structure of Pyrococcus furiosus have been performed to elucidate the effects of the environment on the structure of rubredoxin, and proteins in general. Three 1-ns simulations are reported here: two crystalline state simulations at 123 and 300 K, and a solution state simulation at 300 K. These simulations show that temperature has a greater impact on the protein structure than the close molecular contacts of the crystal matrix in rubredoxin, although both have an effect on its dynamic properties. These results indicate that differences between NMR solution structures and X-ray crystal structures will be relatively minor if they are done at similar temperatures. In addition, the crystal simulations appears to mimic previous crystallographic experiments on the effects of cryo-temperature on temperature factors, and might provide a useful tool in the structural analysis of protein structures solved at cryo-temperatures.     

Solution Structures of Rat Amylin Peptide: Simulation, Theory, and Experiment

Amyloid deposits of amylin in the pancreas are an important characteristic feature found in patients with Type-2 diabetes. The aggregate has been considered important in the disease pathology and has been studied extensively. However, the secondary structures of the individual peptide have not been clearly identified. In this work, we present detailed solution structures of rat amylin using a combination of Monte Carlo and molecular dynamics simulations. A new Monte Carlo method is presented to determine the free energy of distinct biomolecular conformations. Both folded and random-coil conformations of rat amylin are observed in water and their relative stability is examined in detail. The former contains an α-helical segment comprised of residues 7-17. We find that at room temperature the folded structure is more stable, whereas at higher temperatures the random-coil structure predominates. From the configurations and weights we calculate the α-carbon NMR chemical shifts, with results that are in reasonable agreement with experiments of others. We also calculate the infrared spectrum in the amide I stretch regime, and the results are in fair agreement with the experimental line shape presented herein.     

Kinetic Partitioning Modulates Human Telomere DNA G-Quadruplex Structural Polymorphism

Telomeres are specialized chromatin structures found at the end of chromosomes and are crucial to the maintenance of eukaryotic genome stability. Human telomere DNA is comprised of the repeating sequence (T₂AG₃)_n, which is predominantly double-stranded but terminates with a 3’ single-stranded tail. The guanine-rich tail can fold into secondary structures known as a G-quadruplexes (GQs) that may exist as a polymorphic mixture of anti-parallel, parallel, and several hybrid topological isomers. Using single-molecule Förster resonance energy transfer (smFRET), we have reconstructed distributions of telomere DNA GQ conformations generated by an in situ refolding protocol commonly employed in single-molecule studies of GQ structure, or using a slow cooling DNA annealing protocol typically used in the preparation of GQ samples for ensemble biophysical analyses. We find the choice of GQ folding protocol has a marked impact on the observed distributions of DNA conformations under otherwise identical buffer conditions. A detailed analysis of the kinetics of GQ folding over timescales ranging from minutes to hours revealed the distribution of GQ structures generated by in situ refolding gradually equilibrates to resemble the distribution generated by the slow cooling DNA annealing protocol. Interestingly, conditions of low ionic strength, which promote transient GQ unfolding, permit the fraction of folded DNA molecules to partition into a distribution that more closely approximates the thermodynamic folding equilibrium. Our results are consistent with a model in which kinetic partitioning occurs during in situ folding at room temperature in the presence of K⁺ ions, producing a long-lived non-equilibrium distribution of GQ structures in which the parallel conformation predominates on the timescale of minutes. These results suggest that telomere DNA GQ folding kinetics, and not just thermodynamic stability, likely contributes to the physiological ensemble GQ structures.     

Validation of X-ray Crystal Structure Ensemble Representations of SARS-CoV-2 Main Protease by Solution NMR Residual Dipolar Couplings

Considerable debate has focused on whether sampling of molecular dynamics trajectories restrained by crystallographic data can be used to develop realistic ensemble models for proteins in their natural, solution state. For the SARS-CoV-2 main protease, M^pro, we evaluated agreement between solution residual dipolar couplings (RDCs) and various recently reported multi-conformer and dynamic-ensemble crystallographic models. Although Phenix-derived ensemble models showed only small improvements in crystallographic R_free, substantially improved RDC agreement over fits to a conventionally refined 1.2-Å X-ray structure was observed, in particular for residues with above average disorder in the ensemble. For a set of six lower resolution (1.55–2.19 Å) M^pro X-ray ensembles, obtained at temperatures ranging from 100 to 310 K, no significant improvement over conventional two-conformer representations was found. At the residue level, large differences in motions were observed among these ensembles, suggesting high uncertainties in the X-ray derived dynamics. Indeed, combining the six ensembles from the temperature series with the two 1.2-Å X-ray ensembles into a single 381-member “super ensemble” averaged these uncertainties and substantially improved agreement with RDCs. However, all ensembles showed excursions that were too large for the most dynamic fraction of residues. Our results suggest that further improvements to X-ray ensemble refinement are feasible, and that RDCs provide a sensitive benchmark in such endeavors. Remarkably, a weighted ensemble of 350 PDB M^pro X-ray structures provided slightly better cross-validated agreement with RDCs than any individual ensemble refinement, implying that differences in lattice confinement also limit the fit of RDCs to X-ray coordinates.     

Enhancing the quality of protein conformation ensembles with relative populations

The function and dynamics of many proteins are best understood not from a single structure but from an ensemble. A high quality ensemble is necessary for accurately delineating protein dynamics. However, conformations in an ensemble are generally given equal weights. Few attempts were made to assign relative populations to the conformations, mainly due to the lack of right experimental data. Here we propose a method for assigning relative populations to ensembles using experimental residue dipolar couplings (RDC) as constraints, and show that relative populations can significantly enhance an ensemble’s ability in representing the native states and dynamics. The method works by identifying conformation states within an ensemble and assigning appropriate relative populations to them. Each of these conformation states is represented by a sub-ensemble consisting of a subset of the conformations. Application to the ubiquitin X-ray ensemble clearly identifies two key conformation states, with relative populations in excellent agreement with previous work. We then apply the method to a reprotonated ERNST ensemble that is enhanced with a switched conformation, and show that as a result of population reweighting, not only the reproduction of RDCs is significantly improved, but common conformational features (particularly the dihedral angle distributions of ? ₅₃ and ψ ₅₂) also emerge for both the X-ray ensemble and the reprotonated ERNST ensemble.     

Highly Disordered Amyloid-β Monomer Probed by Single-Molecule FRET and MD Simulation

Monomers of amyloid-β (Aβ) protein are known to be disordered, but there is considerable controversy over the existence of residual or transient conformations that can potentially promote oligomerization and fibril formation. We employed single-molecule Förster resonance energy transfer (FRET) spectroscopy with site-specific dye labeling using an unnatural amino acid and molecular dynamics simulations to investigate conformations and dynamics of Aβ isoforms with 40 (Aβ40) and 42 residues (Aβ42). The FRET efficiency distributions of both proteins measured in phosphate-buffered saline at room temperature show a single peak with very similar FRET efficiencies, indicating there is apparently only one state. 2D FRET efficiency-donor lifetime analysis reveals, however, that there is a broad distribution of rapidly interconverting conformations. Using nanosecond fluorescence correlation spectroscopy, we measured the timescale of the fluctuations between these conformations to be ～35 ns, similar to that of disordered proteins. These results suggest that both Aβ40 and Aβ42 populate an ensemble of rapidly reconfiguring unfolded states, with no long-lived conformational state distinguishable from that of the disordered ensemble. To gain molecular-level insights into these observations, we performed molecular dynamics simulations with a force field optimized to describe disordered proteins. We find, as in experiments, that both peptides populate configurations consistent with random polymer chains, with the vast majority of conformations lacking significant secondary structure, giving rise to very similar ensemble-averaged FRET efficiencies.     

A Correspondence Between Solution-State Dynamics of an Individual Protein and the Sequence and Conformational Diversity of its Family

Conformational ensembles are increasingly recognized as a useful representation to describe fundamental relationships between protein structure, dynamics and function. Here we present an ensemble of ubiquitin in solution that is created by sampling conformational space without experimental information using “Backrub” motions inspired by alternative conformations observed in sub-Angstrom resolution crystal structures. Backrub-generated structures are then selected to produce an ensemble that optimizes agreement with nuclear magnetic resonance (NMR) Residual Dipolar Couplings (RDCs). Using this ensemble, we probe two proposed relationships between properties of protein ensembles: (i) a link between native-state dynamics and the conformational heterogeneity observed in crystal structures, and (ii) a relation between dynamics of an individual protein and the conformational variability explored by its natural family. We show that the Backrub motional mechanism can simultaneously explore protein native-state dynamics measured by RDCs, encompass the conformational variability present in ubiquitin complex structures and facilitate sampling of conformational and sequence variability matching those occurring in the ubiquitin protein family. Our results thus support an overall relation between protein dynamics and conformational changes enabling sequence changes in evolution. More practically, the presented method can be applied to improve protein design predictions by accounting for intrinsic native-state dynamics.     

Conformational ensemble of an intrinsically flexible loop in mitochondrial import protein Tim21 studied by modeling and molecular dynamics simulations

BackgroundTim21, a subunit of a highly dynamic translocase of the inner mitochondrial membrane (TIM23) complex, translocates proteins by interacting with subunits in the translocase of the outer membrane (TOM) complex and Tim23 channel in the TIM23 complex. A loop segment in Tim21, which is in close proximity of the binding site of Tim23, has different conformations in X-ray, NMR and new crystal contact-free space (CCFS) structures. MD simulations can provide information on the structure and dynamics of the loop in solution.MethodsThe conformational ensemble of the loop was characterized using loop modeling and molecular dynamics (MD) simulations.ResultsMD simulations confirmed mobility of the loop. Multidimensional scaling and clustering were used to characterize the dynamic conformational ensemble of the loop. Free energy landscape showed that the CCFS crystal structure occupied a low energy region as compared to the conventional X-ray crystal structure. Analysis of crystal packing indicates that the CCFS provides larger conformational space for the motions of the loop.ConclusionsOur work reported the conformational ensemble of the loop in solution, which is in agreement with the structure obtained from CCFS approach. The combination of the experimental techniques and computational methods is beneficial for studying highly flexible regions of proteins.General significanceComputational methods, such as loop modeling and MD simulations, have proved to be useful for studying conformational flexibility of proteins. These methods in integration with experimental techniques such as CCFS has the potential to transform the studies on flexible regions of proteins.     

Influence of the g− conformation of Ser and Thr on the structure of transmembrane helices

In order to study the influence of Ser and Thr on the structure of transmembrane helices we have analyzed a database of helix stretches extracted from crystal structures of membrane proteins and an ensemble of model helices generated by molecular dynamics simulations. Both complementary analyses show that Ser and Thr in the g? conformation induce and/or stabilize a structural distortion in the helix backbone. Using quantum mechanical calculations, we have attributed this effect to the electrostatic repulsion between the side chain Oγ atom of Ser and Thr and the backbone carbonyl oxygen at position i ? 3. In order to minimize the repulsive force between these negatively charged oxygens, there is a modest increase of the helix bend angle as well as a local opening of the helix turn preceding Ser/Thr. This small distortion can be amplified through the helix, resulting in a significant displacement of the residues located at the other side of the helix. The crystal structures of aquaporin Z and the β₂-adrenergic receptor are used to illustrate these effects. Ser/Thr-induced structural distortions can be implicated in processes as diverse as ligand recognition, protein function and protein folding.     

Confirming the revised C-terminal domain of the MscL crystal structure

The structure of the C-terminal domain of the mechanosensitive channel of large conductance (MscL) has generated significant controversy. As a result, several structures have been proposed for this region: the original crystal structure (1MSL) of the Mycobacterium tuberculosis homolog (Tb), a model of the Escherichia coli homolog, and, most recently, a revised crystal structure of Tb-MscL (2OAR). To understand which of these structures represents a physiological conformation, we measured the impact of mutations to the C-terminal domain on the thermal stability of Tb-MscL using circular dichroism and performed molecular dynamics simulations of the original and the revised crystal structures of Tb-MscL. Our results imply that this region is helical and adopts an α-helical bundle conformation similar to that observed in the E. coli MscL model and the revised Tb-MscL crystal structure.     

Backbone Dynamics and Global Effects of an Activating Mutation in Minimized Mtu RecA Inteins

Inteins mediate protein splicing, which has found many applications in biotechnology and protein engineering. A single valine-to-leucine mutation (V67L) can globally enhance splicing and related cleavage reactions in minimized Mycobacterium tuberculosis RecA inteins. However, V67L mutation causes little change in crystal structures. To test whether protein dynamics contribute to activity enhancement in the V67L mutation, we have studied the conformations and dynamics of the minimized and engineered intein ΔΔIhh-V67CM and a single V67L mutant, ΔΔIhh-L67CM, by solution NMR. Chemical shift perturbations established that the V67L mutation causes global changes, including changes at the N-terminus and C-terminus of the intein, which are active sites for protein splicing. The single V67L mutation significantly slows hydrogen-exchange rates globally, indicating a shift to more stable conformations and reduction in ensemble distribution. Whereas the V67L mutation causes little change for motions on the picosecond-to-nanosecond timescale, motions on the microsecond-to-millisecond timescale affect a region involving the conserved F-block histidine and C-terminal asparagine, which are residues important for C-terminal cleavage. The V67L mutation is proposed to activate splicing by reducing the ensemble distribution of the intein structure and by modifying the active sites.     

Biochemical Studies Provide Insights into the Necessity for Multiple Arabidopsis thaliana Protein-Only RNase P Isoenzymes

RNase P catalyzes removal of the 5′ leader from precursor tRNAs (pre-tRNAs) in all three domains of life. Some eukaryotic cells contain multiple forms of the protein-only RNase P (PRORP) variant, prompting efforts to unravel this seeming redundancy. Previous studies concluded that there were only modest differences in the processing of typical pre-tRNAs by the three isoforms in Arabidopsis thaliana [AtPRORP1 (organellar), AtPRORP2 and AtPRORP3 (nuclear)]. Here, we investigated if different physical attributes of the three isoforms might engender payoffs under specific conditions. Our temperature–activity profiling studies revealed that AtPRORPs display substrate-identity dependent behavior at elevated temperatures (37–45 °C), with the organellar variant outperforming the nuclear counterparts. Echoing these findings, molecular dynamics simulations revealed that AtPRORP2 relative to AtPRORP1 samples a wider conformational ensemble that deviates from the crystal structure. Results from our biochemical studies and molecular dynamics simulations support the idea that AtPRORPs have overlapping but not necessarily redundant attributes and inspire new perspectives on the suitability of each variant to perform its function(s) in a specific cellular locale.     

Active Site Flexibility as a Hallmark for Efficient PET Degradation by I. sakaiensis PETase

Polyethylene terephthalate (PET) is one of the most-consumed synthetic polymers, with an annual production of 50 million tons. Unfortunately, PET accumulates as waste and is highly resistant to biodegradation. Recently, fungal and bacterial thermophilic hydrolases were found to catalyze PET hydrolysis with optimal activities at high temperatures. Strikingly, an enzyme from Ideonella sakaiensis, termed PETase, was described to efficiently degrade PET at room temperature, but the molecular basis of its activity is not currently understood. Here, a crystal structure of PETase was determined at 2.02 Å resolution and employed in molecular dynamics simulations showing that the active site of PETase has higher flexibility at room temperature than its thermophilic counterparts. This flexibility is controlled by a novel disulfide bond in its active site, with its removal leading to destabilization of the catalytic triad and reduction of the hydrolase activity. Molecular docking of a model substrate predicts that PET binds to PETase in a unique and energetically favorable conformation facilitated by several residue substitutions within its active site when compared to other enzymes. These computational predictions are in excellent agreement with recent mutagenesis and PET film degradation analyses. Finally, we rationalize the increased catalytic activity of PETase at room temperature through molecular dynamics simulations of enzyme-ligand complexes for PETase and other thermophilic PET-degrading enzymes at 298, 323, and 353 K. Our results reveal that both the binding pose and residue substitutions within PETase favor proximity between the catalytic residues and the labile carbonyl of the substrate at room temperature, suggesting a more favorable hydrolytic reaction. These results are valuable for enabling detailed evolutionary analysis of PET-degrading enzymes and for rational design endeavors aiming at increasing the efficiency of PETase and similar enzymes toward plastic degradation.     

Three-dimensional structures of enzyme-substrate complexes of the hydroxynitrile lyase from Hevea brasiliensis.

The 3D structures of complexes between the hydroxynitrile lyase from Hevea brasiliensis (Hb-HNL) and several substrate and/or inhibitor molecules, including trichloracetaldehyde, hexafluoracetone, acetone, and rhodanide, were determined by X-ray crystallography. The complex with trichloracetaldehyde showed a covalent linkage between the protein and the inhibitor, which had apparently resulted from nucleophilic attack of the catalytic Ser80-Ogamma. All other complexes showed the substrate or inhibitor molecule merely hydrogen bonded to the protein. In addition, the native crystal structure of Hb-HNL was redetermined at cryo-temperature and at room temperature, eliminating previous uncertainties concerning residual electron density within the active site, and leading to the observation of two conserved water molecules. One of them was found to be conserved in all complex structures and appears to have mainly structural significance. The other water molecule is conserved in all structures except for the complex with rhodanide; it is hydrogen bonded to the imidazole of the catalytic His235 and appears to affect the Hb-HNL catalyzed reaction. The observed 3D structural data suggest implications for the enzyme mechanism. It appears that the enzyme-catalyzed cyanohydrin formation is unlikely to proceed via a hemiacetal or hemiketal intermediate covalently attached to the enzyme, despite the observation of such an intermediate for the complex with trichloracetaldehyde. Instead, the data are consistent with a mechanism where the incoming substrate is activated by hydrogen bonding with its carbonyl oxygen to the Ser80 and Thr11 hydroxy groups. A hydrogen cyanide molecule subsequently replaces a water molecule and is deprotonated presumably by the His235 base. Deprotonation is facilitated by the proximity of the positive charge of the Lys236 side chain.     

Exploring the dynamic information content of a protein NMR structure: Comparison of a molecular dynamics simulation with the NMR and X‐ray structures of Escherichia coli ribonuclease HI

The multiconformer nature of solution nuclear magnetic resonance (NMR) structures of proteins results from the effects of intramolecular dynamics, spin diffusion and an uneven distribution of structural restraints throughout the molecule. A delineation of the former from the latter two contributions is attempted in this work for an ensemble of 15 NMR structures of the protein Escherichia coli ribonuclease HI (RNase HI). Exploration of the dynamic information content of the NMR ensemble is carried out through correlation with data from two crystal structures and a 1.7‐ns molecular dynamics (MD) trajectory of RNase HI in explicit solvent. Assessment of the consistency of the crystal and mean MD structures with nuclear Overhauser effect (NOE) data showed that the NMR ensemble is overall more compatible with the high‐resolution (1.48 Å) crystal structure than with either the lower‐resolution (2.05 Å) crystal structure or the MD simulation. Furthermore, the NMR ensemble is found to span more conformational space than the MD simulation for both the backbone and the sidechains of RNase HI. Nonetheless, the backbone conformational variability of both the NMR ensemble and the simulation is especially consistent with NMR relaxation measurements of two loop regions that are putative sites of substrate recognition. Plausible side‐chain dynamic information is extracted from the NMR ensemble on the basis of (i) rotamericity and syn‐pentane character of variable torsion angles, (ii) comparison of the magnitude of atomic mean‐square fluctuations (msf) with those deduced from crystallographic thermal factors, and (iii) comparison of torsion angle conformational behavior in the NMR ensemble and the simulation. Several heterogeneous torsion angles, while adopting non‐rotameric/syn‐pentane conformations in the NMR ensemble, exist in a unique conformation in the simulation and display low X‐ray thermal factors. These torsions are identified as sites whose variability is likely to be an artifact of the NMR structure determination procedure. A number of other torsions show a close correspondence between the conformations sampled in the NMR and MD ensembles, as well as significant correlations among crystallographic thermal factors and atomic msf calculated from the NMR ensemble and the simulation. These results indicate that a significant amount of dynamic information is contained in the NMR ensemble. The relevance of the present findings for the biological function of RNase HI, protein recognition studies, and previous investigations of the motional content of protein NMR structures are discussed. Proteins 1999;36:87–110. © 1999 Wiley‐Liss, Inc.     

Dynameomics: A Consensus View of the Protein Unfolding/Folding Transition State Ensemble across a Diverse Set of Protein Folds

The Dynameomics project aims to simulate a representative sample of all globular protein metafolds under both native and unfolding conditions. We have identified protein unfolding transition state (TS) ensembles from multiple molecular dynamics simulations of high-temperature unfolding in 183 structurally distinct proteins. These data can be used to study individual proteins and individual protein metafolds and to mine for TS structural features common across all proteins. Separating the TS structures into four different fold classes (all proteins, all-α, all-β, and mixed α/β and α + β) resulted in no significant difference in the overall protein properties. The residues with the most contacts in the native state lost the most contacts in the TS ensemble. On average, residues beginning in an α-helix maintained more structure in the TS ensemble than did residues starting in β-strands or any other conformation. The metafolds studied here represent 67% of all known protein structures, and this is, to our knowledge, the largest, most comprehensive study of the protein folding/unfolding TS ensemble to date. One might have expected broad distributions in the average global properties of the TS relative to the native state, indicating variability in the amount of structure present in the TS. Instead, the average global properties converged with low standard deviations across metafolds, suggesting that there are general rules governing the structure and properties of the TS.     

Hierarchical Conformational Analysis of Native Lysozyme Based on Sub-Millisecond Molecular Dynamics Simulations

Hierarchical organization of free energy landscape (FEL) for native globular proteins has been widely accepted by the biophysics community. However, FEL of native proteins is usually projected onto one or a few dimensions. Here we generated collectively 0.2 milli-second molecular dynamics simulation trajectories in explicit solvent for hen egg white lysozyme (HEWL), and carried out detailed conformational analysis based on backbone torsional degrees of freedom (DOF). Our results demonstrated that at micro-second and coarser temporal resolutions, FEL of HEWL exhibits hub-like topology with crystal structures occupying the dominant structural ensemble that serves as the hub of conformational transitions. However, at 100ns and finer temporal resolutions, conformational substates of HEWL exhibit network-like topology, crystal structures are associated with kinetic traps that are important but not dominant ensembles. Backbone torsional state transitions on time scales ranging from nanoseconds to beyond microseconds were found to be associated with various types of molecular interactions. Even at nanoseconds temporal resolution, the number of conformational substates that are of statistical significance is quite limited. These observations suggest that detailed analysis of conformational substates at multiple temporal resolutions is both important and feasible. Transition state ensembles among various conformational substates at microsecond temporal resolution were observed to be considerably disordered. Life times of these transition state ensembles are found to be nearly independent of the time scales of the participating torsional DOFs.     

Ensemble MD simulations restrained via crystallographic data: Accurate structure leads to accurate dynamics

Currently, the best existing molecular dynamics (MD) force fields cannot accurately reproduce the global free‐energy minimum which realizes the experimental protein structure. As a result, long MD trajectories tend to drift away from the starting coordinates (e.g., crystallographic structures). To address this problem, we have devised a new simulation strategy aimed at protein crystals. An MD simulation of protein crystal is essentially an ensemble simulation involving multiple protein molecules in a crystal unit cell (or a block of unit cells). To ensure that average protein coordinates remain correct during the simulation, we introduced crystallography‐based restraints into the MD protocol. Because these restraints are aimed at the ensemble‐average structure, they have only minimal impact on conformational dynamics of the individual protein molecules. So long as the average structure remains reasonable, the proteins move in a native‐like fashion as dictated by the original force field. To validate this approach, we have used the data from solid‐state NMR spectroscopy, which is the orthogonal experimental technique uniquely sensitive to protein local dynamics. The new method has been tested on the well‐established model protein, ubiquitin. The ensemble‐restrained MD simulations produced lower crystallographic R factors than conventional simulations; they also led to more accurate predictions for crystallographic temperature factors, solid‐state chemical shifts, and backbone order parameters. The predictions for ¹⁵N relaxation rates are at least as accurate as those obtained from conventional simulations. Taken together, these results suggest that the presented trajectories may be among the most realistic protein MD simulations ever reported. In this context, the ensemble restraints based on high‐resolution crystallographic data can be viewed as protein‐specific empirical corrections to the standard force fields.     

Conformational dynamics of human protein kinase CK2α and its effect on function and inhibition

Protein kinase, casein kinase II (CK2), is ubiquitously expressed and highly conserved protein kinase that shows constitutive activity. It phosphorylates a diverse set of proteins and plays crucial role in several cellular processes. The catalytic subunit of this enzyme (CK2α) shows remarkable flexibility as evidenced in numerous crystal structures determined till now. Here, using analysis of multiple crystal structures and long timescale molecular dynamics simulations, we explore the conformational flexibility of CK2α. The enzyme shows considerably higher flexibility in the solution as compared to that observed in crystal structure ensemble. Multiple conformations of hinge region, located near the active site, were observed during the dynamics. We further observed that among these multiple conformations, the most populated conformational state was inadequately represented in the crystal structure ensemble. The catalytic spine, was found to be less dismantled in this state as compared to the “open” hinge/αD state crystal structures. The comparison of dynamics in unbound (Apo) state and inhibitor (CX4945) bound state exhibits inhibitor induced suppression in the overall dynamics of the enzyme. This is especially true for functionally important glycine‐rich loop above the active site. Together, this work gives novel insights into the dynamics of CK2α in solution and relates it to the function. This work also explains the effect of inhibitor on the dynamics of CK2α and paves way for development of better inhibitors.