Evidence for domain motion in proteins affecting global diffusion properties: a nuclear magnetic resonance study

HullRad: Fast Calculations of Folded and Disordered Protein and Nucleic Acid Hydrodynamic Properties

Hydrodynamic properties are useful parameters for estimating the size and shape of proteins and nucleic acids in solution. The calculation of such properties from structural models informs on the solution properties of these molecules and complements corresponding structural studies. Here we report, to our knowledge, a new method to accurately predict the hydrodynamic properties of molecular structures. This method uses a convex hull model to estimate the hydrodynamic volume of the molecule and is orders of magnitude faster than common methods. It works well for both folded proteins and ensembles of conformationally heterogeneous proteins and for nucleic acids. Because of its simplicity and speed, the method should be useful for the modification of computer-generated, intrinsically disordered protein ensembles and ensembles of flexible, but folded, molecules in which rapid calculation of experimental parameters is needed. The convex hull method is implemented in a Python script called HullRad. The use of the method is facilitated by a web server and the code is freely available for batch applications.

The underappreciated role of allostery in the cellular network

Allosteric propagation results in communication between distinct sites in the protein structure; it also encodes specific effects on cellular pathways, and in this way it shapes cellular response. One example of long-range effects is binding of morphogens to cell surface receptors, which initiates a cascade of protein interactions that leads to genome activation and specific cellular action. Allosteric propagation results from combinations of multiple factors, takes place through dynamic shifts of conformational ensembles, and affects the equilibria of macromolecular interactions. Here, we (a) emphasize the well-known yet still underappreciated role of allostery in conveying explicit signals across large multimolecular assemblies and distances to specify cellular action; (b) stress the need for quantitation of the allosteric effects; and finally, (c) propose that each specific combination of allosteric effectors along the pathway spells a distinct function. The challenges are colossal; the inspiring reward will be predicting function, misfunction, and outcomes of drug regimes.

SRLS analysis of 15N relaxation from bacteriophage T4 lysozyme: a tensorial perspective that features domain motion

Bacteriophage T4L lysozyme (T4L) comprises two domains connected by a helical linker. Several methods detected ns domain motion associated with the binding of the peptidoglycan substrate. An ESR study of nitroxide-labeled T4L, based on the slowly relaxing local structure (SRLS) approach, detected ns local motion involving the nitroxide and the helix housing it. (15)N−H spin relaxation data from T4L acquired at magnetic fields of 11.7 and 18.8 T, and 298 K, were analyzed previously with the model-free (MF) method. The results did not detect domain motion. SRLS is the generalization of MF. Here, we apply it to the same data analyzed previously with MF. The restricted local N−H motion is described in terms of tilted axial local ordering (S) and local diffusion (D(2)) tensors; dynamical coupling to the global tumbling is accounted for. We find that D(2,⊥) is 1.62 × 10(7) (1.56 × 10(7)) s(−1) for the N-terminal (C-terminal) domain. This dynamic mode represents domain motion. For the linker D(2,⊥) is the same as the rate of global tumbling, given by (1.46 ± 0.04) × 10(7) s(−1). D(2,∥) is 1.3 × 10(9), 1.8 × 10(9) and 5.3 × 10(9) s(−1) for the N-terminal domain, the C-terminal domain, and the linker, respectively. This dynamic mode represents N−H bond vector fluctuations. The principal axis of D(2) is virtually parallel to the N−H bond. The order parameter, S(0)(2), is 0.910 ± 0.046 for most N−H bonds. The principal axis of S is tilted from the C(i−1)(α) −C(i)(α) axis by −2° to 6° for the N-, and C-terminal domains, and by 2.5° for the linker. The tensorial-perspective-based and mode-coupling-based SRLS picture provides new insights into the structural dynamics of bacteriophage T4 lysozyme.

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.

Electron-nuclear interactions as probes of domain motion in proteins

Long range interactions between nuclear spins and paramagnetic ions can serve as a sensitive monitor of internal motion of various parts of proteins, including functional loops and separate domains. In the case of interdomain motion, the interactions between the ion and NMR-observable nuclei are modulated in direction and magnitude mainly by a combination of overall and interdomain motions. The effects on observable parameters such as paramagnetic relaxation enhancement (PRE) and pseudocontact shift (PCS) can, in principle, be used to characterize motion. These parameters are frequently used for the purpose of structural refinements. However, their use to probe actual domain motions is less common and is lacking a proper theoretical treatment from a motional perspective. In this work, a suitable spin Hamiltonian is incorporated in a two body diffusion model to produce the time correlation function for the nuclear spin-paramagnetic ion interactions. Simulated observables for nuclei in different positions with respect to the paramagnetic ion are produced. Based on these simulations, it demonstrated that both the PRE and the PCS can be very sensitive probes of domain motion. Results for different nuclei within the protein sense different aspects of the motions. Some are more sensitive to the amplitude of the internal motion, others are more sensitive to overall diffusion rates, allowing separation of these contributions. Experimentally, the interaction strength can also be tuned by substitution of different paramagnetic ions or by varying magnetic field strength (in the case of lanthanides) to allow the use of more detailed diffusion models without reducing the reliability of data fitting.

The physical basis of model-free analysis of NMR relaxation data from proteins and complex fluids

NMR relaxation experiments have provided a wealth of information about molecular motions in macromolecules and ordered fluids. Even though a rigorous theory of spin relaxation is available, the complexity of the investigated systems often makes the interpretation of limited datasets challenging and ambiguous. To allow physically meaningful information to be extracted from the data without commitment to detailed dynamical models, several versions of a model-free (MF) approach to data analysis have been developed. During the past 2 decades, the MF approach has been used in the vast majority of all NMR relaxation studies of internal motions in proteins and other macromolecules, and it has also played an important role in studies of colloidal systems. Although the MF approach has been almost universally adopted, substantial disagreement remains about its physical foundations and range of validity. It is our aim here to clarify these issues. To this end, we first present rigorous derivations of the three well-known MF formulas for the time correlation function relevant for isotropic solutions. These derivations are more general than the original ones, thereby substantially extending the range of validity of the MF approach. We point out several common misconceptions and explain the physical significance of the approximations involved. In particular, we discuss symmetry requirements and the dynamical decoupling approximation that plays a key role in the MF approach. We also derive a new MF formula, applicable to anisotropic fluids and solids, including microcrystalline protein samples. The so-called slowly relaxing local structure (SRLS) model has been advanced as an alternative to the MF approach that does not require dynamical decoupling of internal and global motions. To resolve the existing controversy about the relative merits of the SRLS model and the MF approach, we formulate and solve a planar version of the SRLS model. The analytical solution of this model reveals the unphysical consequences of the symmetrical two-body Smoluchowski equation as applied to protein dynamics, thus refuting the widely held belief that the SRLS model is more accurate than the MF approach. The different results obtained by analyzing data with these two approaches therefore do not indicate the importance of dynamical coupling between internal and global motions. Finally, we explore the two principal mechanisms of dynamical coupling in proteins: torque-mediated and friction-mediated coupling. We argue by way of specific analytically solvable models that torque-mediated coupling (which the SRLS model attempts to capture) is unimportant because the relatively slow internal motions that might couple to the global motion tend to be intermittent (jumplike) in character, whereas friction-mediated coupling (which neither the SRLS model nor the MF approach incorporates) may be important for proteins with unstructured parts or flexibly connected domains.