From Steroid and Drug Metabolism to Glycobiology, Using Sulfotransferase Structures to Understand and Tailor Function

Sulfotransferases are ubiquitous enzymes that transfer a sulfo group from the universal cofactor donor 3'-phosphoadenosine 5'-phosphosulfate to a broad range of acceptor substrates. In humans, the cytosolic sulfotransferases are involved in the sulfation of endogenous compounds such as steroids, neurotransmitters, hormones, and bile acids as well as xenobiotics including drugs, toxins, and environmental chemicals. The Golgi associated membrane-bound sulfotransferases are involved in post-translational modification of macromolecules from glycosaminoglycans to proteins. The sulfation of small molecules can have profound biologic effects on the functionality of the acceptor, including activation, deactivation, or enhanced metabolism and elimination. Sulfation of macromolecules has been shown to regulate a number of physiologic and pathophysiological pathways by enhancing binding affinity to regulatory proteins or binding partners. Over the last 25 years, crystal structures of these enzymes have provided a wealth of information on the mechanisms of this process and the specificity of these enzymes. This review will focus on the general commonalities of the sulfotransferases, from enzyme structure to catalytic mechanism as well as providing examples into how structural information is being used to either design drugs that inhibit sulfotransferases or to modify the enzymes to improve drug synthesis. SIGNIFICANCE STATEMENT: This manuscript honors Dr. Masahiko Negishi's contribution to the understanding of sulfotransferase mechanism, specificity, and roles in biology by analyzing the crystal structures that have been solved over the last 25 years.

Comprehensive Test of the Brink-Axel Hypothesis in the Energy Region of the Pygmy Dipole Resonance

The validity of the Brink-Axel hypothesis, which is especially important for numerous astrophysical calculations, is addressed for ^{116,120,124}Sn below the neutron separation energy by means of three independent experimental methods. The γ-ray strength functions (GSFs) extracted from primary γ-decay spectra following charged-particle reactions with the Oslo method and with the shape method demonstrate excellent agreement with those deduced from forward-angle inelastic proton scattering at relativistic beam energies. In addition, the GSFs are shown to be independent of excitation energies and spins of the initial and final states. The results provide a critical test of the generalized Brink-Axel hypothesis in heavy nuclei, demonstrating its applicability in the energy region of the pygmy dipole resonance.

Preferential DNA Polymerase β Reverse Reaction with Imidodiphosphate

DNA replication and repair reactions involve the addition of a deoxynucleoside monophosphate onto a growing DNA strand with the loss of pyrophosphate. This chemical reaction is also reversible; the addition of pyrophosphate generates a deoxynucleoside triphosphate, thereby shortening the DNA by one nucleotide. The forward DNA synthesis and reverse pyrophosphorolysis reactions strictly require the presence of divalent metals, usually magnesium, at the reactive center as cofactors. The overall equilibrium enzymatic reaction strongly favors DNA synthesis over pyrophosphorolysis with natural substrates. The DNA polymerase β chemical reaction has been structurally and kinetically characterized, employing natural and chemically modified substrates. Substituting an imido-moiety (NH) for the bridging oxygen between Pβ and Pγ of dGTP dramatically decreased the overall enzymatic activity and resulted in a chemical equilibrium that strongly favors the reverse reaction (i.e., K ≪ 1). Using QM/MM calculations in conjunction with the utilization of parameters such as quantum mechanically derived atomic charges, we have examined the chemical foundation for the altered equilibrium with this central biological reaction. The calculations indicate that the rapid reverse reaction is likely due, in part, to the increased nucleophilicity of the reactive oxygen on the tautomeric form of imidodiphosphate.

Evolution of Octupole Deformation in Radium Nuclei from Coulomb Excitation of Radioactive ^{222}Ra and ^{228}Ra Beams

There is sparse direct experimental evidence that atomic nuclei can exhibit stable "pear" shapes arising from strong octupole correlations. In order to investigate the nature of octupole collectivity in radium isotopes, electric octupole (E3) matrix elements have been determined for transitions in ^{222,228}Ra nuclei using the method of sub-barrier, multistep Coulomb excitation. Beams of the radioactive radium isotopes were provided by the HIE-ISOLDE facility at CERN. The observed pattern of E3 matrix elements for different nuclear transitions is explained by describing ^{222}Ra as pear shaped with stable octupole deformation, while ^{228}Ra behaves like an octupole vibrator.

The observation of vibrating pear-shapes in radon nuclei

There is a large body of evidence that atomic nuclei can undergo octupole distortion and assume the shape of a pear. This phenomenon is important for measurements of electric-dipole moments of atoms, which would indicate CP violation and hence probe physics beyond the Standard Model of particle physics. Isotopes of both radon and radium have been identified as candidates for such measurements. Here, we observed the low-lying quantum states in ²²⁴Rn and ²²⁶Rn by accelerating beams of these radioactive nuclei. We show that radon isotopes undergo octupole vibrations but do not possess static pear-shapes in their ground states. We conclude that radon atoms provide less favourable conditions for the enhancement of a measurable atomic electric-dipole moment.

A Kinetic Model Describing the Interaction of Bovine Prothrombin Fragment 1 with Calcium Ions

A kinetic model is derived for the interaction of bovine prothrombin fragment 1 with calcium ions. The model requires binding of a minimum of two calcium ions for induction of the observed biphasic fluorescence decrease as a function of time. The model is shown to be consistent with experimental kinetic and equilibrium data by fitting theoretical curves for the biphasic fluorescence change to the data through exact solution of the nonlinear differential rate equations derived from the model. The rate constants for the binding of these two required calcium ions are calculated from the solutions as best fit parameters. The thermodynamic equilibrium constants, K1 and K2, for the binding of these two calcium ions are calculated from ratios of the forward and reverse rate constants as 0.6 × 104 and 5.4 × 104, respectively. Thus, the model correctly predicts positively cooperative calcium ion binding for at least the two calcium ions required to induce fluorescence quenching.

Revealing the role of the product metal in DNA polymerase β catalysis

DNA polymerases catalyze a metal-dependent nucleotidyl transferase reaction during extension of a DNA strand using the complementary strand as a template. The reaction has long been considered to require two magnesium ions. Recently, a third active site magnesium ion was identified in some DNA polymerase product crystallographic structures, but its role is not known. Using quantum mechanical/ molecular mechanical calculations of polymerase β, we find that a third magnesium ion positioned near the newly identified product metal site does not alter the activation barrier for the chemical reaction indicating that it does not have a role in the forward reaction. This is consistent with time-lapse crystallographic structures following insertion of Sp-dCTPαS. Although sulfur substitution deters product metal binding, this has only a minimal effect on the rate of the forward reaction. Surprisingly, monovalent sodium or ammonium ions, positioned in the product metal site, lowered the activation barrier. These calculations highlight the impact that an active site water network can have on the energetics of the forward reaction and how metals or enzyme side chains may interact with the network to modulate the reaction barrier. These results also are discussed in the context of earlier findings indicating that magnesium at the product metal position blocks the reverse pyrophosphorolysis reaction.

Hiding in Plain Sight: The Bimetallic Magnesium Covalent Bond in Enzyme Active Sites

The transfer of phosphate groups is an essential function of many intracellular biological enzymes. The transfer is in many cases facilitated by a protein scaffold involving two closely spaced magnesium "ions". It has long been a mystery how these "ions" can retain their closely spaced positions throughout enzymatic phosphate transfer: Coulomb's law would dictate large repulsive forces between these ions at the observed distances. Here we show, however, that the electron density can be borrowed from nearby electron-rich oxygens to populate a bonding molecular orbital that is largely localized between the magnesium "ions". The result is that the Mg-Mg core of these phosphate transfer enzymes is surprisingly similar to a metastable [Mg₂]²⁺ ion in the gas phase, an ion that has been identified experimentally and studied with high-level quantum-mechanical calculations. This similarity is confirmed by comparative computations of the electron densities of [Mg₂]²⁺ in the gas phase and the Mg-Mg core in the structures derived from QM/MM studies of high-resolution X-ray crystal structures. That there is a level of covalent bonding between the two Mg "ions" at the core of these enzymes is a novel concept that enables an improved vision of how these enzymes function at the molecular level. The concept is broader than magnesium-other biologically relevant metals (e.g., Mn and Zn) can also form similar stabilizing covalent Me-Me bonds in both organometallic and inorganic crystals.

Do the crystallographic forms of prethrombin-2 revert to a single form in solution?

It has been earlier established (Pozzi et al. Biochemistry 50 (2011) 10195-10202) that prethrombin-2 crystallizes into two similar but distinct forms: a collapsed form and an alternative form. We employed long molecular dynamics (MD) simulations for these two forms to obtain solvent-equilibrated forms. We find that, at 200ns, the simulated solution collapsed form is quite similar to the X-ray crystal collapsed form, while the simulated solution alternative form deviates from the X-ray crystal alternative form as well as from the solution collapsed form. A detailed structural analysis suggests that the fluctuation of the 140s-loop, in cross-talk with the 220s-loop, may alter the conformation of the W215-E217 segment near the nascent thrombin active site. A rationale is provided for the manner in which interactions of prethrombin-2 with FVa may affect the equilibrium between the two forms of prethrombin-2.

Applications of quantum mechanical/molecular mechanical methods to the chemical insertion step of DNA and RNA polymerization

We review theoretical attempts to model the chemical insertion reactions of nucleoside triphosphates catalyzed by the nucleic acid polymerases using combined quantum mechanical/molecular mechanical methodology. Due to an existing excellent database of high-resolution X-ray crystal structures, the DNA polymerase β system serves as a useful template for discussion and comparison. The convergence of structures of high-quality complexes and continued developments of theoretical techniques suggest a bright future for understanding the global features of nucleic acid polymerization.

Analysis on long-range residue-residue communication using molecular dynamics

We investigated the possibility of inter-residue communication of side chains in barstar, an 89 residue protein, using mutual information theory. The normalized mutual information (NMI) of the dihedral angles of the side chains was obtained from all-atom molecular dynamics simulations. The accumulated NMI from an explicit solvent equilibrated trajectory (600 ns) with free backbone exhibits a parabola-shaped distribution over the inter-residue distances (0-36 Å): smaller at the end regimes but larger in the middle regime. This analysis, plus several other measures, does not find unusual long-range communication for free backbone in explicit solvent simulations.

Amino acid substitution in the active site of DNA polymerase β explains the energy barrier of the nucleotidyl transfer reaction

DNA polymerase β (pol β) is a bifunctional enzyme widely studied for its roles in base excision DNA repair, where one key function is gap-filling DNA synthesis. In spite of significant progress in recent years, the atomic level mechanism of the DNA synthesis reaction has remained poorly understood. Based on crystal structures of pol β in complex with its substrates and theoretical considerations of amino acids and metals in the active site, we have proposed that a nearby carboxylate group of Asp256 enables the reaction by accepting a proton from the primer O3'group, thus activating O3'as the nucleophile in the reaction path. Here, we tested this proposal by altering the side chain of Asp256 to Glu and then exploring the impact of this conservative change on the reaction. The D256E enzyme is more than 1000-fold less active than the wild-type enzyme, and the crystal structures are subtly different in the active sites of the D256E and wild-type enzymes. Theoretical analysis of DNA synthesis by the D256E enzyme shows that the O3'proton still transfers to the nearby carboxylate of residue 256. However, the electrostatic stabilization and location of the O3' proton transfer during the reaction path are dramatically altered compared with wild-type. Surprisingly, this is due to repositioning of the Arg254 side chain in the Glu256 enzyme active site, such that Arg254 is not in position to stabilize the proton transfer from O3'. The theoretical results with the wild-type enzyme indicate an early charge reorganization associated with the O3' proton transfer, and this does not occur in the D256E enzyme. The charge reorganization is mediated by the catalytic magnesium ion in the active site.

Weakly Antiferromagentic Coupling Via Superexchange Interaction Between Mn(II)-Mn(II) Atoms: A QM/MM Study of the Active Site of Human Cytosolic X-Propyl Aminopeptidase P

We investigate the dinuclear manganese, Mn(II)-Mn(II), active site of human cytosolic X-propyl aminopeptidase (XPNPEP1) employing the QM/MM method. The optimized structure supports two manganese atoms at the active site and excludes the possibility of a single Mn(II) atom or other combination of divalent metal ions: Ca(II), Fe(II), Mg(II). A broken symmetry solution verifies an antiferromagnetically coupled state between the Mn(II)-Mn(II) pair, which is the ground state. From the energy difference between the high spin state (HS) and the broken symmetry state (BS), we estimate the exchange coupling constant, J, to be 5.15 cm(-1). Also, we observe multiple bridges (p orbitals) from solvent and two carboxylate linking to the Mn(II)-Mn(II), which leads to the weakly antiferromagnetic interaction of d(5)-d(5) electrons through superexchange coupling.

HPAM: Hirshfeld Partitioned Atomic Multipoles

An implementation of the Hirshfeld (HD) and Hirshfeld-Iterated (HD-I) atomic charge density partitioning schemes is described. Atomic charges and atomic multipoles are calculated from the HD and HD-I atomic charge densities for arbitrary atomic multipole rank l(max) on molecules of arbitrary shape and size. The HD and HD-I atomic charges/multipoles are tested by comparing molecular multipole moments and the electrostatic potential (ESP) surrounding a molecule with their reference ab initio values. In general, the HD-I atomic charges/multipoles are found to better reproduce ab initio electrostatic properties over HD atomic charges/multipoles. A systematic increase in precision for reproducing ab initio electrostatic properties is demonstrated by increasing the atomic multipole rank from l(max) = 0 (atomic charges) to l(max) = 4 (atomic hexadecapoles). Both HD and HD-I atomic multipoles up to rank l(max) are shown to exactly reproduce ab initio molecular multipole moments of rank L for L ≤ l(max). In addition, molecular dipole moments calculated by HD, HD-I, and ChelpG atomic charges only (l(max) = 0) are compared with reference ab initio values. Significant errors in reproducing ab initio molecular dipole moments are found if only HD or HD-I atomic charges used.

A finite field method for calculating molecular polarizability tensors for arbitrary multipole rank

A finite field method for calculating spherical tensor molecular polarizability tensors α(lm;l'm') = ∂Δ(lm)/∂ϕ(l'm')* by numerical derivatives of induced molecular multipole Δ(lm) with respect to gradients of electrostatic potential ϕ(l'm')* is described for arbitrary multipole ranks l and l'. Interconversion formulae for transforming multipole moments and polarizability tensors between spherical and traceless Cartesian tensor conventions are derived. As an example, molecular polarizability tensors up to the hexadecapole-hexadecapole level are calculated for water using the following ab initio methods: Hartree-Fock (HF), Becke three-parameter Lee-Yang-Parr exchange-correlation functional (B3LYP), Møller-Plesset perturbation theory up to second order (MP2), and Coupled Cluster theory with single and double excitations (CCSD). In addition, intermolecular electrostatic and polarization energies calculated by molecular multipoles and polarizability tensors are compared with ab initio reference values calculated by the Reduced Variation Space method for several randomly oriented small molecule dimers separated by a large distance. It is discussed how higher order molecular polarizability tensors can be used as a tool for testing and developing new polarization models for future force fields.

Molecular insights into DNA polymerase deterrents for ribonucleotide insertion

DNA polymerases can misinsert ribonucleotides that lead to genomic instability. DNA polymerase β discourages ribonucleotide insertion with the backbone carbonyl of Tyr-271; alanine substitution of Tyr-271, but not Phe-272, resulted in a >10-fold loss in discrimination. The Y271A mutant also inserted ribonucleotides more efficiently than wild type on a variety of ribonucleoside (rNMP)-containing DNA substrates. Substituting Mn(2+) for Mg(2+) decreased sugar discrimination for both wild-type and mutant enzymes primarily by increasing the affinity for rCTP. This facilitated crystallization of ternary substrate complexes of both the wild-type and Y271A mutant enzymes. Crystallographic structures of Y271A- and wild type-substrate complexes indicated that rCTP is well accommodated in the active site but that O2' of rCTP and the carbonyl oxygen of Tyr-271 or Ala-271 are unusually close (∼2.5 and 2.6 Å, respectively). Structure-based modeling indicates that the local energetic cost of positioning these closely spaced oxygens is ∼2.2 kcal/mol for the wild-type enzyme. Because the side chain of Tyr-271 also hydrogen bonds with the primer terminus, loss of this interaction affects its catalytic positioning. Our results support a model where DNA polymerase β utilizes two strategies, steric and geometric, with a single protein residue to deter ribonucleotide insertion.

Modeling of the DNA-binding site of yeast Pms1 by mass spectrometry

Mismatch repair (MMR) corrects replication errors that would otherwise lead to mutations and, potentially, various forms of cancer. Among several proteins required for eukaryotic MMR, MutLα is a heterodimer comprised of Mlh1 and Pms1. The two proteins dimerize along their C-terminal domains (CTDs), and the CTD of Pms1 houses a latent endonuclease that is required for MMR. The highly conserved N-terminal domains (NTDs) independently bind DNA and possess ATPase active sites. Here we use two protein footprinting techniques, limited proteolysis and oxidative surface mapping, coupled with mass spectrometry to identify amino acids involved along the DNA-binding surface of the Pms1-NTD. Limited proteolysis experiments elucidated several basic residues that were protected in the presence of DNA, while oxidative surface mapping revealed one residue that is uniquely protected from oxidation. Furthermore, additional amino acids distributed throughout the Pms1-NTD were protected from oxidation either in the presence of a non-hydrolyzable analog of ATP or DNA, indicating that each ligand stabilizes the protein in a similar conformation. Based on the recently published X-ray crystal structure of yeast Pms1-NTD, a model of the Pms1-NTD/DNA complex was generated using the mass spectrometric data as constraints. The proposed model defines the DNA-binding interface along a positively charged groove of the Pms1-NTD and complements prior mutagenesis studies of Escherichia coli and eukaryotic MutL.

A hetero-dimer model for concerted action of vitamin K carboxylase and vitamin K reductase in vitamin K cycle

Vitamin K carboxylase (VKC) is believed to convert vitamin K, in the vitamin K cycle, to an alkoxide-epoxide form which then reacts with CO(2) and glutamate to generate γ-carboxyglutamic acid (Gla). Subsequently, vitamin K epoxide reductase (VKOR) is thought to convert the alkoxide-epoxide to a hydroquinone form. By recycling vitamin K, the two integral-membrane proteins, VKC and VKOR, maintain vitamin K levels and sustain the blood coagulation cascade. Unfortunately, NMR or X-ray crystal structures of the two proteins have not been characterized. Thus, our understanding of the vitamin K cycle is only partial at the molecular level. In this study, based on prior biochemical experiments on VKC and VKOR, we propose a hetero-dimeric form of VKC and VKOR that may explain the efficient oxidation and reduction of vitamin K during the vitamin K cycle.

Lee Pedersen’s work in theoretical and computational chemistry and biochemistry

Nature at the lab level in biology and chemistry can be described by the application of quantum mechanics. In many cases, a reasonable approximation to quantum mechanics is classical mechanics realized through Newton’s equations of motion. Dr. Pedersen began his career using quantum mechanics to describe the properties of small molecular complexes that could serve as models for biochemical systems. To describe large molecular systems required a drop-back to classical means and this led surprisingly to a major improvement in the classical treatment of electrostatics for all molecules, not just biological molecules. Recent work has involved the application of quantum mechanics for the putative active sites of enzymes to gain greater insight into the key steps in enzyme catalysis.

Atomic forces for geometry-dependent point multipole and gaussian multipole models

In standard treatments of atomic multipole models, interaction energies, total molecular forces, and total molecular torques are given for multipolar interactions between rigid molecules. However, if the molecules are assumed to be flexible, two additional multipolar atomic forces arise because of (1) the transfer of torque between neighboring atoms and (2) the dependence of multipole moment on internal geometry (bond lengths, bond angles, etc.) for geometry-dependent multipole models. In this study, atomic force expressions for geometry-dependent multipoles are presented for use in simulations of flexible molecules. The atomic forces are derived by first proposing a new general expression for Wigner function derivatives partial derivative D(m'm)(l)/partial derivative Omega. The force equations can be applied to electrostatic models based on atomic point multipoles or gaussian multipole charge density. Hydrogen-bonded dimers are used to test the intermolecular electrostatic energies and atomic forces calculated by geometry-dependent multipoles fit to the ab initio electrostatic potential. The electrostatic energies and forces are compared with their reference ab initio values. It is shown that both static and geometry-dependent multipole models are able to reproduce total molecular forces and torques with respect to ab initio, whereas geometry-dependent multipoles are needed to reproduce ab initio atomic forces. The expressions for atomic force can be used in simulations of flexible molecules with atomic multipoles. In addition, the results presented in this work should lead to further development of next generation force fields composed of geometry-dependent multipole models.

Quantum Chemical Study of the Mechanism of Action of Vitamin K Carboxylase in Solvent

We investigate the post-translational generation of Gla (γ-carboxy glutamic acid) from Glu (glutamic acid) by vitamin K carboxylase (VKC) in solvent. VKC is thought to convert vitamin K, in the vitamin K cycle, to an alkoxide-epoxide form, which then reacts with CO(2) to generate an essential ingredient in blood coagulation, γ-carboxyglutamic acid (Gla). The generation of Gla from Glu is found to be exergenic (-15 kcal/mol) in aqueous solution with the SM6 method. We also produced the free energy profile for this model biochemical process with other solvent methods (polarizable continuum model, dielectric polarizable continuum model) and different dielectric constants. The biological implications are discussed.

Steric, quantum, and electrostatic effects on S(N)2 reaction barriers in gas phase

Biomolecular nucleophilic substitution reactions, S(N)2, are fundamental and commonplace in chemistry. It is the well-documented experimental finding in the literature that vicinal substitution with bulkier groups near the reaction center significantly slows the reaction due to steric hindrance, but theoretical understanding in the quantitative manner about factors dictating the S(N)2 reaction barrier height is still controversial. In this work, employing the new quantification approach that we recently proposed for the steric effect from the density functional theory framework, we investigate the relative contribution of three independent effects-steric, electrostatic, and quantum-to the S(N)2 barrier heights in gas phase for substituted methyl halide systems, R(1)R(2)R(3)CX, reacting with the fluorine anion, where R(1), R(2), and R(3) denote substituting groups and X = F or Cl. We found that in accordance with the experimental finding, for these systems, the steric effect dominates the transition state barrier, contributing positively to barrier heights, but this contribution is largely compensated by the negative, stabilizing contribution from the quantum effect due to the exchange-correlation interactions. Moreover, we find that it is the component from the electrostatic effect that is linearly correlated with the S(N)2 barrier height for the systems investigated in the present study. In addition, we compared our approach with the conventional method of energy decomposition in density functional theory as well as examined the steric effect from the wave function theory for these systems via natural bond orbital analysis.

Recent estimates of the structure of the factor VIIa (FVIIa)/tissue factor (TF) and factor Xa (FXa) ternary complex

The putative structure of the Tissue Factor/Factor VIIa/Factor Xa (TF/FVIIa/FXa) ternary complex is reconsidered. Two independently derived docking models proposed in 2003 (one for our laboratory: CHeA and one from the Scripps laboratory: Ss) are dynamically equilibrated for over 10 ns in an electrically neutral solution using all-atom molecular dynamics. Although the dynamical models (CHeB and Se) differ in atomic detail, there are similarities in that TF is found to interact with the gamma-carboxyglutamic acid (Gla) and Epidermal Growth Factor-like 1 (EGF-1) domains of FXa, and FVIIa is found to interact with the Gla, EGF-2 and serine protease (SP) domains of FXa in both models. FVIIa does not interact with the FXa EGF-1 domain in Se and the EGF domains of FVIIa do not interact with FXa in the CHeB. Both models are consistent with experimentally suggested contacts between the SP domain of FVIIa with the EGF-2 and SP domains of FXa.

Molecular acidity: A quantitative conceptual density functional theory description

Accurate predictions of molecular acidity using ab initio and density functional approaches are still a daunting task. Using electronic and reactivity properties, one can quantitatively estimate pKa values of acids. In a recent paper [S. B. Liu and L. G. Pedersen, J. Phys. Chem. A 113, 3648 (2009)], we employed the molecular electrostatic potential (MEP) on the nucleus and the sum of valence natural atomic orbital (NAO) energies for the purpose. In this work, we reformulate these relationships on the basis of conceptual density functional theory and compare the results with those from the thermodynamic cycle method. We show that MEP and NAO properties of the dissociating proton of an acid should satisfy the same relationships with experimental pKa data. We employ 27 main groups and first to third row transition metal-water complexes as illustrative examples to numerically verify the validity of these strong linear correlations. Results also show that the accuracy of our approach and that of the conventional method through the thermodynamic cycle are statistically similar.

Gaussian Multipole Model (GMM)

An electrostatic model based on charge density is proposed as a model for future force fields. The model is composed of a nucleus and a single Slater-type contracted Gaussian multipole charge density on each atom. The Gaussian multipoles are fit to the electrostatic potential (ESP) calculated at the B3LYP/6-31G* and HF/aug-cc-pVTZ levels of theory and tested by comparing electrostatic dimer energies, inter-molecular density overlap integrals, and permanent molecular multipole moments with their respective ab initio values. For the case of water, the atomic Gaussian multipole moments Q(lm) are shown to be a smooth function of internal geometry (bond length and bond angle), which can be approximated by a truncated linear Taylor series. In addition, results are given when the Gaussian multipole charge density is applied to a model for exchange-repulsion energy based on the inter-molecular density overlap.

Estimation of molecular acidity via electrostatic potential at the nucleus and valence natural atomic orbitals

An effective approach of estimating molecular pK(a) values from simple density functional calculations is proposed in this work. Both the molecular electrostatic potential (MEP) at the nucleus of the acidic atom and the sum of valence natural atomic orbitals are employed for three categories of compounds, amines and anilines, carbonyl acids and alcohols, and sulfonic acids and thiols. A strong correlation between experimental pK(a) values and each of these two quantities for each of the three categories has been discovered. Moreover, if the MEP is subtracted by the isolated atomic MEP for each category of compounds, we observe a single unique linear relationship between the resultant MEP difference and experimental pK(a) data of amines, anilines, carbonyl acids, alcohols, sulfonic acids, thiols, and their substituents. These results can generally be utilized to simultaneously estimate pK(a) values at multiple sites with a single calculation for either relatively small molecules in drug design or amino acids in proteins and macromolecules.

Conformational change path between closed and open forms of C2 domain of coagulation factor V on a two-dimensional free-energy surface

We test a hypothesis that the closed form of the C2 domain of coagulation factor V is more stable than the open form in an aqueous environment using a two-dimensional free-energy calculation with a simple dielectric solvent model. Our result shows that while the free-energy difference between two forms is small, favoring the closed form, a two-dimensional free-energy surface (FES) reveals that a transition state (1.53 kcal/mol) exists between the two conformations. By mapping the one-dimensional order parameter DeltaQ onto the two-dimensional FES, we search the conformational change path with the highest Boltzmann weighting factor between the closed and open form of the factor V C2 domain. The predicted transition path from the closed to open form is not that of simple side chain movements, but instead concerted movements of several loops. We also present a one-dimensional free-energy profile using a collective order parameter, which in a coarse manner locates the energy barriers found on the two-dimensional FES.

Reaction mechanism of the epsilon subunit of E. coli DNA polymerase III: insights into active site metal coordination and catalytically significant residues

The 28 kDa epsilon subunit of Escherichia coli DNA polymerase III is the exonucleotidic proofreader responsible for editing polymerase insertion errors. Here, we study the mechanism by which epsilon carries out the exonuclease activity. We performed quantum mechanics/molecular mechanics calculations on the N-terminal domain containing the exonuclease activity. Both the free-epsilon and a complex epsilon bound to a theta homologue (HOT) were studied. For the epsilon-HOT complex Mg(2+) or Mn(2+) were investigated as the essential divalent metal cofactors, while only Mg(2+) was used for free-epsilon. In all calculations a water molecule bound to the catalytic metal acts as the nucleophile for hydrolysis of the phosphate bond. Initially, a direct proton transfer to H162 is observed. Subsequently, the nucleophilic attack takes place followed by a second proton transfer to E14. Our results show that the reaction catalyzed with Mn(2+) is faster than that with Mg(2+), in agreement with experiment. In addition, the epsilon-HOT complex shows a slightly lower energy barrier compared to free-epsilon. In all cases the catalytic metal is observed to be pentacoordinated. Charge and frontier orbital analyses suggest that charge transfer may stabilize the pentacoordination. Energy decomposition analysis to study the contribution of each residue to catalysis suggests that there are several important residues. Among these, H98, D103, D129, and D146 have been implicated in catalysis by mutagenesis studies. Some of these residues were found to be structurally conserved on human TREX1, the exonuclease domains from E. coli DNA-Pol I, and the DNA polymerase of bacteriophage RB69.

A computational modeling and molecular dynamics study of the Michaelis complex of human protein Z-dependent protease inhibitor (ZPI) and factor Xa (FXa)

Protein Z-dependent protease inhibitor (ZPI) and antithrombin III (AT3) are members of the serpin superfamily of protease inhibitors that inhibit factor Xa (FXa) and other proteases in the coagulation pathway. While experimental structural information is available for the interaction of AT3 with FXa, at present there is no structural data regarding the interaction of ZPI with FXa, and the precise role of this interaction in the blood coagulation pathway is poorly understood. In an effort to gain a structural understanding of this system, we have built a solvent equilibrated three-dimensional structural model of the Michaelis complex of human ZPI/FXa using homology modeling, protein-protein docking and molecular dynamics simulation methods. Preliminary analysis of interactions at the complex interface from our simulations suggests that the interactions of the reactive center loop (RCL) and the exosite surface of ZPI with FXa are similar to those observed from X-ray crystal structure-based simulations of AT3/FXa. However, detailed comparison of our modeled structure of ZPI/FXa with that of AT3/FXa points to differences in interaction specificity at the reactive center and in the stability of the inhibitory complex, due to the presence of a tyrosine residue at the P1 position in ZPI, instead of the P1 arginine residue in AT3. The modeled structure also shows specific structural differences between AT3 and ZPI in the heparin-binding and flexible N-terminal tail regions. Our structural model of ZPI/FXa is also compatible with available experimental information regarding the importance for the inhibitory action of certain basic residues in FXa.

Exploring the origin of the internal rotational barrier for molecules with one rotatable dihedral angle

Continuing our recent endeavor, we systematically investigate in this work the origin of internal rotational barriers for small molecules using the new energy partition scheme proposed recently by one of the authors [S. B. Liu, J. Chem. Phys. 126, 244103 (2007)], where the total electronic energy is decomposed into three independent components, steric, electrostatic, and fermionic quantum. Specifically, we focus in this work on six carbon, nitrogen, and oxygen containing hydrides, CH(3)CH(3), CH(3)NH(2), CH(3)OH, NH(2)NH(2), NH(2)OH, and H(2)O(2), with only one rotatable dihedral angle [angle]H-X-Y-H (X,Y=C,N,O). The relative contributions of the different energy components to the total energy difference as a function of the internal dihedral rotation will be considered. Both optimized-geometry (adiabatic) and fixed-geometry (vertical) differences are examined, as are the results from the conventional energy partition and natural bond orbital analysis. A wealth of strong linear relationships among the total energy difference and energy component differences for different systems have been observed but no universal relationship applicable to all systems for both cases has been discovered, indicating that even for simple systems such as these, there exists no omnipresent, unique interpretation on the nature and origin of the internal rotation barrier. Different energy components can be employed for different systems in the rationalization of the barrier height. Confirming that the two differences, adiabatic and vertical, are disparate in nature, we find that for the vertical case there is a unique linear relationship applicable to all the six molecules between the total energy difference and the sum of the kinetic and electrostatic energy differences. For the adiabatic case, it is the total potential energy difference that has been found to correlate well with the total energy difference except for ethane whose rotation barrier is dominated by the quantum effect.

Computational study of the putative active form of protein Z (PZa): sequence design and structural modeling

Although protein Z (PZ) has a domain arrangement similar to the essential coagulation proteins FVII, FIX, FX, and protein C, its serine protease (SP)-like domain is incomplete and does not exhibit proteolytic activity. We have generated a trial sequence of putative activated protein Z (PZa) by identifying amino acid mutations in the SP-like domain that might reasonably resurrect the serine protease catalytic activity of PZ. The structure of the activated form was then modeled based on the proposed sequence using homology modeling and solvent-equilibrated molecular dynamics simulations. In silico docking of inhibitors of FVIIa and FXa to the putative active site of equilibrated PZa, along with structural comparison with its homologous proteins, suggest that the designed PZa can possibly act as a serine protease.

Proposed structural models of human factor Va and prothrombinase

BACKGROUND

The prothrombinase complex consists of factor Xa, FVa, calcium ions, and phospholipid membrane. The prothrombinase complex plays a key role in the blood coagulation process.

OBJECTIVE

To derive solvent-equilibrated models of human FVa and the prothrombinase complex.

METHODS

Several modeling techniques have been employed, including homology modeling, protein-protein docking, and molecular dynamics simulation methods, to build the structural models.

RESULTS AND CONCLUSIONS

We found, upon simulation, a possibly significant shift towards planarity of the five FVa domains. To estimate a prothrombinase structure, we docked an FXa model to the equilibrated FVa model using experimental data as docking filters. We found that simulation of the docked complex led to some changes in the protein-protein contacts, but not buried surface area, as compared to the initial docking model. Possible locations of prothrombin binding to prothrombinase are indicated.

A quantum chemical study of the mechanism of action of Vitamin K epoxide reductase (VKOR)

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.

A quantum chemical study of the mechanism of action of Vitamin K carboxylase (VKC)

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.

Transcription levels of endothelin-1 and endothelin receptors are associated with age and leaflet location in porcine mitral valves

The aim of the study was to investigate the expression levels of endothelin-1 (ET-1) and ET(A) and ET(B) receptors (ET(A)-R and ET(B)-R) in porcine mitral valves and associate the transcription levels to age, leaflet location and deposition of mucopolysaccharides (MPS). Tissue samples from the chordal and inter-chordal insertion area of the anterior mitral valve leaflet from 11 sows (> or = 2 years of age) and 10 slaughter pigs (approximately 6 months old) were obtained and the relative gene expression levels of ET-1, ET(A)-R and ET(B)-R measured by semi-quantitative real-time PCR. A separate tissue sample was taken for histopathological grading of MPS deposition. The transcription levels of ET-1 (P < 0.0001) and ET(A)-R (P < 0.0004) were significantly higher in leaflets from the sows compared with slaughter pigs. The gene expression of ET(B)-R was not associated to age (P = 0.38), but increased in chordal insertion areas compared with inter-chordal areas (P = 0.01). The expression of ET-1 and ET(A)-R mRNA did not differ significantly between the two leaflet locations. The valve leaflets from sows had a significantly increased degree of MPS deposition compared with slaughter pigs upon histological examination (P = 0.04). In conclusion, an age-related valvular degeneration is observed in porcine mitral valve leaflets and ET-1 is suggested to be involved through action of both ET(A) and ET(B) receptors.

Direct measurements of nitric oxide release in relation to expression of endothelial nitric oxide synthase in isolated porcine mitral valves

The aim of this study was to measure the direct release of nitric oxide (NO) from the porcine mitral valve using a NO microelectrode. Furthermore, the expression and localization of endothelial nitric oxide synthase (eNOS) in the mitral valve was studied using immunohistochemistry, Western blotting and RT-PCR. Results show that bradykinin increases NO release from mitral valves (DeltaBradykinin: 33.71 +/- 10.41 nm NO, P < 0.001, n = 10), whereas N-nitro-l-arginine methyl esther (l-NAME) decreases NO release when compared with basal level (Deltal-NAME: 82.69 +/- 15.66 nm NO, P < 0.005, n = 4). Both protein and mRNA expression of eNOS in mitral valves and in isolated valvular endothelial cells suggest that the NO release is mainly associated with the mitral valve endothelium. It is concluded that direct NO release from porcine mitral valves coincides with eNOS expression. This study documents useful techniques for investigations into the role of local NO release in mitral valve diseases.

Quantum chemical study of the mechanism of action of vitamin K carboxylase (VKC). IV. Intermediates and transition states

We studied proposed steps for the enzymatic formation of gamma-carboxyglutamic acid by density functional theory (DFT) quantum chemistry. Our results for one potentially feasible mechanism show that a vitamin K alkoxide intermediate can abstract a proton from glutamic acid at the gamma-carbon to form a carbanion and vitamin K epoxide. The hydrated carbanion can then react with CO2 to form gamma-carboxyglutamic acid. Computations at the B3LYP/6-311G** level were used to determine the intermediates and transition states for the overall process. The activation free energy for the gas-phase path is 22 kcal/mol, with the rate-limiting step for the reaction being the attack of the carbanion on CO2. Additional solvation studies, however, indicate that the formation of the carbanion step can be competitive with the CO2 attack step in high-dielectric systems. We relate these computations to the entire vitamin K cycle in the blood coagulation cascade, which is essential for viability of vertebrates.

Increased expression of endothelin B receptor in static stretch exposed porcine mitral valve leaflets

The aim of this study was to evaluate the effect of mechanical stretch on the expression of ET-1 and ET(A)- and ET(B)-receptors in porcine mitral valve leaflets. Leaflet segments from 10 porcine mitral valves were exposed to a static stretch load of 1.5 N for 3.5h in buffer at 37 degrees C together with matching control segments. Subsequently, the mRNA expression of ET-1, ET(A)-R and ET(B)-R was measured by real-time RT-PCR in the chordal insertion areas. The analyses showed an increased transcription of ET(B)-receptors in stretch-exposed leaflet segments compared to unstretched segments median 2.23 (quartiles 1.37 and 2.70) vs. median 1.56 (quartiles 1.38 and 2.17, P=0.03) whereas the mRNA expression of ET(A)-receptors (P=0.90) and ET-1 (P=0.51) remained unchanged. Stretch increased the expression of ET(B)-receptors in porcine mitral valve leaflets. The finding could lead to a better understanding of the pathogenesis of myxomatous mitral valve disease.

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.