Transition path sampling (TPS) has been applied to the chemical step of human purine nucleoside phosphorylase (PNP). The transition path ensemble provides insight into the detailed mechanistic dynamics and atomic motion involved in transition state passage. The reaction mechanism involves early loss of the ribosidic bond to form a transition state with substantial ribooxacarbenium ion character, followed by dynamic motion from the enzyme and a conformational change in the ribosyl group leading to migration of the anomeric carbon toward phosphate, to form the product ribose 1-phosphate. Calculations of the commitment probability along reactive paths demonstrated the presence of a broad energy barrier at the transition state. TPS identified (i) compression of the O4⅐⅐⅐O5 vibrational motion, (ii) optimized leaving group interactions, and (iii) activation of the phosphate nucleophile as the reaction proceeds through the transition state region. Dynamic motions on the femtosecond timescale provide the simultaneous optimization of these effects and coincide with transition state formation.dynamics in catalysis ͉ promoting vibrations ͉ purine nucleoside phosphorylase ͉ transition path sampling ͉ transition state lifetime P rotein motions are widely accepted to be intimately connected to the catalytic function of enzymes. Important motions on different timescales have been experimentally and computationally observed, and a direct link between microsecond to millisecond domain motions and enzymatic function has been established (1-3). Although gating modes have been identified that link dynamics of the protein environment to hydrogen transfer (4) the role of dynamics on the timescale of bond vibrations (femtoseconds to picoseconds) in enzyme catalytic efficiency is less well understood. We have proposed that these sub-ps motions, which we termed ''protein promoting vibrations,'' may be coupled to the reaction coordinate and influence catalysis (5-8).Molecular dynamics (MD) simulations can provide details at an atomistic level, however, the accessible simulation time is usually in the picosecond to nanosecond timescale. It is currently not possible for a conventional MD simulation to study a rare chemical event, such as an enzymatic reaction, whose overall turnover rate is on the millisecond timescale. Techniques such as umbrella sampling, which impose a bias potential on a guess for the reaction coordinates, have been widely used in studies of enzymatic reaction mechanisms. Chandler and coworkers developed the Transition Path Sampling (TPS) method to overcome both the long timescale problem and the lack of knowledge about reaction mechanisms (9, 10). This method performs Monte Carlo (MC) walks through either trajectory space (MD) or stochastic (MC) space and is used in conjunction with either MD or MC. The principal advantage of TPS is to allow the study of reaction mechanisms in atomic detail, without prior knowledge of the reaction coordinate and transition state.The TPS algorithm has been applied in studies of lactase...

It has been found that with mutation of two surface residues (Lys(22) --> Glu and His(104) --> Arg) in human purine nucleoside phosphorylase (hPNP), there is an enhancement of catalytic activity in the chemical step. This is true although the mutations are quite remote from the active site, and there are no significant changes in crystallographic structure between the wild-type and mutant active sites. We propose that dynamic coupling from the remote residues to the catalytic site may play a role in catalysis, and it is this alteration in dynamics that causes an increase in the chemical step rate. Computational results indicate that the mutant exhibits stronger coupling between promotion of vibrations and the reaction coordinate than that found in native hPNP. Power spectra comparing native and mutant proteins show a correlation between the vibrations of Immucillin-G (ImmG):O5'...ImmG:N4' and H257:Ndelta...ImmG:O5' consistent with a coupling of these motions. These modes are linked to the protein promoting vibrations. Stronger coupling of motions to the reaction coordinate increases the probability of reaching the transition state and thus lowers the activation free energy. This motion has been shown to contribute to catalysis. Coincident with the approach to the transition state, the sum of the distances of ImmG:O4'...ImmG:O5'...H257:Ndelta became smaller, stabilizing the oxacarbenium ion formed at the transition state. Combined results from crystallography, mutational analysis, chemical kinetics, and computational analysis are consistent with dynamic compression playing a significant role in forming the transition state. Stronger coupling of these pairs is observed in the catalytically enhanced mutant enzyme. That motion and catalysis are enhanced by mutations remote from the catalytic site implicates dynamic coupling through the protein architecture as a component of catalysis in hPNP.

The spread of acquired immune deficiency syndrome has increasingly become a great concern owing largely to the failure of chemotherapies. The G48V is considered the key signature residue mutation of HIV-1 protease developing with saquinavir therapy. Molecular dynamics simulations of the wild-type and the G48V HIV-1 protease complexed with saquinavir were carried out to explore structure and interactions of the drug resistance. The molecular dynamics results combined with the quantum-based and molecular mechanics Poisson-Boltzmann surface area calculations indicated a monoprotonation took place on D25, one of the triad active site residues. The inhibitor binding of the triad residues and its interaction energy in the mutant were similar to those in the wild-type. The overall structure of both complexes is almost identical. However, the steric conflict of the substituted valine results in the conformational change of the P2 subsite and the disruption of hydrogen bonding between the -NH of the P2 subsite and the backbone -CO of the mutated residue. The magnitude of interaction energy changes was comparable to the experimental K(i) data. The designing for a new drug should consider a reduction of steric repulsion on P2 to enhance the activity toward this mutant strain.

No abstract

Human purine nucleoside phosphorylase (PNP) is a homotrimer, containing three nonconserved tryptophan residues at positions 16, 94, and 178, all remote from the catalytic site. The Trp residues were replaced with Tyr to produce Trp-free PNP (Leuko-PNP). Leuko-PNP showed near-normal kinetic properties. It was used (1) to determine the tautomeric form of guanine that produces strong fluorescence when bound to PNP, (2) for thermodynamic binding analysis of binary and ternary complexes with substrates, (3) in temperature-jump perturbation of complexes for evidence of multiple conformational complexes, and (4) to establish the ionization state of a catalytic site tyrosine involved in phosphate nucleophile activation. The (13)C NMR spectrum of guanine bound to Leuko-PNP, its fluorescent properties, and molecular orbital electronic transition analysis establish that its fluorescence originates from the lowest singlet excited state of the N1H, 6-keto, N7H guanine tautomer. Binding of guanine and phosphate to PNP and Leuko-PNP are random, with decreased affinity for formation of ternary complexes. Pre-steady-state kinetics and temperature-jump studies indicate that the ternary complex (enzyme-substrate-phosphate) forms in single binding steps without kinetically significant protein conformational changes as monitored by guanine fluorescence. Spectral changes of Leuko-PNP upon phosphate binding establish that the hydroxyl of Tyr88 is not ionized to the phenolate anion when phosphate is bound. A loop region (residues 243-266) near the purine base becomes highly ordered upon substrate/inhibitor binding. A single Trp residue was introduced into the catalytic loop of Leuko-PNP (Y249W-Leuko-PNP) to determine effects on catalysis and to introduce a fluorescence catalytic site probe. Although Y249W-Leuko-PNP is highly fluorescent and catalytically active, substrate binding did not perturb the fluorescence. Thermodynamic boxes, constructed to characterize the binding of phosphate, guanine, and hypoxanthine to native, Leuko-, and Y249W-Leuko-PNPs, establish that Leuko-PNP provides a versatile protein scaffold for introduction of specific Trp catalytic site probes.