Functional roles of enzyme dynamics in accelerating active site chemistry: Emerging techniques and changing concepts

Vibrational Energy Landscapes and Energy Flow in GPCRs

We construct and analyze disconnectivity graphs to provide the first graphical representation of the vibrational energy landscape of a protein, in this study β2AR, a G-protein coupled receptor (GPCR), in active and inactive states. The graphs, which indicate the relative free energy of each residue and the minimum free energy barriers for energy transfer between them, reveal important composition, structural and dynamic properties that mediate the flow of energy. Prolines and glycines, which contribute to GPCR plasticity and function, are identified as bottlenecks to energy transport along the backbone from which alternative pathways for energy transport via nearby noncovalent contacts emerge, seen also in the analysis of first passage time (FPT) distributions presented here. Striking differences between the disconnectivity graphs and FPT distributions for the inactive and active states of β2AR are found where structural and dynamic changes occur upon activation, contributing to allosteric regulation.

Experimental Evidence for the Role of Dynamics in pH-Dependent Enzymatic Activity

Enzymatic activity is heavily influenced by pH, but the rationale for the dynamical mechanism of pH-dependent enzymatic activity has not been fully understood. In this work, combined neutron scattering techniques, including quasielastic neutron scattering (QENS) and small angle neutron scattering (SANS), are used to study the structural and dynamic changes of a model enzyme, xylanase, under different pH and temperature environments. The QENS results reveal that xylanase at optimal pH exhibits faster relaxational dynamics and a lower energy barrier between conformational substates. The SANS results demonstrate that pH affects both xylanase's stability and monodispersity. Our findings indicate that enzymes have optimized stability and function under their optimal pH conditions, with both structure and dynamics being affected. The current study offers valuable insights into enzymatic functionality mechanisms, allowing for broad industrial applications.

Cryo-EM and Molecular Dynamics Simulations Reveal Hidden Conformational Dynamics Controlling Ammonia Transport in Human Asparagine Synthetase

How motions in enzymes might be linked to catalytic function is of considerable general interest. Recent advances in X-ray crystallography and cryogenic electron microscopy offer the promise of elucidating functionally relevant motions in proteins that are not easily amenable to study by other biophysical methods. Here we use 3D variability analysis (3DVA) on cryo-EM maps for wild type (WT) human asparagine synthetase (ASNS) and the R142I ASNS variant to identify conformational changes in the Arg-142 side chain, which mediates the formation of a catalytically relevant intramolecular tunnel. Our 3DVA results for WT ASNS are consistent with independent molecular dynamics (MD) simulations on a model generated from the X-ray structure of human ASNS. Moreover, MD simulations of computational models for the ASNS/β-aspartyl-AMP/MgPPi and R142I/β-aspartyl-AMP/MgPPi ternary complexes, suggest that the structural integrity of the tunnel is impaired in the R142I variant when β-aspartyl-AMP is present in the synthetase active site. The kinetic properties of the R142I ASNS variant support the proposed function of Arg-142. These studies illustrate the power of cryo-EM to identify localized motions and dissect the conformational landscape of large proteins. When combined with MD simulations, 3DVA is a powerful approach to understanding how conformational dynamics might regulate function in multi-domain enzymes possessing multiple active sites.

Anisotropic dynamics of an interfacial enzyme active site observed using tethered substrate analogs and ultrafast 2D IR spectroscopy

Enzymes accelerate the rates of biomolecular reactions by many orders of magnitude compared to bulk solution, and it is widely understood that this catalytic effect arises from a combination of polar pre-organization and electrostatic transition state stabilization. A number of recent reports have also implicated ultrafast (femtosecond-picosecond) timescale motions in enzymatic activity. However, complications arising from spatially-distributed disorder, the occurrence of multiple substrate binding modes, and the influence of hydration dynamics on solvent-exposed active sites still confound many experimental studies. Here we use ultrafast two-dimensional infrared (2D IR) spectroscopy and covalently-tethered substrate analogs to examine dynamical properties of the promiscuous Pyrococcus horikoshii ene-reductase (PhENR) active site in two binding configurations mimicking proposed "inactive" and "reactive" Michaelis complexes. Spectral diffusion measurements of aryl-nitrile substrate analogs reveal an end-to-end tradeoff between fast (sub-ps) and slow (>5 ps) motions. Fermi resonant aryl-azide analogs that sense interactions of coupled oscillators are described. Lineshape and quantum beat analyses of these probes reveal characteristics that correlate with aryl-nitrile frequency fluctuation correlation functions parameters, demonstrating that this anisotropy is an intrinsic property of the water-exposed active site, where countervailing gradients of fast dynamics and disorder in the reactant ground state are maintained near the hydration interface. Our results suggest several plausible factors leading to state-selective rate enhancement and promiscuity in PhENR. This study also highlights a strategy to detect perturbations to vibrational modes outside the transparent window of the mid-IR spectrum, which may be extended to other macromolecular systems.

Effects of Hydrogen/Deuterium Exchange on Protein Stability in Solution and in the Gas Phase

Mass spectrometry (MS)-based techniques are widely used for probing protein structure and dynamics in solution. H/D exchange (HDX)-MS is one of the most common approaches in this context. HDX is often considered to be a "benign" labeling method, in that it does not perturb protein behavior in solution. However, several studies have reported that D₂O pushes unfolding equilibria toward the native state. The origin, and even the existence of this protein stabilization remain controversial. Here we conducted thermal unfolding assays in solution to confirm that deuterated proteins in D₂O are more stable, with 2-4 K higher melting temperatures than unlabeled proteins in H₂O. Previous studies tentatively attributed this phenomenon to strengthened H-bonds after deuteration, an effect that may arise from the lower zero-point vibrational energy of the deuterated species. Specifically, it was proposed that strengthened water-water bonds (W···W) in D₂O lower the solubility of nonpolar side chains. The current work takes a broader view by noting that protein stability in solution also depends on water-protein (W···P) and protein-protein (P···P) H-bonds. To help unravel these contributions, we performed collision-induced unfolding (CIU) experiments on gaseous proteins generated by native electrospray ionization. CIU profiles of deuterated and unlabeled proteins were indistinguishable, implying that P···P contacts are insensitive to deuteration. Thus, protein stabilization in D₂O is attributable to solvent effects, rather than alterations of intraprotein H-bonds. Strengthening of W···W contacts represents one possible explanation, but the stabilizing effect of D₂O can also originate from weakened W···P bonds. Future work will be required to elucidate which of these two scenarios is correct, or if both contribute to protein stabilization in D₂O. In any case, the often-repeated adage that "D-bonds are more stable than H-bonds" does not apply to intramolecular contacts in native proteins.

Systematic mapping of the conformational landscape and dynamism of soluble fibrinogen

BACKGROUND

Fibrinogen is a soluble, multisubunit, and multidomain dimeric protein, which, upon its proteolytic cleavage by thrombin, is converted to insoluble fibrin, initiating polymerization that substantially contributes to clot growth. Fibrinogen contains numerous, transiently accessible "cryptic" epitopes for hemostatic and immunologic proteins, suggesting that fibrinogen exhibits conformational flexibility, which may play functional roles in its temporal and spatial interactions. Hitherto, there have been limited integrative approaches characterizing the solution structure and internal flexibility of fibrinogen.

METHODS

Here, utilizing a multipronged, biophysical approach involving 2 solution-based techniques, temperature-dependent hydrogen-deuterium exchange mass spectrometry and small angle X-ray scattering, corroborated by negative stain electron microscopy, we present a holistic, conformationally dynamic model of human fibrinogen in solution.

RESULTS

Our data reveal 4 major and distinct conformations of fibrinogen accommodated by a high degree of internal protein flexibility along its central scaffold. We propose that the fibrinogen structure in the solution consists of a complex, conformational landscape with multiple local minima. This is further supported by the location of numerous point mutations that are linked to dysfibrinogenemia and posttranslational modifications, residing near the identified fibrinogen flexions.

CONCLUSION

This work provides a molecular basis for the structural "dynamism" of fibrinogen that is expected to influence the broad swath of its functionally diverse macromolecular interactions and fine-tune the structural and mechanical properties of blood clots.

Enzyme dynamics-a brief review

Debate has been simmering for some years regarding the importance of internal thermal motions of enzymes to catalysis. Recent developments in protein design may bring resolution of the more contentious points a little closer.

Atomistic description of the relationship between protein dynamics and catalysis with transition path sampling

Despite initial resistance, it has been increasingly accepted that protein dynamics plays a role in enzymatic catalysis. There have been two lines of research. Some works study slow conformational motions that are not coupled to the reaction coordinate, but guide the system towards catalytically competent conformations. Understanding at the atomistic level how this is accomplished has remained elusive except for a few systems. In this review we focus on fast sub-picosecond motions that are coupled to the reaction coordinate. The use of Transition Path Sampling has allowed us an atomistic description of how these rate-promoting vibrational motions are incorporated in the reaction mechanism. We will also show how we used insights from rate-promoting motions in protein design.

Engineering the Active Site Lid Dynamics to Improve the Catalytic Efficiency of Yeast Cytosine Deaminase

Conformational dynamics is important for enzyme catalysis. However, engineering dynamics to achieve a higher catalytic efficiency is still challenging. In this work, we develop a new strategy to improve the activity of yeast cytosine deaminase (yCD) by engineering its conformational dynamics. Specifically, we increase the dynamics of the yCD C-terminal helix, an active site lid that controls the product release. The C-terminal is extended by a dynamical single α-helix (SAH), which improves the product release rate by up to ~8-fold, and the overall catalytic rate kcat by up to ~2-fold. It is also shown that the kcat increase is due to the favorable activation entropy change. The NMR H/D exchange data indicate that the conformational dynamics of the transition state analog complex increases as the helix is extended, elucidating the origin of the enhanced catalytic entropy. This study highlights a novel dynamics engineering strategy that can accelerate the overall catalysis through the entropy-driven mechanism.

Nonreciprocal forces enable cold-to-hot heat transfer between nanoparticles

We study the heat transfer between two nanoparticles held at different temperatures that interact through nonreciprocal forces, by combining molecular dynamics simulations with stochastic thermodynamics. Our simulations reveal that it is possible to construct nano refrigerators that generate a net heat transfer from a cold to a hot reservoir at the expense of power exerted by the nonreciprocal forces. Applying concepts from stochastic thermodynamics to a minimal underdamped Langevin model, we derive exact analytical expressions predictions for the fluctuations of work, heat, and efficiency, which reproduce thermodynamic quantities extracted from the molecular dynamics simulations. The theory only involves a single unknown parameter, namely an effective friction coefficient, which we estimate fitting the results of the molecular dynamics simulation to our theoretical predictions. Using this framework, we also establish design principles which identify the minimal amount of entropy production that is needed to achieve a certain amount of uncertainty in the power fluctuations of our nano refrigerator. Taken together, our results shed light on how the direction and fluctuations of heat flows in natural and artificial nano machines can be accurately quantified and controlled by using nonreciprocal forces.

Dynamical activation of function in metalloenzymes

Formulations of hydrogen tunneling in enzyme-catalysed C-H activation reactions indicate enthalpic barriers to reaction that are independent of chemical steps and dependent on the protein scaffold. A tool to identify catalytically relevant site-specific protein thermal networks has emerged from temperature-dependent hydrogen deuterium exchange (TDHDX). Focusing on mutant enzyme forms with altered activation energies for catalysis, TDHDX provides a comparative analysis of the impact of mutation on Ea for local protein unfolding. Identified thermal networks appear unrelated to protein scaffold conservation and track to the dictates of the catalysed reaction, including sites for metal binding. The positions of thermal networks provide a framework for further understanding of time-dependent, functionally relevant protein motions. Measurement of nanosecond Stokes shifts at the surface of the thermal network in soybean lipoxygenase yields activation energies that are identical to Ea values measured for kcat . This finding identifies a rapid (> nanosecond), long-range and cooperative structural reorganization as the thermal barrier to catalysis. A model for protein dynamics is put forward that integrates broadly distributed protein conformational sampling with protein embedded thermal networks.