Computational studies of the domain movement and the catalytic mechanism of thymidine phosphorylase.

Thymidine phosphorylase (TP) is a dual substrate enzyme with two domains. Each domain binds a substrate. In the crystal structure of Escherichia coli TP, the two domains are arranged so that the two substrate binding sites are too far away for the two substrates to directly react. Molecular dynamics simulations reveal a different structure of the enzyme in which the two domains have moved to place the two substrates in close contact. This structure has a root-mean-square deviation from the crystal structure of 4.1 A. Quantum mechanical calculations using this structure find that the reaction can proceed by a direct nucleophilic attack with a low barrier. This mechanism is not feasible in the crystal structure environment and is consistent with the mechanism observed for other N-glycosidic enzymes. Important catalytic roles are found for the three highly conserved residues His 85, Arg 171, and Lys 190.

Sequence diversity of SIV(Mne) Nef in vivo and in vitro.

We have compared nef gene sequences isolated by PCR from peripheral blood lymphocyte DNA of macaques which had been inoculated with either biologically or molecularly cloned SIV(Mne). Two samples from each animal obtained either early after infection (week 2-8) or after significant CD4+ depletion (week 21-137) were analyzed. Three substitutions in the predicted Nef amino acid sequence were seen in all animals at the late time point, and two more in all but one. Two of the common exchanges are located about 40 residues apart in the Nef core sequence, but are in proximity on the tertiary structure as judged by computer modelling using the structure of the HIV Nef core protein as a guide. Most recurring in vivo changes replaced a residue found in the cloned Nef sequence with one present in a consensus derived by aligning the Nef sequences of the SIVsm/HIV-2 groups. Animals inoculated with virus already containing the "late version" nef gene developed a more aggressive disease. The macaque adapted (MA)nef conferred a threefold higher infectivity to the cloned virus, but had no effects on CD4 downregulation. Propagation of virus with MAnef in tissue culture resulted in the rapid emergence of variants with newly attenuated nef. These findings suggest that the selective pressure on nef in vivo and in vitro are different.

Reaction path and free energy calculations of the transition between alternate conformations of HIV-1 protease.

Two different structures of ligand-free HIV protease have been determined by X-ray crystallography. These structures differ in the position of two 12 residue, beta-hairpin regions (or "flaps") which cap the active site. The movements of the flaps must be involved in the binding of substrates since, in either conformation, the flaps block the binding site. One of these structures is similar to structures of the ligand-bound enzyme; however, the importance of both structures to enzyme function is unclear. This transformation takes place on a time scale too long for conventional molecular dynamics simulations, so the process was studied by first identifying a reaction path between the two structures and then calculating the free energy along this path using umbrella sampling. For the ligand-free enzyme, it is found that the two structures are nearly equally stable, with the ligand-bound-type structure being less stable, consistent with X-ray crystallography data. The more stable open structure does not have a lower potential energy, but is stabilized by entropy. The transition occurs through a collapse and reformation of the beta-sheet structure of the conformationally flexible, glycine-rich flap ends. Additionally, some problems in studying conformational changes in proteins through the use of a single reaction path are addressed.

Molecular mechanisms of resistance: free energy calculations of mutation effects on inhibitor binding to HIV-1 protease.

The changes in the inhibitor binding constants due to the mutation of isoleucine to valine at position 84 of HIV-1 protease are calculated using molecular dynamics simulations. The calculations are done for three potent inhibitors--KNI-272, L-735,524 (indinavir or MK-639), and Ro 31-8959 (saquinavir). The calculations agree with the experimental data both in terms of an overall trend and in the magnitude of the resulting free energy change. HIV-1 protease is a homodimer, so each mutation causes two changes in the enzyme. The decrease in the binding free energy from each mutated side chain differs among the three inhibitors and correlates well with the size of the cavities induced in the protein interior near the mutated residue. The cavities are created as a result of a mutation to a smaller side chain, but the cavities are less than would be predicted from the wild-type structures, indicating that there is significant relaxation to partially fill the cavities.