Use of quantitative <Superscript>1</Superscript>H NMR chemical shift changes for ligand docking into barnase

¹H NMR complexation-induced changes in chemical shift (CIS) of HN protons have been used to characterize the complexes of barnase with the deoxyoligonucleotides d(GC) and d(CGAC). Quantitative shift changes are used not only to locate the most probable binding site (using ring-current shifts), but also to determine the orientation of the ligand within the binding site, based on a more complete shift calculation including bond magnetic anisotropies and electric field effects. For both ligands, the guanine is in the same binding site cleft, in the same position as identified in the crystal structure of the d(CGAC) complex. By contrast, a previous X-ray crystal structure of the d(GC) complex showed the ligand in the mouth of the active site, rather than at the guanyl-specific site, implying that the location may be an artifact of the crystallisation process.     

Flap opening and dimer-interface flexibility in the free and inhibitor-bound HIV protease, and their implications for function.

BACKGROUND: (1)H and (15)N transverse relaxation measurements on perdeuterated proteins are ideally suited for detecting backbone conformational fluctuations on the millisecond-microsecond timescale. The identification of conformational exchange on this timescale by measuring the relaxation of both (1)H and (15)N holds great promise for the elucidation of functionally relevant conformational changes in proteins. RESULTS: We measured the transverse (1)H and (15)N relaxation rates of backbone amides of HIV-1 protease in its free and inhibitor-bound forms. An analysis of these rates, obtained as a function of the effective rotating frame field, provided information about the timescale of structural fluctuations in several regions of the protein. The flaps that cover the active site of the inhibitor-bound protein undergo significant changes of backbone (φ,psi) angles, on the 100 micros timescale, in the free protein. In addition, the intermonomer beta-sheet interface of the bound form, which from protease structure studies appears to be rigid, was found to fluctuate on the millisecond timescale. CONCLUSIONS: We present a working model of the flap-opening mechanism in free HIV-1 protease which involves a transition from a semi-open to an open conformation that is facilitated by interaction of the Phe53 ring with the substrate. We also identify a surprising fluctuation of the beta-sheet intermonomer interface that suggests a structural requirement for maturation of the protease. Thus, slow conformational fluctuations identified by (1)H and (15)N transverse relaxation measurements can be related to the biological functions of proteins.     

Effect of nucleotide substrate binding on the pKa of catalytic residues in barnase

A thermodynamic cycle is used to describe barnase catalysis, which considers explicitly the presence of different ionic states of the catalytic residues Glu-73 and His-102 in barnase during the enzyme-substrate recognition process. Reinterpretation of published experimental data using rate equations derived from this cycle provides estimates of the ionization constants of these catalytic side chains, in the free enzyme and in the barnase-GpA complex. In addition, the electrostatic properties of the barnase-d(CGAC) crystal complex and of a barnase-5′3′(AAGAAp)-O-methyl ester modeled complex are investigated by means of a continuum approach to account for solvent polarization effects. Taking GpA as a reference substrate, it is shown that Increasing the length of the bound nucleotide induces pK_a shifts in the catalytic side chains, which modulate the fraction of enzyme in the correct ionic form for achieving the transesterification reaction. The computed results are in good agreement with the experimental variation of the optimum pH of barnase activity. The present analysis underscores the influence of pH effects on the k_cat and K_M kinetic constants of barnase and provides the basic formalism for linking the effective kinetic parameters, which usually depend on the pH, to the theoretical estimates of the true kinetic constants. © 1996 Wiley-Liss, Inc.     

Characterization of phosphate binding in the active site of barnase by site-directed mutagenesis and NMR

Phosphate is a competitive inhibitor of transesterification of GpC by the ribonuclease barnase. Barnase is significantly stabilized in the presence of phosphate against urea denaturation. The data are consistent with the existence of a single phosphate binding site in barnase with a dissociation constant, Kd, of 1.3 mM. The 2D 1H NMR spectrum of wild-type barnase with bound phosphate is assigned. Changes in chemical shifts and NOEs for wild type with bound phosphate compared with free wild type indicate that phosphate binds in the active site and that only small conformational changes occur on binding. Site-directed mutagenesis of the active site residues His-102, Lys-27, and Arg-87 to Ala increases the magnitude of Kd for phosphate by more than 20-fold. The 2D 1H NMR spectra of the mutants His-102----Ala, Lys-27----Ala, and Arg-87----Ala are assigned. Comparison with the spectra of wild-type barnase reveals that His-102----Ala and Lys-27----Ala have essentially the same structure as weild type, while some structural changes occur in Arg-87----Ala. It appears that phosphate binding by barnase is effected mainly by positively charge residues including His-102, Lys-27, and Arg-87. This may have applications for the design of phosphate binding sites in other proteins.     

Molecular dynamics simulation shows large volume fluctuations of proteins

In this paper we present a new approach to study the volume fluctuations of proteins. From a 1 ns molecular dynamics simulation, the volume fluctuation of human lysozyme has been calculated. We used two different ways for the calculation. In the first one, the volume fluctuation is extracted directly from the trajectory. For the second one, a newly developed formalism based on principal component analysis is used. The r.m.s. volume fluctuations obtained from the two analyses agree well with each other. The isothermal intrinsic compressibility was found to be larger than the one reported by experiment. The difference is discussed and suggested to exist in the assumed uncertainty of the compressibility of hydrated water to deduce the isothermal intrinsic compressibility from the experimental value. Spectral analysis shows that low-frequency dynamics dominate the total volume fluctuation. The same aspect is found in the study using principal component analysis. This low-frequency region is related to large and slow motions of proteins. Therefore a long time dynamics simulation is necessary to describe the volume fluctuations of proteins.     

Downhill binding energy surface of the barnase–barstar complex

We have employed biased molecular dynamics simulations in explicit solvent to characterize the one‐dimensional potential of mean force for the dissociation process of the barnase–barstar protein–protein complex. Unbinding of barstar from wild‐type barnase was compared with dissociation from four charge‐deletion mutants of barnase. Interestingly, we find in all cases that unbinding of barnase and barstar is an uphill process on a smooth, tilted energy landscape. The total free energy difference between the dissociated and bound state was similar for wild‐type barnase–barstar and for the R87A mutant of barnase. The values for the three other mutant barnase mutants K27A, R59A, and R83Q were only about half as much. Besides, we have analyzed the conformational dynamics of important residues at the barnase–barstar interface. In the bound state, their conformational fluctuations are reduced relatively to the free state because of the formation of intermolecular contacts. Interestingly, we find that some residues also show decreased mobility at intermediate stages of the unbinding process suggesting that these residues may be involved in the first contacts being formed on binding. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 977–985, 2010.     

Compressibility-structure relationship of globular proteins

The adiabatic compressibility, -beta s, of 11 globular proteins in water was determined by means of sound velocity measurements at 25 degrees C. All the proteins studied except for subtilisin showed positive -beta s values, indicating the large internal compressibility of the protein molecules. The intrinsic compressibility of proteins free from the hydration effect appeared to be comparable to that of normal ice. The compressibility data for 25 proteins, including 14 reported previously [Gekko, K., & Noguchi, H. (1979) J. Phys. Chem. 83, 2706-2714], were statistically analyzed to examine the correlation of the compressibility with some structural parameters and the amino acid compositions of proteins. It was found that -beta s increases with increasing partial specific volume and hydrophobicity of proteins. The helix element also seemed to be a dynamic domain to increase -beta s. Four amino acid residues (Leu, Glu, Phe, and His) greatly increased -beta s, and another four (Asn, Gly, Ser, and Thr) decreased it. Some empirical equations were derived for the estimation of the -beta s values of unknown proteins on the basis of their amino acid compositions. The volume fluctuations of proteins revealed by the compressibility data were in the range of 30-200 mL/mol, which corresponded to about 0.3% of the total protein volume. The conformational fluctuation seemed to enhance the thermal stability of proteins.     

Pressure-dependent changes in the solution structure of hen egg-white lysozyme

The "rules" governing protein structure and stability are still poorly understood. Important clues have come from proteins that operate under extreme conditions, because these clarify the physical constraints on proteins. One obvious extreme is pressure, but so far little is known of the behavior of proteins under pressure, largely for technical reasons. We have therefore developed new methodology for calculating structure change in solution with pressure, using NMR chemical shift changes, and we report the change in structure of lysozyme on going from 30 bar to 2000 bar, this being the first solution structure of a globular protein under pressure. The alpha-helical domain is compressed by approximately 1%, due to tighter packing between helices. The interdomain region is also compressed. By contrast, the beta-sheet domain displays very little overall compression, but undergoes more structural distortion than the alpha-domain. The largest volume changes tend to occur close to hydrated cavities. Because isothermal compressibility is related to volume fluctuation, this suggests that buried water molecules play an important role in conformational fluctuation at normal pressures, and are implicated as the nucleation sites for structural changes leading to pressure denaturation or channel opening.     

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.     

Conformational deformation in deoxymyoglobin by hydrostatic pressure

Pressure effect on the equilibrium conformation in sperm whale deoxymyoglobin and its volume fluctuation are studied by the normal mode analysis and strain tensor analysis. The pressure-induced deformation of interhelix regions are found to be remarkably more compressed than the other parts of the molecule. The intrahelix compressibility is shown to be relatively small. We also calculate the compressibility of the three hydrophobic clusters, located at the bottom, distal, and proximal side of the heme. Its value is found to decrease in the indicated order. The average compressibility of these hydrophobic clusters is less than the average interhelix compressibility, even though there are large cavities in these clusters. In order to study overall deformation, we define a linear compressibility and calculate it for all pairs of C^αatoms. The pressure-induced deformation is observed to be heterogeneous also in this analysis. The calculated root-mean-square displacement of the constituent atoms in the equilibrium conformation at 1,000 atm from those at 0 atm is 0.12 Å, which is much smaller in magnitude than the average value of the atomic fluctuations at room temperature. In the water solvent, the volume excluded by the protein molecule in the equilibrium conformation is reduced by 0.9%, when the pressure is raised from 0 to 1,000 atm. The calculated magnitude of the root-mean-square volume fluctuation is 0.3% of the excluded volume at room temperature. The square of the volume fluctuation is given as a sum of contributions from individual normal modes. Contributions from low frequency normal modes are found to dominate over those from higher frequency normal modes. The estimated value of the isothermal compressibility of deoxymyoglobin is 9.37 × 10^?12 cm²/dyn. © 1993 Wiley-Liss, Inc.     

Observation volumes and {gamma}-factors in two-photon fluorescence fluctuation spectroscopy

Fluorescence fluctuation spectroscopy has become an important measurement tool for investigating molecular dynamics, molecular interactions, and chemical kinetics in biological systems. Although the basic theory of fluctuation spectroscopy is well established, it is not widely recognized that saturation of the fluorescence excitation can dramatically alter the size and profile of the fluorescence observation volume from which fluorescence fluctuations are measured, even at relatively modest excitation levels. A precise model for these changes is needed for accurate analysis and interpretation of fluctuation spectroscopy data. We here introduce a combined analytical and computational approach to characterize the observation volume under saturating conditions and demonstrate how the variation in the volume is important in two-photon fluorescence correlation spectroscopy. We introduce a simple approach for analysis of fluorescence correlation spectroscopy data that can fully account for the effects of saturation, and demonstrate its success for characterizing the observed changes in both the amplitude and relaxation timescale of measured correlation curves. We also discuss how a quantitative model for the observed phenomena may be of broader importance in fluorescence fluctuation spectroscopy.     

Insight derived from molecular dynamics simulation into dynamics and molecular motions of cuticle-degrading serine protease Ver112

Cuticle-degrading serine protease Ver112, which derived from a nematophagous fungus Lecanicillium psalliotae, has been exhibited to have high cuticle-degrading and nematicidal activities. We have performed molecular dynamics (MD) simulation based on the crystal structure of Ver112 to investigate its dynamic properties and large-scale concerted motions. The results indicate that the structural core of Ver112 shows a small fluctuation amplitude, whereas the substrate binding sites, and the regions close to and opposite the substrate binding sites experience significant conformational fluctuations. The large concerted motions obtained from essential dynamics (ED) analysis of MD trajectory can lead to open or close of the substrate binding sites, which are proposed to be linked to the functional properties of Ver112, such as substrate binding, orientation, catalytic, and release. The significant motion in the loop regions that is located opposite the binding sites are considered to play an important role in modulating the dynamics of the substrate binding sites. Furthermore, the bottom of free energy landscape (FEL) of Ver112 are rugged, which is mainly caused by the fluctuations of substrate binding regions and loops located opposite the binding site. In addition, the mechanism underlying the high flexibility and catalytic activity of Ver112 was also discussed. Our simulation study complements the biochemical and structural studies, and provides insight into the dynamics-function relationship of cuticle-degrading serine protease Ver112.     

High pressure effect on the main transition from the ripple gel P′β phase to the liquid crystal (Lα) phase in dipalmitoylphosphatidylcholine. Microcalorimetric study

Scanning microcalorimetry has been used to study the high pressure effect on the main transition from the ripple gel P′_β phase to the liquid crystal (L_α) phase in DPPC (dipalmitoylphosphatidylcholine). It has been demonstrated that an increase of the pressure by 200 MPa shifts the transition to higher temperatures by 36.4 degrees. The pressure increase does not affect the cooperativity of transition but reduces noticeably its enthalpy. The changes of the molar partial volume, isothermal compressibility as well as volume thermal expansibility during transition in DPPC suspension have been estimated. It has been shown that monovalent ions (Na⁺, Cl⁻) in solution slightly affect the main thermodynamic parameters of the transition. Calcium ions significantly decrease distinction in compressibility and thermal expansibility between liquid-crystal and ripple gel phases of lipid suspension, which in its turn reflects less difference in their volume fluctuations.     

The influence of correlated protein-water volume fluctuations on the apparent compressibility of proteins determined by ultrasonic velocimetry

The elasticity of proteins, expressed by the compressibility, is potentially one of the most important properties of proteins because of the close relationship with its functionality. The compressibility of solutions can be determined by measurements of sound velocity and density. These quantities are related by the Newton-Laplace equation. In order to interpret the apparent compressibility of solutes in highly dilute solutions, it is required to consider the relation between compressibility and sound velocity of the solution using an appropriate reference system. The classical approach usually gives too small values for the apparent compressibility when compared with other methods. We show that the difference can partially be explained if the correlated volume fluctuations of the solvent are taken into consideration. A special attention is given to the compressibility of proteins. Finally, the present paper is not intended to replace established approaches, but it wants to create awareness that the classical mixing rules refer to ideal gasses assuming uncorrelated volume fluctuations and that a considerable part of the hydration effects could be explained by correlated volume fluctuations.     

Binding of bovine pancreatic trypsin inhibitor to trypsinogen: spectroscopic and volumetric studies

We have investigated the binding of bovine pancreatic trypsin inhibitor (BPTI) to bovine trypsinogen by combining ultrasonic velocimetry, high precision densimetry, and fluorescence spectroscopy. We report the changes in volume, adiabatic compressibility, van't Hoff enthalpy, entropy, and free energy that accompany the association of the two proteins at 25 degrees C and pH 8.0. We have used the measured changes in volume and compressibility in conjunction with available structural data to characterize the binding-induced changes in the hydration properties and intrinsic packing of the two proteins. Our estimate reveals that 110 +/- 40 water molecules become released to the bulk from the hydration shells of BPTI and trypsinogen. Furthermore, we find that the intrinsic coefficient of adiabatic compressibility of the two proteins decreases by 14 +/- 2%, which is suggestive of the binding-induced rigidification of the proteins' interior. BPTI-trypsinogen association is an entropy-driven event which proceeds with an unfavorable change in enthalpy. The favorable change in entropy results from partial compensation between two predominant terms. Namely, a large favorable change in hydrational entropy slightly prevails over a close in magnitude but opposite in sign change in configurational entropy. The reduction in configurational entropy and, consequently, protein dynamics is consistent with the observed decrease in intrinsic compressibility. In general, results of this work emphasize the vital role that water plays in modulating protein recognition events.     

Undulation Contributions to the Area Compressibility in Lipid Bilayer Simulations

It is here shown that there is a considerable system size-dependence in the area compressibility calculated from area fluctuations in lipid bilayers. This is caused by the contributions to the area fluctuations from undulations. This is also the case in experiments. At present, such a contribution, in most cases, is subtracted from the experimental values to obtain a true area compressibility. This should also be done with the simulation values. Here, this is done by extrapolating area compressibility versus system size, down to very small (zero) system size, where undulations no longer exist. The area compressibility moduli obtained from such simulations do not agree with experimental true area compressibility moduli as well as the uncorrected ones from contemporary or earlier simulations, but tend, instead, to be ∼50% too large. As a byproduct, the bending modulus can be calculated from the slope of the compressibility modulus versus system-size. The values obtained in this way for the bending modulus are then in good agreement with experiment.     

Ligand-dependent dynamics and intramolecular signaling in a PDZ domain

Allosteric communication is a fundamental process that proteins use to propagate signals from one site to functionally important distal sites. Although allostery is usually associated with multimeric proteins and enzymes, “long-range” communication may be a fundamental property of proteins. In some cases, communication occurs with minimal structural change. PDZ (post-synaptic density-95/discs large/zonula occludens-1) domains are small, protein-protein binding modules that can use multiple surfaces for docking diverse molecules. Furthermore, these domains have long-range energetic couplings that link the ligand-binding site to distal regions of the structure. Here, we show that allosteric behavior in a representative member of the PDZ domain family may be directly detected using side-chain methyl dynamics measurements. The changes in side-chain dynamics parameters in the second PDZ domain from the human tyrosine phosphatase 1E (hPTP1E) were determined upon binding a peptide target. Long-range dynamic effects were detected that correspond to previously observed pair-wise energetic couplings. These results provide one of the first experimental examples for the potential role of ps-ns timescale dynamics in propagating long-range signals within a protein, and reinforce the idea that dynamic fluctuations in proteins contribute to allosteric signal transduction.     

Contribution of the Cytoskeleton to the Compressive Properties and Recovery Behavior of Single Cells

The cytoskeleton is known to play an important role in the biomechanical nature and structure of cells, but its particular function in compressive characteristics has not yet been fully examined. This study focused on the contribution of the main three cytoskeletal elements to the bulk compressive stiffness (as measured by the compressive modulus), volumetric or apparent compressibility changes (as further indicated by apparent Poisson's ratio), and recovery behavior of individual chondrocytes. Before mechanical testing, cytochalasin D, acrylamide, or colchicine was used to disrupt actin microfilaments, intermediate filaments, or microtubules, respectively. Cells were subjected to a range of compressive strains and allowed to recover to equilibrium. Analysis of the video recording for each mechanical event yielded relevant compressive properties and recovery characteristics related to the specific cytoskeletal disrupting agent and as a function of applied axial strain. Inhibition of actin microfilaments had the greatest effect on bulk compressive stiffness (∼50% decrease compared to control). Meanwhile, intermediate filaments and microtubules were each found to play an integral role in either the diminution (compressibility) or retention (incompressibility) of original cell volume during compression. In addition, microtubule disruption had the largest effect on the “critical strain threshold” in cellular mechanical behavior (33% decrease compared to control), as well as the characteristic time for recovery (∼100% increase compared to control). Elucidating the role of the cytoskeleton in the compressive biomechanical behavior of single cells is an important step toward understanding the basis of mechanotransduction and the etiology of cellular disease processes.     

Fluctuations of the red blood cell membrane: relation to mechanical properties and lack of ATP dependence

We have analyzed the fluctuations of the red blood cell membrane in both the temporal ((ω(s⁻¹)) and spatial (q(m⁻¹)) frequency domains. The cells were examined over a range of osmolarities leading to cell volumes from 50% to 170% of that in the isotonic state. The fluctuations of the isotonic cell showed an ∼q⁻³-dependence, indicative of a motion dominated by bending, with an inferred bending modulus of ∼9 × 10⁻¹⁹J. When the cells were osmotically swollen to just below the point of lysis (166% of physiological volume), a q⁻¹-dependence of the fluctuations supervened, implying that the motion was now dominated by membrane tension; estimated as ∼1.3 × 10⁻⁴ nm⁻¹. When, on the other hand, the cells were osmotically dehydrated, the fluctuation amplitude progressively decreased. This was caused by a rise in internal viscosity, as shown by measurements on resealed ghosts containing a reduced hemoglobin concentration, which displayed no such effect. We examined, in addition, cells depleted of ATP, before the onset of echinocytosis, and could observe no change in fluctuation amplitude. We conclude that the membrane fluctuations of the red cell are governed by bending modulus, membrane tension, and cytosolic viscosity, with little or no dependence on the presence or absence of ATP.