Shift in nucleotide conformational equilibrium contributes to increased rate of catalysis of GpAp versus GpA in barnase

Why the Energy Landscape of Barnase Is Hierarchical

We have used NMR and computational methods to characterize the dynamics of the ribonuclease barnase over a wide range of timescales in free and inhibitor-bound states. Using temperature- and denaturant-dependent measurements of chemical shift, we show that barnase undergoes frequent and highly populated hinge bending. Using relaxation dispersion, we characterize a slower and less populated motion with a rate of 750 ± 200 s⁻¹, involving residues around the lip of the active site, which occurs in both free and bound states and therefore suggests conformational selection. Normal mode calculations characterize correlated hinge bending motions on a very rapid timescale. These three measurements are combined with previous measurements and molecular dynamics calculations on barnase to characterize its dynamic landscape on timescales from picoseconds to milliseconds and length scales from 0.1 to 2.5 nm. We show that barnase has two different large-scale fluctuations: one on a timescale of 10⁻⁹−10⁻⁶ s that has no free energy barrier and is a hinge bending that is determined by the architecture of the protein; and one on a timescale of milliseconds (i.e., 750 s⁻¹) that has a significant free energy barrier and starts from a partially hinge-bent conformation. These two motions can be described as hierarchical, in that the more highly populated faster motion provides a platform for the slower (less probable) motion. The implications are discussed. The use of temperature and denaturant is suggested as a simple and general way to characterize motions on the intermediate ns-μs timescale.

Pressure-Dependent Conformation and Fluctuation in Folded Protein Molecules

Hydrostatic pressure leads to nonuniform compression of proteins. The structural change is on average only about 0.1 Å kbar(-1), and is therefore within the range of fluctuations at ambient pressure. The largest changes are around cavities and buried water molecules. Sheets distort much more than helices. Hydrogen bonds compress about 0.012 Å kbar(-1), although there is a wide range, including some hydrogen bonds that lengthen. In the presence of ligands and inhibitors, structural changes are smaller. Pressure has little effect on rapid fluctuations, but with larger scale slower motions, pressure increases the population of excited states (if they have smaller overall volume), and slows the fluctuations. In barnase, pressure is shown to be a useful way to characterise fluctuations on the timescale of microseconds, and helps to show that fluctuations in barnase are hierarchical, with the faster fluctuations providing a platform for the slower ones. The excited states populated at high pressure are probably functionally important.

Pressure-dependent structure changes in barnase on ligand binding reveal intermediate rate fluctuations

In this work we measured 1H NMR chemical shifts for the ribonuclease barnase at pressures from 3 MPa to 200 MPa, both free and bound to d(CGAC). Shift changes with pressure were used as restraints to determine the change in structure with pressure. Free barnase is compressed by approximately 0.7%. The largest changes are on the ligand-binding face close to Lys-27, which is the recognition site for the cleaved phosphate bond. This part of the protein also contains the buried water molecules. In the presence of d(CGAC), the compressibility is reduced by approximately 70% and the region of structural change is altered: the ligand-binding face is now almost incompressible, whereas changes occur at the opposite face. Because compressibility is proportional to mean square volume fluctuation, we conclude that in free barnase, volume fluctuation is largest close to the active site, but when the inhibitor is bound, the fluctuations become much smaller and are located mainly on the opposite face. The timescale of the fluctuations is nanoseconds to microseconds, consistent with the degree of ordering required for the fluctuations, which are intermediate between rapid uncorrelated side-chain dynamics and slow conformational transitions. The high-pressure technique is therefore useful for characterizing motions on this relatively inaccessible timescale.

Coupling of the guanosine glycosidic bond conformation and the ribonucleotide cleavage reaction: implications for barnase catalysis

To examine the possible relationship of guanine-dependent GpA conformations with ribonucleotide cleavage, two potential of mean force (PMF) calculations were performed in aqueous solution. In the first calculation, the guanosine glycosidic (Gchi) angle was used as the reaction coordinate, and computations were performed on two GpA ionic species: protonated (neutral) or deprotonated (negatively charged) guanosine ribose O2 '. Similar energetic profiles featuring two minima corresponding to the anti and syn Gchi regions were obtained for both ionic forms. For both forms the anti conformation was more stable than the syn, and barriers of approximately 4 kcal/mol were obtained for the anti --> syn transition. Structural analysis showed a remarkable sensitivity of the phosphate moiety to the conformation of the Gchi angle, suggesting a possible connection between this conformation and the mechanism of ribonucleotide cleavage. This hypothesis was confirmed by the second PMF calculations, for which the O2 '--P distance for the deprotonated GpA was used as reaction coordinate. The computations were performed from two selected starting points: the anti and syn minima determined in the first PMF study of the deprotonated guanosine ribose O2'. The simulations revealed that the O2 ' attack along the syn Gchi was more favorable than that along the anti Gchi: energetically, significantly lower barriers were obtained in the syn than in the anti conformation for the O--P bond formation; structurally, a lesser O2 '--P initial distance, and a better suited orientation for an in-line attack was observed in the syn relative to the anti conformation. These results are consistent with the catalytically competent conformation of barnase-ribonucleotide complex, which requires a guanine syn conformation of the substrate to enable abstraction of the ribose H2 ' proton by the general base Glu73, thereby suggesting a coupling between the reactive substrate conformation and enzyme structure and mechanism.

Propagation of Dynamic Changes in Barnase Upon Binding of Barstar: An NMR and Computational Study

NMR spectroscopy and computer simulations were used to examine changes in chemical shifts and in dynamics of the ribonuclease barnase that result upon binding to its natural inhibitor barstar. Although the spatial structures of free and bound barnase are very similar, binding results in changes of the dynamics of both fast side-chains, as revealed by (2)H relaxation measurements, and NMR chemical shifts in an extended beta-sheet that is located far from the binding interface. Both side-chain dynamics and chemical shifts are sensitive to variations in the ensemble populations of the inter-converting molecular states, which can escape direct structural observation. Molecular dynamics simulations of free barnase and barnase in complex with barstar, as well as a normal mode analysis of barnase using a Gaussian network model, reveal relatively rigid domains that are separated by the extended beta-sheet mentioned above. The observed changes in NMR parameters upon ligation can thus be rationalized in terms of changes in inter-domain dynamics and in populations of exchanging states, without measurable structural changes. This provides an alternative model for the propagation of a molecular response to ligand binding across a protein that is based exclusively on changes in dynamics.

The catalytic power of enzymes: Conformational selection or transition state stabilization?

The mechanism by which enzymes produce enormous rate enhancements in the reactions they catalyze remains unknown. Two viewpoints, selection of ground state conformations and stabilization of the transition state, are present in the literature in apparent opposition. To provide more insight into current discussion about enzyme efficiency, a two-state model of enzyme catalysis was developed. The model was designed to include both the pre-chemical (ground state conformations) and the chemical (transition state) components of the process for the substrate both in water and in the enzyme. Although the model is of general applicability, the chorismate to prephenate reaction catalyzed by chorismate mutase was chosen for illustrative purposes. The resulting kinetic equations show that the catalytic power of enzymes, quantified as the k(cat)/k(uncat) ratio, is the product of two terms: one including the equilibrium constants for the substrate conformational states and the other including the rate constants for the uncatalyzed and catalyzed chemical reactions. The model shows that these components are not mutually exclusive and can be simultaneously present in an enzymic system, being their relative contribution a property of the enzyme. The developed mathematical expressions reveal that the conformational and reaction components of the process perform differently for the translation of molecular efficiency (changes in energy levels) into observed enzymic efficiency (changes in k(cat)), being, in general, more productive the component involving the transition state.