Dynamic structures of horse liver alcohol dehydrogenase (HLADH): results of molecular dynamics simulations of HLADH-NAD(+)-PhCH(2)OH, HLADH-NAD(+)-PhCH(2)O(-), and HLADH-NADH-PhCHO.

Molecular dynamics simulations of the oxidation of benzyl alcohol by horse liver alcohol dehydrogenase (HLADH) have been carried out. The following three states have been studied: HLADH.PhCH(2)OH.NAD(+) (MD1), HLADH.PhCH(2)O(-).NAD(+) (MD2), and HLADH.PhCHO.NADH (MD3). MD1, MD2, and MD3 simulations were carried out on one of the subunits of the dimeric enzyme covered in a 32-A-radius sphere of TIP3P water centered on the active site. The proton produced on ionization of the alcohol when HLADH.PhCH(2)OH.NAD(+) --> HLADH.PhCH(2)O(-).NAD(+) is transferred from the active site to solvent water via a hydrogen bonding network consisting of serine48 hydroxyl, ribose 2'- and 3'-hydroxyl groups, and Hist51. Hydrogen bonding of the 3'OH of ribose to Ile269 carbonyl maintains this proton in position to be transferred to water. Molecular dynamic simulations have been employed to track water1287 from the TIP3 water pool to the active site, thus exhibiting the mode of entrance of water to the active site. With time the water1287 accumulates in two different positions in order to accept the proton from the ribose 3'-OH and from His51. There can be identified two structural substates for proton passage. In the first substate the imidazole Ne2 of His51 is adjacent to the nicotinamide ribose C2'-OH and hydrogen bonding distances for proton transfer through the hydrogen bonded relay series PhCH(2)OH...Ser48-OH...Ribose2'-OH...His51...OH(2) (path 1) average 2.0, 2.0, and 2.1 A and (for His51...OH(2)) minimal distances less or equal to 2.5 A. The structure for path 1 is present 20% of the time span. And in the second substate, there are two possible proton passages: path 1 as before and path 2. Path 2 involves the hydrogen-bonded relay series PhCH(2)OH...Ser48-OH...Ribose2'-OH...Ribose3'-OH...His51.OH(2) with the average bonding distances being 2.0, 2.0, 2.1, and 2.0 A and (for His51...OH(2)) minimal distances less or equal to 2.5 A (20% probability of the time span), respectively. During the molecular dynamics simulation the NAD(+) ribose conformations have stabilized at the C2'-endo-C3'-exo or the C2'-endo conformations. With the C2'-endo conformation the first and second substates are able to persist for different time spans, while with the C2'-endo-C3'-exo conformation the only possible pathway involves the first substate. For both first and second substates the fluctuation of the distances between the ribose-OH protons and N epsilon 2 of His51 imidazole ring is partially contributed by the "windshield wiper" motion of the His51 imidazole ring. Since the imidazole of His-51 contributes only about 10-fold to activity, as estimated from the decrease in activity upon substitution with a Gln, there must be an alternate route for the proton to pass to solvent without going through this histidine. A third pathway involves ribose C3'-OH and Ile-269. In MD2, near attack conformers (NACs) for hydride transfer from PhCH(2)O(-) to NAD(+) represent approximately 60% of E.S conformers. The molecular dynamic study of MD3 at mildly basic pH reveals that reactive ground state conformers (NACs) for hydride transfer from NADH to PhCHO amount to 12 mol % of conformers. In MD3, anisotropic bending of the dihydronicotinamide ring of NADH (average value of alpha(c) = 4.0 degrees and alpha(n) = 0.5 degrees, respectively) is observed.     

The importance of being r: greater oxidative stability of RNA compared with DNA

The 2'-hydroxyl group of ribose imparts hydrolytic lability on RNA, which provides a mechanism for numerous biological functions. Recent evidence from chemical cleavage studies shows that this hydroxyl group also stabilizes the sugar moiety in RNA towards oxidation relative to DNA. Is this just because RNA needs to be distinguishable from DNA or does it have other evolutionary significance?     

Detection of scalar couplings involving 2'-hydroxyl protons across hydrogen bonds in a frameshifting mRNA pseudoknot

The -1 frameshift-stimulating mRNA pseudoknot from pea enation mosaic virus-1 (PEMV-1) is composed nearly entirely of RNA triple helix.(4) The 2'-OH hydroxyl protons of riboses C15 and C16 are hydrogen bond donors to the N1 atoms of adenosines A27 and A25, respectively, positioned in the minor groove of pseudoknot stem S1. In this paper, a nonrefocused (1)H,(15)N CPMG HSQC of uniformly (13)C,(15)N-labeled 33-mer PEMV-1 RNA has been tailored to reveal a correlation of the 2'-OH hydroxyl proton of C15 to the N1 nitrogen resonance of A27 mediated by a cross hydrogen bond scalar coupling. The (1h)J(2'OH,N) cross hydrogen bond scalar coupling constant determined from a quantitative 1D (15N) spin-echo difference experiment for the C15/A27 interaction is 1.7 +/- 0.1 Hz, while that for the C16/A25 interaction appears larger, 3.5 +/- 0.3 Hz, despite the fact that the corresponding direct correlation between the 2'-OH hydroxyl proton of C16 and the N1 of A25 is missing due to unfavorable solvent exchange properties. These findings reveal a detailed picture of critical noncanonical O-H.N hydrogen-bonding loop-stem interactions in an RNA triple helical structure.     

2'-hydroxyl proton positions in helical RNA from simultaneously measured heteronuclear scalar couplings and NOEs

The 2'-hydroxyl group in RNA plays an important structural role; it defines hydration in the minor groove, impacts thermodynamic stability of RNA, and often participates in RNA catalysis. To better study this important functional group in RNA, we describe a constant-time HMQC-IPAP-NOESY 3D NMR experiment. It simultaneously yields highly resolved 13C-separated NOEs from ribose protons to the 2'OH proton, as well as E.COSY-type measurement of JC-OH couplings, thereby permitting a quantitative study of the orientation of the 2'OH proton in RNA. The observed NOE patterns indicate that the 2'OH bonds in A-form helical RNA are primarily oriented toward the base domain, as further supported by small (1.3 +/- 0.7 Hz) 3JC1'-2'OH and relatively large (4.2 +/- 0.7 Hz) 3JC3'-2'OH couplings. The constant-time HMQC-IPAP-NOESY is suitable for measurement of interactions of rapidly exchanging protons in proteins and nucleic acids.     

The hydrogen bonding and hydration of 2'-OH in adenosine and adenosine 3'-ethyl phosphate

The 2'-OH group has major structural implications in the recognition, processing, and catalytic properties of RNA. We report here intra- and intermolecular H-bonding of 2'-OH in adenosine 3'-ethyl phosphate (1), 3'-deoxyadenosine (2), and adenosine (3) by both temperature- and concentration-dependent NMR studies, as well as by detailed endo ((3)J(H,H)) and exocyclic ((3)J(H,OH)) coupling constant analyses. We have also examined the nature of hydration and exchange processes of 2'-OH with water by a combination of NOESY and ROESY experiments in DMSO-d(6) containing 2 mol % HOD. The NMR-constrained molecular modeling (by molecular mechanics as well as by ab initio methods both in the gas and solution phase) has been used to characterize the energy minima among the four alternative dihedrals possible from the solution of the Karplus equation for (3)J(H2',OH) and (3)J(H3',OH) to delineate the preferred orientation of 2'-O-H proton in 1 and 2 as well as for 2'/3'-O-H protons in 3. The NMR line shape analysis of 2'-OH gave the DeltaG(H-bond)(298K) of 7.5 kJ mol(-1) for 1 and 8.4 kJ mol(-1) for 3; similar analyses of the methylene protons of 3'-ethyl phosphate moiety in 1 also gave comparable DeltaG(H-bond)(298K) of 7.3 kJ mol(-1). The donor nature of the 2'-OH in the intramolecular H-bonding in 3 is evident from its relatively reduced flexibility [-TDeltaS++](2'-OH) = -17.9(+/-0.5) kJ mol(-1)] because of the loss of conformational freedom owing to the intramolecular 2'O-H...O3' H-bonding, compared to the acceptor 3'-OH in 3 [-TDeltaS++](3'-OH) = -19.8 (+/- 0.6) kJ mol(-1)] at 298 K. The presence of intramolecular 2'-OH...O3' H-bonding in 3 is also corroborated by the existence of weak long-range (4)J(H2',OH3') in 3 (i.e., W conformation of H2'-C2'-C3'-O3'-H) as well as by (3)J(H,OH) dependent orientation of the 2'- and 3'-OH groups. The ROESY spectra for 1 and 3 at 308 K, in DMSO-d(6), show a clear positive ROE contact of both 2'- and 3'-OH with water. The presence of a hydrophilic 3'-phosphate group in 1 causes a much higher water activity in the vicinity of its 2'-OH, which in turn causes the 2'-OH to exchange faster, culminating in a shorter exchange lifetime (tau) for 2'-OH proton with HOD in 1 (tau2'-OH: 489 ms) compared to that in 3 (tau2'-OH: 6897 ms). The activation energy (E(a)) of the exchange with the bound-water for 2'- and 3'-OH in 3 (48.3 and 45.0 kJ mol(-1), respectively) is higher compared to that of 2'-OH in 1 (31.9 kJ mol(-1)), thereby showing that the kinetic availability of hydrated 2'-OH in 1 for any inter- and intramolecular interactions, in general, is owing to the vicinal 3'-phosphate residue. It also suggests that 2'-OH in native RNA can mediate other inter- or intramolecular interactions only in competition with the bound-water, depending upon the specific chemical nature and spatial orientation of other functions with potential for hydrogen bonding in the neighborhood. This availability of the bound water around 2'-OH in RNA would, however, be dictated by whether the vicinal phosphate is exposed to the bulk water or not. This implies that relatively poor hydration around a specific 2'-OH across a polyribonucleotide chain, owing to some hydrophobic microenvironmental pocket around that hydroxyl, may make it more accessible to interact with other donor or acceptor functions for H-bonding interactions, which might then cause the RNA to fold in a specific manner generating a new motif leading to specific recognition and function. Alternatively, a differential hydration of a specific 2'-OH may modulate its nucleophilicity to undergo stereospecific transesterification reaction as encountered in ubiquitous splicing of pre-mRNA to processed RNA or RNA catalysis, in general.     

Non-Watson-Crick base pairing in RNA. quantum chemical analysis of the cis Watson-Crick/sugar edge base pair family

Large RNA molecules exhibit an astonishing variability of base-pairing patterns, while many of the RNA base-pairing families have no counterparts in DNA. The cis Watson-Crick/sugar edge (cis WC/SE) RNA base pairing is investigated by ab initio quantum chemical calculations. A detailed structural and energetic characterization of all 13 crystallographically detected members of this family is provided by means of B3LYP/6-31G and RIMP2/aug-cc-pVDZ calculations. Further, a prediction is made for the remaining 3 cis WC/SE base pairs which are yet to be seen in the experiments. The interaction energy calculations point at the key role of the 2'-OH group in stabilizing the sugar-base contact and predict all 16 cis WC/SE base-pairing patterns to be nearly isoenergetic. The perfect correlation of the main geometrical parameters in the gas-phase optimized and X-ray structures shows that the principle of isosteric substitutions in RNA is rooted from the intrinsic structural similarity of the isolated base pairs. The present quantum chemical calculations for the first time analyze base pairs involving the ribose 2'-OH group and unambiguously correlate the structural information known from experiments with the energetics of interactions. The calculations further show that the relative importance and absolute value of the dispersion energy in the cis WC/SE base pairs are enhanced compared to the standard base pairs. This may by an important factor contributing to the strength of such interactions when RNA folds in its polar environment. The calculations further demonstrate that the Cornell et al. force field commonly used in molecular modeling and simulations provides satisfactory performance for this type of RNA interactions.     

2'-Fluoro substituents can mimic native 2'-hydroxyls within structured RNA

The ability of fluorine in a C-F bond to act as?a hydrogen bond acceptor is controversial. To test such ability in complex RNA macromolecules, we have replaced native 2'-OH groups with 2'-F and 2'-H groups in two related systems, the Tetrahymena group I ribozyme and the ΔC209 P4-P6 RNA domain. In three cases the introduced 2'-F mimics the native 2'-OH group, suggesting that the fluorine atom can accept a hydrogen bond. In each of these cases the native hydroxyl group interacts with a purine exocyclic amine. Our results give insight about the properties of organofluorine and suggest a possible general biochemical signature for tertiary interactions between 2'-hydroxyl groups and exocyclic amino groups within RNA.     

Synthesis of RNA using 2'-O-DTM protection

tert-Butyldithiomethyl (DTM), a novel hydroxyl protecting group, cleavable under reductive conditions, was developed and applied for the protection of 2'-OH during solid-phase RNA synthesis. This function is compatible with all standard protecting groups used in oligonucleotide synthesis, and allows for fast and high-yield synthesis of RNA. Oligonucleotides containing the 2'-O-DTM groups can be easily deprotected under the mildest possible aqueous and homogeneous conditions. The preserved 5'-O-DMTr function can be used for high-throughput cartridge RNA purification.     

The hydration of monosaccharides—an NMR study

At temperatures close to 0°C proton exchange between sugar hydroxyl groups and water is slow, and separate proton resonance peaks can be detected for the hydroxyl protons. All are shifted downfield of the water resonance, the anomeric hydroxyl proton shift being the greatest. Axial anomeric hydroxyl protons are shifted less than corresponding equatorial protons. Proton exchange with water is strongly acid and base catalyzed, but, at least in some cases, there seems to be an additional pH-independent mechanism involved. From the temperature effect on the shifts, and the effect of added dimethyl sulfoxide, we conclude that each hydroxyl group is bonded on average to two water molecules. This estimate of the hydration number for monosaccharides is far greater than those previously deduced from relaxation studies. It is suggested that the source of this difference lies in the residence times of the bound water molecules. Shifts of the hydroxyl proton resonances for sugars in methanol are compared with those for aqueous solutions and are found to be very similar. Hence it is concluded that these shifts do not reveal any special effects due to water structure. There are quite marked differences in the shifts for different sugars, and, in particular, the anomeric hydroxyl proton shifts for ketoses are smaller than those for aldoses.Taken as solvation spectra, Part 56.     

Principles of RNA base pairing: structures and energies of the trans Watson-Crick/sugar edge base pairs

Due to the presence of the 2'-OH hydroxyl group of ribose, RNA molecules utilize an astonishing variability of base pairing patterns to build up their structures and perform the biological functions. Many of the key RNA base pairing families have no counterparts in DNA. In this study, the trans Watson-Crick/sugar edge (trans WC/SE) RNA base pair family has been characterized using quantum chemical and molecular mechanics calculations. Gas-phase optimized geometries from density functional theory (DFT) calculations and RIMP2 interaction energies are reported for the 10 crystallographically identified trans WC/SE base pairing patterns. Further, stable structures are predicted for all of the remaining six possible members of this family not seen in RNAs so far. Among these novel six base pairs, the computations substantially refine two structures suggested earlier based on simple isosteric considerations. For two additional trans WC/SE base pairs predicted in this study, no arrangement was suggested before. Thus, our study brings a complete set of trans WC/SE base pairing patterns. The present results are also contrasted with calculations reported recently for the cis WC/SE base pair family. The computed base pair sizes are in sound correlation with the X-ray data for all WC/SE pairing patterns including both their cis and trans isomers. This confirms that the isostericity of RNA base pairs, which is one of the key factors determining the RNA sequence conservation patterns, originates in the properties of the isolated base pairs. In contrast to the cis structures, however, the isosteric subgroups of the trans WC/SE family differ not only in their H-bonding patterns and steric dimensions but also in the intrinsic strength of the intermolecular interactions. The distribution of the total interaction energy over the sugar-base and base-base contributions is controlled by the cis-trans isomerism.     

Oligodeoxyribonucleotide analogues with bridging dimethylene sulfide,sulfoxide, and sulfone groups. Toward a second-generation model of nucleic acid structure

Short DNA analogues with bridging dimethylene sulfide, sulfoxide, and sulfone groups replacing the phosphate diesters (S-DNAs) were synthesized from building blocks prepared via two routes, both starting from D-glucose. Building blocks for RNA analogues were prepared by stereoselective introduction of nucleobase into a 2'-acylated ribose analogue. The ribose analogues were converted to deoxyribose analogues by replacement of a 3'-OH group by a thioacetyl unit, followed by photolytic deoxygenation or radical-based 2'-deoxygenation. DNA analogues joined via CH(2)(-)S-CH(2) units were prepared by S(N)2 displacement of a 6'-mesyl group on one building block using a thiolate nucleophile of another. 4,4'-Dimethoxytrityl protection and deprotection schemes were established for both the thiol and hydroxyl groups. The corresponding sulfoxide DNA analogues were obtained by oxidation with hydrogen peroxide. Sulfone DNA analogues were obtained by oxidation of the sulfide DNA with persulfate or hydrogen peroxide in the presence of a titanium silicate catalyst. The physical properties of several representative oligonucleotide analogues were examined, and interpreted in light of a "second-generation" model for DNA strand-strand recognition, a model that emphasizes the role of the polyanionic backbone in diminishing unwanted tendencies of highly functionalized molecules to form "structure" in solution. Even short sulfide-linked DNA analogues displayed association properties different from those displayed by standard DNA molecules. Complex formation observed with sulfide-linked tetramers by HPLC study in different solvents suggested that the complex is formed using hydrogen bonding. Sulfone-linked dinucleotides display Watson-Crick behavior; the tetramer, however, displayed self-structure. Self-structure and self-aggregation become more prominent as the length of the oligonucleotide analogues increases. The tendency to self-aggregate can be decreased by adding a charged sulfonate group to the 3'-end of the DNA analogue. Features of the second-generation model are important for many areas of nucleic acid chemistry, from the design of nucleic acid therapeutic agents to the search for life on other planets.     

Theoretical study of the scalar coupling constants across the noncovalent contacts in RNA base pairs: the cis- and trans-watson-crick/sugar edge base pair family

The structure and function of RNA molecules are substantially affected by non-Watson-Crick base pairs actively utilizing the 2'-hydroxyl group of ribose. Here we correlate scalar coupling constants across the noncovalent contacts calculated for the cis- and trans-WC/SE (Watson-Crick/sugar edge) RNA base pairs with the geometry of base to base and sugar to base hydrogen bond(s). 23 RNA base pairs from the 32 investigated were found in RNA crystal structures, and the calculated scalar couplings are therefore experimentally relevant with regard to the binding patterns occurring in this class of RNA base pairs. The intermolecular scalar couplings 1hJ(N,H), 2hJ(N,N), 2hJ(C,H), and 3hJ(C,N) were calculated for the N-H...N and N-H...O=C base to base contacts and various noncovalent links between the sugar hydroxyl and RNA base. Also, the intramolecular 1J(N,H) and 2J(C,H) couplings were calculated for the amino or imino group of RNA base and the ribose 2'-hydroxyl group involved in the noncovalent interactions. The calculated scalar couplings have implications for validation of local geometry, show specificity for the amino and imino groups of RNA base involved in the linkage, and can be used for discrimination between the cis- and trans-WC/SE base pairs. The RNA base pairs within an isosteric subclass of the WC/SE binding patterns can be further sorted according to the scalar couplings calculated across different local noncovalent contacts. The effect of explicit water inserted in the RNA base pairs on the magnitude of the scalar couplings was calculated, and the data for discrimination between the water-inserted and direct RNA base pairs are presented. The calculated NMR data are significant for structural interpretation of the scalar couplings in the noncanonical RNA base pairs.     

Structural rationalization of a large difference in RNA affinity despite a small difference in chemistry between two 2'-O-modified nucleic acid analogues

Chemical modification of nucleic acids at the 2'-position of ribose has generated antisense oligonucleotides (AONs) with a range of desirable properties. Electron-withdrawing substituents such as 2'-O-[2-(methoxy)ethyl] (MOE) confer enhanced RNA affinity relative to that of DNA by conformationally preorganizing an AON for pairing with the RNA target and by improving backbone hydration. 2'-Substitution of the ribose has also been shown to increase nuclease resistance and cellular uptake via changes in lipophilicity. Interestingly, incorporation of either 2'-O-[2-(methylamino)-2-oxoethyl]- (NMA) or 2'-O-(N-methylcarbamate)-modified (NMC) residues into AONs has divergent effects on RNA affinity. Incorporation of 2'-O-NMA-T considerably improves RNA affinity while incorporation of 2'-O-NMC-T drastically reduces RNA affinity. Crystal structures at high resolution of A-form DNA duplexes containing either 2'-O-NMA-T or 2'-O-NMC-T shed light on the structural origins of the surprisingly large difference in stability given the relatively minor difference in chemistry between NMA and NMC. NMA substituents adopt an extended conformation and use either their carbonyl oxygen or amino nitrogen to trap water molecules between phosphate group and sugar. The conformational properties of NMA and the observed hydration patterns are reminiscent of those found in the structures of 2'-O-MOE-modified RNA. Conversely, the carbonyl oxygen of NMC and O2 of T are in close contact, providing evidence that an unfavorable electrostatic interaction and the absence of a stable water structure are the main reasons for the loss in thermodynamic stability as a result of incorporation of 2'-O-NMC-modified residues.     

A ribose sugar conformational switch in the LTR-retrotransposon Ty3 polypurine tract-containing RNA/DNA hybrid

To probe structural features of a polypurine tract (PPT) that mediate its specific recognition and processing, a model 20 bp RNA/DNA hybrid duplex, which includes the full PPT sequence of the Saccharomyces cerevisiae LTR-retrotransposon Ty3, has been investigated using solution NMR spectroscopy. While homonuclear NOESY and DQF-COSY analyses indicate that this PPT-containing RNA/DNA hybrid adopts an overall A-form-like helical geometry, an unexpected sugar pucker switch has been detected for the ribose at position +1, relative to the cleavage site, on the RNA strand. A model of the conformational changes induced by the A- to B-type sugar pucker switch shows a significant change in the backbone trajectory of the RNA strand, which alters the presentation of backbone phosphate and 2' hydroxyl groups 3' of this residue. This observation implies that part of the mechanism governing RNase H fidelity may be through distortion of the RNA/DNA helix one base ahead of the scissile bond.     

The role of 23S ribosomal RNA residue A2451 in peptide bond synthesis revealed by atomic mutagenesis

Peptide bond formation is a fundamental reaction in biology, catalyzed by the ribosomal peptidyl-transferase ribozyme. Although all active-site 23S ribosomal RNA nucleotides are universally conserved, atomic mutagenesis suggests that these nucleobases do not carry functional groups directly involved in peptide bond formation. Instead, a single ribose 2'-hydroxyl group at A2451 was identified to be of pivotal importance. Here, we altered the chemical characteristics by replacing its 2'-hydroxyl with selected functional groups and demonstrate that hydrogen donor capability is essential for transpeptidation. We propose that the A2451-2'-hydroxyl directly hydrogen bonds to the P-site tRNA-A76 ribose. This promotes an effective A76 ribose C2'-endo conformation to support amide synthesis via a proton shuttle mechanism. Simultaneously, the direct interaction of A2451 with A76 renders the intramolecular transesterification of the peptide from the 3'- to 2'-oxygen unfeasible, thus promoting effective peptide bond synthesis.     

13C-NMR quantification of proton exchange at LewisX hydroxyl groups in water

NMR-based analysis of glycans by directly observing hydroxyl protons has been difficult because of their inherently fast exchange with water. We observed hydroxyl proton exchanges in a LewisX-LewisX interaction by using deuterium isotope shifts on (13)C-NMR. This strategy is suitable for analyzing weak interactions by identifying involved protons.     

Multinuclear NMR investigations of the oxygen, water, and hydroxyl environments in sodium hexaniobate

Solid-state (1)H, (17)O MAS NMR, (1)H-(93)Nb TRAPDOR NMR, and (1)H double quantum 2D MAS NMR experiments were used to characterize the oxygen, water, and hydroxyl environments in the monoprotonated hexaniobate material, Na(7)[HNb(6)O(19)].15H(2)O. These solid-state NMR experiments demonstrate that the proton is located on the bridging oxygen of the [Nb(6)O(19)](8-) cluster. The solid-state NMR results also show that the NbOH protons are spatially isolated from similar protons, but undergo proton exchange with the water species located in the crystal lattice. On the basis of double quantum (1)H MAS NMR measurements, it was determined that the water species in the crystal lattice have restricted motional dynamics. Two-dimensional (1)H-(17)O MAS NMR correlation experiments show that these restricted waters are preferentially associated with the bridging oxygen. Solution (17)O NMR experiments show that the hydroxyl proton is also attached to the bridging oxygen for the compound in solution. In addition, solution (17)O NMR kinetic studies for the hexaniobate allowed the measurement of relative oxygen exchange rates between the bridging, terminal, and hydroxyl oxygen and the oxygen of the solvent as a function of pH and temperature. These NMR experiments are some of the first investigations into the proton location, oxygen and proton exchange processes, and water dynamics for a base stable polyoxoniobate material, and they provide insight into the chemistry and reactivity of these materials.     

Computer simulation of the chemical catalysis of DNA polymerases: discriminating between alternative nucleotide insertion mechanisms for T7 DNA polymerase

Understanding the chemical step in the catalytic reaction of DNA polymerases is essential for elucidating the molecular basis of the fidelity of DNA replication. The present work evaluates the free energy surface for the nucleotide transfer reaction of T7 polymerase by free energy perturbation/empirical valence bond (FEP/EVB) calculations. A key aspect of the enzyme simulation is a comparison of enzymatic free energy profiles with the corresponding reference reactions in water using the same computational methodology, thereby enabling a quantitative estimate for the free energy of the nucleotide insertion reaction. The reaction is driven by the FEP/EVB methodology between valence bond structures representing the reactant, pentacovalent intermediate, and the product states. This pathway corresponds to three microscopic chemical steps, deprotonation of the attacking group, a nucleophilic attack on the P(alpha) atom of the dNTP substrate, and departure of the leaving group. Three different mechanisms for the first microscopic step, the generation of the RO(-) nucleophile from the 3'-OH hydroxyl of the primer, are examined: (i) proton transfer to the bulk solvent, (ii) proton transfer to one of the ionic oxygens of the P(alpha) phosphate group, and (iii) proton transfer to the ionized Asp654 residue. The most favorable reaction mechanism in T7 pol is predicted to involve the proton transfer to Asp654. This finding sheds light on the long standing issue of the actual role of conserved aspartates. The structural preorganization that helps to catalyze the reaction is also considered and analyzed. The overall calculated mechanism consists of three subsequent steps with a similar activation free energy of about 12 kcal/mol. The similarity of the activation barriers of the three microscopic chemical steps indicates that the T7 polymerase may select against the incorrect dNTP substrate by raising any of these barriers. The relative height of these barriers comparing right and wrong dNTP substrates should therefore be a primary focus of future computational studies of the fidelity of DNA polymerases.     

Density functional study of hydrogen bond formation between methanol and organic molecules containing Cl, F, NH2, OH, and COOH functional groups

Various hydrogen-bonded complexes of methanol with different proton accepting and proton donating molecules containing Cl, F, NH(2), OH, OR, and COOH functional groups have been modeled using DFT with hybrid B3LYP and M05-2X functionals. The latter functional was found to provide more accurate estimates of the structural and thermodynamic parameters of the complexes of halides, amines, and alcohols. The characteristics of these complexes are influenced not only by the principle hydrogen bond of the methanol OH with the proton acceptor heteroatom, but also by additional hydrogen bonds of a C-H moiety with methanol oxygen as a proton acceptor. The contribution of the former hydrogen bond in the total binding enthalpy increases in the order chlorides < fluorides < alcohols < amines, while the contribution of the second type of hydrogen bond increases in the reverse order. A general correlation was found between the binding enthalpy of the complex and the electrostatic potential at the hydrogen center participating in the formation of the hydrogen bond. The calculated binding enthalpies of different complexes were used to clarify which functional groups can potentially form a hydrogen bond to the 2'-OH hydroxyl group in ribose, which is strong enough to block it from participation in the intramolecular catalytic activation of the peptide bond synthesis. Such blocking could result in inhibition of the protein biosynthesis in the living cell if the corresponding group is delivered as a part of a drug molecule in the vicinity of the active site in the ribosome. According to our results, such activity can be accomplished by secondary or tertiary amines, alkoxy groups, deprotonated carboxyl groups, and aliphatic fluorides, but not by the other modeled functional groups.