Estimation of Molecular Acidity via Electrostatic Potential at the Nucleus and Valence Natural Atomic Orbitals

An effective approach of estimating molecular pKa 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 pKa 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 pKa data of amines, anilines, carbonyl acids, alcohols, sulfonic acids, thiols, and their substituents. These results can generally be utilized to simultaneously estimate pKa values at multiple sites with a single calculation for either relatively small molecules in drug design or amino acids in proteins and macromolecules

Clar Sextet Analysis of Triangular, Rectangular and Honeycomb Graphene Antidot Lattices

Pristine graphene is a semimetal and thus does not have a band gap. By making a nanometer scale periodic array of holes in the graphene sheet a band gap may form; the size of the gap is controllable by adjusting the parameters of the lattice. The hole diameter, hole geometry, lattice geometry and the separation of the holes are parameters that all play an important role in determining the size of the band gap, which, for technological applications, should be at least of the order of tenths of an eV. We investigate four different hole configurations: the rectangular, the triangular, the rotated triangular and the honeycomb lattice. It is found that the lattice geometry plays a crucial role for size of the band gap: the triangular arrangement displays always a sizable gap, while for the other types only particular hole separations lead to a large gap. This observation is explained using Clar sextet theory, and we find that a sufficient condition for a large gap is that the number of sextets exceeds one third of the total number of hexagons in the unit cell. Furthermore, we investigate non-isosceles triangular structures to probe the sensitivity of the gap in triangular lattices to small changes in geometry

Analysis on long-range residue-residue communication using molecular dynamics: Analysis of Long-Range Residue-Residue Communication

We investigated the possibility of inter-residue communication of side chains in barstar, an 89 residue protein, employing 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

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 ΔQ 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

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

A revisit to the one-form kinetic model of prothrombinase: A comment on the rebuttal

Thrombin is generated enzymatically from prothrombin by two pathways with the intermediates of meizothrombin and prethrombin-2. Experimental concentration profiles from two independent groups for these two pathways have been re-analyzed. By rationally combining the independent data sets, a simple mechanism can be established and rate constants determined. A structural model is consistent with the data-derived finding that mechanisms that feature channeling or ratcheting are not necessary to describe thrombin production

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 lmax 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 lmax = 0 (atomic charges) to lmax = 4 (atomic hexadecapoles). Both HD and HD-I atomic multipoles up to rank lmax are shown to exactly reproduce ab initio molecular multipole moments of rank L for L ≤ lmax. In addition, molecular dipole moments calculated by HD, HD-I, and ChelpG atomic charges only (lmax = 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

Spin-orbit interaction and the 'metal-insulator' transition observed in two-dimensional hole systems

We present calculations of the spin and phase relaxation rates in GaAs/AlGaAs $p$-type quantum wells. These rates are used to derive the temperature dependence of the weak-localization correction to the conductivity. In $p$-type quantum wells both weak localization and weak anti-localization are present due to the strong spin-orbit interaction. When determining the total conductivity correction one also have to include the term due to hole-hole interaction. The magnitude of the latter depends on the ratio between the thermal energy and the Fermi energy, $k_{\rm B}T/E_{\rm F}$ and whether the system can be considered as ballistic $(k_{\rm B}T \tau_{\rm tr} / \hbar>1)$ or diffusive ($k_{\rm B}T \tau_{\rm tr}/\hbar<1$). We argue that due to the relatively low Fermi energy and the moderate mobilities, in the $p$-type systems in question, the conductivity correction arising from hole-hole interactions is negligible at the highest temperatures accessible in the experiments. Hence the 'metal-insulator' transition observed at these relatively high temperatures could be caused by interference effects. We compare our calculations of the weak anti-localization correction with the experimental results from different independent groups with special emphasis on the experiments by Simmons et al. We find good agreement between predicted and observed transistion density $p_{c}$.Comment: 6 pages, 4 figures. Accepted to PRB (15 June, 2002

Weak Antiferromagnetic Coupling via a Superexchange Interaction between Mn(II)–Mn(II) Ions: 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 d5-d5 electrons through superexchange coupling

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 γ-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