Trimethylamine N-oxide stabilizes proteins via a distinct mechanism compared with betaine and glycine.

We report experimental and computational studies investigating the effects of three osmolytes, trimethylamine N-oxide (TMAO), betaine, and glycine, on the hydrophobic collapse of an elastin-like polypeptide (ELP). All three osmolytes stabilize collapsed conformations of the ELP and reduce the lower critical solution temperature (LSCT) linearly with osmolyte concentration. As expected from conventional preferential solvation arguments, betaine and glycine both increase the surface tension at the air-water interface. TMAO, however, reduces the surface tension. Atomically detailed molecular dynamics (MD) simulations suggest that TMAO also slightly accumulates at the polymer-water interface, whereas glycine and betaine are strongly depleted. To investigate alternative mechanisms for osmolyte effects, we performed FTIR experiments that characterized the impact of each cosolvent on the bulk water structure. These experiments showed that TMAO red-shifts the OH stretch of the IR spectrum via a mechanism that was very sensitive to the protonation state of the NO moiety. Glycine also caused a red shift in the OH stretch region, whereas betaine minimally impacted this region. Thus, the effects of osmolytes on the OH spectrum appear uncorrelated with their effects upon hydrophobic collapse. Similarly, MD simulations suggested that TMAO disrupts the water structure to the least extent, whereas glycine exerts the greatest influence on the water structure. These results suggest that TMAO stabilizes collapsed conformations via a mechanism that is distinct from glycine and betaine. In particular, we propose that TMAO stabilizes proteins by acting as a surfactant for the heterogeneous surfaces of folded proteins.

Overdamped dynamics of folded protein domains within a locally harmonic basin using coarse graining based on a partition of compact flexible clusters.

A coarse-graining method based on the partitioning of atoms into compact flexible clusters is used to formulate the dynamics of the nonequilibrium response of a protein to ligand dissociation. The α-carbon positions are used as the degrees of freedom. The net stiffness between each pair of neighboring α-carbons is calculated for the quasi-static, overdamped regime within the harmonic (quadratic potential energy surface) using the equivalent stiffness matrix of the network of atoms occupying the intervening space within the locally interacting region. This localized approach realizes a divide and conquer strategy that results in a substantial reduction in computational complexity while accurately predicting relaxations under general loading conditions. A close correlation between the shapes and time scales of the relaxation curves of the coarse-grained and all-atom instances of two medium-sized proteins, T4 lysozyme and ferric binding protein (each of which having known apo and holo structures), was observed for the holo to the apo transitions. Furthermore, for both proteins the dominant modes of motion and the decay rates of the temporal relaxation profiles monitoring the separation distance between select amino acid pairs were found to be nearly identical when calculated on the coarse-grained and all-atom scales.

Response of rotation-translation blocked proteins using Langevin dynamics on a locally harmonic landscape.

Langevin dynamics is used to compute the time evolution of the nonequilibrium motion of the atomic coordinates of a protein in response to ligand dissociation. The protein potential energy surface (PES) is approximated by a harmonic basin about the minimum of the unliganded state. Upon ligand dissociation, the protein undergoes relaxation from the bound to the unbound state. A coarse graining scheme based on rotation translation blocks (RTB) is applied to the relaxation of the two domain iron transport protein, ferric binding protein. This scheme provides a natural and efficient way to freeze out the small amplitude, high frequency motions within each rigid fragment, thereby allowing for the number of dynamical degrees of freedom to be reduced. The results obtained from all flexible atom (constraint free) dynamics are compared to those obtained using RTB-Langevin dynamics. To assess the impact of the assumed rigid fragment clustering on the temporal relaxation dynamics of the protein molecule, three distinct rigid block decompositions were generated and their responses compared. Each of the decompositions was a variant of the one-block-per-residue grouping, with their force and friction matrices being derived from their fully flexible counterpart. Monitoring the time evolution of the distance separating a selected pair of amino acids, the response curves of the blocked decompositions were similar in shape to each other and to the control system in which all atomic degrees of freedom are fully independent. The similar shape of the blocked responses showed that the variations in grouping had only a minor impact on the kinematics. Compared with the all atom responses, however, the blocked responses were faster as a result of the instantaneous transmission of force throughout each rigid block. This occurred because rigid blocking does not permit any intrablock deformation that could store or divert energy. It was found, however, that this accelerated response could be successfully corrected by scaling each eigenvalue in the appropriate propagation matrix by the least-squares fitted slope of the blocked vs nonblocked eigenvalue spectra. The RTB responses for each test system were dominated by small eigenvalue overdamped Langevin modes. The large eigenvalue members of each response dissipated within the first 5 ps, after which the long time response was dominated by a modest set of low energy, overdamped normal modes, that were characterized by highly cooperative, functionally relevant displacements. The response assuming that the system is in the overdamped limit was compared to the full phase space Langevin dynamics results. The responses after the first 5 ps were nearly identical, confirming that the inertial components were significant only in the initial stages of the relaxation. Since the propagator matrix in the overdamped formulation is real-symmetric and does not require the inertial component in the propagator, the computation time and memory footprint was reduced by 1 order of magnitude.