Biomedical subjects

Lee G Pedersen

Publications and source records attributed to Lee G Pedersen.

12 recordsLinked to original sources

PubMed2006

A quantum chemical study of the mechanism of action of Vitamin K epoxide reductase (VKOR) II. Transition states.

A reaction path including transition states is generated for the Silverman mechanism [R.B. Silverman, Chemical model studies for the mechanism of Vitamin K epoxide reductase, J. Am. Chem. Soc. 103 (1981) 5939-5941] of action for Vitamin K epoxide reductase (VKOR) using quantum mechanical methods (B3LYP/6-311G**). VKOR, an essential enzyme in mammalian systems, acts to convert Vitamin K epoxide, formed by Vitamin K carboxylase, to its (initial) quinone form for cellular reuse. This study elaborates on a prior work that focused on the thermodynamics of VKOR [D.W. Deerfield II, C.H. Davis, T. Wymore, D.W. Stafford, L.G. Pedersen, Int. J. Quant. Chem. 106 (2006) 2944-2952]. The geometries of proposed model intermediates and transition states in the mechanism are energy optimized. We find that once a key disulfide bond is broken, the reaction proceeds largely downhill. An important step in the conversion of the epoxide back to the quinone form involves initial protonation of the epoxide oxygen. We find that the source of this proton is likely a free mercapto group rather than a water molecule. The results are consistent with the current view that the widely used drug Warfarin likely acts by blocking binding of Vitamin K at the VKOR active site and thereby effectively blocking the initiating step. These results will be useful for designing more complete QM/MM studies of the enzymatic pathway once three-dimensional structural data is determined and available for VKOR.

Epoxy Compounds↗

PubMed2006

A quantum chemical study of the mechanism of action of Vitamin K carboxylase (VKC) III. Intermediates and transition states.

A reaction path including transition states is generated for the Dowd mechanism [P. Dowd, R. Hershlne, S.W. Ham, S. Naganathan. Vitamin K and energy transduction: a base strength amplification mechanism. Science 269 (2005) 1684-1691] of action for Vitamin K carboxylase (VKC) using quantum chemical methods (B3LYP/6-311G**). VKC, an essential enzyme in mammalian systems, catalyzes the conversion of hydroquinone form of Vitamin K to the epoxide form in the presence of oxygen. An intermediate species of the oxidation of Vitamin K, an alkoxide, acts apparently to abstract the gamma hydrogen from specifically located glutamate residues. We are able to follow the Dowd proposed path to generate this alkoxide species. The geometries of the proposed model intermediates and transition states in the mechanism are energy optimized. We find that the most energetic step in the mechanism is the uni-deprotonation of the hydroquinone - once this occurs, there is only a small barrier of 3.5kcal/mol for the interaction of oxygen with the carbon to be attacked - and then the reaction proceeds downhill in free energy to form the critical alkoxide species. The results are consistent with the idea that the enzyme probably acts to facilitate the formation of the epoxide by reducing the energy required to deprotonate the hydroquinone form.

Carbon-Carbon Ligases↗

PubMed2006

Energy analysis of chemistry for correct insertion by DNA polymerase beta.

X-ray crystallographic structures of human DNA polymerase beta with nonhydrolyzable analogs containing all atoms in the active site required for catalysis provide a secure starting point for a theoretical analysis (quantum mechanics/molecular mechanics) of the mechanism of chemistry without biasing of modeling assumptions as required in previous studies. These structures provide the basis for a detailed quantum mechanics/molecular mechanics study of the path for the complete transfer of a monophosphate nucleoside donor to the sugar acceptor in the active site. The reaction is largely associative with the main energetic step preceded by proton transfer from the terminal primer deoxyribose O3' to Asp-256. The key residues that provide electrostatic stabilization of the transition state are identified and compared with those identified by mutational studies.

Binding Sites↗

PubMed2006

Towards accurate solvation dynamics of divalent cations in water using the polarizable amoeba force field: From energetics to structure.

Molecular dynamics simulations were performed using a modified amoeba force field to determine hydration and dynamical properties of the divalent cations Ca2+ and Mg2+. The extension of amoeba to divalent cations required the introduction of a cation specific parametrization. To accomplish this, the Thole polarization damping model parametrization was modified based on the ab initio polarization energy computed by a constrained space orbital variation energy decomposition scheme. Excellent agreement has been found with condensed phase experimental results using parameters derived from gas phase ab initio calculations. Additionally, we have observed that the coordination of the calcium cation is influenced by the size of the periodic water box, a recurrent issue in first principles molecular dynamics studies.

Amoeba↗

PubMed2005

Deuterium and its role in the machinery of evolution.

An argument is presented that the spontaneous mutation rate, the core of evolution theory, may be dictated by the deuterium/hydrogen (D/H) abundance ratio. This argument is supported by quantum mechanical calculations of the zero-point energy reduction for DNA base pairs upon deuterium substitution for hydrogen and recent experiments that show that the rate of catalytic dsDNA unwinding is dependent on the stability of the dsDNA.

Biological Evolution↗

PubMed2004

Towards an accurate representation of electrostatics in classical force fields: efficient implementation of multipolar interactions in biomolecular simulations.

The accurate simulation of biologically active macromolecules faces serious limitations that originate in the treatment of electrostatics in the empirical force fields. The current use of "partial charges" is a significant source of errors, since these vary widely with different conformations. By contrast, the molecular electrostatic potential (MEP) obtained through the use of a distributed multipole moment description, has been shown to converge to the quantum MEP outside the van der Waals surface, when higher order multipoles are used. However, in spite of the considerable improvement to the representation of the electronic cloud, higher order multipoles are not part of current classical biomolecular force fields due to the excessive computational cost. In this paper we present an efficient formalism for the treatment of higher order multipoles in Cartesian tensor formalism. The Ewald "direct sum" is evaluated through a McMurchie-Davidson formalism [L. McMurchie and E. Davidson, J. Comput. Phys. 26, 218 (1978)]. The "reciprocal sum" has been implemented in three different ways: using an Ewald scheme, a particle mesh Ewald (PME) method, and a multigrid-based approach. We find that even though the use of the McMurchie-Davidson formalism considerably reduces the cost of the calculation with respect to the standard matrix implementation of multipole interactions, the calculation in direct space remains expensive. When most of the calculation is moved to reciprocal space via the PME method, the cost of a calculation where all multipolar interactions (up to hexadecapole-hexadecapole) are included is only about 8.5 times more expensive than a regular AMBER 7 [D. A. Pearlman et al., Comput. Phys. Commun. 91, 1 (1995)] implementation with only charge-charge interactions. The multigrid implementation is slower but shows very promising results for parallelization. It provides a natural way to interface with continuous, Gaussian-based electrostatics in the future. It is hoped that this new formalism will facilitate the systematic implementation of higher order multipoles in classical biomolecular force fields.

Algorithms↗

PubMed2003

Glucosaminylglycan biosynthesis: what we can learn from the X-ray crystal structures of glycosyltransferases GlcAT1 and EXTL2.

The X-ray crystal structures of two glycosyltransferases (GTs)--beta 1,3-glucuronyltransferase 1 (GlcAT1) and alpha 1,4-N-acetylhexosaminyltransferase (EXTL2)--have now been determined in the presence of both donor and acceptor substrates. These enzymes are involved in glucosaminylglycan (GAG) synthesis where they catalyze inverting and retaining transfer reactions, respectively. As members of a large family of enzymes that transfer sugar groups from donor nucleotide-sugars to acceptor substrates, GlcAT1 and EXTL2 retain conserved GT folds. Comparative analysis of these structures reveals signature features for selecting specific donor sugars. Adaptive binding of the disaccharide moiety of the acceptor sugars enables the enzymes to catalyze either an inverting S(N)2-type displacement reaction or a retaining S(N)i-like transfer reaction.

Amino Acid Sequence↗

PubMed2003

Crystal structure of an alpha 1,4-N-acetylhexosaminyltransferase (EXTL2), a member of the exostosin gene family involved in heparan sulfate biosynthesis.

EXTL2, an alpha1,4-N-acetylhexosaminyltransferase, catalyzes the transfer reaction of N-acetylglucosamine and N-acetylgalactosamine from the respective UDP-sugars to the non-reducing end of [glucuronic acid]beta1-3[galactose]beta1-O-naphthalenemethanol, an acceptor substrate analog of the natural common linker of various glycosylaminoglycans. We have solved the x-ray crystal structure of the catalytic domain of mouse EXTL2 in the apo-form and with donor substrates UDP-N-acetylglucosamine and UDP-N-acetylgalactosamine. In addition, a structure of the ternary complex with UDP and the acceptor substrate analog [glucuronic acid]beta1-3[galactose]beta1-O-naphthalenemethanol has been determined. These structures reveal three highly conserved residues, Asn-243, Asp-246, and Arg-293, located at the active site. Mutation of these residues greatly decreases the activity. In the ternary complex, an interaction exists between the beta-phosphate of the UDP leaving group and the acceptor hydroxyl of the substrate that may play a functional role in catalysis. These structures represent the first structures from the exostosin gene family and provide important insight into the mechanisms of alpha1,4-N-acetylhexosaminyl transfer in heparan biosynthesis.

Amino Acid Sequence↗

PubMed2003

Explicit water near the catalytic I helix Thr in the predicted solution structure of CYP2A4.

The solution structure of mouse cytochrome P450 2A4 (CYP2A4), a monooxygenase of deoxysteroids, was obtained using homology modeling and molecular dynamics. The solvent-equilibrated CYP2A4 preserves the essential features of CYP450s. A comparison of the models CYP2A4 and CYP2A4 with testosterone bound CYP2A4/T illustrates the changes induced by the binding of the substrate. Experimental evidence links four amino acid residues to the catalytic activity, substrate specificity, and regioselectivity of this enzyme. Three of the four amino acids are found within contact distance of the testosterone substrate, and therefore may control the binding of the substrate through direct interaction. Remarkably, a water complex previously observed in x-ray crystal structure forms near the bulge in the central I helix that contains a conserved Thr. The properties of the I helix are computed in the context of the presence or absence of ligand.

Amino Acid Sequence↗

PubMed2002

Three-dimensional solution structure of Tropidechis carinatus venom extract trocarin: a structural homologue of Xa and prothrombin activator.

Trocarin belongs to group D of prothrombin activators derived from snake venom of Tropidechis carinatus and is a rich non-hepatic source of Xa, the only known hepatic prothrombin activator. The structural and functional similarity with Xa makes trocarin an interesting target for exploring the structure-functional relationship with Xa. Herein we report a predicted complete three-dimensional all-atom structural model of trocarin equilibrated in explicit water using 4 ns of molecular dynamics simulation. The tertiary structure was modeled using the structure of human blood coagulation factor Xa. The conformational and structural features of trocarin are then compared with the X-ray crystal and solution simulation structures of human factor Xa. The modeled structure of trocarin has four individual domains (Gla, EGF1, EGF2 and SP) connected along the long axis with similar secondary structural elements to Xa. The simulations suggest that sodium ion binding in the serine protease domain is impaired in trocarin as compared to Xa. In contrast to Xa, for which the sodium ion forms an octahedral coordination network that brings two loop regions connecting four anti-parallel beta-sheets together, we do not find a similar pattern of network in trocarin. We observe that the difference in the binding pattern of sodium ion leads to a approximately 2-A "shrinkage" of the beta2 strand (B2), in comparison to human Xa, as marked by a shorter distance between 189Asp373 (S1-site residue) and 195Ser379 (active-site residue) in the B2 strand. We propose that these differences may be linked to the experimentally observed lower amidolytic activity of trocarin as compared to Xa.

Amino Acid Sequence↗

PubMed2002

Predicted solution structure of zymogen human coagulation FVII.

A model solution structure for the complete tissue factor-free calcium ion-bound human zymogen FVII (residues 1-406) (FVII) has been constructed to study possible conformational changes associated with the activation process and tissue factor (TF) binding. The initial structure for the present model was constructed using the X-ray crystallographic structure of human coagulation FVIIa/TF complex bound with calcium ions (Banner et al., Nature 1996, 380, 41-46). This model was subsequently subjected to lengthy molecular dynamics simulations. The Amber force field in conjunction with the PME electrostatic summation method was employed. The estimated TF free solution structure was then compared with the currently available X-ray crystal structures of FVIIa (with or without TF, variable inhibitor bound) to estimate the restructuring of FVII due to TF binding and activation. The solution structure of the zymogen FVII in the absence of TF is predicted to be an extended domain structure similar to that of the TF-bound X-ray crystal structure. An additional extension of the serine protease (SP) domain of the zymogen above a reference lipid surface by approximately 7 A was in agreement with experiment. Significant Gla-EGF1 and EGF1-EGF2 interdomain motions in the zymogen were observed. Carbohydrate dimers attached to Ser-52 and Ser-60 did not cause restructuring in this domain. Minimal restructuring of the SP domain is found upon inference of the zymogen from the activated form. The catalytic triad residues maintain the H-bonded network while Lys-341 occupies the S1 specific site in the zymogen.

Binding Sites↗

PubMed2002

Structure and dynamics of zymogen human blood coagulation factor X.

The solution structure and dynamics of the human coagulation factor X (FX) have been investigated to understand the key structural elements in the zymogenic form that participates in the activation process. The model was constructed based on the 2.3-A-resolution x-ray crystallographic structure of active-site inhibited human FXa (PDB:1XKA). The missing gamma-carboxyglutamic acid (GLA) and part of epidermal growth factor 1 (EGF1) domains of the light chain were modeled based on the template of GLA-EGF1 domains of the tissue factor (TF)-bound FVIIa structure (PDB:1DAN). The activation peptide and other missing segments of FX were introduced using homology modeling. The full calcium-bound model of FX was subjected to 6.2 ns of molecular dynamics simulation in aqueous medium using the AMBER6.0 package. We observed significant reorientation of the serine-protease (SP) domain upon activation leading to a compact multi-domain structure. The solution structure of zymogen appears to be in a well-extended conformation with the distance between the calcium ions in the GLA domain and the catalytic residues estimated to be approximately 95 A in contrast to approximately 83 A in the activated form. The latter is in close agreement with fluorescence studies on FXa. The S1-specificity residues near the catalytic triad show significant differences between the zymogen and activated structures.

1-Carboxyglutamic Acid↗