Internal Protein Motions Research Articles

A Highly Ordered Nitroxide Side Chain for Distance Mapping and Monitoring Slow Structural Fluctuations in Proteins.

The online version contains supplementary material available at 10.1007/s00723-023-01618-8.

Local and Large-Scale Conformational Dynamics in Unfolded Proteins and IDPs. I. Effect of Solvent Viscosity and Macromolecular Crowding.

Protein/solvent interactions largely influence protein dynamics, particularly motions in unfolded and intrinsically disordered proteins (IDPs). Here, we apply triplet-triplet energy transfer (TTET) to investigate the coupling of internal protein motions to solvent motions by determining the effect of solvent viscosity (η) and macromolecular crowding on the rate constants of loop formation (kc) in several unfolded polypeptide chains including IDPs. The results show that the viscosity dependence of loop formation depends on amino acid sequence, loop length, and co-solute size. Below a critical size (rc), co-solutes exert a maximum effect, indicating that under these conditions microviscosity experienced by chain motions matches macroviscosity of the solvent. rc depends on chain stiffness and reflects the length scale of the chain motions, i.e., it is related to the persistence length. Above rc, the effect of solvent viscosity decreases with increasing co-solute size. For co-solutes typically used to mimic cellular environments, a scaling of kc ∝ η-0.1 is observed, suggesting that dynamics in unfolded proteins are only marginally modulated in cells. The effect of solvent viscosity on kc in the small co-solute limit (below rc) increases with increasing chain length and chain flexibility. Formation of long and very flexible loops exhibits a kc ∝ η-1 viscosity dependence, indicating full solvent coupling. Shorter and less flexible loops show weaker solvent coupling with values as low as kc ∝ η-0.75±0.02. Coupling of formation of short loops to solvent motions is very little affected by amino acid sequence, but solvent coupling of long-range loop formation is decreased by side chain sterics.

Perspective: How Fast Dynamics Affect Slow Function in Protein Machines.

Internal motions in proteins take place on a broad range of time- and space-scales. The potential roles of these dynamics in the biochemical functions of proteins have intrigued biophysicists for many years, and multiple mechanisms to couple motions to function have been proposed. Some of these mechanisms have relied on equilibrium concepts. For example, the modulation of dynamics was proposed to change the entropy of a protein, hence affecting processes such as binding. This so-called dynamic allostery scenario has been demonstrated in several recent experiments. Perhaps even more intriguing may be models that involve out-of-equilibrium operation, which by necessity require the input of energy. We discuss several recent experimental studies that expose such potential mechanisms for coupling dynamics and function. In Brownian ratchets, for example, directional motion is promoted by switching a protein between two free energy surfaces. An additional example involves the effect of microsecond domain-closure dynamics of an enzyme on its much slower chemical cycle. These observations lead us to propose a novel two-time-scale paradigm for the activity of protein machines: fast equilibrium fluctuations take place on the microsecond-millisecond time scale, while on a slower time scale, free energy is invested in order to push the system out of equilibrium and drive functional transitions. Motions on the two time scales affect each other and are essential for the overall function of these machines.

Study of Protein Dynamics via Neutron Spin Echo Spectroscopy.

Most human body proteins' activity and functionality are related to configurational changes of entire subdomains within the protein crystal structure. The crystal structures build the basis for any calculation that describes the structure or dynamics of a protein, most of the time with strong geometrical restrictions. However, these restrictions from the crystal structure are not present in the solution. The structure of the proteins in the solution may differ from the crystal due to rearrangements of loops or subdomains on the pico to nanosecond time scale (i.e., the internal protein dynamics time regime). The present work describes how slow motions on timescales of several tens of nanoseconds can be accessed using neutron scattering. In particular, the dynamical characterization of two major human proteins, an intrinsically disordered protein that lacks a well-defined secondary structure and a classical antibody protein, is addressed by neutron spin echo spectroscopy (NSE) combined with a wide range of laboratory characterization methods. Further insights into protein domain dynamics were achieved using mathematical modeling to describe the experimental neutron data and determine the crossover between combined diffusive and internal protein motions. The extraction of the internal dynamic contribution to the intermediate scattering function obtained from NSE, including the timescale of the various movements, allows further vision into the mechanical properties of single proteins and the softness of proteins in their nearly natural environment in the crowded protein solution.

Time- and ensemble-average statistical mechanics of the Gaussian network model

We present analytical results for a set of time- and ensemble-averaged physical observables in the non-Hookean Gaussian network model (GNM)—a generalization of the Rouse model to elastic networks with links that display a certain degree of extensional and rotational stiffness. We focus on a set of coarse-grained observables that may be of interest in the analysis of GNM in the context of internal motions in proteins and mechanical frames in contact with a heat bath. A C++ computer code is made available that implements all analytical results.

Intermediate scattering functions of a rigid body monoclonal antibody protein in solution studied by dissipative particle dynamic simulation.

In the past decade, there was increased research interest in studying internal motions of flexible proteins in solution using Neutron Spin Echo (NSE) as NSE can simultaneously probe the dynamics at the length and time scales comparable to protein domain motions. However, the collective intermediate scattering function (ISF) measured by NSE has the contributions from translational, rotational, and internal motions, which are rather complicated to be separated. Widely used NSE theories to interpret experimental data usually assume that the translational and rotational motions of a rigid particle are decoupled and independent to each other. To evaluate the accuracy of this approximation for monoclonal antibody (mAb) proteins in solution, dissipative particle dynamic computer simulation is used here to simulate a rigid-body mAb for up to about 200 ns. The total ISF together with the ISFs due to only the translational and rotational motions as well as their corresponding effective diffusion coefficients is calculated. The aforementioned approximation introduces appreciable errors to the calculated effective diffusion coefficients and the ISFs. For the effective diffusion coefficient, the error introduced by this approximation can be as large as about 10% even though the overall agreement is considered reasonable. Thus, we need to be cautious when interpreting the data with a small signal change. In addition, the accuracy of the calculated ISFs due to the finite computer simulation time is also discussed.

Membrane Proteins Have Distinct Fast Internal Motion and Residual Conformational Entropy.

The internal motions of integral membrane proteins have largely eluded comprehensive experimental characterization. Here the fast side‐chain dynamics of the α‐helical sensory rhodopsin II and the β‐barrel outer membrane protein W have been investigated in lipid bilayers and detergent micelles by solution NMR relaxation techniques. Despite their differing topologies, both proteins have a similar distribution of methyl‐bearing side‐chain motion that is largely independent of membrane mimetic. The methyl‐bearing side chains of both proteins are, on average, more dynamic in the ps–ns timescale than any soluble protein characterized to date. Accordingly, both proteins retain an extraordinary residual conformational entropy in the folded state, which provides a counterbalance to the absence of the hydrophobic effect. Furthermore, the high conformational entropy could greatly influence the thermodynamics underlying membrane‐protein functions, including ligand binding, allostery, and signaling.

Membrane Proteins Have Distinct Fast Internal Motion and Residual Conformational Entropy

AbstractThe internal motions of integral membrane proteins have largely eluded comprehensive experimental characterization. Here the fast side‐chain dynamics of the α‐helical sensory rhodopsin II and the β‐barrel outer membrane protein W have been investigated in lipid bilayers and detergent micelles by solution NMR relaxation techniques. Despite their differing topologies, both proteins have a similar distribution of methyl‐bearing side‐chain motion that is largely independent of membrane mimetic. The methyl‐bearing side chains of both proteins are, on average, more dynamic in the ps–ns timescale than any soluble protein characterized to date. Accordingly, both proteins retain an extraordinary residual conformational entropy in the folded state, which provides a counterbalance to the absence of the hydrophobic effect. Furthermore, the high conformational entropy could greatly influence the thermodynamics underlying membrane‐protein functions, including ligand binding, allostery, and signaling.

Reducing bias in the analysis of solution-state NMR data with dynamics detectors.

Nuclear magnetic resonance (NMR) is sensitive to dynamics on a wide range of correlation times. Recently, we have shown that analysis of relaxation rates via fitting to a correlation function with a small number of exponential terms could yield a biased characterization of molecular motion in solid-state NMR due to limited sensitivity of experimental data to certain ranges of correlation times. We introduced an alternative approach based on "detectors" in solid-state NMR, for which detector responses characterize motion for a range of correlation times and reduce potential bias resulting from the use of simple models for the motional correlation functions. Here, we show that similar bias can occur in the analysis of solution-state NMR relaxation data. We have thus adapted the detector approach to solution-state NMR, specifically separating overall tumbling motion from internal motions and accounting for contributions of chemical exchange to transverse relaxation. We demonstrate that internal protein motions can be described with detectors when the overall motion and the internal motions are statistically independent. We illustrate the detector analysis on ubiquitin with typical relaxation data sets recorded at a single high magnetic field or at multiple high magnetic fields and compare with results of model-free analysis. We also compare our methodology to LeMaster's method of dynamics analysis.

Introduction to Enzymes

Microbes utilize enzymes to perform a variety of functions. Enzymes are biocatalysts working as highly efficient machines at the molecular level. In the past, enzymes have been viewed as static entities and their function has been explained on the basis of direct structural interactions between the enzyme and the substrate. A variety of experimental and computational techniques, however, continue to reveal that proteins are dynamically active machines, with various parts exhibiting internal motions at a wide range of time-scales. Increasing evidence also indicates that these internal protein motions play a role in promoting protein function such as enzyme catalysis. Moreover, the thermodynamical fluctuations of the solvent, surrounding the protein, have an impact on internal protein motions and, therefore, on enzyme function. In this review, we describe recent biochemical and theoretical investigations of internal protein dynamics linked to enzyme catalysis. In the enzyme cyclophilin A, investigations have lead to the discovery of a network of protein vibrations promoting catalysis. Cyclophilin A catalyzes peptidyl-prolyl cis/trans isomerization in a variety of peptide and protein substrates. Recent studies of cyclophilin A are discussed in detail and other enzymes (dihydrofolate reductase and liver alcohol dehydrogenase) where similar discoveries have been reported are also briefly discussed. The detailed characterization of the discovered networks indicates that protein dynamics plays a role in rate-enhancement achieved by enzymes. An integrated view of enzyme structure, dynamics and function have wide implications in understanding allosteric and co-operative effects, as well as protein engineering of more efficient enzymes and novel drug design.

Frataxin Structure and Function.

Mammalian frataxin is a small mitochondrial protein involved in iron sulfur cluster assembly. Frataxin deficiency causes the neurodegenerative disease Friedreich's Ataxia. Valuable knowledge has been gained on the structural dynamics of frataxin, metal-ion-protein interactions, as well as on the effect of mutations on protein conformation, stability and internal motions. Additionally, laborious studies concerning the enzymatic reactions involved have allowed for understanding the capability of frataxin to modulate Fe-S cluster assembly function. Remarkably, frataxin biological function depends on its interaction with some proteins to form a supercomplex, among them NFS1 desulfurase and ISCU, the scaffolding protein. By combining multiple experimental tools including high resolution techniques like NMR and X-ray, but also SAXS, crosslinking and mass-spectrometry, it was possible to build a reliable model of the structure of the desulfurase supercomplex NFS1/ACP-ISD11/ISCU/frataxin. In this chapter, we explore these issues showing how the scientific view concerning frataxin structure-function relationships has evolved over the last years.

A Biophysical Perspective on Enzyme Catalysis.

Even after a century of investigation, our understanding of how enzymes work remains far from complete. In particular, several factors that enable enzymes to achieve high catalytic efficiencies remain only poorly understood. A number of theories have been developed, which propose or reaffirm that enzymes work as structural scaffolds, serving to bring together and properly orient the participants so that the reaction can proceed; therefore, leading to enzymes being viewed as only passive participants in the catalyzed reaction. A growing body of evidence shows that enzymes are not rigid structures but are constantly undergoing a wide range of internal motions and conformational fluctuations. In this Perspective, on the basis of studies from our group, we discuss the emerging biophysical model of enzyme catalysis that provides a detailed understanding of the interconnection among internal protein motions, conformational substates, enzyme mechanisms, and the catalytic efficiency of enzymes. For a number of enzymes, networks of conserved residues that extend from the surface of the enzyme all the way to the active site have been discovered. These networks are hypothesized to serve as pathways of energy transfer that enables thermodynamical coupling of the surrounding solvent with enzyme catalysis and play a role in promoting enzyme function. Additionally, the role of enzyme structure and electrostatic effects has been well acknowledged for quite some time. Collectively, the recent knowledge gained about enzyme mechanisms suggests that the conventional paradigm of enzyme structure encoding function is incomplete and needs to be extended to structure encodes dynamics, and together these enzyme features encode function including catalytic rate acceleration.

Predicting X-ray solution scattering from flexible macromolecules.

Wide-angle X-ray solution scattering (WAXS) patterns contain substantial information about the structure and dynamics of a protein. Solution scattering from a rigid protein can be predicted from atomic coordinate sets to within experimental error. However, structural fluctuations of proteins in solution can lead to significant changes in the observed intensities. The magnitude and form of these changes contain information about the nature and spatial extent of structural fluctuations in the protein. Molecular dynamics (MD) simulations based on a crystal structure and selected force field generate models for protein internal motions, and here we demonstrate that they can be used to predict the impact of structural fluctuations on solution scattering data. In cases where the observed and calculated intensities correspond, we can conclude that the X-ray scattering provides direct experimental validation of the structural and MD results. In cases where calculated and observed intensities are at odds, the inconsistencies can be used to determine the origins of these discrepancies. They may be because of overestimates or underestimates of structural fluctuations in MD simulations, under-sampling of the structural ensemble in the simulations, errors in the structural model, or a mismatch between the experimental conditions and the parameters used in carrying out the MD simulation.

Functional Role of Solvent Entropy and Conformational Entropy of Metal Binding in a Dynamically Driven Allosteric System.

Allostery is a regulatory phenomenon whereby ligand binding to one site influences the binding of the same or a different ligand to another site on a macromolecule. The physical origins of allosteric regulation remain under intense investigation. In general terms, ligand-induced structural changes, perturbations of residue-specific dynamics, and surrounding solvent molecules all potentially contribute to the global energetics of allostery. While the role of solvent is generally well understood in regulatory events associated with major protein structural rearrangements, the degree to which protein dynamics impact solvent degrees of freedom is unclear, particularly in cases of dynamically driven allostery. With the aid of new crystal structures, extensive calorimetric and residue-specific dynamics studies over a range of time scales and temperatures, we dissect for the first time the relative degree to which changes in solvent entropy and residue-specific dynamics impact dynamically driven, allosteric inhibition of DNA binding by Zn in the zinc efflux repressor, CzrA (chromosomal zinc-regulated repressor). We show that non-native residue-specific dynamics in allosterically impaired CzrA mutants are accompanied by significant perturbations in solvent entropy that cannot be predicted from crystal structures. We conclude that functional dynamics are not necessarily restricted to protein residues but involve surface water molecules that may be responding to ligand (Zn)-mediated perturbations in protein internal motions that define the conformational ensemble, rather than major structural rearrangements.

Internal protein motions in molecular-dynamics simulations of Bragg and diffuse X-ray scattering.

Molecular-dynamics (MD) simulations of Bragg and diffuse X-ray scattering provide a means of obtaining experimentally validated models of protein conformational ensembles. This paper shows that compared with a single periodic unit-cell model, the accuracy of simulating diffuse scattering is increased when the crystal is modeled as a periodic supercell consisting of a 2×2× 2 layout of eight unit cells. The MD simulations capture the general dependence of correlations on the separation of atoms. There is substantial agreement between the simulated Bragg reflections and the crystal structure; there are local deviations, however, indicating both the limitation of using a single structure to model disordered regions of the protein and local deviations of the average structure away from the crystal structure. Although it was anticipated that a simulation of longer duration might be required to achieve maximal agreement of the diffuse scattering calculation with the data using the supercell model, only a microsecond is required, the same as for the unit cell. Rigid protein motions only account for a minority fraction of the variation in atom positions from the simulation. The results indicate that protein crystal dynamics may be dominated by internal motions rather than packing interactions, and that MD simulations can be combined with Bragg and diffuse X-ray scattering to model the protein conformational ensemble.

Slow Dynamics of Tryptophan-Water Networks in Proteins.

Water has a profound effect on the dynamics of biomolecules and governs many biological processes, leading to the concept that function is slaved to solvent dynamics within and surrounding the biomolecule. Protein conformational changes on μs-ms time scales are frequently associated with protein function, but little is known about the behavior of protein-bound water on these time scales. Here we have used NMR relaxation dispersion measurements to probe the tryptophan indoles in the enzyme dihydrofolate reductase (DHFR). We find that during structural changes on the μs-ms time scale, large chemical shift changes are often observed for the NH proton on the indole ring, while relatively smaller chemical shift changes are observed for the ring nitrogen atom. Comparison with experimental chemical shifts and density functional theory-based chemical shift predictions show that during the structural change the tryptophan indole NHs remain bound to water, but the geometry of the protein-bound water networks changes. These results establish that relaxation dispersion measurements can indirectly probe water dynamics and indicate that water can influence, or be influenced by, protein conformational changes on the μs-ms time scale. Our data show that structurally conserved bound water molecules can play a critical role in transmitting information between functionally important regions of the protein and provide evidence that internal protein motions can be coupled through the mediation of hydrogen-bonded water bound in the protein structure.

Homology Modeling and Molecular Dynamics Simulation Combined with X-ray Solution Scattering Defining Protein Structures of Thromboxane and Prostacyclin Synthases.

A combination of molecular dynamics (MD) simulations and X-ray scattering (SAXS) has emerged as the approach of choice for studying protein structures and dynamics in solution. This approach has potential applications for membrane proteins that neither are soluble nor form crystals easily. We explore the water-coupled dynamic structures of thromboxane synthase (TXAS) and prostacyclin synthase (PGIS) from scanning HPLC-SAXS measurements combined with MD ensemble analyses. Both proteins are heme-containing enzymes in the cytochrome P450 family, known as prostaglandin H2 (PGH2) isomerase, with counter-functions in regulation of platelet aggregation. Currently, the X-ray crystallographic structures of PGIS are available, but those for TXAS are not. The use of homology modeling of the TXAS structure with ns-μs explicit water solvation MD simulations allows much more accurate estimation of the configuration space with loop motion and origin of the protein behaviors in solution. In contrast to the stability of the conserved PGIS structure in solution, the pronounced TXAS flexibility has been revealed to have unstructured loop regions in connection with the characteristic P450 structural elements. The MD-derived and experimental-solution SAXS results are in excellent agreement. The significant protein internal motions, whole-molecule structures, and potential problems with protein folding, crystallization, and functionality are examined.

Haloalkane dehalogenases as a subject for crystallographic studies

An integrated view of protein structure, solvation, dynamics, and function is emerging, where proteins are considered as dynamically active machines and internal protein motions are closely linked to function such as enzyme catalysis. It is widely believed that hydration is necessary for enzyme catalytic function.

15N CEST data and traditional model-free analysis capture fast internal dynamics of DJ-1

Previous studies have shown that relaxation parameters and fast protein dynamics can be quickly elucidated from 15N-CEST experiments [1]. Longitudinal R1 and transverse R2 values were reliably derived from fitting of CEST profiles. Herein we show that 15N-CEST experiments and traditional modelfree analysis provide the internal dynamics of three states of human protein DJ-1 at physiological temperature. The chemical exchange profiles show the absence of a minor state conformation and, in conjunction with 1H-15N NOEs, show increased mobility. R1 and R2 values remained relatively unchanged at the three naturally occurring oxidation states of DJ-1, but exhibit striking NOE differences. The NOE data was, therefore, essential in determining the internal motions of the DJ-1 proteins. To the authors' knowledge, we present the first study that combines 15N CEST data with traditional model-free analyses in the study of a biological system and affirm that more ‘lean’ model-free approaches should be used cautiously.

Determination of Protein ps-ns Motions by High-Resolution Relaxometry.

Many of the functions of biomacromolecules can be rationalized by the characterization of their conformational energy landscapes: the structures of the dominant states, transitions between states and motions within states. Nuclear magnetic resonance (NMR) spectroscopy is the technique of choice to study internal motions in proteins. The determination of motions on picosecond to nanosecond timescales requires the measurement of nuclear spin relaxation rates at multiple magnetic fields. High sensitivity and resolution are obtained only at high magnetic fields, so that, until recently, site-specific relaxation rates in biomolecules were only measured over a narrow range of high magnetic fields. This limitation was particularly striking for the quantification of motions on nanosecond timescales, close to the correlation time for overall rotational diffusion. High-resolution relaxometry is an emerging technique to investigate picosecond-nanosecond motions of proteins. This approach uses a high-field NMR spectrometer equipped with a sample shuttle device, which allows for the measurement of the relaxation rate constants at low magnetic fields, while preserving the sensitivity and resolution of a high-field NMR spectrometer. The combined analysis of high-resolution relaxometry and standard high-field relaxation data provides a more accurate description of the dynamics of proteins, in particular in the nanosecond range. The purpose of this chapter is to describe how to perform high-resolution relaxometry experiments and how to analyze the rates measured with this technique.