Role Of Protein Dynamics Research Articles

The plastoquinol-plastoquinone exchange mechanism in photosystem II: insight from molecular dynamics simulations.

In the photosystem II (PSII) of oxygenic photosynthetic organisms, the reaction center (RC) core mediates the light-induced electron transfer leading to water splitting and production of reduced plastoquinone molecules. The reduction of plastoquinone to plastoquinol lowers PSII affinity for the latter and leads to its release. However, little is known about the role of protein dynamics in this process. Here, molecular dynamics simulations of the complete PSII complex embedded in a lipid bilayer have been used to investigate the plastoquinol release mechanism. A distinct dynamic behavior of PSII in the presence of plastoquinol is observed which, coupled to changes in charge distribution and electrostatic interactions, causes disruption of the interactions seen in the PSII-plastoquinone complex and leads to the "squeezing out" of plastoquinol from the binding pocket. Displacement of plastoquinol closes the second water channel, recently described in a 2.9Å resolution PSII structure (Guskov et al. in Nat Struct Mol Biol 16:334-342, 2009), allowing to rule out the proposed "alternating" mechanism of plastoquinol-plastoquinone exchange, while giving support to the "single-channel" one. The performed simulations indicated a pivotal role of D1-Ser264 in modulating the dynamics of the plastoquinone binding pocket and plastoquinol-plastoquinone exchange via its interaction with D1-His252 residue. The effects of the disruption of this hydrogen bond network on the PSII redox reactions were experimentally assessed in the D1 site-directed mutant Ser264Lys.

Dynamics and mechanism of ultrafast water–protein interactions

Protein hydration is essential to its structure, dynamics, and function, but water-protein interactions have not been directly observed in real time at physiological temperature to our awareness. By using a tryptophan scan with femtosecond spectroscopy, we simultaneously measured the hydration water dynamics and protein side-chain motions with temperature dependence. We observed the heterogeneous hydration dynamics around the global protein surface with two types of coupled motions, collective water/side-chain reorientation in a few picoseconds and cooperative water/side-chain restructuring in tens of picoseconds. The ultrafast dynamics in hundreds of femtoseconds is from the outer-layer, bulk-type mobile water molecules in the hydration shell. We also found that the hydration water dynamics are always faster than protein side-chain relaxations but with the same energy barriers, indicating hydration shell fluctuations driving protein side-chain motions on the picosecond time scales and thus elucidating their ultimate relationship.

Atomistic details of protein dynamics and the role of hydration water

BackgroundThe importance of protein dynamics for their biological activity is now well recognized. Different experimental and computational techniques have been employed to study protein dynamics, hierarchy of different processes and the coupling between protein and hydration water dynamics. Yet, understanding the atomistic details of protein dynamics and the role of hydration water remains rather limited. Scoop of reviewBased on overview of neutron scattering, molecular dynamic simulations, NMR and dielectric spectroscopy results we present a general picture of protein dynamics covering time scales from faster than ps to microseconds and the influence of hydration water on different relaxation processes. Major conclusionsInternal protein dynamics spread over a wide time range from faster than picosecond to longer than microseconds. We suggest that the structural relaxation in hydrated proteins appears on the microsecond time scale, while faster processes present mostly motion of side groups and some domains. Hydration water plays a crucial role in protein dynamics on all time scales. It controls the coupled protein-hydration water relaxation on 10–100ps time scale. This process defines the friction for slower protein dynamics. Analysis suggests that changes in amount of hydration water affect not only general friction, but also influence significantly the protein's energy landscape. General significanceThe proposed atomistic picture of protein dynamics provides deeper understanding of various relaxation processes and their hierarchy, similarity and differences between various biological macromolecules, including proteins, DNA and RNA.This article is part of a Special Issue entitled “Science for Life” Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo”.

Structural and Kinetic Studies of Formate Dehydrogenase from Candida boidinii.

The structure of formate dehydrogenase from Candida boidinii (CbFDH) is of both academic and practical interests. First, this enzyme represents a unique model system for studies on the role of protein dynamics in catalysis, but so far these studies have been limited by the availability of structural information. Second, CbFDH and its mutants can be used in various industrial applications (e.g., CO2 fixation or nicotinamide recycling systems), and the lack of structural information has been a limiting factor in commercial development. Here, we report the crystallization and structural determination of both holo- and apo-CbFDH. The free-energy barrier for the catalyzed reaction was computed and indicates that this structure indeed represents a catalytically competent form of the enzyme. Complementing kinetic examinations demonstrate that the recombinant CbFDH has a well-organized reactive state. Finally, a fortuitous observation has been made: the apoenzyme crystal was obtained under cocrystallization conditions with a saturating concentration of both the cofactor (NAD(+)) and inhibitor (azide), which has a nanomolar dissociation constant. It was found that the fraction of the apoenzyme present in the solution is less than 1.7 × 10(-7) (i.e., the solution is 99.9999% holoenzyme). This is an extreme case where the crystal structure represents an insignificant fraction of the enzyme in solution, and a mechanism rationalizing this phenomenon is presented.

Role of protein dynamics in ion selectivity and allosteric coupling in the NaK channel

Flux-dependent inactivation that arises from functional coupling between the inner gate and the selectivity filter is widespread in ion channels. The structural basis of this coupling has only been well characterized in KcsA. Here we present NMR data demonstrating structural and dynamic coupling between the selectivity filter and intracellular constriction point in the bacterial nonselective cation channel, NaK. This transmembrane allosteric communication must be structurally different from KcsA because the NaK selectivity filter does not collapse under low-cation conditions. Comparison of NMR spectra of the nonselective NaK and potassium-selective NaK2K indicates that the number of ion binding sites in the selectivity filter shifts the equilibrium distribution of structural states throughout the channel. This finding was unexpected given the nearly identical crystal structure of NaK and NaK2K outside the immediate vicinity of the selectivity filter. Our results highlight the tight structural and dynamic coupling between the selectivity filter and the channel scaffold, which has significant implications for channel function. NaK offers a distinct model to study the physiologically essential connection between ion conduction and channel gating.

Role of Protein Dynamics in Allosteric Control of the Catalytic Phosphoryl Transfer of Insulin Receptor Kinase.

The catalytic and allosteric mechanisms of insulin receptor kinase (IRK) are investigated by a combination of ab initio and semiempirical quantum mechanical and molecular mechanical (QM/MM) methods and classical molecular dynamics (MD) simulations. The simulations reveal that the catalytic reaction proceeds in two steps, starting with the transfer of a proton from substrate Tyr to the catalytic Asp1132, followed by the phosphoryl transfer from ATP to substrate Tyr. The enhancement of the catalytic rate of IRK upon phosphorylations in the enzyme's activation loop is found to occur mainly via changes to the free energy landscape of the proton transfer step, favoring the proton transfer in the fully phosphorylated enzyme. In contrast, the effects of the phosphorylations on the phosphoryl transfer are smaller. Equilibrium MD simulations show that IRK phosphorylations affect the protein dynamics of the enzyme before the proton transfer to Asp1132 with only a minor effect after the proton transfer. This finding is consistent with the large change in the proton transfer free energy and the smaller change in the free energy barrier of phosphoryl transfer found by QM/MM simulations. Taken together, the present results provide details on how IRK phosphorylation exerts allosteric control of the catalytic activity via modifications of protein dynamics and free energy landscape of catalytic reaction. The results also highlight the importance of protein dynamics in connecting protein allostery and catalysis to control catalytic activity of enzymes.

Wide‐angle X‐ray Solution Scattering Studies of Slow Correlated Motions in Proteins

Slow correlated motions in the submicrosecond to second range are particularly important for protein function but remain challenging to study. Wide‐angle x‐ray solution scattering (WAXS) generates information about the length scale of intramolecular motions difficult to obtain by other methods; is highly sensitive to these motions, being able to detect changes due to single amino acid replacements, ligand binding or alterations in solution conditions; and is sufficiently fast as to make practical the screening of hundreds of conditions/ligands/mutations to identify those that impact protein motions.Direct comparison of the results of WAXS and molecular dynamics (MD) is accomplished through introduction of the sigma‐r plot, a plot of the standard deviation of interatomic vector lengths as a function of that length ‐ derivable from WAXS data and MD trajectories. MD simulations generate extraordinary detail about protein trajectories; WAXS provides experimental measures of length scale of protein motion under a large number of conditions. Direct comparison of WAXS with MD will enable detailed molecular interpretations of WAXS data using the power of MD; while making possible the extrapolation of MD results to understand molecular behavior of a large number of mutants in multiple ligation states and solution conditions. Application of this approach to HIV protease provides new insights into the role of protein dynamics in its enzymatic function.

Dynamics of protein kinases: insights from nuclear magnetic resonance.

Protein kinases are ubiquitous enzymes with critical roles in cellular processes and pathology. As a result, researchers have studied their activity and regulatory mechanisms extensively. Thousands of X-ray structures give snapshots of the architectures of protein kinases in various states of activation and ligand binding. However, the extent of and manner by which protein motions and conformational dynamics underlie the function and regulation of these important enzymes is not well understood. Nuclear magnetic resonance (NMR) methods provide complementary information about protein conformation and dynamics in solution. However, until recently, the large size of these enzymes prevented researchers from using these methods with kinases. Developments in transverse relaxation-optimized spectroscopy (TROSY)-based techniques and more efficient isotope labeling strategies are now allowing researchers to carry out NMR studies on full-length protein kinases. In this Account, we describe recent insights into the role of dynamics in protein kinase regulation and catalysis that have been gained from NMR measurements of chemical shift changes and line broadening, residual dipolar couplings, and relaxation. These findings show strong associations between protein motion and events that control kinase activity. Dynamic and conformational changes occurring at ligand binding sites and other regulatory domains of these proteins propagate to conserved kinase core regions that mediate catalytic function. NMR measurements of slow time scale (microsecond to millisecond) motions also reveal that kinases carry out global exchange processes that synchronize multiple residues and allosteric interconversion between conformational states. Activating covalent modifications or ligand binding to form the Michaelis complex can induce these global processes. Inhibitors can also exploit the exchange properties of kinases by using conformational selection to form dynamically quenched states. These investigations have revealed that kinases are highly dynamic enzymes, whose regulation by interdomain interactions, ligand binding, and covalent modifications involve changes in motion and conformational equilibrium in a manner that can be correlated with function. Thus, NMR provides a unique window into the role of protein dynamics in kinase regulation and catalysis with important implications for drug design.

Toward fluorescent probes for G-protein-coupled receptors (GPCRs).

G-protein-coupled receptors (GPCRs), a superfamily of cell-surface receptors that are the targets of about 40% of prescription drugs on the market, can sense numerous critical extracellular signals. Recent breakthroughs in structural biology, especially in holo-form X-ray crystal structures, have contributed to our understanding of GPCR signaling. However, actions of GPCRs at the cellular and molecular level, interactions between GPCRs, and the role of protein dynamics in receptor activities still remain controversial. To overcome these dilemmas, fluorescent probes of GPCRs have been employed, which have advantages of in vivo safety and real-time monitoring. Various probes that depend on specific mechanisms and/or technologies have been used to study GPCRs. The present review focuses on surveying the design and applications of fluorescent probes for GPCRs that are derived from small molecules or using protein-labeling techniques, as well as discussing some design strategies for new probes.

The role of protein dynamics and thermal fluctuations in regulating cytochrome c/cytochrome c oxidase electron transfer

The role of protein dynamics and thermal fluctuations in regulating cytochrome c/cytochrome c oxidase electron transfer

Different dynamical effects in mesophilic and hyperthermophilic dihydrofolate reductases.

The role of protein dynamics in the reaction catalyzed by dihydrofolate reductase from the hyperthermophile Thermotoga maritima (TmDHFR) has been examined by enzyme isotope substitution (15N, 13C, 2H). In contrast to all other enzyme reactions investigated previously, including DHFR from Escherichia coli (EcDHFR), for which isotopic substitution led to decreased reactivity, the rate constant for the hydride transfer step is not affected by isotopic substitution of TmDHFR. TmDHFR therefore appears to lack the coupling of protein motions to the reaction coordinate that have been identified for EcDHFR catalysis. Clearly, dynamical coupling is not a universal phenomenon that affects the efficiency of enzyme catalysis.

Mechanistic studies of different thymidylate synthases and the role of protein dynamics in catalysis (110.1)

Thymidylate synthases (TSases) catalyze the last committed step in the de novo biosynthesis of thymidylate, one of the four DNA bases (T). Classical TSases catalyze a complex cascade, and we developed methods to study the physical nature of its mechanism, and to test the role of enzyme dynamics in catalysis. In a related study, a newly discovered family of TSases was examined. We used a combination of isotopic labeling, mutations, quench-flow, and X-ray crystallography (Fig. 1), and found an unusual chemical mechanism for these flavin-dependent thymidylate synthases (FDTS). Figure 1. An active site view of crystal structures of TmFDTS in complex with FAD, dUMP, and the co-factor CH2H4folate.

Contraction-Induced Changes in Hydrogen Bonding of Muscle Hydration Water.

Protein–water interaction plays a crucial role in protein dynamics and hence function. To study the chemical environment of water and proteins with high spatial resolution, synchrotron radiation-Fourier transform infrared (SR-FTIR) spectromicroscopy was used to probe skeletal muscle myofibrils. Observing the OH stretch band showed that water inside of relaxed myofibrils is extensively hydrogen-bonded with little or no free OH. In higher-resolution measurements obtained with single isolated myofibrils, the water absorption peaks were relatively higher within the center region of the sarcomere compared to those in the I-band region, implying higher hydration capacity of thick filaments compared to the thin filaments. When specimens were activated, changes in the OH stretch band showed significant dehydrogen bonding of muscle water; this was indicated by increased absorption at ∼3480 cm–1 compared to relaxed myofibrils. These contraction-induced changes in water were accompanied by splitting of the amide I (C=O) peak, implying that muscle proteins transition from α-helix to β-sheet-rich structures. Hence, muscle contraction can be characterized by a loss of order in the muscle–protein complex, accompanied by a destructuring of hydration water. The findings shed fresh light on the molecular mechanism of muscle contraction and motor protein dynamics.

Effects of Cavities at the Nicotinamide Binding Site of Liver Alcohol Dehydrogenase on Structure, Dynamics and Catalysis

A role for protein dynamics in enzymatic catalysis of hydrogen transfer has received substantial scientific support, but the connections between protein structure and catalysis remain to be established. Valine residues 203 and 207 are at the binding site for the nicotinamide ring of the coenzyme in liver alcohol dehydrogenase and have been suggested to facilitate catalysis with “protein-promoting vibrations” (PPV). We find that the V207A substitution has small effects on steady-state kinetic constants and the rate of hydrogen transfer; the introduced cavity is empty and is tolerated with minimal effects on structure (determined at 1.2 Å for the complex with NAD+ and 2,3,4,5,6-pentafluorobenzyl alcohol). Thus, no evidence is found to support a role for Val-207 in the dynamics of catalysis. The protein structures and ligand geometries (including donor–acceptor distances) in the V203A enzyme complexed with NAD+ and 2,3,4,5,6-pentafluorobenzyl alcohol or 2,2,2-trifluoroethanol (determined at 1.1 Å) are very similar to those for the wild-type enzyme, except that the introduced cavity accommodates a new water molecule that contacts the nicotinamide ring. The structures of the V203A enzyme complexes suggest, in contrast to previous studies, that the diminished tunneling and decreased rate of hydride transfer (16-fold, relative to that of the wild-type enzyme) are not due to differences in ground-state ligand geometries. The V203A substitution may alter the PPV and the reorganization energy for hydrogen transfer, but the protein scaffold and equilibrium thermal motions within the Michaelis complex may be more significant for enzyme catalysis.

Escherichia coli SufE Sulfur Transfer Protein Modulates the SufS Cysteine Desulfurase through Allosteric Conformational Dynamics

Fe-S clusters are critical metallocofactors required for cell function. Fe-S cluster biogenesis is carried out by assembly machinery consisting of multiple proteins. Fe-S cluster biogenesis proteins work together to mobilize sulfide and iron, form the nascent cluster, traffic the cluster to target metalloproteins, and regulate the assembly machinery in response to cellular Fe-S cluster demand. A complex series of protein-protein interactions is required for the assembly machinery to function properly. Despite considerable progress in obtaining static three-dimensional structures of the assembly proteins, little is known about transient protein-protein interactions during cluster assembly or the role of protein dynamics in the cluster assembly process. The Escherichia coli cysteine desulfurase SufS (EC 2.8.1.7) and its accessory protein SufE work together to mobilize persulfide from L-cysteine, which is then donated to the SufB Fe-S cluster scaffold. Here we use amide hydrogen/deuterium exchange mass spectrometry (HDX-MS) to characterize SufS-SufE interactions and protein dynamics in solution. HDX-MS analysis shows that SufE binds near the SufS active site to accept persulfide from Cys-364. Furthermore, SufE binding initiates allosteric changes in other parts of the SufS structure that likely affect SufS catalysis and alter SufS monomer-monomer interactions. SufE enhances the initial l-cysteine substrate binding to SufS and formation of the external aldimine with pyridoxal phosphate required for early steps in SufS catalysis. Together, these results provide a new picture of the SufS-SufE sulfur transferase pathway and suggest a more active role for SufE in promoting the SufS cysteine desulfurase reaction for Fe-S cluster assembly.

Structural and dynamics characterization of norovirus protease

Norovirus protease is an essential enzyme for proteolytic maturation of norovirus nonstructural proteins and has been implicated as a potential target for antiviral drug development. Although X-ray structural studies of the protease give us wealth of structural information including interactions of the protease with its substrate and dimeric overall structure, the role of protein dynamics in the substrate recognition and the biological relevance of the protease dimer remain unclear. Here we determined the solution NMR structure of the 3C-like protease from Norwalk virus (NV 3CLpro), a prototype strain of norovirus, and analyzed its backbone dynamics and hydrodynamic behavior in solution. ¹⁵N spin relaxation and analytical ultracentrifugation analyses demonstrate that NV 3CLpro is predominantly a monomer in solution. Solution structure of NV 3CLpro shows significant structural variation in C-terminal domain compared with crystal structures and among lower energy structure ensembles. Also, ¹⁵N spin relaxation and Carr-Purcell-Meiboom-Gill (CPMG)-based relaxation dispersion analyses reveal the dynamic properties of residues in the C-terminal domain over a wide range of timescales. In particular, the long loop spanning residues T123-G133 show fast motion (ps-ns), and the residues in the bII-cII region forming the large hydrophobic pocket (S2 site) undergo conformational exchanges on slower timescales (μs-ms), suggesting their important role in substrate recognition.

Mass Modulation of Protein Dynamics Associated with Barrier Crossing in Purine Nucleoside Phosphorylase

The role of protein dynamics on different time scales in enzyme catalysis remains an area of active debate. The connection between enzyme dynamics on the femtosecond time scale and transition state formation has been demonstrated in human purine nucleoside phosphorylase (PNP) through the study of a mass-altered enzyme. Isotopic substitution in human PNP (heavy PNP) decreased the rate of on-enzyme chemistry but did not alter either the transition state structure or steady-state kinetic parameters. Here we investigate the underlying atomic motions associated with altered barrier crossing probability for heavy PNP. Transition path sampling was employed to illuminate the molecular differences between barrier crossing in light and heavy enzymes. The mass effect is apparent in promoting vibrations that polarize the N-ribosidic bond, and that promote the stability of the purine leaving group. These motions facilitate barrier crossing.

Protein Motions Are Coupled to the Reaction Chemistry in Coenzyme B12‐Dependent Ethanolamine Ammonia Lyase

The role of protein dynamics in promoting catalysis is hotly debated. Infrared data from both ultrafast flash photolysis and stopped-flow studies show that not only does there appear to be vibrational coupling between the cofactor and protein in B12-dependent ethanolamine ammonia lyase, but also that there are significant protein motions coupled to the reaction that follows substrate binding (see picture).

Changes in Protein Dynamics of the DNA Repair Dioxygenase AlkB upon Binding of Fe2+ and 2-Oxoglutarate

The Escherichia coli DNA repair enzyme AlkB is a 2-oxoglutarate (2OG)-dependent Fe(2+) binding dioxygenase that removes methyl lesions from DNA and RNA. To date, nine human AlkB homologues are known: ABH1 to ABH8 and the obesity-related FTO. Similar to AlkB, these homologues exert their activity on nucleic acids, although for some homologues the biological substrate remains to be identified. 2OG dioxygenases require binding of the cofactors Fe(2+) and 2OG in the active site to form a catalytically competent complex. We present a structural analysis of AlkB using NMR, fluorescence, and CD spectroscopy to show that AlkB is a dynamic protein exhibiting different folding states. In the absence of the cofactors Fe(2+) and 2OG, apoAlkB is a highly dynamic protein. Binding of either Fe(2+) or 2OG alone does not significantly affect the protein dynamics. Formation of a fully folded and catalytically competent holoAlkB complex only occurs when both 2OG and Fe(2+) are bound. These findings provide the first insights into protein folding of 2OG-dependent dioxygenases. A role for protein dynamics in the incorporation of the metal cofactor is discussed.

Taking Ockham's razor to enzyme dynamics and catalysis

The role of protein dynamics in enzyme catalysis is a matter of intense current debate. Enzyme-catalysed reactions that involve significant quantum tunnelling can give rise to experimental kinetic isotope effects with complex temperature dependences, and it has been suggested that standard statistical rate theories, such as transition-state theory, are inadequate for their explanation. Here we introduce aspects of transition-state theory relevant to the study of enzyme reactivity, taking cues from chemical kinetics and dynamics studies of small molecules in the gas phase and in solution--where breakdowns of statistical theories have received significant attention and their origins are relatively better understood. We discuss recent theoretical approaches to understanding enzyme activity and then show how experimental observations for a number of enzymes may be reproduced using a transition-state-theory framework with physically reasonable parameters. Essential to this simple model is the inclusion of multiple conformations with different reactivity.