Response of small-scale, methyl rotors to protein-ligand association: a simulation analysis of calmodulin-peptide binding

Halogen bonding as a supramolecular dynamics catalyst

Dynamic processes have many implications in functional molecules, including catalysts, enzymes, host-guest complexes, and molecular machines. Here, we demonstrate via deuterium NMR relaxation experiments how halogen bonding directly impacts the dynamics in solid 2,3,5,6-tetramethylpyrazine cocrystals, catalyzing the methyl group rotation. On average, we observe a reduction of 56% in the rotational activation energy of the methyl groups in the halogen bonded cocrystals, contrasting the reduction of 36% in the hydrogen bonded cocrystals, with respect to pure 2,3,5,6-tetramethylpyrazine. Density functional theory calculations attribute this superior catalytic ability of the halogen bond to the simultaneous destabilization of the staggered conformation and stabilization of the gauche conformation, overall reducing the rotational energy barrier. Furthermore, the calculations suggest that the catalytic ability of the halogen bond may be tuneable, with stronger halogen bond donors acting as superior dynamics catalysts. Thus, halogen bonding may play a role in both assembly and promoting dynamical processes.

The halogen bond is well known for its ability to assemble supramolecules. Here, using NMR experiments, the authors reveal the role of these bonds in dynamic processes, finding that the halogen bond directly catalyzes dynamical rotation in solid cocrystals by reducing the associated energy barrier.

Direct Determination of Site-Specific Noncovalent Interaction Strengths of Proteins from NMR-Derived Fast Side Chain Motional Parameters

A novel approach to accurately determine residue-specific noncovalent interaction strengths (ξ) of proteins from NMR-measured fast side chain motional parameters (O_axis²) is presented. By probing the environmental sensitivity of side chain conformational energy surfaces of individual residues of a diverse set of proteins, the microscopic connections between ξ, O_axis², conformational entropy (S_conf), conformational barriers, and rotamer stabilities established here are found to be universal among proteins. The results reveal that side chain flexibility and conformational entropy of each residue decrease with increasing ξ and that for each residue type there exists a critical range of ξ, determined primarily by the mean side chain conformational barriers, within which flexibility of any residue can be reversibly tuned from highly flexible (with O_axis² ∼ 0) to highly restricted (with O_axis² ∼ 1) by increasing ξ by ∼3 kcal/mol. Beyond this critical range of ξ, both side chain flexibility and conformational entropy are insensitive to ξ. The interrelationships between conformational dynamics, conformational entropy, and noncovalent interactions of protein side chains established here open up new avenues to probe perturbation-induced (for example, ligand-binding, temperature, pressure) changes in fast side chain dynamics and thermodynamics of proteins by comparing their conformational energy surfaces in the native and perturbed states.

Entropic Enhancement of Protein-DNA Affinity by Oxygen-to-Sulfur Substitution in DNA Phosphate

Dithioation of DNA phosphate is known to enhance binding affinities, at least for some proteins. We mechanistically characterized this phenomenon for the Antennapedia homeodomain-DNA complex by integrated use of fluorescence, isothermal titration calorimetry, NMR spectroscopy, and x-ray crystallography. By fluorescence and isothermal titration calorimetry, we found that this affinity enhancement is entropy driven. By NMR, we investigated the ionic hydrogen bonds and internal motions of lysine side-chain NH3(+) groups involved in ion pairs with DNA. By x-ray crystallography, we compared the structures of the complexes with and without dithioation of the phosphate. Our NMR and x-ray data show that the lysine side chain in contact with the DNA phosphate becomes more dynamic upon dithioation. Our thermodynamic, structural, and dynamic investigations collectively suggest that the affinity enhancement by the oxygen-to-sulfur substitution in DNA phosphate is largely due to an entropic gain arising from mobilization of the intermolecular ion pair at the protein-DNA interface.

Distinguishing unfolding and functional conformational transitions of calmodulin using ultraviolet resonance Raman spectroscopy

Calmodulin (CaM) is a ubiquitous moderator protein for calcium signaling in all eukaryotic cells. This small calcium-binding protein exhibits a broad range of structural transitions, including domain opening and folding-unfolding, that allow it to recognize a wide variety of binding partners in vivo. While the static structures of CaM associated with its various binding activities are fairly well-known, it has been challenging to examine the dynamics of transition between these structures in real-time, due to a lack of suitable spectroscopic probes of CaM structure. In this article, we examine the potential of ultraviolet resonance Raman (UVRR) spectroscopy for clarifying the nature of structural transitions in CaM. We find that the UVRR spectral change (with 229 nm excitation) due to thermal unfolding of CaM is qualitatively different from that associated with opening of the C-terminal domain in response to Ca(2+) binding. This spectral difference is entirely due to differences in tertiary contacts at the interdomain tyrosine residue Tyr138, toward which other spectroscopic methods are not sensitive. We conclude that UVRR is ideally suited to identifying the different types of structural transitions in CaM and other proteins with conformation-sensitive tyrosine residues, opening a path to time-resolved studies of CaM dynamics using Raman spectroscopy.

Three entropic classes of side chain in a globular protein

The relationship between the NMR methyl group axial order parameter and the side chain conformational entropy is investigated in inhibitor-bound and apo human HIV protease using molecular dynamics simulation. Three distinct entropic classes of methyl-bearing side chains, determined by the topological distance of the methyl group from the protein backbone (i.e., the number of χ-bonds between the Cα and the carbon of the CH3 group), are revealed by atomistic trajectory analyses performed in the local frame of reference of individual methyl probes. The results demonstrate that topologically equivalent methyl groups experience similar nonbonded microenvironments regardless of the type of residues to which they are attached. Similarly, methyl groups that belong to the same side chain but that are not topologically equivalent exhibit different thermodynamic and dynamic properties. The two-parameter classification (based upon entropy and methyl axial order parameter) of side chains described here permits improved estimates of the conformational entropies of proteins from NMR motional parameters.

Direct observation of the ion-pair dynamics at a protein-DNA interface by NMR spectroscopy

Ion pairing is one of the most fundamental chemical interactions and is essential for molecular recognition by biological macromolecules. From an experimental standpoint, very little is known to date about ion-pair dynamics in biological macromolecular systems. Absorption, infrared, and Raman spectroscopic methods were previously used to characterize dynamic properties of ion pairs, but these methods can be applied only to small compounds. Here, using NMR (15)N relaxation and hydrogen-bond scalar (15)N-(31)P J-couplings ((h3)J(NP)), we have investigated the dynamics of the ion pairs between lysine side-chain NH3(+) amino groups and DNA phosphate groups at the molecular interface of the HoxD9 homeodomain-DNA complex. We have determined the order parameters and the correlation times for C-N bond rotation and reorientation of the lysine NH3(+) groups. Our data indicate that the NH3(+) groups in the intermolecular ion pairs are highly dynamic at the protein-DNA interface, which should lower the entropic costs for protein-DNA association. Judging from the C-N bond-rotation correlation times along with experimental and quantum-chemically derived (h3)J(NP) hydrogen-bond scalar couplings, it seems that breakage of hydrogen bonds in the ion pairs occurs on a sub-nanosecond time scale. Interestingly, the oxygen-to-sulfur substitution in a DNA phosphate group was found to enhance the mobility of the NH3(+) group in the intermolecular ion pair. This can partially account for the affinity enhancement of the protein-DNA association by the oxygen-to-sulfur substitution, which is a previously observed but poorly understood phenomenon.

Reconstruction of protein side-chain conformational free energy surfaces from NMR-derived methyl axis order parameters

An analytical approach is developed for reconstructing site-specific methyl-bearing protein side-chain conformational energy surfaces from NMR methyl axis order parameters (O(axis)(2)). Application of an enhanced sampling algorithm (adaptive biasing force) to molecular dynamics simulation of a protein, calcium-bound calmodulin, reveals a nonlinear correlation between O(axis)(2) and the populations of rotamer states of protein side-chains, permitting the rotamer populations to be extracted directly from O(axis)(2). The analytical approach yields side-chain conformational distributions that are in excellent agreement with those obtained from the enhanced-sampling MD results.

Estimation of conformational entropy in protein-ligand interactions: a computational perspective

Conformational entropy is an important component of the change in free energy upon binding of a ligand to its target protein. As a consequence, development of computational techniques for reliable estimation of conformational entropies is currently receiving an increased level of attention in the context of computational drug design. Here, we review the most commonly used techniques for conformational entropy estimation from classical molecular dynamics simulations. Although by-and-large still not directly used in practical drug design, these techniques provide a golden standard for developing other, computationally less-demanding methods for such applications, in addition to furthering our understanding of protein-ligand interactions in general. In particular, we focus on the quasi-harmonic approximation and discuss different approaches that can be used to go beyond it, most notably, when it comes to treating anharmonic and/or correlated motions. In addition to reviewing basic theoretical formalisms, we provide a concrete set of steps required to successfully calculate conformational entropy from molecular dynamics simulations, as well as discuss a number of practical issues that may arise in such calculations.

Three-dimensional mapping of microenvironmental control of methyl rotational barriers

Sterical (van der Waals-induced) rotational barriers of methyl groups are investigated theoretically, using ab initio and empirical force field calculations, for various three-dimensional microenvironmental conditions around the methyl group rotator of a model neopentane molecule. The destabilization (reducing methyl rotational barriers) or stabilization (increasing methyl rotational barriers) of the staggered conformation of the methyl rotator depends on a combination of microenvironmental contributions from (i) the number of atoms around the rotator, (ii) the distance between the rotator and the microenvironmental atoms, and (iii) the dihedral angle between the stator, rotator, and molecular environment around the rotator. These geometrical criteria combine their respective effects in a linearly additive fashion, with no apparent cooperative effects, and their combination in space around a rotator may increase, decrease, or leave the rotator's rotational barrier unmodified. This is exemplified in a geometrical analysis of the alanine dipeptide crystal where microenvironmental effects on methyl rotators' barrier of rotation fit the geometrical mapping described in the neopentane model.

Nuclear magnetic resonance study of the role of M42 in the solution dynamics of Escherichia coli dihydrofolate reductase

It is widely recognized that key positions throughout a protein's structure contribute unequally to function. In light of recent studies that suggest protein dynamics are required for function, a number of these residues may serve to promote motions required for ligand binding and catalysis. In this nuclear magnetic resonance (NMR) study, the conformational dynamics of the dihydrofolate reductase (DHFR) mutant M42W, in the presence of methotrexate and NADPH, are characterized and compared to those of the wild-type enzyme. M42 is distal to the active site, yet the M42W substitution regulates catalysis and ligand affinity and is therefore analogous to an allosteric modulator of DHFR function. To gain understanding of how this mutation regulates activity, we employ a "pandynamic" strategy by measuring conformational fluctuations of backbone amide and side-chain methyl groups on multiple time scales. Changes in pico- to nanosecond dynamics indicate that the mutational effects are propagated throughout a network of interacting residues within DHFR, consistent with a role for M42 as a dynamic communication hub. On the micro- to millisecond time scale, mutation increases the rate of switching in the catalytic core. Mutation also introduces switching in the adenosine binding subdomain that occurs at a higher frequency than in the catalytic core and which correlates with the rate of product release for M42W-DHFR. Finally, a structurally inferred analysis of side-chain dynamics suggests that the M42W mutation dampens motional contributions from nonlocal sources. These data show that the M42W mutation alters the dynamics of DHFR and are consistent with theoretical analysis that suggests this mutation disrupts motion that promotes catalysis.