PubMed Health⌕ Search

Biomedical subjects

Darrin M York

Publications and source records attributed to Darrin M York.

At least 19 recordsLinked to original sources
PubMed2007

CHARMM force field parameters for simulation of reactive intermediates in native and thio-substituted ribozymes.

Force field parameters specifically optimized for residues important in the study of RNA catalysis are derived from density-functional calculations, in a fashion consistent with the CHARMM27 all-atom empirical force field. Parameters are presented for residues that model reactive RNA intermediates and transition state analogs, thio-substituted phosphates and phosphoranes, and bound Mg(2+) and di-metal bridge complexes. Target data was generated via density-functional calculations at the B3LYP/6-311++G(3df,2p)// B3LYP/6-31++G(d,p) level. Partial atomic charges were initially derived from CHelpG electrostatic potential fitting and subsequently adjusted to be consistent with the CHARMM27 charges. Lennard-Jones parameters were determined to reproduce interaction energies with water molecules. Bond, angle, and torsion parameters were derived from the density-functional calculations and renormalized to maintain compatibility with the existing CHARMM27 parameters for standard residues. The extension of the CHARMM27 force field parameters for the nonstandard biological residues presented here will have considerable use in simulations of ribozymes, including the study of freeze-trapped catalytic intermediates, metal ion binding and occupation, and thio effects.

Catalysis↗
PubMed2006

Enzyme dynamics and tunneling enhanced by compression in the hydrogen abstraction catalyzed by soybean lipoxygenase-1.

A fully microscopical simulation of the rate-limiting hydrogen abstraction catalyzed by soybean lipoxygenase-1 (SLO-1) has been carried out. This enzyme exhibits the largest, and weakly temperature dependent, experimental H/D kinetic isotope effect (KIE) reported for a biological system. The theoretical model used here includes the complete enzyme with a solvation shell of water molecules, the Fe(III)-OH- cofactor, and the linoleic acid substrate. We have used a hybrid QM(PM3/d-SRP)/MM method to describe the potential energy surface of the whole system, and the ensemble-averaged variational transition-state theory with multidimensional tunneling (EA-VTST/MT) to calculate the rate constant and the primary KIE. The computational results show that the compression of the wild-type active site enzyme results in the huge contribution of tunneling (99%) to the rate of the hydrogen abstraction. Importantly, the active site becomes more flexible in the Ile553Ala mutant reactant complex simulation (for which a markedly temperature dependent KIE has been experimentally determined), thus justifying the proposed key role of the gating promoting mode in the reaction catalyzed by SLO-1. Finally, the results indicate that the calculated KIE for the wild-type enzyme has an important dependence on the barrier width.

Binding Sites↗
PubMed2006

Transesterification thio effects of phosphate diesters: free energy barriers and kinetic and equilibrium isotope effects from density-functional theory.

Primary and secondary kinetic and equilibrium isotope effects are calculated with density-functional methods for the in-line dianionic methanolysis of the native (unsubstituted) and thio-substituted cyclic phosphates. These reactions represent reverse reaction models for RNA transesterification under alkaline conditions. The effect of solvent is treated with explicit (single and double) water molecules and self-consistently with an implicit (continuum) solvation model. Singly substituted reactions at the nonbridging O(P1) position and bridging O(2)('), O(3)('), and O(5)(') positions and a doubly substituted reaction at the O(P1) and O(P2) positions were considered. Aqueous free energy barriers are calculated, and the structures and bond orders of the rate-controlling transition states are characterized. The results are consistent with available experimental data and provide useful information for the interpretation of measured isotope and thio effects used to probe mechanism in phosphoryl transfer reactions catalyzed by enzymes and ribozymes.

Catalysis↗
PubMed2006

Nucleophilic attack on phosphate diesters: a density functional study of in-line reactivity in dianionic, monoanionic, and neutral systems.

A density functional study of the hydrolysis reaction of phosphodiesters with a series of attacking nucleophiles in the gas phase and in solution is presented. The nucleophiles HOH, HO-, CH3OH, and CH3O- were studied in reactions with ethylene phosphate, 2'3'-ribose cyclic phosphate and in their neutral (protonated) and monoanionic forms. Stationary-point geometries for the reactions were determined at the density functional B3LYP/6-31++G(d,p) level followed by energy refinement at the B3LYP/6-311++G(3df,2p) level. Solvation effects were estimated by using a dielectric approximation with the polarizable continuum model (PCM) at the gas-phase optimized geometries. This series of reactions characterizes factors that influence the intrinsic reactivity of the model phosphate compounds, including the effect of nucleophile, protonation state, cyclic structure, and solvent. The present study of the in-line mechanism for phosphodiester hydrolysis, a reaction of considerable biological importance, has implications for enzymatic mechanisms. The analysis generally supports the associative mechanism for phosphate ester hydrolysis. The results highlight the importance for the reaction barrier of charge neutralization resulting from the protonation of the nonbridging phosphoryl oxygens and the role of internal hydrogen transfer in the gas-phase mechanism. It also shows that solvent stabilization has a profound influence on the relative barrier heights for the dianionic, monoanionic, and neutral reactions. The calculations provide a comprehensive data set for the in-line hydrolysis mechanisms that can be used for the development of improved semiempirical quantum models for phosphate hydrolysis reactions.

Anions↗
PubMed2006

QCRNA 1.0: a database of quantum calculations for RNA catalysis.

This work outlines a new on-line database of quantum calculations for RNA catalysis (QCRNA) available via the worldwide web at http://theory.chem.umn.edu/QCRNA. The database contains high-level density functional calculations for a large range of molecules, complexes and chemical mechanisms important to phosphoryl transfer reactions and RNA catalysis. Calculations are performed using a strict, consistent protocol such that a wealth of cross-comparisons can be made to elucidate meaningful trends in biological phosphate reactivity. Currently, around 2000 molecules have been collected in varying charge states in the gas phase and in solution. Solvation was treated with both the PCM and COSMO continuum solvation models. The data can be used to study important trends in reactivity of biological phosphates, or used as benchmark data for the design of new semiempirical quantum models for hybrid quantum mechanical/molecular mechanical simulations.

Computer Graphics↗
PubMed2006

Multilevel and density functional electronic structure calculations of proton affinities and gas-phase basicities involved in biological phosphoryl transfer.

Five multilevel model chemistries (CBS-QB3, G3B3, G3MP2B3, MCG3/3, and MC-QCISD/3) and seven hybrid density functional methods (PBE0, B1B95, B3LYP, MPW1KCIS, PBE1KCIS, and MPW1B95) have been applied to the calculation of gas-phase basicity and proton affinity values for a series of 17 molecules relevant to the study of biological phosphoryl transfer. In addition, W1 calculations were performed on a subset of molecules. The accuracy of the methods was assessed and the nature of systematic errors was explored, leading to the introduction of a set of effective bond enthalpy and entropy correction terms. The multicoefficient correlation methods (MCG3/3 and MC-QCISD), with inclusion of specific zero-point scale factors, slightly outperform the other multilevel methods tested (CBS-QB3, G3B3, and G3MP2B3), with significantly less computational cost, and in the case of MC-QCISD, slightly less severe scaling. Four density functional methods, PBE1KCIS, MPW1B95, PBE0, and B1B95 perform nearly as well as the multilevel methods. These results provide an important set of benchmarks relevant to biological phosphoryl transfer reactions.

Electrons↗
PubMed2006

A charge-scaling implementation of the variational electrostatic projection method.

Two new charge-scaling methods for efficient modeling of the solvated macromolecular environment in hybrid QM/MM calculations of biological reactions are presented. The methods are extensions of the variational electrostatic projection (VEP) method, and allows a subset of atomic charges in the external environment to be adjusted to mimic, in the active dynamical region, the electrostatic potential and field due to the large surrounding macromolecule and solvent. The method has the advantages that it offers improved accuracy, does not require the use of a three-dimensional grid or auxiliary set of fitting points, and requires only minor molecular simulation code modifications. The VEP-cs and VEP-RVM+cs methods are able to attain very high accuracy (relative force errors of 10(-7) or better with appropriate choice of control parameters), and take advantage of a recently introduced set of high-order discretization schemes and Gaussian exponents for boundary element solvation and VEP methods. The methods developed here serve as potentially powerful tools in the arsenal of computational techniques used in multiscale computational modeling problems.

Biopolymers↗
PubMed2005

Solvent polarization and kinetic isotope effects in nitroethane deprotonation and implications to the nitroalkane oxidase reaction.

We have carried out a mixed molecular dynamics and centroid path integral simulation using a combined quantum mechanical and molecular mechanical (QM/MM) potential to study the anomalous Brønsted relationship between rates and equilibria for deprotonation of nitroalkanes in water, which is known as the nitroalkane anomaly. The deprotonation process is catalyzed by nitroalkane oxidase. Our results show that the difference in solvent polarization effects for the TS and products is a major factor for the differential solvent effects on rate and equilibrium of nitroalkane deprotonation. This is due to poor charge delocalization as a result of slow rehybridization compared to bond breaking. Although solvent effects do not affect significantly the computed kinetic isotope effects in comparison with the gas-phase value, there is slight solvent-induced increase in tunneling. The present results suggest that an effective means by which the transition state can be stabilized in the enzyme nitroalkane oxidase is to facilitate the Calpha rehybridization.

Computer Simulation↗
PubMed2005

Density functional study of the in-line mechanism of methanolysis of cyclic phosphate and thiophosphate esters in solution: insight into thio effects in RNA transesterification.

Density functional calculations of thio effects on the in-line mechanism of methanolysis of ethylene phosphate (a reverse reaction model for RNA phosphate transesterification) are presented. A total of 12 reaction mechanisms are examined using the B3LYP functional with large basis sets, and the effects of solvation were treated using the PCM, CPCM, and SM5 solvation models. Single thio substitutions at all of the distinct phosphoryl oxygen positions (2', 3', 5', pro-R) and a double thio substitution at the nonbridging (pro-R/pro-S) positions were considered. Profiles for each reaction were calculated in the dianionic and monoanionic/monoprotic states, corresponding to reaction models under alkaline and nonalkaline conditions, respectively. These models provide insight into the mechanisms of RNA transesterification thio effects and serve as a set of high-level quantum data that can be used in the design of new semiempirical quantum models for hybrid quantum mechanical/molecular mechanical simulations and linear-scaling electronic structure calculations.

Cyclization↗
PubMed2005

Improvement of semiempirical response properties with charge-dependent response density.

The present work outlines a new method for treatment of charge-dependent polarizability in semiempirical quantum models for use in combined quantum-mechanical/molecular mechanical simulations of biological reactions. The method addresses a major shortcoming in the performance of conventional semiempirical models for these simulations that is tied to the use of a localized minimal atomic-orbital basis set. The present approach has the advantages that it uses a density basis that retains a set of linear-response equations, does not increase the atomic-orbital basis, and avoids the problem of artificial charge transfer and scaling of the polarizability seen in related models that allow atomic charges to fluctuate. The model introduces four new atom-based parameters and has been tested with the modified neglect of differential overlap d-orbital Hamiltonian against 1132 molecules and ions and shown to decrease the dipole moment and polarizability errors by factors of 2 and 10, respectively, with respect to density-functional results. The method performs impressively for a variety of charge states (from 2+ to 2-), and offers a potentially powerful extension in the design of next generation semiempirical quantum models for accurate simulations of highly charged biological reactions.

Anions↗
PubMed2005

Theoretical studies of dissociative phosphoryl transfer in interconversion of phosphoenolpyruvate to phosphonopyruvate: solvent effects, thio effects, and implications for enzymatic reactions.

The conversion of phosphoenolpyruvate (PEP) to phosphonopyruvate (P-pyr) is catalyzed by PEP mutase via a dissociative mechanism. In this work, we investigate the uncatalyzed reaction using ab initio methods, density functional theory, and the semiempirical MNDO/d method. Comparisons of geometries and relative energies of stationary points (minima and transition states) with density functional results indicate that the semiempirical method is reasonably accurate. Solvent effects are examined using implicit solvent models, including the recently extended smooth conductor-like screening model. Due to the large negative charge carried by the system, solvation is found to drastically alter the location and energy of stationary points along the dissociative reaction pathways. The influence of substituting a nonbridging phosphoryl oxygen by sulfur (thio effects) was also investigated. Implications of these results for the enzymatic reaction are discussed.

Calibration↗
PubMed2005

Benchmark calculations of proton affinities and gas-phase basicities of molecules important in the study of biological phosphoryl transfer.

Benchmark calculations of proton affinities and gas-phase basicities of molecules most relevant to biological phosphoryl transfer reactions are presented and compared with available experimental results. The accuracy of proton affinity and gas-phase basicity results obtained from several multi-level model chemistries (CBS-QB3, G3B3, and G3MP2B3) and density-functional quantum models (PBE0, B1B95, and B3LYP) are assessed and compared. From these data, a set of empirical bond enthalpy, entropy, and free energy corrections are introduced that considerably improve the accuracy and predictive capability of the methods. These corrections are applied to the prediction of proton affinity and gas-phase basicity values of important biological phosphates and phosphoranes for which experimental data does not currently exist. Comparison is made with results from semiempirical quantum models that are commonly employed in hybrid quantum mechanical/molecular mechanical simulations. Data suggest that the design of improved semiempirical quantum models with increased accuracy for relative proton affinity values is necessary to obtain quantitative accuracy for phosphoryl transfer reactions in solution, enzymes, and ribozymes.

Algorithms↗
PubMed2005

Smooth solvation method for d-orbital semiempirical calculations of biological reactions. 2. Application to transphosphorylation thio effects in solution.

Density-functional and semiempirical quantum methods and continuum dielectric and explicit solvation models are applied to study the role of solvation on the stabilization of native and thio-substituted transphosphorylation reactions. Extensive comparison is made between results obtained from the different methods. For the semiempirical methods, explicit solvation was treated using a hybrid quantum mechanical/molecular mechanical (QM/MM) approach and the implicit solvation was treated using a recently developed smooth solvation model implemented into a d-orbital semiempirical framework (MNDO/d-SCOSMO) within CHARMM. The different quantum and solvation methods were applied to the transesterification of a 3'-ribose,5'-methyl phosphodiester that serves as a nonenzymatic model for the self-cleavage reaction catalyzed by the hammerhead and hairpin ribozymes. Thio effects were studied for a double sulfur substitution at the nonbridging phosphoryl oxygen positions. The reaction profiles of both the native and double sulfur-substituted reactions from the MNDO/d-SCOSMO calculations were similar to the QM/MM results and consistent with the experimentally observed trends. These results underscore the need for a d-orbital semiempirical representation for phosphorus and sulfur for the study of experimentally observed thio effects in enzymatic and nonenzymatic phosphoryl transfer reactions. One of the major advantages of the present approach is that it can be applied to model chemical reactions at a significantly lower computational cost than either the density-functional calculations with implicit solvation or the semiempirical QM/MM simulations with explicit solvent.

Algorithms↗
PubMed2005

Smooth solvation method for d-orbital semiempirical calculations of biological reactions. 1. Implementation.

The present paper describes the extension of a recently developed smooth conductor-like screening model for solvation to a d-orbital semiempirical framework (MNDO/d-SCOSMO) with analytic gradients that can be used for geometry optimizations, transition state searches, and molecular dynamics simulations. The methodology is tested on the potential energy surfaces for separating ions and the dissociative phosphoryl transfer mechanism of methyl phosphate. The convergence behavior of the smooth COSMO method with respect to discretization level is examined and the numerical stability of the energy and gradient are compared to that from conventional COSMO calculations. The present method is further tested in applications to energy minimum and transition state geometry optimizations of neutral and charged metaphosphates, phosphates, and phosphoranes that are models for stationary points in transphosphorylation reaction pathways of enzymes and ribozymes. The results indicate that the smooth COSMO method greatly enhances the stability of quantum mechanical geometry optimization and transition state search calculations that would routinely fail with conventional solvation methods. The present MNDO/d-SCOSMO method has considerable computational advantages over hybrid quantum mechanical/molecular mechanical methods with explicit solvation, and represents a potentially useful tool in the arsenal of multi-scale quantum models used to study biochemical reactions.

Algorithms↗
PubMed2005

High-order discretization schemes for biochemical applications of boundary element solvation and variational electrostatic projection methods.

A series of high-order surface element discretization schemes for variational boundary element methods are introduced. The surface elements are chosen in accord with angular quadrature rules for integration of spherical harmonics. Surface element interactions are modeled by Coulomb integrals between spherical Gaussian functions with exponents chosen to reproduce the exact variational energy and Gauss's law for a point charge in a spherical cavity. The present work allows high-order surface element expansions to be made for variational methods such as the conductorlike screening model for solvation and the variational electrostatic projection method for generalized solvent boundary potentials in molecular simulations.

Algorithms↗
PubMed2005

Pseudorotation barriers of biological oxyphosphoranes: a challenge for simulations of ribozyme catalysis.

Pseudorotation reactions of biologically relevant oxyphosphoranes were studied by using density functional and continuum solvation methods. A series of 16 pseudorotation reactions involving acyclic and cyclic oxyphosphoranes in neutral and monoanionic (singly deprotonated) forms were studied, in addition to pseudorotation of PF5. The effect of solvent was treated by using three different solvation models for comparison. The barriers to pseudorotation ranged from 1.5 to 8.1 kcal mol(-1) and were influenced systematically by charge state, apicophilicity of ligands, intramolecular hydrogen bonding, cyclic structure and solvation. Barriers to pseudorotation for monoanionic phosphoranes occur with the anionic oxo ligand as the pivotal atom, and are generally lower than for neutral phosphoranes. The OCH3 groups were observed to be more apicophilic than OH groups, and hence pseudorotations that involve axial OCH3/equatorial OH exchange had higher reaction and activation free energy values. Solvent generally lowered barriers relative to the gas-phase reactions. These results, together with isotope 18O exchange experiments, support the assertion that dianionic phosphoranes are not sufficiently long-lived to undergo pseudorotation. Comparison of the density functional results with those from several semiempirical quantum models highlight a challenge for new-generation hybrid quantum mechanical/molecular mechanical potentials for non-enzymatic and enzymatic phosphoryl transfer reactions: the reliable modeling of pseudorotation processes.

Catalysis↗
PubMed2005

Ellipticity: a convenient tool to characterize electrocyclic reactions.

The ellipticity of the electron density is proposed as a convenient tool to characterize electrocyclic reactions. The study of the electron density offers several advantages since it is an experimental observable. This simple topological index has proven powerful and robust for the characterization of a wide variety of electrocyclic and pseudoelectrocyclic processes. This property is sensitive to changes in the anisotropy of the electron density in the bond-forming region and provides an insightful description of the events occurring along the reaction coordinate.

Cyclization↗