Kinetics of formation of hypoxanthine containing base pairs by HIV-RT: RNA template effects on the base substitution frequencies

Hypoxanthine (H), the deamination product of adenine, has been implicated in the high frequency of A to G transitions observed in retroviral and other RNA genomes. Although H·C base pairs are thermodynamically more stable than other H·N pairs, polymerase selection may be determined in part by kinetic factors. Therefore, the hypoxanthine induced substitution pattern resulting from replication by viral polymerases may be more complex than that predicted from thermodynamics. We have examined the steady-state kinetics of formation of base pairs opposite template H in RNA by HIV-RT, and for the incorporation of dITP during first- and second-strand synthesis. Hypoxanthine in an RNA template enhances the k2app for pairing with standard dNTPs by factors of 10–1000 relative to adenine at the same sequence position. The order of base pairing preferences for H in RNA was observed to be H·C >> H·T > H·A > H·G. Steady-state kinetics of insertion for all possible mispairs formed with dITP were examined on RNA and DNA templates of identical sequence. Insertion of dITP opposite all bases occurs 2–20 times more frequently on RNA templates. This bias for higher insertion frequencies on RNA relative to DNA templates is also observed for formation of mispairs at template A. This kinetic advantage afforded by RNA templates for mismatches and pairing involving H suggests a higher induction of mutations at adenines during first-strand synthesis by HIV-RT.

Energetics of the interaction between water and the helical peptide group and its role in determining helix propensities

The alanine helix provides a model system for studying the energetics of interaction between water and the helical peptide group, a possible major factor in the energetics of protein folding. Helix formation is enthalpy-driven (−1.0 kcal/mol per residue). Experimental transfer data (vapor phase to aqueous) for amides give the enthalpy of interaction with water of the amide group as ≈−11.5 kcal/mol. The enthalpy of the helical peptide hydrogen bond, computed for the gas phase by quantum mechanics, is −4.9 kcal/mol. These numbers give an enthalpy deficit for helix formation of −7.6 kcal/mol. To study this problem, we calculate the electrostatic solvation free energy (ESF) of the peptide groups in the helical and β-strand conformations, by using the delphi program and parse parameter set. Experimental data show that the ESF values of amides are almost entirely enthalpic. Two key results are: in the β-strand conformation, the ESF value of an interior alanine peptide group is −7.9 kcal/mol, substantially less than that of N-methylacetamide (−12.2 kcal/mol), and the helical peptide group is solvated with an ESF of −2.5 kcal/mol. These results reduce the enthalpy deficit to −1.5 kcal/mol, and desolvation of peptide groups through partial burial in the random coil may account for the remainder. Mutant peptides in the helical conformation show ESF differences among nonpolar amino acids that are comparable to observed helix propensity differences, but the ESF differences in the random coil conformation still must be subtracted.

Prediction of the stability of DNA triplexes.

We present rules that allow one to predict the stability of DNA pyrimidine.purine.pyrimidine (Y.R.Y) triple helices on the basis of the sequence. The rules were derived from van't Hoff analysis of 23 oligonucleotide triplexes tested at a variety of pH values. To predict the enthalpy of triplex formation (delta H degrees), a simple nearest-neighbor model was found to be sufficient. However, to accurately predict the free energy of the triplex (delta G degrees), a combination model consisting of five parameters was needed. These parameters were (i) the delta G degrees for helix initiation, (ii) the delta G degrees for adding a T-A.T triple, (iii) the delta G degrees for adding a C(+)-G.C triple, (iv) the penalty for adjacent C bases, and (v) the pH dependence of the C(+)-G.C triple's stability. The fitted parameters are highly consistent with thermodynamic data from the basis set, generally predicting both delta H degrees and delta G degrees to within the experimental error. Examination of the parameters points out several interesting features. The combination model predicts that C(+) -G.C. triples are much more stabilizing than T-A.T triples below pH 7.0 and that the stability of the former increases approximately equal to 1 kcal/mol per pH unit as the pH is decreased. Surprisingly though, the most stable sequence is predicted to be a CT repeat, as adjacent C bases partially cancel the stability of one another. The parameters successfully predict tm values from other laboratories, with some interesting exceptions.

The invariant system of coordinates of antibody molecules: prediction of the "standard" C alpha framework of VL and VH domains.

A new approach of comparing protein structures that does not involve the procedure of superposition is suggested. An invariant system of coordinates for immunoglobulin molecules that is based on the geometrical symmetry inherent to the variable domain light-chain (VL)-heavy-chain (VH) complex is described. The coordinates of the Calpha atoms in 22 immunoglobulin structures are calculated in the invariant system of coordinates. We found that 76 identical positions in this Calpha framework are symmetrical about the twofold axis. Comparison of the identical positions in these molecules allows us to select 96 positions in the light chains and 87 positions in the heavy chains whose Calpha atom coordinates are approximately the same. To check whether the average coordinates of Calpha atoms in these positions complies with the stereochemical requirements, we calculated Calpha-Calpha distances. Seventy-three positions of the light chains and 72 positions of the heavy chains satisfy the Calpha-Calpha distance criterion. The Calpha atoms in these positions are used for constructing the "standard" Calpha framework of VL and VH complexes. The average coordinates of Calpha atoms are presented.

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.

Isolation of HLA-DR1.(staphylococcal enterotoxin A)2 trimers in solution.

Mutational studies indicate that the superantigen staphylococcal enterotoxin A (SEA) has two separate binding sites for major histocompatibility complex (MHC) class II molecules. Direct evidence is provided here for the formation of SEA-MHC class II trimers in solution. Isoelectric focusing separated SEA-HLA-DR1 complexes into both dimers and HLA-DR1.SEA2 trimers. The molar ratio of components was determined by dual isotope labeling. The SEA mutant SEA-F47S, L48S, Y92A, which is deficient in MHC class II alpha-chain binding, formed only dimers with HLA-DR1, whereas a second SEA mutant, SEA-H225A, which lacks high-affinity MHC class II beta-chain binding was incapable of forming any complexes. Thus SEA binding to its MHC receptor is a two-step process involving initial beta-chain binding followed by cooperative binding of a second SEA molecule to the class II alpha chain.

Changing nucleosome positions through modification of the DNA rotational information.

The effects of the rotational information of DNA in determining the in vitro localization of nucleosomal core particles (ncps) have been studied in the Saccharomyces cerevisiae 5S rRNA repeat gene. We have altered the distribution of the phased series of flexibility signals present on this DNA by inserting a 25-bp tract, and we have analyzed the effects of this mutation on the distribution and on the frequencies of ncps, as compared with the wild type and a reference 21-bp insertion mutant. The variation of the standard free energy of nucleosome reconstitution was determined. The results show that the DNA rotational information is a major determinant of ncps positioning, define how many rotationally phased signals are required for the formation of a stable particle, and teach how to modify their distribution through the alteration of the rotational signals.