Exploring the role of structure and dynamics in the function of chymotrypsin inhibitor 2

Double Mutant of Chymotrypsin Inhibitor 2 Stabilized through Increased Conformational Entropy

The conformational heterogeneity of a folded protein can affect not only its function but also stability and folding. We recently discovered and characterized a stabilized double mutant (L49I/I57V) of the protein CI2 and showed that state-of-the-art prediction methods could not predict the increased stability relative to the wild-type protein. Here, we have examined whether changed native-state dynamics, and resulting entropy changes, can explain the stability changes in the double mutant protein, as well as the two single mutant forms. We have combined NMR relaxation measurements of the ps-ns dynamics of amide groups in the backbone and the methyl groups in the side chains with molecular dynamics simulations to quantify the native-state dynamics. The NMR experiments reveal that the mutations have different effects on the conformational flexibility of CI2: a reduction in conformational dynamics (and entropy estimated from this) of the native state of the L49I variant correlates with its decreased stability, while increased dynamics of the I57V and L49I/I57V variants correlates with their increased stability. These findings suggest that explicitly accounting for changes in native-state entropy might be needed to improve the predictions of the effect of mutations on protein stability.

Molecular insight into chymotrypsin inhibitor 2 resisting proteolytic degradation

Chymotrypsin inhibitor 2 (CI2) is a special serine protease inhibitor which can resist hydrolysis for several days with a rapid equilibrium between the Michaelis complex and acyl-enzyme intermediate. The energies and conformational changes for subtilisin-catalyzed proteolysis of CI2 were examined in this paper for the first time by employing pseudo bond ab initio QM/MM MD simulations. In the acylation reaction, a low-barrier hydrogen bond between His64 and Asp32 in the transition state together with the lack of covalent backbone constraints makes the peptide bonds of CI2 break more easily than in other serine protease inhibitors. After acyl-enzyme formation, molecular dynamics simulations showed that the access of hydrolytic water to the active site requires partial dissociation of the leaving group. However, retention of the leaving group mainly by the intra- and inter-molecular H-bonding networks hinders the access of water and retards the deacylation reaction. Instead of the dissociation constant of inhibitors, we suggest employing the free energy at the acyl-enzyme state to predict the relative hydrolysis rates of CI2 mutants, which are testified by the experimental relative hydrolysis rates.

Widespread Perturbation of Function, Structure, and Dynamics by a Conservative Single-Atom Substitution in Thymidylate Synthase

Thymidylate synthase (TSase) is responsible for synthesizing the sole de novo source of dTMP in all organisms. TSase is a drug target, and as such, it has been well studied in terms of both structure and reaction mechanism. Cysteine 146 in Escherichia coli TSase is universally conserved because it serves as the nucleophile in the enzyme mechanism. Here we use the C146S mutation to probe the role of the sulfur atom in early events in the catalytic cycle beyond serving as the nucleophile. Surprisingly, the single-atom substitution severely decreases substrate binding affinity, and the unfavorable ΔΔG°_bind is comprised of roughly equal enthalpic and entropic components at 25 °C. Chemical shifts in the free and dUMP-bound states show the mutation causes perturbations throughout TSase, including regions important for complex stability, in agreement with a less favorable enthalpy change. We measured the nuclear magnetic resonance methyl symmetry axis order parameter (S²_axis), a proxy for conformational entropy, for TSase at all vertices of the dUMP binding/C146S mutation thermodynamic cycle and found that the calculated TΔΔS°_conf is similar in sign and magnitude to the calorimetric TΔΔS°. Further, we ascribed minor resonances in wild-type-dUMP spectra to a state with a covalent bond between Sγ of C146 and C6 of dUMP and find S²_axis values are unaffected by covalent bond formation, indicating this reaction step is neutral with respect to ΔS°_conf. Lastly, the C146S mutation allowed us to measure cofactor analog binding by isothermal titration calorimetry without the confounding heat signature of covalent bond formation. Raltitrexed binds free and singly bound TSase with similar affinities, yet the two binding events have different enthalpy changes, providing further evidence of communication between the two active sites.

Advances in NMR Methods To Map Allosteric Sites: From Models to Translation

The last five years have witnessed major developments in the understanding of the allosteric phenomenon, broadly defined as coupling between remote molecular sites. Such advances have been driven not only by new theoretical models and pharmacological applications of allostery, but also by progress in the experimental approaches designed to map allosteric sites and transitions. Among these techniques, NMR spectroscopy has played a major role given its unique near-atomic resolution and sensitivity to the dynamics that underlie allosteric couplings. Here, we highlight recent progress in the NMR methods tailored to investigate allostery with the goal of offering an overview of which NMR approaches are best suited for which allosterically relevant questions. The picture of the allosteric "NMR toolbox" is provided starting from one of the simplest models of allostery (i.e., the four-state thermodynamic cycle) and continuing to more complex multistate mechanisms. We also review how such an "NMR toolbox" has assisted the elucidation of the allosteric molecular basis for disease-related mutations and the discovery of novel leads for allosteric drugs. From this overview, it is clear that NMR plays a central role not only in experimentally validating transformative theories of allostery, but also in tapping the full translational potential of allosteric systems.

Temperature-accelerated molecular dynamics gives insights into globular conformations sampled in the free state of the AC catalytic domain

The catalytic domain of the adenyl cyclase (AC) toxin from Bordetella pertussis is activated by interaction with calmodulin (CaM), resulting in cAMP overproduction in the infected cell. In the X-ray crystallographic structure of the complex between AC and the C terminal lobe of CaM, the toxin displays a markedly elongated shape. As for the structure of the isolated protein, experimental results support the hypothesis that more globular conformations are sampled, but information at atomic resolution is still lacking. Here, we use temperature-accelerated molecular dynamics (TAMD) simulations to generate putative all-atom models of globular conformations sampled by CaM-free AC. As collective variables, we use centers of mass coordinates of groups of residues selected from the analysis of standard molecular dynamics (MD) simulations. Results show that TAMD allows extended conformational sampling and generates AC conformations that are more globular than in the complexed state. These structures are then refined via energy minimization and further unrestrained MD simulations to optimize inter-domain packing interactions, thus resulting in the identification of a set of hydrogen bonds present in the globular conformations.