Experimental, Hartree−Fock, and Density Functional Theory Investigations of the Charge Density, Dipole Moment, Electrostatic Potential, and Electric Field Gradients inl-Asparagine Monohydrate

Using general computational chemistry strategy to unravel the reactivity of emerging pollutants: An example of sulfonamide chlorination

Increasing number of emerging pollutants in environments requires an effective approach which can facilitate the prediction of reactivity and provide insights into the reaction mechanisms. Computational chemistry is exactly the tool to fulfill this demand with its good performance in theoretical investigation of chemical reactions at molecular level. In this study, chlorination of sulfonamide antibiotics is used as an illustration to present a systematic strategy demonstrating how computational chemistry can be applied to investigate the reaction behavior of emerging pollutants. Sulfonamides is a class of micropollutants that contain the common structure of 4-aminobenzenesulfonmaide while differ in their heterocycles. Based on the calculated conceptual DFT indices, the reactive sites of sulfonamide are successfully predicted, which locate on their common structure of 4-aminobenzenesulfonmaide. Therefore, all sulfonamides follow the similar reaction pathway. Product identification by LTQ-Orbitrap MS further verifies the in silico prediction. Three critical pathways are discovered, i.e., S-N bond cleavage, Cl-substitution onto aniline-N, and the following rearrangement to lose -SO₂- group, among which Cl-substitution is the key step due to its lowest free energy barrier. Heterocycles impact the reaction rate by affecting the electronic density of aniline group. In general, the more electron-donating the heterocycle is, the more readily sulfonamides to be chlorinated.

Mechanism of Catalysis by l-Asparaginase

Two bacterial type II l-asparaginases, from Escherichia coli and Dickeya chrysanthemi, have played a critical role for more than 40 years as therapeutic agents against juvenile leukemias and lymphomas. Despite a long history of successful pharmacological applications and the apparent simplicity of the catalytic reaction, controversies still exist regarding major steps of the mechanism. In this report, we provide a detailed description of the reaction catalyzed by E. coli type II l-asparaginase (EcAII). Our model was developed on the basis of new structural and biochemical experiments combined with previously published data. The proposed mechanism is supported by quantum chemistry calculations based on density functional theory. We provide strong evidence that EcAII catalyzes the reaction according to the double-displacement (ping-pong) mechanism, with formation of a covalent intermediate. Several steps of catalysis by EcAII are unique when compared to reactions catalyzed by other known hydrolytic enzymes. Here, the reaction is initiated by a weak nucleophile, threonine, without direct assistance of a general base, although a distant general base is identified. Furthermore, tetrahedral intermediates formed during the catalytic process are stabilized by a never previously described motif. Although the scheme of the catalytic mechanism was developed only on the basis of data obtained from EcAII and its variants, this novel mechanism of enzymatic hydrolysis could potentially apply to most (and possibly all) l-asparaginases.

Electric Dipole Descriptor for Machine Learning Prediction of Catalyst Surface-Molecular Adsorbate Interactions

The challenge of evaluating catalyst surface-molecular adsorbate interactions holds the key for rational design of catalysts. Finding an experimentally measurable and theoretically computable descriptor for evaluating surface-adsorbate interactions is a significant step toward achieving this goal. Here we show that the electric dipole moment can serve as a convenient yet accurate descriptor for establishing structure-property relationships for molecular adsorbates on metal catalyst surfaces. By training a machine learning neural network with a large data set of first-principles calculations, we achieve quick and accurate predictions of molecular adsorption energy and transferred charge. The training model using NO/CO@Au(111) can be extended to study additional substrates such as Au(001) or Ag(111), thus exhibiting extraordinary transferability. These findings validate the effectiveness of the electric dipole descriptor, providing an efficient modality for future catalyst design.

Molecular Dynamics model of peptide-protein conjugation: case study of covalent complex between Sos1 peptide and N-terminal SH3 domain from Grb2

We have investigated covalent conjugation of VPPPVPPRRRX′ peptide (where X′ denotes N^ε-chloroacetyl lysine) to N-terminal SH3 domain from adapter protein Grb2. Our experimental results confirmed that the peptide first binds to the SH3 domain noncovalently before establishing a covalent linkage through reaction of X′ with the target cysteine residue C32. We have also confirmed that this reaction involves a thiolate-anion form of C32 and follows the S_N2 mechanism. For this system, we have developed a new MD-based protocol to model the formation of covalent conjugate. The simulation starts with the known coordinates of the noncovalent complex. When two reactive groups come into contact during the course of the simulation, the reaction is initiated. The reaction is modeled via gradual interpolation between the two sets of force field parameters that are representative of the noncovalent and covalent complexes. The simulation proceeds smoothly, with no appreciable perturbations to temperature, pressure or volume, and results in a high-quality MD model of the covalent complex. The validity of this model is confirmed using the experimental chemical shift data. The new MD-based approach offers a valuable tool to explore the mechanics of protein-peptide conjugation and build accurate models of covalent complexes.

Geometric considerations support the double-displacement catalytic mechanism of l-asparaginase

Twenty crystal structures of the complexes of l-asparaginase with l-Asn, l-Asp, and succinic acid that are currently available in the Protein Data Bank, as well as 11 additional structures determined in the course of this project, were analyzed in order to establish the level of conservation of the geometric parameters describing interactions between the substrates and the active site of the enzymes. We found that such stereochemical relationships are highly conserved, regardless of the organism from which the enzyme was isolated, specific crystallization conditions, or the nature of the ligands. Analysis of the geometry of the interactions, including Bürgi-Dunitz and Flippin-Lodge angles, indicated that Thr12 (Escherichia coli asparaginase II numbering) is optimally placed to be the primary nucleophile in the most likely scenario utilizing a double-displacement mechanism, whereas catalysis through a single-displacement mechanism appears to be the least likely.

Assessment on third order non linearity and other optical analyses of L-Asparagine Monohydrate single crystal: An efficient candidate for harmonic conversions

Single crystal of l-Asparagine Monohydrate, an organic material has been successfully grown by slow evaporation solution growth technique at ambient condition. The lattice parameters and its strain of the grown crystal have been evaluated from powder X-ray diffraction and found that it belongs to orthorhombic crystal system. The polarizability has been measured by using the Clausius-Mossotti relation. The crystalline perfection of grown single crystal has been examined by high resolution X-ray diffraction and its imperfection in the diffraction plane was clearly visible by recording topographical image of the plane. From the high resolution XRD, it confirms that the crystal contained high crystalline perfection. The optical behavior was analyzed by photoluminescence and birefringence methods. In the photoluminescence, a broad peak has been observed at 475 nm which suggest that it emits blue light. The decay tendency of the material has also been observed by calculating decay constant. The optical homogeneity has been determined by the dispersion pattern of the material. The two photon absorption coefficient was further calculated by Z-scan, which gives the information about the third order non linear optical behavior of the material. The value of two-photon absorption coefficient is 4.25 × 10(-12)m/W. The thermal parameters like thermal effusivity, thermal diffusivity, specific heat and thermal conductivity was obtained by using photopyroelectric technique. The ferroelectric behavior of the grown specimen was analyzed from PE (polarization VS electric field) loop. The loop suggests that the material was a nearly equivalent to ideal capacitor.

Unraveling the Molecular Mechanism of Enthalpy Driven Peptide Folding by Polyol Osmolytes

Many polyols and carbohydrates serve in different organisms as protective osmolytes that help to stabilize proteins in their native, functional state, even under a variety of environmental stresses. However, despite their important role, much of the molecular mechanism by which these osmolytes exert their action remains elusive. We have recently shown experimentally that, although polyols and carbohydrates are excluded from protein and peptide interfaces, as also expected for the known entropic "crowding" mechanism, the osmolyte folding action can in fact primarily be enthalpic in nature. To follow this newly resolved enthalpically driven stabilization mechanism, we report here on molecular dynamics simulations of a model peptide that can fold in solution into a β-hairpin. In agreement with experiments, our simulations indicate that sorbitol, a representative polyol, promotes peptide folding by preferential exclusion. At the molecular level, simulations further show that peptide stabilization can be explained by sorbitol's perturbation of the solution hydrogen bonding network in the peptide first hydration shells. Consequently, fewer hydrogen bonds between peptide and solvating water are lost upon folding, and additional internal peptide hydrogen bonds are formed in the presence of sorbitol, while internal peptide and water-associated hydrogen bonds are strengthened, resulting in stabilization of the peptide folded state. We further find that changes in water orientational entropy are reduced upon folding in sorbitol solution, reflecting the struggle of water molecules to maintain optimal hydrogen bonding in the presence of competing polyols. By providing first molecular underpinnings for enthalpically driven osmolyte stabilization of peptides and proteins, this mechanism should allow a better understanding of the variety of physical forces by which protective osmolytes act in biologically realistic solutions.

Experimental charge density of an L-phenylalanine formic acid complex with a short hydrogen bond determined at 25 K

The experimental charge density and related atomic and bond topological properties of an L-phenylalanine formic acid complex were derived from a high resolution X-ray data set (sin θ/λ = 1.18 Å–1/d = 0.42 Å) measured at 25 K. The complex consists of a zwitterionic and a cationic phenylalanine molecule with formate as counterion. Special focus was directed on the density distribution in the region of a strong O—H ·· O hydrogen bond (O ·· O = 2.491(1) Å) which is formed between the two phenylalanine units. The obtained results are compared with the 15 previously derived experimental amino acid charge density data, with various theoretical calculations at experimental geometries and with the complete set of topological descriptors based on ab initio calculations of the neutral forms of all 20 amino acids published recently in the literature. A comparison of all available data in this biologically important class of compounds gives an impression about the significance of the quantitative results from experimental and theoretical charge density determinations.

Dipole Moment Enhancement in Molecular Crystals from X-ray Diffraction Data

Although reliable determination of the molecular dipole moment from experimental charge density analyses on molecular crystals is a challenging undertaking, these values are becoming increasingly common experimental results. We collate all known experimental determinations and use this database to identify broad trends in the dipole moment enhancements implied by these measurements as well as outliers for which enhancements are pronounced. Compelling evidence emerges that molecular dipole moments from X-ray diffraction data can provide a wealth of information on the change in the molecular charge distribution that results from crystal formation. Most importantly, these experiments are unrivalled in their potential to provide this information in such detail and deserve to be exploited to a much greater extent. The considerable number of experimental determinations now available has enabled us to pinpoint those studies that merit further attention, either because they point unequivocally to a considerable enhancement in the crystal (of 50 % or more), or because the experimental determinations suggest enhancements of 100 % or more--much larger than independent theoretical estimates. In both cases further detailed experimental and theoretical studies are indicated.

A Solid State 13C NMR, Crystallographic, and Quantum Chemical Investigation of Phenylalanine and Tyrosine Residues in Dipeptides and Proteins

We report the results of a solid-state NMR and quantum chemical investigation of the 13C gamma NMR chemical shifts in phenylalanine and tyrosine in dipeptides and proteins. Accurate computation of the experimental shifts is shown to require a good description of local electrostatic field effects, and we find the best results (R2=0.94, rmsd=1.6 ppm, range = 17.1 ppm, N=14) by using a self-consistent reaction field continuum model. There are no obvious correlations with phi, psi, chi 1, or chi2 torsion angles, unlike the results seen with other amino acids. There is, however, a linear relation between computed C gamma atomic charges and shifts for the 14 peptide as well as 18 protein residues investigated. This result is similar to the correlation reported in the 1960s between pi-electron density and 13C shifts for classical 4n + 2 (n=0, 1, 2) pi-electron aromatic species, such as cyclopentadienide and the tropylium cation, and in fact, we found that the shielding/atomic charge correlation seen in the peptides and proteins is virtually identical to that seen with a broad range of aromatic carbocations/carbanions. These results suggest the dominance of an electrostatic field polarization model in which increasing pi electron density results in an increase in C gamma atomic charge and increased shielding (of sigma 11 and sigma 22, perpendicular to the pi orbital) in Phe and Tyr, as well as in the other aromatic species. These results are of general interest since they demonstrate the importance of electrostatic field effects on Phe and Tyr C gamma chemical shifts in peptides and proteins and imply that inclusion of these effects will be necessary in order to interpret the shifts of other aromatic species, such as drug molecules, bound to proteins.

Reproducibility and transferability of topological data: experimental charge density study of two modifications ofl-alanyl-l-tyrosyl-l-alanine

Two crystalline modifications of the tripeptide L-Ala-L-Tyr-L-Ala, which have different solvent molecules in the crystal structure (water and ethanol for modifications 1 and 2), were the subject of experimental charge density studies based on high resolution X-ray data collected at ultra-low temperatures of 9 K (1) and 20 K (2), respectively. The molecular structures and the intermolecular interactions were found to be rather similar in the two crystal lattices, so that this study allowed the reproducibility of the charge density of a given molecule in different (but widely comparable) crystalline environments to be examined. With respect to bond topological and atomic properties, the agreement between the two modifications of the title tripeptide was in the same range as found from the comparison with the previously reported results of tri-L-alanine. It follows that the reproducibility and transferability of quantitative topological data are comparable and that within the accuracy of experimental charge density work the replacement of the central amino acid residue L-Ala by L-Tyr has no significant influence, neither on bond nor on the atomic properties of the oligopeptide main chain. Intermolecular interactions in the form of hydrogen bonds were characterized quantitatively and qualitatively by topological criteria and by mapping the charge density distribution on the Hirshfeld surface.

31P NMR Chemical Shifts in Hypervalent Oxyphosphoranes and Polymeric Orthophosphates

We report the first quantum chemical investigation of the solid- and solution-state 31P NMR chemical shifts in models for phosphoryl transfer enzyme reaction intermediates and in polymeric inorganic phosphates. The 31P NMR chemical shifts of five- and six-coordinate oxyphosphoranes containing a variety of substitutions at phosphorus, as well as four-coordinate polymeric orthophosphates and four-coordinate phosphonates, are predicted with a slope of 1.00 and an R2= 0.993 (N = 34), corresponding to a 3.8 ppm (or 2.1%) error over the entire 178.3 ppm experimental chemical shift range, using Hartree-Fock methods. For the oxyphosphoranes, we used either experimental crystallographic structures or, when these were not available, fully geometry optimized molecular structures. For the four-coordinate phosphonates we used X-ray structures together with charge field perturbation, to represent lattice interactions. For the three-dimensional orthophosphates (BPO4, AlPO4, GaPO4 we again used X-ray structures, but for these inorganic systems we employed a self-consistent charge field perturbation approach on large clusters, to deduce peripheral atom charges. For pentaoxyphosphoranes, the solvent effect on 31P NMR chemical shieldings was found to be very small (<0.5 ppm). The 31P NMR chemical shielding tensors in the pentaoxyphosphoranes were in most cases found to be close to axially symmetric and were dominated by changes in the shielding tensor components in the equatorial plane (sigma22 and sigma33). The isotropic shifts were highly correlated (R2= 0.923) with phosphorus natural bonding orbital charges, with the larger charges being associated with shorter axial P-O bond lengths and, hence, more shielding. Overall, these results should facilitate the use of 31P NMR techniques in investigating the structures of more complex systems, such as phosphoryl transfer enzymes, as well as in investigating other, complex oxide structures.

A Solid State 13C NMR, Crystallographic, and Quantum Chemical Investigation of Chemical Shifts and Hydrogen Bonding in Histidine Dipeptides

We report the first solid-state NMR, crystallographic, and quantum chemical investigation of the origins of the 13C NMR chemical shifts of the imidazole group in histidine-containing dipeptides. The chemical shift ranges for Cgamma and Cdelta2 seen in eight crystalline dipeptides were very large (12.7-13.8 ppm); the shifts were highly correlated (R2= 0.90) and were dominated by ring tautomer effects and intermolecular interactions. A similar correlation was found in proteins, but only for buried residues. The imidazole 13C NMR chemical shifts were predicted with an overall rms error of 1.6-1.9 ppm over a 26 ppm range, by using quantum chemical methods. Incorporation of hydrogen bond partner molecules was found to be essential in order to reproduce the chemical shifts seen experimentally. Using AIM (atoms in molecules) theory we found that essentially all interactions were of a closed shell nature and the hydrogen bond critical point properties were highly correlated with the N...H...O (average R2= 0.93) and Nepsilon2...H...N (average R2= 0.98) hydrogen bond lengths. For Cepsilon1, the 13C chemical shifts were also highly correlated with each of these properties (at the Nepsilon2 site), indicating the dominance of intermolecular interactions for Cepsilon1. These results open up the way to analyzing 13C NMR chemical shifts, tautomer states (from Cdelta2, Cepsilon1 shifts), and hydrogen bond properties (from Cepsilon1 shifts) of histidine residue in proteins and should be applicable to imidazole-containing drug molecules bound to proteins, as well.

Quantum chemical studies of protein structure

Quantum chemical methods now permit the prediction of many spectroscopic observables in proteins and related model systems, in addition to electrostatic properties, which are found to be in excellent accord with those determined from experiment. I discuss the developments over the past decade in these areas, including predictions of nuclear magnetic resonance chemical shifts, chemical shielding tensors, scalar couplings and hyperfine (contact) shifts, the isomer shifts and quadrupole splittings in Mössbauer spectroscopy, molecular energies and conformations, as well as a range of electrostatic properties, such as charge densities, the curvatures, Laplacians and Hessians of the charge density, electrostatic potentials, electric field gradients and electrostatic field effects. The availability of structure/spectroscopic correlations from quantum chemistry provides a basis for using numerous spectroscopic observables in determining aspects of protein structure, in determining electrostatic properties which are not readily accessible from experiment, as well as giving additional confidence in the use of these techniques to investigate questions about chemical bonding and chemical reactions.

Multiple-quantum magic-angle spinning spectroscopy using nonlinear sampling

NMR spectroscopy is a relatively insensitive technique and many biomolecular applications operate near the limits of sensitivity and resolution. A particularly challenging example is detection of the quadrupolar nucleus 17O, due to its low natural abundance, large quadrupole couplings, and low gyromagnetic ratio. Yet the chemical shift of 17O spans almost 1000 ppm in organic molecules and it serves as a potentially unique reporter of hydrogen bonding in peptides, nucleic acids, and water, and as a valuable complement to 13C and 15N NMR. Recent developments including the multiple-quantum magic-angle spinning (MQMAS) experiment have enabled the detection of 17O in biological solids, but very long data acquisitions are required to achieve sufficient sensitivity and resolution. Here, we perform nonlinear sampling in the indirect dimension of MQMAS experiments to substantially reduce the total acquisition time and improve sensitivity and resolution. Nonlinear sampling prevents the use of the discrete Fourier transform; instead, we employ maximum entropy (MaxEnt) reconstruction. Nonlinearly sampled MQMAS spectra are shown to provide high resolution and sensitivity in several systems, including lithium sulfate monohydrate (LiSO(4)-H(2)17O) and L-asparagine monohydrate (H(2)17O). The combination of nonlinear sampling and MaxEnt reconstruction promises to make the application of 17O MQMAS practical in a wider range of biological systems.

Chemical shifts in amino acids, peptides, and proteins: from quantum chemistry to drug design

This chapter discusses recent progress in the investigation and use of (13)C, (15)N, and (19)F nuclear magnetic resonance (NMR) chemical shifts and chemical shift tensors in proteins and model systems primarily using quantum chemical (ab initio Hartree-Fock and density functional theory) techniques. Correlations between spectra and structure are made and the techniques applied to other spectroscopic and electrostatic properties as well, including hydrogen bonding, ligand binding to heme proteins, J-couplings, electric field gradients, and atoms-in-molecules theory, together with a brief review of the use of NMR chemical shifts in drug design.

Ab initio conformational studies on diols and binary diol-water systems using DFT methods. Intramolecular hydrogen bonding and 1:1 complex formation with water

Studies on the conformational equilibrium for the following diols, ethane-1,2-diol (12EG, CAS 107-21-1), 2R-D-(-)-propane-1,2-diol (12PG, CAS 4254-14-2), (2S,3S)-L-(+)-butane-2,3-diol (L23BD, CAS 19132-06-0), and (2S,3R)-meso-butane-2,3-diol (m23BD, CAS 5341-95-7), are described using Gaussian ab initio calculations involving density functional theory (DFT) methods. We also report in this article results on the stability and conformation for the 1:1 water-diol complex formed by ethane-1,2-diol, propane-1,2-diol, and L- and meso-butane-2,3-diol. The relative stability of the intramolecular (internal) hydrogen bond in a range of diols (n = 2 to 6), based on ab initio geometry optimization and determination of the -O...H- distance, dOH, and -O-H...O- angle, theta, increases through the sequence 1,2 approximately equals 2,3 < 1,3 < 1,4 approximately equals 1,5 approximately equals 1,6, as judged from the bond linearity and -O...H- separation. Quantum mechanical and topological analysis of possible intramolecular hydrogen bonding in this complete series of diols provides convincing evidence for this in diols in which the hydroxyl groups are separated by three or more carbon atoms, that is, in (n, n+m) diols for m > or = 2, but not for ethane-1,2-diol or other vicinal diols, which do not satisfy Popelier's topological and electron density criteria based on the AIM theory of Bader. Based on these criteria it is unlikely that vicinal diols are in fact capable of forming an intramolecular hydrogen bond, in spite of geometric and spectroscopic data in the literature suggesting otherwise.

The hydrogen bond in the solid state

The hydrogen bond is the most important of all directional intermolecular interactions. It is operative in determining molecular conformation, molecular aggregation, and the function of a vast number of chemical systems ranging from inorganic to biological. Research into hydrogen bonds experienced a stagnant period in the 1980s, but re-opened around 1990, and has been in rapid development since then. In terms of modern concepts, the hydrogen bond is understood as a very broad phenomenon, and it is accepted that there are open borders to other effects. There are dozens of different types of X-H.A hydrogen bonds that occur commonly in the condensed phases, and in addition there are innumerable less common ones. Dissociation energies span more than two orders of magnitude (about 0.2-40 kcal mol(-1)). Within this range, the nature of the interaction is not constant, but its electrostatic, covalent, and dispersion contributions vary in their relative weights. The hydrogen bond has broad transition regions that merge continuously with the covalent bond, the van der Waals interaction, the ionic interaction, and also the cation-pi interaction. All hydrogen bonds can be considered as incipient proton transfer reactions, and for strong hydrogen bonds, this reaction can be in a very advanced state. In this review, a coherent survey is given on all these matters.

New developments in applying quantum mechanics to proteins

Algorithmic improvements of quantum mechanical methodologies have increased our ability to study the electronic structure of fragments of a biomolecule (e.g. an enzyme active site) or entire biomolecules. Three main strategies have emerged as ways in which quantum mechanics can be applied to biomolecules. The supermolecule approach continues to be utilized, but it is slowly being replaced by the so-called coupled quantum mechanical/molecular mechanical methodologies. An exciting new direction is the continued development and application of linear-scaling quantum mechanical approaches to biomolecular systems.