Polar Versus Non-polar Local Ordering at Mobile Sites in Proteins: Slowly Relaxing Local Structure Analysis of (15)N Relaxation in the Third Immunoglobulin-Binding Domain of Streptococcal Protein G

Local Structures in Proteins from Microsecond Molecular Dynamics Simulations: 2. The Role of Symmetry in GTPase Binding and Dimer Formation

The Rho GTPase binding domain of Plexin-B1 (RBD) prevails in solution as dimer. Under appropriate circumstances, it binds the small GTPase Rac1 to yield the complex RBD-Rac1. Here, we study RBD dimerization and complex formation from a symmetry-based perspective using data derived from 1 μs long MD simulations. The quantities investigated are the local potentials, u(MD), prevailing at the N-H sites of the protein. These potentials are statistical in character providing an empirical description of the local structure. To establish more methodical description, a method for approximating them by explicit functions, u(simulated), was developed in the preceding article in this journal issue. These functions are combinations of analytical Wigner functions, DL,K, belonging to the D2h point group. The D2h subgroups Ag and B2u are found to dominate u(simulated); the B1u subgroup contributes in some cases. The Ag (B2u) functions have axial or rhombic symmetry. For the first time, local potentials in proteins can be quantitatively characterized in terms of their strength (rhombicity) evaluated by axial Ag (rhombic Ag and B2u) contributions. Until now, the chain-segment [β3-L3-β4] and to some extent the α2-helix have been associated with GTPase binding. Here, we find that this process causes an increase (decrease) in the potential strength of β3 and β4 (the preceding L2 loop and the remote chain-segment [(α2-helix)-(α2/β5-turn)-(β5-strand)]), suggesting effects of counterbalancing and allostery. There is evidence for the L2 loop being associated with RBD-GTPase binding. Until now only the L4 loop has been associated with RBD dimerization. The latter process is found to cause an increase (decrease) in the potential strength and rhombicity of the L4 loop (the adjacent chain-segment [(α2-helix)-(α2/β5-turn)-(β5-strand)]), suggesting counterbalancing activity. On average, the RBD dimer features stronger local potentials than RBD-Rac1. The novel information inherent in these findings is mesoscopic in character. Prospects of interest include exploring relation to atomistic force-field parameters.

Local Structures in Proteins from Microsecond Molecular Dynamics Simulations: A Symmetry-Based Perspective

We report on a new method for the characterization of local structures in proteins based on extensive molecular dynamics (MD) simulations, here, 1 μs in length. The N-H bond of the Rho GTPase binding domain of plexin-B1 (RBD) serves as a probe and the potential, u(MD), which restricts its internal motion, as a qualifier of the local dynamic structure. u(MD) is derived from the MD trajectory as a function of the polar angles, (θ, φ), which specify the N-H orientation in the protein. u(MD) is statistical in character yielding empirical descriptions. To establish more insightful methodical descriptions, we develop a comprehensive method which approximates u(MD) by combinations of analytical Wigner functions that belong to the D2h point group. These combinations, called u(simulated), make it possible to gain a new perspective of local dynamic structures in proteins based on explicit potentials/free energy surfaces and associated probability densities, entropy, and ordering. A simpler method was developed previously using 100 ns MD simulations. In that case, the traditional "perpendicular N-H ordering" setting centered at Cα-Cα with (θ, φ) = (90, 90) and generally, featuring positive φ, prevailed. u(MD) derived from 1 μs MD simulations is considerably more complex requiring substantial model enhancement. The enhanced method applies to the well-structured sections of the RBD. It only applies partly to its loops where u(MD) extends into the negative-φ region where we detect nonperpendicular N-H ordering. This arrangement requires devising new reference structures and making substantial algorithmic changes, to be performed in future work. Here, we focus on developing the comprehensive method and using it to investigate perpendicular ordering settings. We find that secondary structures (loops) exhibit varying (virtually invariant) potentials with Ag, B2u, and B1u (Ag and B2u) D2h symmetry. Application to RBD dimerization and RBD binding to the GTPase Rac1 is described in the subsequent article. Applications to other probes, proteins, and biological functions, based on explicit local potentials, probability densities, entropy, and ordering, are possible.

Conformational Entropy from Restricted Bond-Vector Motion in Proteins: The Symmetry of the Local Restrictions and Relation to NMR Relaxation

Locally mobile bond-vectors contribute to the conformational entropy of the protein, given by S_k ≡ S/k = -∫(P_eq ln P_eq)dΩ - ln∫dΩ. The quantity P_eq = exp(-u)/Z is the orientational probability density, where Z is the partition function and u is the spatially restricting potential exerted by the immediate internal protein surroundings at the site of the motion of the bond-vector. It is appropriate to expand the potential, u, which restricts local rotational reorientation, in the basis set of the real combinations of the Wigner rotation matrix elements, D_0K^L. For small molecules dissolved in anisotropic media, one typically keeps the lowest even L, L = 2, nonpolar potential in axial or rhombic form. For bond-vectors anchored at the protein, the lowest odd L, L = 1, polar potential is to be used in axial or rhombic form. Here, we investigate the effect of the symmetry and polarity of these potentials on S_k. For L = 1 (L = 2), S_k is the same (differs) for parallel and perpendicular ordering. The plots of S_k as a function of the coefficients of the rhombic L = 1 (L = 2) potential exhibit high-symmetry (specific low-symmetry) patterns with parameter-range-dependent sensitivity. Similar statements apply to analogous plots of the potential minima. S_k is also examined as a function of the order parameters defined in terms of u. Graphs displaying these correlations, and applications illustrating their usage, are provided. The features delineated above are generally useful for devising orienting potentials that best suit given physical circumstances. They are particularly useful for bond-vectors acting as NMR relaxation probes in proteins, when their restricted local motion is analyzed with stochastic models featuring Wigner-function-made potentials. The relaxation probes could also be molecules adsorbed at surfaces, inserted into membranes, or interlocked within metal-organic frameworks.

Local ordering and dynamics in anisotropic media by magnetic resonance: from liquid crystals to proteins

Magnetic resonance methods have been used extensively for over 50 years to elucidate molecular structure and dynamics of liquid crystals (LCs), providing information quite unique in its rigour and extent. The ESR- or NMR-active probe is often a solute molecule reporting on characteristics associated with the surrounding (LC) medium, which exerts the spatial restrictions on the probe. The theoretical approaches developed for LCs are applicable to anisotropic media in general. Of particular interest is the interior space of a globular protein labelled, e.g. with a nitroxide moiety or a ¹⁵N-¹H bond. The ESR or NMR label plays the role of the probe and the internal protein surroundings the role of the anisotropic medium. A general feature of the restricted motions is the local ordering, i.e. the nature, magnitude and symmetry of the spatial restraints exerted at the site of the moving probe. This property is the main theme of the present review article. We outline its treatment in our work from both the theoretical and the experimental points of view, highlighting the new physical insights gained. Our illustrations include studies on thermotropic (nematic and smectic) and lyotropic liquid crystals formed by phospholipids, in addition to studies of proteins.

15N-H-Related Conformational Entropy Changes Entailed By Plexin-B1 RBD Dimerization: Combined Molecular Dynamics/NMR Relaxation Approach

We report on a new method for determining function-related conformational entropy changes in proteins. Plexin-B1 RBD dimerization serves as example, and internally mobile N-H bonds serve as probes. S_k (entropy in units of k_BT) is given by -∫(P_eqlnP_eq)dΩ, where P_eq = exp(-u) is the probability density for probe orientation, and u the local potential. Previous slowly relaxing local structure (SRLS) analyses of ¹⁵N-H relaxation in proteins determined linear combinations of D₀₀²(Ω) and (D₀₂²(Ω) + D_0-2²(Ω)) (D_0K^L(Ω) represents a Wigner rotation matrix element in uniaxial local medium) as "best-fit" form of u. SRLS also determined the "best-fit" orientation of the related ordering tensor. On the basis of this information the coefficients (in the linear combination) of the terms specified above are determined with molecular dynamics (MD) simulations. With the explicit expression for u thus in hand, S_k is calculated. We find that in general S_k decreases, i.e., the local order increases, upon plexin-B1 RBD dimerization. The largest decrease in S_k occurs in the helices α₁ and α₂, followed by the α₂/β₆ turn. Only the relatively small peripheral β₂ strand, β₂/α₁ turn, and L3 loop become more disordered. That α-helices dominate ΔS_k = S_k(dimer) - S_k(monomer), a few peripheral outliers partly counterbalance the overall decrease in S_k, and the probability density function, P_eq, has rhombic symmetry given that the underlying potential function, u, has rhombic symmetry, are interesting features. We also derive S² (the proxy of u in the simple "model-free (MF)" limit of SRLS) with MD. Its conversion into a potential requires assumptions and adopting a simple axial form of u. Ensuing ΔS_k(MF) profiles are u-dependent and differ from ΔS_k(SRLS). A method that provides consistent, general, and accurate S_k, atomistic/mesoscopic in nature, has been developed. Its ability to provide new insights in protein research has been illustrated.