Cooperative symmetric to asymmetric conformational transition of the apo-form of scavenger decapping enzyme revealed by simulations

Decapping is a central step in eukaryotic mRNA turnover and in gene expression regulation. The human scavenger decapping enzyme, DcpS, catalyses cap hydrolysis following mRNA degradation. DcpS is a dimeric enzyme, with two active sites. Crystal structures suggest that DcpS must undergo significant conformational changes upon ligand binding, but the mechanism of this transition is unknown. Here, we report two long timescale (20 ns) molecular dynamics simulations of the apo-form of DcpS. The dimer is observed to undergo a strikingly cooperative motion, with one active site closing while the other opens. The amplitude of the conformational change is 6-21 A and the apparent timescale is 4-13 ns. These findings indicate that the crystallographically observed symmetric conformation of apo-form of DcpS is only a minor conformation in solution. The simulations also show that active sites are structurally connected via the domain-swapped dimer structure of the N-terminal domain, even in the absence of a bound ligand. These findings suggest a functional reason for the enzyme existing as a dimer, and may be widely relevant, also for other dimeric proteins.     

Classification and annotation of the relationship between protein structural change and ligand binding

The causal relationship between protein structural change and ligand binding was classified and annotated for 839 nonredundant pairs of crystal structures in the Protein Data Bank—one with and the other without a bound low-molecular-weight ligand molecule. Protein structural changes were first classified into either domain or local motions depending on the size of the moving protein segments. Whether the protein motion was coupled with ligand binding was then evaluated based on the location of the ligand binding site and by application of the linear response theory of protein structural change. Protein motions coupled with ligand binding were further classified into either closure or opening motions. This classification revealed the following: (i) domain motions coupled with ligand binding are dominated by closure motions, which can be described by the linear response theory; (ii) local motions frequently accompany order-disorder or α-helix-coil conformational transitions; and (iii) transferase activity (Enzyme Commission   number 2) is the predominant function among coupled domain closure motions. This could be explained by the closure motion acting to insulate the reaction site of these enzymes from environmental water.     

Coupling between Catalytic Loop Motions and Enzyme Global Dynamics

Catalytic loop motions facilitate substrate recognition and binding in many enzymes. While these motions appear to be highly flexible, their functional significance suggests that structure-encoded preferences may play a role in selecting particular mechanisms of motions. We performed an extensive study on a set of enzymes to assess whether the collective/global dynamics, as predicted by elastic network models (ENMs), facilitates or even defines the local motions undergone by functional loops. Our dataset includes a total of 117 crystal structures for ten enzymes of different sizes and oligomerization states. Each enzyme contains a specific functional/catalytic loop (10–21 residues long) that closes over the active site during catalysis. Principal component analysis (PCA) of the available crystal structures (including apo and ligand-bound forms) for each enzyme revealed the dominant conformational changes taking place in these loops upon substrate binding. These experimentally observed loop reconfigurations are shown to be predominantly driven by energetically favored modes of motion intrinsically accessible to the enzyme in the absence of its substrate. The analysis suggests that robust global modes cooperatively defined by the overall enzyme architecture also entail local components that assist in suitable opening/closure of the catalytic loop over the active site.     

Filtering molecular dynamics trajectories to reveal low-frequency collective motions: phospholipase A2

A novel method for analysing molecular dynamics trajectories has been developed, which filters out high frequencies using digital signal processing techniques and facilitates focusing on the low-frequency collective motions of proteins. These motions involve low energy slow motions, which lead to important biological phenomena such as domain closure and allosteric effects in enzymes. The filtering method treats each of the atomic trajectories obtained from the molecular dynamics simulation as a "signal". The trajectories of each of the atoms in the system (or any subset of interest) are Fourier transformed to the frequency domain, a filtering function is applied and then an inverse transformation back to the time domain yields the filtered trajectory. The filtering method has been used to study the dynamics of the enzyme phospholipase A2. In the filtered trajectory, all the high frequency bond and valence angle vibrations were eliminated, leaving only low-frequency motion, mainly fluctuations in torsions and conformational transitions. Analysis of this trajectory revealed interesting motions of the protein, including concerted movements of helices, and changes in shape of the active site cavity. Unlike normal mode analysis, which has been used to study the motion of proteins, this method does not require converged minimizations or diagonalization of a matrix of second derivatives. In addition, anharmonicity, multiple minima and conformational transitions are treated explicitly. Thus, the filtering method avoids most of the approximations implicit in other investigations of the dynamic behaviour of large systems.     

Molecular dynamics simulations of arachidonic acid complexes with COX-1 and COX-2: insights into equilibrium behavior

The cyclooxygenase (COX) enzymes are responsible for the committed step in prostaglandin biosynthesis, the generation of prostaglandin H(2). As a result, these enzymes are pharmacologically important targets for nonsteroidal antiinflammatory drugs, such as aspirin and newer COX-2 selective inhibitors. The cyclooxygenases are functional homodimers, and each subunit contains both a cyclooxygenase and a peroxidase active site. These enzymes are quite interesting mechanistically, as the conversion of arachidonic acid to prostaglandin H(2) requires two oxygenation and two cyclization reactions, resulting in the formation of five new chiral centers with nearly absolute regio- and stereochemical fidelity. We have used molecular dynamics (MD) simulations to investigate the equilibrium behavior of both COX-1 and COX-2 enzyme isoforms with bound arachidonate. These simulations were compared with reference simulations of arachidonate in solution to explore the effect of enzyme on substrate conformation and positioning in the active site. The simulations suggest that the substrate has greater conformational freedom in the COX-2 active site, consistent with the larger COX-2 active site volume observed in X-ray crystal structures. The simulations reveal different conformational behavior for arachidonate in each subunit over the course of extended equilibrium MD simulations. The simulations also provide detailed information for several protein channels that might be important for oxygen and water transport to or from active sites or for intermediate trafficking between the cyclooxygenase and peroxidase active sites. The detailed comparisons for COX-1 versus COX-2 active site structural fluctuations may also provide useful information for design of new isozyme-selective inhibitors.     

Insights into correlated motions and long-range interactions in CheY derived from molecular dynamics simulations

CheY is a response regulator protein involved in bacterial chemotaxis. Much is known about its active and inactive conformations, but little is known about the mechanisms underlying long-range interactions or correlated motions. To investigate these events, molecular dynamics simulations were performed on the unphosphorylated, inactive structure from Salmonella typhimurium and the CheY-BeF(3)(-) active mimic structure (with BeF(3)(-) removed) from Escherichia coli. Simulations utilized both sequences in each conformation to discriminate sequence- and structure-specific behavior. The previously identified conformational differences between the inactive and active conformations of the strand-4-helix-4 loop, which are present in these simulations, arise from the structural, and not the sequence, differences. The simulations identify previously unreported structure-specific flexibility features in this loop and sequence-specific flexibility features in other regions of the protein. Both structure- and sequence-specific long-range interactions are observed in the active and inactive ensembles. In the inactive ensemble, two distinct mechanisms based on Thr-87 or Ile-95 rotameric forms, are observed for the previously identified g+ and g- rotamer sampling by Tyr-106. These molecular dynamics simulations have thus identified both sequence- and structure-specific differences in flexibility, long-range interactions, and rotameric form of key residues. Potential biological consequences of differential flexibility and long-range correlated motion are discussed.     

Structures of mesophilic and extremophilic citrate synthases reveal rigidity and flexibility for function

Citrate synthase (CS) catalyses the entry of carbon into the citric acid cycle and is highly‐conserved structurally across the tree of life. Crystal structures of dimeric CSs are known in both “open” and “closed” forms, which differ by a substantial domain motion that closes the substrate‐binding clefts. We explore both the static rigidity and the dynamic flexibility of CS structures from mesophilic and extremophilic organisms from all three evolutionary domains. The computational expense of this wide‐ranging exploration is kept to a minimum by the use of rigidity analysis and rapid all‐atom simulations of flexible motion, combining geometric simulation and elastic network modeling. CS structures from thermophiles display increased structural rigidity compared with the mesophilic enzyme. A CS structure from a psychrophile, stabilized by strong ionic interactions, appears to display likewise increased rigidity in conventional rigidity analysis; however, a novel modified analysis, taking into account the weakening of the hydrophobic effect at low temperatures, shows a more appropriate decreased rigidity. These rigidity variations do not, however, affect the character of the flexible dynamics, which are well conserved across all the structures studied. Simulation trajectories not only duplicate the crystallographically observed symmetric open‐to‐closed transitions, but also identify motions describing a previously unidentified antisymmetric functional motion. This antisymmetric motion would not be directly observed in crystallography but is revealed as an intrinsic property of the CS structure by modeling of flexible motion. This suggests that the functional motion closing the binding clefts in CS may be independent rather than symmetric and cooperative. Proteins 2014; 82:2657–2670. © 2014 Wiley Periodicals, Inc.     

Projection of Monte Carlo and molecular dynamics trajectories onto the normal mode axes: human lysozyme

A method is presented to describe the internal motions of proteins obtained from molecular dynamics or Monte Carlo simulations as motions of normal mode variables. This method calculates normal mode variables by projecting trajectories of these simulations onto the axes of normal modes and expresses the trajectories as a linear combination of normal mode variables. This method is applied to the result of the molecular dynamics and the Monte Carlo simulations of human lysozyme. The motion of the lowest frequency mode extracted from the simulations represents the hinge bending motion very faithfully. Analysis of the obtained motions of the normal mode variables provides an explanation of the anharmonic aspects of protein dynamics as due first to the anharmonicity of the actual potential energy surface near a minimum and second to trans-minimum conformational changes.     

Protein dynamics and enzyme catalysis: insights from simulations

The role of protein dynamics in enzyme catalysis is one of the most active and controversial areas in enzymology today. Some researchers claim that protein dynamics are at the heart of enzyme catalytic efficiency, while others state that dynamics make no significant contribution to catalysis. What is the biochemist - or student - to make of the ferocious arguments in this area? Protein dynamics are complex and fascinating, as molecular dynamics simulations and experiments have shown. The essential question is: do these complex motions have functional significance? In particular, how do they affect or relate to chemical reactions within enzymes, and how are chemical and conformational changes coupled together? Biomolecular simulations can analyse enzyme reactions and dynamics in atomic detail, beyond that achievable in experiments: accurate atomistic modelling has an essential part to play in clarifying these issues. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.     

Conformational regulation of the EGFR kinase core by the juxtamembrane and C-terminal tail: a molecular dynamics study

The catalytic domain of epidermal growth factor receptor (EGFR) is activated by dimerization, which requires allosteric coupling between distal dimerization and catalytic sites. Although crystal structures of EGFR kinases, solved in various conformational states, have provided important insights into EGFR activation by dimerization, the atomic details of how dimerization signals are dynamically coupled to catalytic regions of the kinase core are not fully understood. In this study, we have performed unrestrained and targeted molecular dynamics simulations on the active and inactive states of EGFR, followed by principal component analysis on the simulated trajectories, to identify correlated motions in the EGFR kinase domain upon dimerization. Our analysis reveals that the conformational changes associated with the catalytic functions of the kinase core are highly correlated with motions in the juxtamembrane (JM) and C-terminal tail, two flexible structural elements that play an active role in EGFR kinase activation and dimerization. In particular, the opening and closing of the ATP binding lobe relative to the substrate binding lobe is highly correlated with motions in the JM and C-terminal tail, suggesting that ATP and substrate binding can be coordinated with dimerization through conformational changes in the JM and C-terminal tail. Our study pinpoints key residues involved in this conformational coupling, and provides new insights into the role of the JM and C-terminal tail segments in EGFR kinase functions.     

Engineered control of enzyme structural dynamics and function

Enzymes undergo a range of internal motions from local, active site fluctuations to large‐scale, global conformational changes. These motions are often important for enzyme function, including in ligand binding and dissociation and even preparing the active site for chemical catalysis. Protein engineering efforts have been directed towards manipulating enzyme structural dynamics and conformational changes, including targeting specific amino acid interactions and creation of chimeric enzymes with new regulatory functions. Post‐translational covalent modification can provide an additional level of enzyme control. These studies have not only provided insights into the functional role of protein motions, but they offer opportunities to create stimulus‐responsive enzymes. These enzymes can be engineered to respond to a number of external stimuli, including light, pH, and the presence of novel allosteric modulators. Altogether, the ability to engineer and control enzyme structural dynamics can provide new tools for biotechnology and medicine.     

A Direct Coupling between Global and Internal Motions in a Single Domain Protein? MD Investigation of Extreme Scenarios

Proteins are not rigid molecules, but exhibit internal motions on timescales ranging from femto- to milliseconds and beyond. In solution, proteins also experience global translational and rotational motions, sometimes on timescales comparable to those of the internal fluctuations. The possibility that internal and global motions may be directly coupled has intriguing implications, given that enzymes and cell signaling proteins typically associate with binding partners and cellular scaffolds. Such processes alter their global motion and may affect protein function. Here, we present molecular dynamics simulations of extreme case scenarios to examine whether a possible relationship exists. In our model protein, a ubiquitin-like RhoGTPase binding domain of plexin-B1, we removed either internal or global motions. Comparisons with unrestrained simulations show that internal and global motions are not appreciably coupled in this single-domain protein. This lack of coupling is consistent with the observation that the dynamics of water around the protein, which is thought to permit, if not stimulate, internal dynamics, is also largely independent of global motion. We discuss implications of these results for the structure and function of proteins.     

Conformational change in an MFS protein: MD simulations of LacY

Molecular dynamics simulations of lactose permease (LacY) in a phospholipid bilayer reveal the conformational dynamics of the protein. In inhibitor-bound simulations (i.e., those closest to the X-ray structure) the protein was stable, showing little conformational change over a 50 ns timescale. Movement of the bound inhibitor, TDG, to an alternative binding mode was observed, so that it interacted predominantly with the N-terminal domain and with residue E269 from the C-terminal domain. In multiple ligand-free simulations, a degree of domain closure occurred. This switched LacY to a state with a central cavity closed at both the intracellular and periplasmic ends. This may resemble a possible intermediate in the transport mechanism. Domain closure occurs by a combination of rigid-body movements of domains and of intradomain motions of helices, especially TM4, TM5, TM10, and TM11. A degree of intrahelix flexibility appears to be important in the conformational change.     

Nanoseconds molecular dynamics simulation of primary mechanical energy transfer steps in F1-ATP synthase

The mitochondrial membrane protein FoF1-ATP synthase synthesizes adenosine triphosphate (ATP), the universal currency of energy in the cell. This process involves mechanochemical energy transfer from a rotating asymmetric gamma-'stalk' to the three active sites of the F1 unit, which drives the bound ATP out of the binding pocket. Here, the primary structural changes associated with this energy transfer in F1-ATP synthase were studied with multi-nanosecond molecular dynamics simulations. By forced rotation of the gamma-stalk that mimics the effect of proton motive Fo-rotation during ATP synthesis, a time-resolved atomic model for the structural changes in the F1 part in terms of propagating conformational motions is obtained. For these, different time scales are found, which allows the separation of nanosecond from microsecond conformational motions. In the simulations, rotation of the gamma-stalk lowers the ATP affinity of the betaTP binding pocket and triggers fast, spontaneous closure of the empty betaE subunit. The simulations explain several mutation studies and the reduced hydrolysis rate of gamma-depleted F1-ATPase.     

Opening and closing motions in the periplasmic vitamin B12 binding protein BtuF

BtuF is the periplasmic binding protein (PBP) in the vitamin B(12) uptake system in Escherichia coli where it is associated with the ABC transporter BtuCD. When the ligand binds, PBPs generally display large conformational changes, commonly termed the Venus flytrap mechanism. BtuF belongs to a group of PBPs that, on the basis of crystal structures, does not appear to display such behavior. Using 480 ns multicopy molecular dynamics simulations of apo and holo forms of the protein, we investigate the dynamics of BtuF. We find BtuF to be more flexible than previously assumed, displaying clear opening and closing motions which are more pronounced in the apo form. The protein behavior is compatible with a PBP functional model that postulates a closed conformation for the ligand-bound state, whereas the empty form fluctuates between open and closed conformations. Elastic network normal-mode analysis suggests that all BtuF-like PBPs are capable of similar opening and closing motions. It also makes the typical Venus flytrap domain motions a likely common means of how PBP-ABC transporter interaction could occur.     

Binding leverage as a molecular basis for allosteric regulation

Allosteric regulation involves conformational transitions or fluctuations between a few closely related states, caused by the binding of effector molecules. We introduce a quantity called binding leverage that measures the ability of a binding site to couple to the intrinsic motions of a protein. We use Monte Carlo simulations to generate potential binding sites and either normal modes or pairs of crystal structures to describe relevant motions. We analyze single catalytic domains and multimeric allosteric enzymes with complex regulation. For the majority of the analyzed proteins, we find that both catalytic and allosteric sites have high binding leverage. Furthermore, our analysis of the catabolite activator protein, which is allosteric without conformational change, shows that its regulation involves other types of motion than those modulated at sites with high binding leverage. Our results point to the importance of incorporating dynamic information when predicting functional sites. Because it is possible to calculate binding leverage from a single crystal structure it can be used for characterizing proteins of unknown function and predicting latent allosteric sites in any protein, with implications for drug design.     

Functionally important correlated motions in the single-stranded DNA-binding protein encoded by filamentous phage Pf3

To elucidate the interplay between different parts of dimeric single-stranded DNA-binding proteins we have studied the correlated motions in the protein encoded by filamentous phage Pf3 via the combined use of 15N-NMR relaxation experiments, molecular dynamics simulations and essential dynamics calculations. These studies provide insight into the mechanism underlying the protein-DNA binding reaction. The most important motions can be described by a few essential modes. Most outstanding is the correlated symmetric motion of the DNA-binding wings, which are far apart in the structure. This motion determines the access of DNA to the DNA-binding domain. A correlation between the motion of the DNA-binding wing and the complex loop is indicated to play a role in the cooperative binding of the protein to DNA. These motions are in the nanosecond regime in correspondence with the 15N-NMR relaxation experiments.     

Conformational State Distributions and Catalytically Relevant Dynamics of a Hinge-Bending Enzyme Studied by Single-Molecule FRET and a Coarse-Grained Simulation

Over the last few decades, a view has emerged showing that multidomain enzymes are biological machines evolved to harness stochastic kicks of solvent particles into highly directional functional motions. These intrinsic motions are structurally encoded, and Nature makes use of them to catalyze chemical reactions by means of ligand-induced conformational changes and states redistribution. Such mechanisms align reactive groups for efficient chemistry and stabilize conformers most proficient for catalysis. By combining single-molecule Förster resonance energy transfer measurements with normal mode analysis and coarse-grained mesoscopic simulations, we obtained results for a hinge-bending enzyme, namely phosphoglycerate kinase (PGK), which support and extend these ideas. From single-molecule Förster resonance energy transfer, we obtained insight into the distribution of conformational states and the dynamical properties of the domains. The simulations allowed for the characterization of interdomain motions of a compact state of PGK. The data show that PGK is intrinsically a highly dynamic system sampling a wealth of conformations on timescales ranging from nanoseconds to milliseconds and above. Functional motions encoded in the fold are performed by the PGK domains already in its ligand-free form, and substrate binding is not required to enable them. Compared to other multidomain proteins, these motions are rather fast and presumably not rate-limiting in the enzymatic reaction. Ligand binding slightly readjusts the orientation of the domains and feasibly locks the protein motions along a preferential direction. In addition, the functionally relevant compact state is stabilized by the substrates, and acts as a prestate to reach active conformations by means of Brownian motions.     

Ribosomal RNA kink-turn motif--a flexible molecular hinge

Ribosomal RNA K-turn motifs are asymmetric internal loops characterized by a sharp bend in the phosphodiester backbone resulting in "V" shaped structures, recurrently observed in ribosomes and showing a high degree of sequence conservation. We have carried out extended explicit solvent molecular dynamics simulations of selected K-turns, in order to investigate their intrinsic structural and dynamical properties. The simulations reveal an unprecedented dynamical flexibility of the K-turns around their X-ray geometries. The K-turns sample, on the nanosecond timescale, different conformational substates. The overall behavior of the simulations suggests that the sampled geometries are essentially isoenergetic and separated by minimal energy barriers. The nanosecond dynamics of isolated K-turns can be qualitatively considered as motion of two rigid helix stems controlled by a very flexible internal loop which then leads to substantial hinge-like motions between the two stems. This internal dynamics of K-turns is strikingly different for example from the bacterial 5S rRNA Loop E motif or BWYV frameshifting pseudoknot which appear to be rigid in the same type of simulations. Bistability and flexibility of K-turns was also suggested by several recent biochemical studies. Although the results of MD simulations should be considered as a qualitative picture of the K-turn dynamics due to force field and sampling limitations, the main advantage of the MD technique is its ability to investigate the region close to K-turn ribosomal-like geometries. This part of the conformational space is not well characterized by the solution experiments due to large-scale conformational changes seen in the experiments. We suggest that K-turns are well suited to act as flexible structural elements of ribosomal RNA. They can for example be involved in mediation of large-scale motions or they can allow a smooth assembling of the other parts of the ribosome.     

Probing the Periplasmic-Open State of Lactose Permease in Response to Sugar Binding and Proton Translocation

Based on the crystal structure of lactose permease (LacY) open to the cytoplasm, a hybrid molecular simulation approach with self-guided Langevin dynamics is used to describe conformational changes that lead to a periplasmic-open state. This hybrid approach consists of implicit (IM) and explicit (EX) membrane simulations and requires self-guided Langevin dynamics to enhance protein motions during the IM simulations. The pore radius of the lumen increases by 3.5 Å on the periplasmic side and decreases by 2.5 Å on the cytoplasmic side (relative to the crystal structure), suggesting a lumen that is fully open to the periplasm to allow for extracellular sugar transport and closed to the cytoplasm. Based on our simulations, the mechanism that triggers this conformational change to the periplasmic-open state is the protonation of Glu269 and binding of the disaccharide. Then, helix packing is destabilized by breaking of several side chains involved in hydrogen bonding (Asn245, Ser41, Glu374, Lys42, and Gln242). For the periplasmic-open conformations obtained from our simulations, helix-helix distances agree well with experimental measurements using double electron-electron resonance, fluorescence resonance energy transfer, and varying sized cross-linkers. The periplasmic-open conformations are also in compliance with various substrate accessibility/reactivity measurements that indicate an opening of the protein lumen on the periplasmic side on sugar binding. The comparison with these measurements suggests a possible incomplete closure of the cytoplasmic half in our simulations. However, the closure is sufficient to prevent the disaccharide from transporting to the cytoplasm, which is in accordance with the well-established alternating access model. Ser53, Gln60, and Phe354 are determined to be important in sugar transport during the periplasmic-open stage of the sugar transport cycle and the sugar is found to undergo an orientational change in order to escape the protein lumen.