X-ray structure and conformational dynamics of the HIV-1 protease in complex with the inhibitor SDZ283-910: agreement of time-resolved spectroscopy and molecular dynamics simulations

Based on the X-ray structure of the human immunodeficiency virus type-1 (HIV-1) protease in complex with the statine-derived inhibitor SDZ283-910, a 542 ps molecular dynamics trajectory was computed. For comparison with the 805 ps trajectory obtained for the uncomplexed enzyme, the theoretical fluorescence anisotropy decay of the unliganded protease and the inhibitor complex was calculated from the trajectories of the Trp6A/Trp6B and Trp42A/Trp42B transition dipole moments. This enabled us to directly compare the simulated data with the experimental picosecond time-resolved fluorescence data. Fitting both experimental and simulated data to the Kohlrausch-Williams-Watts (KWW) function exp(-t/tauk)beta revealed a very good agreement for the uncomplexed protease as well as for the SDZ283-910 complex. Binding of the inhibitor induced a faster decay of both the experimental and the computed protease fluorescence anisotropy decay. By this integrative approach, the atomic detail of inhibitor-induced changes in the conformational dynamics of the HIV-1 protease was experimentally verified and will be used for further inhibitor optimisation.     

Viscosity dependence of the solute quenching of the tryptophanyl fluorescence of proteins

We have studied the viscosity dependence of the acrylamide quenching of the fluorescence on the internal tryptophan residues in cod parvalbumin and ribonuclease T1, as well as the model systems. N-acetyl-L-tryptophanamide and glucagon. For the latter systems, the apparent rate constant, kq(app), for acrylamide quenching shows a typical diffusion-limited behavior. For parvalbumin and ribonuclease T1, however, the viscosity dependence of kq(app) is quite different. There is little change in the kq(app) values on increasing the bulk viscosity from 1 to 10 cP (by addition of glycerol), but a further increase from 10 to 100 cP results in a significant reduction in the kq(app). Both an unfolding mechanism and a quencher penetration mechanism are considered to explain the results. Only the penetration mechanism is found to be consistent, and our data are interpreted as indicating that the rate-limiting step for quenching goes from being that of diffusion through the protein matrix, at low viscosity, to diffusion through the bulk solvent, at high viscosity. By also considering the Kramers' relationship in fitting our data, we are able to obtain insight regarding the coupling between internal fluctuations in the structure of the protein and motion of the bulk solvent. For parvalbumin and ribonuclease T1, the internal dynamics are found to be very weakly coupled to the bulk.     

Molecular dynamics simulations of ribonuclease T1: comparison of the free enzyme and the 2' GMP-enzyme complex

Molecular dynamics simulations were performed on free RNase T1 and the 2'GMP-RNase T1 complex in vacuum and with water in the active site along with crystallographically identified waters, allowing analysis of both active site and overall structural and dynamics changes due to the presence of 2'GMP. Differences in the active site include a closing in the presence of 2'GMP, which is accompanied by a decrease in mobility of active site residues. The functional relevance of the active site fluctuations is discussed. 2'GMP alters the motion of Tyr-45, suggesting a role for that residue in providing a hydrophobic environment for the protein-nucleic acid interactions responsible for the specificity of RNase T1. The presence of 2'GMP causes a structural change of the C-terminus of the alpha-helix, indicating the transmission of structural changes from the active site through the protein matrix. Overall fluctuations of both the free and 2'GMP enzyme forms are in good agreement with X-ray temperature factors. The motion of Trp-59 is influenced by 2'GMP, indicating differences in enzyme dynamics away from the active site, with the calculated changes following those previously seen in time-resolved fluorescence experiments.     

Diffusional barrier crossing in a two-state protein folding reaction.

There has been some debate as to whether protein folding involves diffusive chain motions and thus depends on solvent viscosity. The interpretation of folding kinetics in viscous solvents has remained difficult and controversial, in that viscogenic agents affect folding rates not only by increasing solvent viscosity but also by increasing protein stability. By carefully choosing experimental conditions, we can now eliminate the effect on stability and show that the folding dynamics of the cold shock protein CspB are viscosity dependent. Thus Kramers' theory of reaction rates rather than transition state theory should be used to describe this folding reaction.     

Protein folding as a diffusional process.

A protein chain must move relative to the solvent molecules and explore many conformations when it folds from the extended unfolded state to the compact native state. Experimental and theoretical approaches suggest that diffusional processes in fact contribute to the kinetics of protein folding. We describe here how variations of the solvent viscosity can be employed to uncover the diffusional contributions to a folding reaction and assess the use of transition state theory and Kramers' rate theory for the analysis of protein folding reactions.     

Determination of the Backbone Mobility of Ribonuclease T1 and its 2′GMP Complex Using Molecular Dynamics Simulations and NMR Relaxation Data

Abstract

The results of 1-nanosecond molecular dynamics simulations of the enzyme ribonuclease T1 and its 2′GMP complex in water are presented. A classification of the angular reorientations of the backbone amide groups is achieved via a transformation of NH-vector trajectories into several coordinate frames, thus unravelling contributions of NH-bond librations and backbone dihedral angle fluctuations.

The former turned out to be similar for all amides, as characterized by correlation times of librational motions in a subpicosecond scale, angular amplitudes of about 10–12° for out-of-peptide-plane displacements of the NH-bond and 3–5° for the in-plane displacements, whereas the contributions of much slower backbone dihedral angle fluctuations strongly depend on the secondary structure. Correlation functions relevant for NMR were obtained and analyzed utilizing the ‘model-free’ approach (Lipari, G. and Szabo, A. (1982) J. Am. Chem. Soc. 104, 4546–4559,4559-4570; Clore et al., (1990) J. Am. Chem. Soc. 112, 4989–4991). The dependence of the amplitude of local motion on the residue location in the backbone is in good agreement with the results of NMR relaxation measurements and X-ray data. The protein dynamics is characterized by a highly restricted local motion of those parts of the backbone with defined secondary structure as well as by a high flexibility in loop regions. The comparison of results derived from different periods of the trajectory (of 50 ps and 1 ns duration, 1000 points sampled) reveals a dependence of the observed dynamic picture on the characteristic time scale of the experimental method used. Comparison of the MD data for the free and liganded enzyme clearly indicates a restriction of the mobility within certain regions of the backbone upon inhibitor binding.     

The role of protein–solvent hydrogen bond dynamics in the structural relaxation of a protein in glycerol versus water

We used MD simulations to investigate the dependence of the dynamics of a soluble protein, RNase A, on temperature and solvent environment. Consistent with neutron scattering data, the simulations predict that the protein undergoes a dynamical transition in both glycerol and aqueous solutions that is absent in the dry protein. The temperature of the transition is higher, while the rate of increase with temperature of the amplitudes of motion on the 100 ps timescale is lower, in glycerol versus water. Analysis of the dynamics of hydrogen bonds revealed that the protein dynamical transition is connected to the relaxation of the protein-solvent hydrogen bond network, which, in turn, is associated with solvent translational diffusion. Thus, it appears that the role of solvent dynamics in affecting the protein dynamical transition is qualitatively similar in water and glycerol.     

Molecular dynamics simulations and the conformational mobility of blood group oligosaccharides

Molecular dynamics simulations were carried out without explicit consideration of solvent to explore the conformational mobility of blood group A and H oligosaccharides. The potential energy force field of Rasmussen and co-workers was used with the CHARMM program on a number of disaccharide and trisaccharide models composed of fucose, galactose, glucose, N-acetyl glucosamine, and N-acetyl galactosamine chosen to represent various fragments of blood group oligosaccharides. In agreement with results of earlier studies, stable chair conformations were found for each pyranoside from which no transitions were detected in simulations as long as 800 ps. Exocyclic dihedral angles, including that of C5-C6, generally show numerous transitions on a time scale of approximately 5-30 ps. The dihedral angles of some but not all glycosidic linkages of blood group oligosaccharides show transitions on the time scale of 30-50 ps, implying that the extent of internal motion in blood group oligosaccharides depends strongly on linkage stereochemistry. For certain blood group A and H oligosaccharides that show limited internal motion in these simulations, we argue that the calculations are consistent with our previous analysis of 1H nuclear Overhauser enhancement (NOE) data that imply single conformations over a wide range of temperature and solvent conditions. While the trajectories are consistent with 13C T1 values that have been interpreted as indicating rigid conformations, measurements of 13C-NOE and T1 as a function of magnetic field strength are proposed as an improved method for experimental detection of the internal motion that is suggested for certain oligosaccharides in these simulations. The results of these simulations differ substantially from those of peptides of a similar molecular weight in that the oligosaccharides show much less internal motion.     

Hydration dependent dynamics in RNA

The essential role played by local and collective motions in RNA function has led to a growing interest in the characterization of RNA dynamics. Recent investigations have revealed that even relatively simple RNAs experience complex motions over multiple time scales covering the entire ms–ps motional range. In this work, we use deuterium solid-state NMR to systematically investigate motions in HIV-1 TAR RNA as a function of hydration. We probe dynamics at three uridine residues in different structural environments ranging from helical to completely unrestrained. We observe distinct and substantial changes in ²H solid-state relaxation times and lineshapes at each site as hydration levels increase. By comparing solid-state and solution state ¹³C relaxation measurements, we establish that ns–μs motions that may be indicative of collective dynamics suddenly arise in the RNA as hydration reaches a critical point coincident with the onset of bulk hydration. Beyond that point, we observe smaller changes in relaxation rates and lineshapes in these highly hydrated solid samples, compared to the dramatic activation of motion occurring at moderate hydration.     

Fluorescence and 13C NMR determination of side-chain and backbone dynamics of synthetic melittin and melittin analogues in isotropic solvents

The dynamics in isotopic solvents of selectively 13C labeled synthetic melittin and three analogues have been investigated by using NMR and fluorescence techniques both separately and in combination. In conjunction with the "model-free" approach to interpretation of NMR relaxation data [Lipari, G., & Szabo, A. (1982) J. Am. Chem. Soc. 104, 4546-4570], the availability of steady-state fluorescence anisotropy and lifetime data augment T1, T2, and NOE data to provide quantitative information about fluorophore dynamics in these peptides. A method is presented for using combined fluorescence and NMR data to obtain technique- and model-independent values for parameters describing local motion of 13C-labeled fluorophores in peptides and proteins. The dynamics of melittin and melittin analogues are found to be consistent with structural characteristics inferred from CD, fluorescence, and NMR spectral information presented in the preceding paper (Weaver et al., 1989). In particular, the mobility of the random coil peptide monomers is shown to be quite similar, while side-chain as well as peptide backbone motion in the aggregated or oligomeric species differs markedly among the analogues. For melittin itself, experimentally determined overall rotational correlation times for the monomer and tetramer agree very well with values predicted on the basis of solvent-accessible protein surface area. The local dynamics of selectively 13C-labeled Trp-19 and Gly-12 residues of melittin are also found to be consistent with peptide structure. In random coil melittin monomer, a specific model for the motion indicates that the Trp side chain moves through an approximate angle of +/- 71 degrees about the beta-gamma bond with a correlation time of 159 +/- 24 ps. In melittin tetramer, the indole moiety is spatially more confined with a flip angle of +/- 37 degrees, yet demonstrates an increased rate of motion with a correlation time of 56 +/- 8 ps. The constrained mobility of the Trp-19 side chain is consistent with motional constraints inferred from the X-ray structure of melittin tetramer. These results show that protein side-chain motion, even of moieties as large as indole, can occur on the picosecond time scale and that these motions are reasonably similar to those inferred from molecular dynamics simulations.     

Probing mechanisms for enzymatic activity enhancement of organophosphorus hydrolase in functionalized mesoporous silica

We have previously reported that organophosphorus hydrolase (OPH) can be spontaneously entrapped in functionalized mesoporous silica (FMS) with HOOC- as the functional groups and the entrapped OPH in HOOC-FMS showed enhanced enzyme specific activity. This work is to study the mechanisms that why OPH entrapped in FMS displayed the enhanced activity in views of OPH-FMS interactions using spectroscopic methods. The circular dichroism (CD) spectra show that, comparing to the secondary structure of OPH free in solution, OPH in HOOC-FMS displayed increased α-helix/β-strand transition of OPH with increased OPH loading density. The fluorescence emission spectra of Trp residues were used to assess the tertiary structural changes of the enzyme. There was a 42% increase in fluorescence. This is in agreement with the fact that the fluorescence intensity of OPH was increased accompanying with the increased OPH activity when decreasing urea concentrations in solution. The steady-state anisotropy was increased after OPH entrapping in HOOC-FMS comparing to the free OPH in solution, indicating that protein mobility was reduced upon entrapment. The solvent accessibility of Trp residues of OPH was probed by using acrylamide as a collisional quencher. Trp residues of OPH-FMS had less solvent exposure comparing with free OPH in solution due to its electrostatical binding to HOOC-FMS thereby displaying the increased fluorescence intensity. These results suggest the interactions of OPH with HOOC-FMS resulted in the protein immobilization and a favorable conformational change for OPH in the crowded confinement space and accordingly the enhanced activity.     

Free-radical label--new approach to the study of dynamics of lipid systems

A new method of EPR-spectroscopy--the recombination of free radicals appearing as a result of indirect radiolysis of biological molecules after a low temperature irradiation--is applied to the study of molecular dynamics of phosphatidylcholine dimyristoyl in mass and in the structure of liposomes above and below the transition temperature. It was shown, that the mobility of lipid molecules in crystalline liposomes is higher than in the structure of liquid-crystalline liposomes. The addition of cholesterol in liposome membranes decreases the lateral molecular motion of lipids in crystalline and liquid-crystalline state, in the latter case the effect of cholesterol addition is more pronounced. The activation energy for the displacement of the fragments of lipid molecules and the lipid molecule as a whole was estimated, and it was shown, that lipid matrix possesses a high degree of heterogeneity.     

Investigations of peptide hydration using NMR and molecular dynamics simulations: A study of effects of water on the conformation and dynamics of antamanide

Summary The influence of water binding on the conformational dynamics of the cyclic decapeptide antamanide dissolved in the model lipophilic environment chloroform is investigated by NMR relaxation measurements. The water-peptide complex has a lifetime of 35 s at 250 K, which is longer than typical lifetimes of water-peptide complexes reported in aqueous solution. In addition, there is a rapid intracomplex mobility that probably involves librational motions of the bound water or water molecules hopping between different binding sites. Water binding restricts the flexibility of antamanide. The experimental findings are compared with GROMOS molecular dynamics simulations of antamanide with up to eight bound water molecules. Within the simulation time of 600 ps, no water molecule leaves the complex. Additionally, the simulations show a reduced flexibility for the complex in comparison with uncomplexed antamanide. Thus, there is a qualitative agreement between the experimental NMR results and the computer simulations.     

Conformational dynamics and solvent viscosity effects in carboxypeptidase-A-catalyzed benzoylglycylphenyllactate hydrolysis

We have used a new approach to the dynamics of hydrolytic metalloenzyme catalysis based on investigations of both external solvent viscosity effects and kinetic 2H isotope effects. The former reflects solvent and protein dynamics, and the nuclear reorganization distribution among damped protein motion and intramolecular friction-free nuclear motion. The isotope effect represents proton tunnelling and reorganization in the hydrogen bond network around the active site. We illustrate the approach by new spectrophotometric and pH-titration data for carboxypeptidase-A-catalyzed benzoylglycyl-L-phenyllactate hydrolysis. This substrate exhibits both a significant inverse fractional power law viscosity dependence over wide ranges controlled by glycerol and sucrose, and a kinetic 2H isotope effect of 1.65. The analogous benzoylglycylphenylalanine hydrolysis has a smaller isotope effect (1.3) and no viscosity dependence. Viscosity variation has no effect on the CD spectra in the 180-240-nm range. In terms of stochastic chemical rate theory, the data correspond to an enzyme-peptide substrate complex with a 'tight' structure protected from the solvent. In comparison, the enzyme-ester substrate complex is 'softer', strongly coupled to the solvent, and the rate-determining step is accompanied by proton transfer or by substantial reorganization in the hydrogen bonds near the active site.     

Thermodynamics of a diffusional protein folding reaction

The folding reactions of several proteins are well described as diffusional barrier crossing processes, which suggests that they should be analyzed by Kramers' rate theory rather than by transition state theory. For the cold shock protein Bc-Csp from Bacillus caldolyticus, we measured stability and folding kinetics, as well as solvent viscosity as a function of temperature and denaturant concentration. Our analysis indicates that diffusional folding reactions can be treated by transition state theory, provided that the temperature and denaturant dependence of the solvent viscosity is properly accounted for, either at the level of the measured rate constants or of the calculated activation parameters. After viscosity correction the activation barriers for folding become less enthalpic and more entropic. The transition from an enthalpic to an entropic folding barrier with increasing temperature is, however, apparent in the data before and after this correction. It is a consequence of the negative activation heat capacity of refolding, which is independent of solvent viscosity. Bc-Csp and its mesophilic homolog Bs-CspB from Bacillus subtilis differ strongly in stability but show identical enthalpic and entropic barriers to refolding. The increased stability of Bc-Csp originates from additional enthalpic interactions that are established after passage through the activated state. As a consequence, the activation enthalpy of unfolding is increased relative to Bs-CspB.     

Structure and dynamics of the G121V dihydrofolate reductase mutant: lessons from a transition-state inhibitor complex

It is well known that enzyme flexibility is critical for function. This is due to the observation that the rates of intramolecular enzyme motions are often matched to the rates of intermolecular events such as substrate binding and product release. Beyond this role in progression through the reaction cycle, it has been suggested that enzyme dynamics may also promote the chemical step itself. Dihydrofolate reductase (DHFR) is a model enzyme for which dynamics have been proposed to aid in both substrate flux and catalysis. The G121V mutant of DHFR is a well studied form that exhibits a severe reduction in the rate of hydride transfer yet there remains dispute as to whether this defect is caused by altered structure, dynamics, or both. Here we address this by presenting an NMR study of the G121V mutant bound to reduced cofactor and the transition state inhibitor, methotrexate. NMR chemical shift markers demonstrate that this form predominantly adopts the closed conformation thereby allowing us to provide the first glimpse into the dynamics of a catalytically relevant complex. Based on (15)N and (2)H NMR spin relaxation, we find that the mutant complex has modest changes in ps-ns flexibility with most affected residues residing in the distal adenosine binding domain rather than the active site. Thus, aberrant ps-ns dynamics are likely not the main contributor to the decreased catalytic rate. The most dramatic effect of the mutation involves changes in μs-ms dynamics of the F-G and Met20 loops. Whereas loop motion is quenched in the wild type transition state inhibitor complex, the F-G and Met20 loops undergo excursions from the closed conformation in the mutant complex. These excursions serve to decrease the population of conformers having the correct active site configuration, thus providing an explanation for the G121V catalytic defect.     

Viscosity dependence of protein dynamics

The influence of solvent viscosity on protein dynamics was investigated with molecular dynamics simulations of factor Xa in two solvents differing only in viscosity, by a factor of 10. We obtained this viscosity change by changing the masses of the solvent atoms by a factor of 100. Equilibrium properties of the protein, that is, the average structure, its fluctuations, and the secondary structure, show no significant dependence on the solvent viscosity. The dynamic properties of the protein, that is, the atom-positional correlation times and torsional angle transitions, however, depend on the solvent viscosity. The protein appears to be much more mobile in the solvent of lower viscosity. It feels the influence of the solvent not only on the surface but even in its core. With increasing solvent viscosity, the positional relaxation times of atoms in the protein core increase as much as those of atoms on the protein surface, and the relative increase in the core is even larger than on the surface.     

The dynamics of protein hydration water: a quantitative comparison of molecular dynamics simulations and neutron-scattering experiments

We present results from an extensive molecular dynamics simulation study of water hydrating the protein Ribonuclease A, at a series of temperatures in cluster, crystal, and powder environments. The dynamics of protein hydration water appear to be very similar in crystal and powder environments at moderate to high hydration levels. Thus, we contend that experiments performed on powder samples are appropriate for discussing hydration water dynamics in native protein environments. Our analysis reveals that simulations performed on cluster models consisting of proteins surrounded by a finite water shell with free boundaries are not appropriate for the study of the solvent dynamics. Detailed comparison to available x-ray diffraction and inelastic neutron-scattering data shows that current generation force fields are capable of accurately reproducing the structural and dynamical observables. On the time scale of tens of picoseconds, at room temperature and high hydration, significant water translational diffusion and rotational motion occur. At low hydration, the water molecules are translationally confined but display appreciable rotational motion. Below the protein dynamical transition temperature, both translational and rotational motions of the water molecules are essentially arrested. Taken together, these results suggest that water translational motion is necessary for the structural relaxation that permits anharmonic and diffusive motions in proteins. Furthermore, it appears that the exchange of protein-water hydrogen bonds by water rotational/librational motion is not sufficient to permit protein structural relaxation. Rather, the complete exchange of protein-bound water molecules by translational displacement seems to be required.