Progress toward chemical accuracy in the computer simulation of condensed phase reactions.

We describe a procedure for the generation of chemically accurate computer-simulation models to study chemical reactions in the condensed phase. The process involves (i) the use of a coupled semiempirical quantum and classical molecular mechanics method to represent solutes and solvent, respectively; (ii) the optimization of semiempirical quantum mechanics (QM) parameters to produce a computationally efficient and chemically accurate QM model; (iii) the calibration of a quantum/classical microsolvation model using ab initio quantum theory; and (iv) the use of statistical mechanical principles and methods to simulate, on massively parallel computers, the thermodynamic properties of chemical reactions in aqueous solution. The utility of this process is demonstrated by the calculation of the enthalpy of reaction in vacuum and free energy change in aqueous solution for a proton transfer involving methanol, methoxide, imidazole, and imidazolium, which are functional groups involved with proton transfers in many biochemical systems. An optimized semiempirical QM model is produced, which results in the calculation of heats of formation of the above chemical species to within 1.0 kcal/mol (1 kcal = 4.18 kJ) of experimental values. The use of the calibrated QM and microsolvation QM/MM (molecular mechanics) models for the simulation of a proton transfer in aqueous solution gives a calculated free energy that is within 1.0 kcal/mol (12.2 calculated vs. 12.8 experimental) of a value estimated from experimental pKa values of the reacting species.

A structural basis for a phosphoramide mustard-induced DNA interstrand cross-link at 5'-d(GAC).

Phosphoramide mustard-induced DNA interstrand cross-links were studied both in vitro and by computer simulation. The local determinants for the formation of phosphoramide mustard-induced DNA interstrand cross-links were defined by using different pairs of synthetic oligonucleotide duplexes, each of which contained a single potentially cross-linkable site. Phosphoramide mustard was found to cross-link dG to dG at a 5'-d(GAC)-3'. The structural basis for the formation of this 1,3 cross-link was studied by molecular dynamics and quantum chemistry. Molecular dynamics indicated that the geometrical proximity of the binding sites also favored a 1,3 dG-to-dG linkage over a 1,2 dG-to-dG linkage in a 5'-d(GCC)-3' sequence. While the enthalpies of 1,2 and 1,3 mustard cross-linked DNA were found to be very close, a 1,3 structure was more flexible and may therefore be in a considerably higher entropic state.

Calicheamicin-DNA complexes: warhead alignment and saccharide recognition of the minor groove.

The solution structures of calicheamicin gamma 1I, its cycloaromatized analog (calicheamicin epsilon), and its aryl tetrasaccharide complexed to a common DNA hairpin duplex have been determined by NMR and distance-refined molecular dynamics computations. Sequence specificity is associated with carbohydrate-DNA recognition that places the aryl tetrasaccharide component of all three ligands in similar orientations in the minor groove at the d(T-C-C-T).d(A-G-G-A) segment. The complementary fit of the ligands and the DNA minor groove binding site creates numerous van der Waals contacts as well as hydrogen bonding interactions. Notable are the iodine and sulfur atoms of calicheamicin that hydrogen bond with the exposed amino proton of the 5'- and 3'-guanines, respectively, of the d(A-G-G-A) segment. The sequence-specific carbohydrate binding orients the enediyne aglycone of calicheamicin gamma 1I such that its C3 and C6 proradical centers are adjacent to the cleavage sites. While the enediyne aglycone of calicheamicin gamma 1I is tilted relative to the helix axis and spans the minor groove, the cycloaromatized aglycone is aligned approximately parallel to the helix axis in the respective complexes. Specific localized conformational perturbations in the DNA have been identified from imino proton complexation shifts and changes in specific sugar pucker patterns on complex formation. The helical parameters for the carbohydrate binding site are comparable with corresponding values in B-DNA fibers while a widening of the groove is observed at the adjacent aglycone binding site.

Crystal structure of Pseudomonas aeruginosa catabolic ornithine transcarbamoylase at 3.0-A resolution: a different oligomeric organization in the transcarbamoylase family.

The crystal structure of the Glu-105-->Gly mutant of catabolic ornithine transcarbamoylase (OTCase; carbamoyl phosphate + L-ornithine = orthophosphate + L-citrulline, EC 2.1.3.3) from Pseudomonas aeruginosa has been determined at 3.0-A resolution. This mutant is blocked in the active R (relaxed) state. The structure was solved by the molecular replacement method, starting from a crude molecular model built from a trimer of the catalytic subunit of another transcarbamoylase, the extensively studied aspartate transcarbamoylase (ATCase) from Escherichia coli. This model was used to generate initial low-resolution phases at 8-A resolution, which were extended to 3-A by noncrystallographic symmetry averaging. Four phase extensions were required to obtain an electron density map of very high quality from which the final model was built. The structure, including 4020 residues, has been refined to 3-A, and the current crystallographic R value is 0.216. No solvent molecules have been added to the model. The catabolic OTCase is a dodecamer composed of four trimers organized in a tetrahedral manner. Each monomer is composed of two domains. The carbamoyl phosphate binding domain shows a strong structural homology with the equivalent ATCase part. In contrast, the other domain, mainly implicated in the binding of the second substrate (ornithine for OTCase and aspartate for ATCase) is poorly conserved. The quaternary structures of these two allosteric transcarbamoylases are quite divergent: the E. coli ATCase has pseudo-32 point-group symmetry, with six catalytic and six regulatory chains; the catabolic OTCase has 23 point-group symmetry and only catalytic chains. However, both enzymes display homotropic and heterotropic cooperativity.