An atoms-in-molecules study of the genetically-encoded amino acids: I. Effects of conformation and of tautomerization on geometric, atomic, and bond properties

The theory of Atoms-In-Molecules (AIM) is a partitioning of the real space of a molecule into disjoint atomic constituents as determined by the topology of the electron density, rho(r). This theory identifies an atom in a molecule with a quantum mechanical open system and, consequently, all of the atom's properties are unambiguously defined. AIM recovers the basic empirical cornerstone of chemistry: that atoms and functional groups possess characteristic and additive properties that in many cases exhibit a remarkable transferability between different molecules. As a result, the theory enables the theoretical synthesis of a large molecule and the prediction of its properties by joining fragments that are predetermined as open systems. The present article is the first of a series (in preparation) that explore this possibility for polypeptides by determining the transferability of the building blocks: the amino acid residues. Transferability of group properties requires transferability of the electron density rho(r), which in turn requires the transferability of the geometric parameters. This article demonstrates that these parameters are conformation-insensitive for a representative amino acid, leucine, and that the atomic and bond properties exhibit a corresponding transferability. The effects of hydrogen bonding are determined and a set of geometrical conditions for the occurrence of such bonding is identified. The effects of transforming neutral leucine into its zwitter-ionic form on its atomic and bond properties are shown to be localized primarily to the sites of ionization.     

Atoms-in-molecules study of the genetically encoded amino acids. III. Bond and atomic properties and their correlations with experiment including mutation-induced changes in protein stability and genetic coding

This article presents a study of the molecular charge distributions of the genetically encoded amino acids (AA), one that builds on the previous determination of their equilibrium geometries and the demonstrated transferability of their common geometrical parameters. The properties of the charge distributions are characterized and given quantitative expression in terms of the bond and atomic properties determined within the quantum theory of atoms-in-molecules (QTAIM) that defines atoms and bonds in terms of the observable charge density. The properties so defined are demonstrated to be remarkably transferable, a reflection of the underlying transferability of the charge distributions of the main chain and other groups common to the AA. The use of the atomic properties in obtaining an understanding of the biological functions of the AA, whether free or bound in a polypeptide, is demonstrated by the excellent statistical correlations they yield with experimental physicochemical properties. A property of the AA side chains of particular importance is the charge separation index (CSI), a quantity previously defined as the sum of the magnitudes of the atomic charges and which measures the degree of separation of positive and negative charges in the side chain of interest. The CSI values provide a correlation with the measured free energies of transfer of capped side chain analogues, from the vapor phase to aqueous solution, yielding a linear regression equation with r2 = 0.94. The atomic volume is defined by the van der Waals isodensity surface and it, together with the CSI, which accounts for the electrostriction of the solvent, yield a linear regression (r2 = 0.98) with the measured partial molar volumes of the AAs. The changes in free energies of transfer from octanol to water upon interchanging 153 pairs of AAs and from cyclohexane to water upon interchanging 190 pairs of AAs, were modeled using only three calculated parameters (representing electrostatic and volume contributions) yielding linear regressions with r2 values of 0.78 and 0.89, respectively. These results are a prelude to the single-site mutation-induced changes in the stabilities of two typical proteins: ubiquitin and staphylococcal nuclease. Strong quadratic correlations (r2 approximately 0.9) were obtained between DeltaCSI upon mutation and each of the two terms DeltaDeltaH and TDeltaDeltaS taken from recent and accurate differential scanning calorimetry experiments on ubiquitin. When the two terms are summed to yield DeltaDeltaG, the quadratic terms nearly cancel, and the result is a simple linear fit between DeltaDeltaG and DeltaCSI with r2 = 0.88. As another example, the change in the stability of staphylococcal nuclease upon mutation has been fitted linearly (r2 = 0.83) to the sum of a DeltaCSI term and a term representing the change in the van der Waals volume of the side chains upon mutation. The suggested correlation of the polarity of the side chain with the second letter of the AA triplet genetic codon is given concrete expression in a classification of the side chains in terms of their CSI values and their group dipole moments. For example, all amino acids with a pyrimidine base as their second letter in mRNA possess side-chain CSI < or = 2.8 (with the exception of Cys), whereas all those with CSI > 2.8 possess an purine base. The article concludes with two proposals for measuring and predicting molecular complementarity: van der Waals complementarity expressed in terms of the van der Waals isodensity surface and Lewis complementarity expressed in terms of the local charge concentrations and depletions defined by the topology of the Laplacian of the electron density. A display of the experimentally accessible Laplacian distribution for a folded protein would offer a clear picture of the operation of the "stereochemical code" proposed as the determinant in the folding process.     

Hybridization-displaced charges for amino-acids: a new model using two point charges per atom along with bond-center charges

A new charge distribution is proposed for the amino acids where each atom is associated with two point charges while each bond center is associated with one point charge. Centroids of charges arising due to atomic orbital hybridization called hybridization-displaced charges (HDC) and those located at the atomic sites and bond centers obtained by a modified form of the Mulliken scheme were combined. The density matrix calculations required for this analysis were performed at the B3LYP/6-31G** level of density functional theory. The combination of HDC centroids with the modified Mulliken charges was found to yield dipole moments and surface molecular electrostatic potentials (MEP) of the amino acids in good agreement with those obtained by rigorous DFT calculations or those obtained using the MEP-fitted CHelpG charges. This study shows that the combination of HDC centroids with the modified Mulliken charges is significantly superior to the conventional Mulliken charges.     

Topological analysis of aromatic halogen/hydrogen bonds by electron charge density and electrostatic potentials

In this work, the intermolecular distribution of the electronic charge density in the aromatic hydrogen/halogen bonds is studied within the framework of the atoms in molecules (AIM) theory and the molecular electrostatic potentials (MEP) analysis. The study is carried out in nine complexes formed between benzene and simple lineal molecules, where hydrogen, fluorine and chlorine atoms act as bridge atoms. All the results are obtained at MP2 level theory using cc-pVTZ basis set. Attention is focused on topological features observed at the intermolecular region such as bond, ring and cage critical points of the electron density, as well as the bond path, the gradient of the density maps, molecular graphs and interatomic surfaces. The strength of the interaction increases in the following order: F⋅⋅⋅π < Cl⋅⋅⋅π < H⋅⋅⋅π. Our results show that the fluorine atom has the capability to interact with the π−cloud to form an aromatic halogen bond, as long as the donor group is highly electron withdrawing. The Laplacian topology allows us to state that the halogen atoms can act as nucleophiles as well as electrophiles, showing clearly their dual character.     

Topological description of the bond-breaking and bond-forming processes of the alkene protonation reaction in zeolite chemistry: an AIM study

Density functional theory and atoms in molecules theory were used to study bond breakage and bond formation in the trans-2-butene protonation reaction in an acidic zeolitic cluster. The progress of this reaction along the intrinsic reaction coordinate, in terms of several topological properties of relevant bond critical points and atomic properties of the key atoms involved in these concerted mechanisms, were analyzed in depth. At B3LYP/6-31++G(d,p)//B3LYP/6-31G(d,p) level, the results explained the electron density redistributions associated with the progressive bond breakage and bond formation of the reaction under study, as well as the profiles of the electronic flow between the different atomic basins involved in these electron reorganization processes. In addition, we found a useful set of topological indicators that are useful to show what is happening in each bond/atom involved in the reaction site as the reaction progresses.     

On contribution of known atomic partial charges of protein backbone in electrostatic potential density maps

Partial charges of atoms in a molecule and electrostatic potential (ESP) density for that molecule are known to bear a strong correlation. In order to generate a set of point‐field force field parameters for molecular dynamics, Kollman and coworkers have extracted atomic partial charges for each of all 20 amino acids using restrained partial charge‐fitting procedures from theoretical ESP density obtained from condensed‐state quantum mechanics. The magnitude of atomic partial charges for neutral peptide backbone they have obtained is similar to that of partial atomic charges for ionized carboxylate side chain atoms. In this study, the effect of these known atomic partial charges on ESP is examined using computer simulations and compared with the experimental ESP density recently obtained for proteins using electron microscopy. It is found that the observed ESP density maps are most consistent with the simulations that include atomic partial charges of protein backbone. Therefore, atomic partial charges are integral part of atomic properties in protein molecules and should be included in model refinement.     

Question 9: Quantum Self-Assembly and Photoinduced Electron Tunneling in Photosynthetic Systems of Artificial Minimal Living Cells

Natural and artificial living cells and their substructures are self-assembling, due to electron correlation interactions among biological and water molecules, which lead to attractive dispersion forces and hydrogen bonds. Dispersion forces are weak intermolecular forces that arise from the attractive force between quantum multipoles. A hydrogen bond is a special type of quantum attractive interaction that exists between an electronegative atom and a hydrogen atom bonded to another electronegative atom; and this hydrogen atom exist in two quantum states. The best method to simulate these dispersion forces and hydrogen bonds is to perform quantum mechanical non-local density functional potential calculations of artificial minimal living cells consisting of around 1,000 atoms. The cell systems studied are based on peptide nucleic acid and are 3.0–4.2 nm in diameter. The electron tunneling and associated light absorption of the most intense transitions, as calculated by the time dependent density functional theory method, differs from spectroscopic experiments by only 0.2–0.3 nm, which is within the value of experiment errors. This agreement implies that the quantum mechanically self-assembled structures of artificial minimal living cells very closely approximate realistic ones.     

Theoretical study of the interaction between 5-methylcytosine and acrylamide

The hydrogen-bonded complexes between 5-methylcytosine and acrylamide have been investigated using the density function theory (DFT) method. Five stable complexes have been found with no imaginary frequencies. Complex C3 is the most stable one with interaction energies of -69.01?kJ?mol(-1) corrected for basis set superposition error (BSSE). The charge change in the process of these complexes formation has also been examined. The atoms in molecules (AIM) theory and natural bond orbital (NBO) method have been performed to investigate the hydrogen bonds involved in all the complexes. The electron density and its corresponding Laplacian at the bond and ring critical points have been analyzed. In C3 complex, there is the largest stabilization energy (18.17?kJ?mol(-1)) between N11-H12 antibonding orbital and lone electron pair of O17. It can be seen that the hydrogen bonds play a crucial role in the stability of all the complexes between 5-methylcytosine and acrylamide. The theoretical results could provide helpful information for other researchers in further work.     

Characteristics of beryllium bonds; a QTAIM study

The nature of beryllium bonds formed between BeX2 (X is H, F and Cl) and some Lewis bases have been investigated. The distribution of the Laplacian of electron density shows that there is a region of charge depletion around the Be atom, which, according to Laplacian complementary principal, can interact with a region of charge concentration of an atom in the base and form a beryllium bond. The molecular graphs of the investigated complexes indicate that beryllium in BeH2 and BeF2 can form “beryllium bonds” with O, N and P atoms but not with halogens. In addition, eight criteria based on QTAIM properties, including the values of electron density and its Laplacian at the BCP, penetration of beryllium and acceptor atom, charge, energy, volume and first atomic moment of beryllium atom, have been considered and compared with the corresponding ones in conventional hydrogen bonds. These bonds share many common features with very strong hydrogen bonds, however,some differences have also been observed.     

Effects of halogen substitution on Watson–Crick base pairing: A possible mechanism for radiosensitivity

The halogen substituent effect on geometries and charge distributions of the A–T base pair derivatives was evaluated using density functional theory at B3LYP/6-31G1 level. The results indicate that all of the substitutions affect geometries and charge distributions of the atoms contributing hydrogen bonds. These changes would be the reason of the radiosensitization of these compounds incorporating DNA.     

Topological analysis of tetraphosphorus oxides (P4O6+n (n = 0–4))

Quantum chemical calculations were used to analyze the chemical bonding and the reactivity of phosphorus oxides (P₄O_6+n (n?=?0–4)). The chemical bonding was studied using topological analysis such as atoms in molecules (AIM), electron localization function (ELF), and the reactivity using the Fukui function. A classification of the P-O bonds formed in all structures was done according to the coordination number in each P and O atoms. It was found that there are five P-O bond types and these are distributed among the five phosphorus oxides structures. Results showed that there is good agreement among the evaluated properties (length, bond order, density at the critical point, and disynaptic population) and each P-O bond type. It was found that regardless of the structure in which a P-O bond type is present the topological and geometric properties do not have a significant variation. The topological parameters electron density and Laplacian of electron density show excellent linear correlation with the average length of P-O bond in each bond type for each structure. From the Fukui function analysis it was possible to predict that from P₄O₆ until P₄O₈ the most reactive regions are basins over the P.     

Crystal structure of Bordetella pertussis BugD solute receptor unveils the basis of ligand binding in a new family of periplasmic binding proteins

Periplasmic binding proteins of a new family particularly well represented in Bordetella pertussis have been called Bug receptors. One B.pertussis Bug protein is part of a tripartite tricarboxylate transporter while the functions of the other 77 are unknown. We report the first structure of a Bug receptor, BugD. It adopts the characteristic Venus flytrap motif observed in other periplasmic binding proteins, with two globular domains bisected by a deep cleft. BugD displays a closed conformation resulting from the fortuitous capture of a ligand, identified from the electron density as an aspartate. The structure reveals a distinctive alpha carboxylate-binding motif, involving two water molecules that bridge the carboxylate oxygen atoms to the protein. Both water molecules are hydrogen bonded to a common carbonyl group from Ala14, and each forms a hydrogen bond with one carboxylate oxygen atom of the ligand. Additional hydrogen bonds are found between the ligand alpha carboxylate oxygen atoms and protein backbone amide groups and with a threonine hydroxyl group. This specific ligand-binding motif is highly conserved in Bug proteins, indicating that they may all be receptors of amino acids or other carboxylated solutes, with a similar binding mode. The present structure thus unveils the bases of ligand binding in this large family of periplasmic binding proteins, several hundred members of which have been identified in various bacterial species.     

The occurrence of C--H...O hydrogen bonds in alpha-helices and helix termini in globular proteins

A comprehensive database analysis of C--H...O hydrogen bonds in 3124 alpha-helices and their corresponding helix termini has been carried out from a nonredundant data set of high-resolution globular protein structures resolved at better than 2.0 A in order to investigate their role in the helix, the important protein secondary structural element. The possible occurrence of 5 --> 1 C--H...O hydrogen bond between the ith residue CH group and (i - 4)th residue C==O with C...O < or = 3.8 A is studied, considering as potential donors the main-chain Calpha and the side-chain carbon atoms Cbeta, Cgamma, Cdelta and Cepsilon. Similar analysis has been carried out for 4 --> 1 C--H...O hydrogen bonds, since the C--H...O hydrogen bonds found in helices are predominantly of type 5 --> 1 or 4 --> 1. A total of 17,367 (9310 of type 5 --> 1 and 8057 of type 4 --> 1) C--H...O hydrogen bonds are found to satisfy the selected criteria. The average stereochemical parameters for the data set suggest that the observed C--H...O hydrogen bonds are attractive interactions. Our analysis reveals that the Cgamma and Cbeta hydrogen atom(s) are frequently involved in such hydrogen bonds. A marked preference is noticed for aliphatic beta-branched residue Ile to participate in 5 --> 1 C--H...O hydrogen bonds involving methylene Cgamma 1 atom as donor in alpha-helices. This may be an enthalpic compensation for the greater loss of side-chain conformational entropy for beta-branched amino acids due to the constraint on side-chain torsion angle, namely, chi1, when they occur in helices. The preference of amino acids for 4 --> 1 C--H...O hydrogen bonds is found to be more for Asp, Cys, and for aromatic residues Trp, Phe, and His. Interestingly, overall propensity for C--H...O hydrogen bonds shows that a majority of the helix favoring residues such as Met, Glu, Arg, Lys, Leu, and Gln, which also have large side-chains, prefer to be involved in such types of weak attractive interactions in helices. The amino acid side-chains that participate in C--H...O interactions are found to shield the acceptor carbonyl oxygen atom from the solvent. In addition, C--H...O hydrogen bonds are present along with helix stabilizing salt bridges. A novel helix terminating interaction motif, X-Gly with Gly at C(cap) position having 5 --> 1 Calpha--H...O, and a chain reversal structural motif having 1 --> 5 Calpha-H...O have been identified and discussed. Our analysis highlights that a multitude of local C--H...O hydrogen bonds formed by a variety of amino acid side-chains and Calpha hydrogen atoms occur in helices and more so at the helix termini. It may be surmised that the main-chain Calpha and the side-chain CH that participate in C--H...O hydrogen bonds collectively augment the cohesive energy and thereby contribute together with the classical N--H...O hydrogen bonds and other interactions to the overall stability of helix and therefore of proteins.     

Patterns in ionizable side chain interactions in protein structures

In a selected set of 44 high-resolution, non-homologous protein structures, the intramolecular hydrogen bonds or salt bridges formed by ionizable amino acid side chains were identified and analyzed. The analysis was based on the investigation of several properties of the involved residues such as their solvent exposure, their belonging to a certain secondary structural element, and their position relative to the N- and C-termini of their respective structural element. It was observed that two-thirds of the interactions made by basic or acidic side chains are hydrogen bonds to polar uncharged groups. In particular, the majority (78%) of the hydrogen bonds between ionizable side chains and main chain polar groups (sch:mch bonds) involved at least one buried atom, and in 42% of the cases both interacting atoms were buried. In α-helices, the sch:mch bonds observed in the proximity of the C- and N-termini show a clear preference for acidic and basic side chains, respectively. This appears to be due to the partial charges of peptide group atoms at the termini of α-helices, which establish energetically favorable electrostatic interactions with side chain carrying opposite charge, at distances even greater than 4.5 Å. The sch:mch interactions involving ionizable side chains that belong either to β-strands or to the central part of α-helices are based almost exclusively on basic residues. This results from the presence of main chain carbonyl oxygen atoms in the protein core which have unsatisfied hydrogen bonding capabilities.     

Theoretical investigations on the hydrogen bonding of nitrile isomers with H2O,HF, NH3 and H2S

The hydrogen bonds formed by the interaction of nitriles with water, hydrogen fluoride, ammonia and hydrogen sulphide have been studied using B3LYP and second-order Møller–Plesset perturbation (MP2) methods and 6-311+ + G(d,p) basis set. The energies and structures of 80 hydrogen-bonded complexes between nitriles and small molecules were examined systematically using B3LYP and MP2 procedure. Categorisation of the hydrogen bonds involved in the various complexes led to an ordering of hydrogen bond donor and acceptor abilities for some functional groups. The interaction energies have been corrected for the basis set superposition error using Boy's counterpoise correction method. The Morokuma energy decomposition analysis reveals that the strong interactions are due to the attractive contributions from the electrostatic (ES), polarisation (PL) and charge transfer (CT) components. The topological parameters, electron density and Laplacian of electron density show excellent correlation with the hydrogen bond length. Natural bond orbital (NBO) analysis has also been performed to study the CT from proton acceptor to the antibonding orbital of the H–Y bond in the proton donor part of complexes. The frequency analysis of C–H…Y bond in the complexes indicates the blue-shifting nature largely in case of sp² hybridised carbon atom.     

A geometry optimization and molecular electrostatic potential mapping study of structure-activity relationship for some anti-Alzheimer agents.

Molecular geometries of some substituted (pyrroloamino)pyridines which possess anti-Alzheimer activity were optimized and potential-derived CHelpG point charges were computed using ab initio SCF molecular orbital approach employing the 3-21G basis set. AM1 molecular orbital calculations were performed using these optimized geometries and thus optimized Hybridization. Displacement Charges (HDC) combined with L?wdin charges continuously distributed in three dimension were obtained. Molecular electrostatic potential (MEP) maps of the molecules were obtained in two ways: (i) using the HDC-based model with the help of which MEP minima near the molecules were located, and (ii) using the CHelpG point charges, MEP values on the van der Waals surfaces of the molecules were computed. The MEP maps computed using both the methods have negative MEP regions near the pyridine nitrogen atom which appears to be the main binding site of the molecules with the appropriate receptor. Both electrostatic interaction and lipophilic association between these molecules and the receptor appear to contribute to biological activity.     

Nearest-neighbor effects and structural preferences in dipeptides are a function of the electronic properties of amino acid side-chains

The electronic properties of amino acid side-chains are emerging as an important factor in the preference for secondary structure in proteins. These properties have not been fully characterized, nor has their role in the behavior of peptides been explored in any detail. The present studies sought to evaluate several possibilities: 1) that hydrophilicity can be expressed solely in electronic terms, 2) that substituent effects of side-chains extend across the peptide bond, and (3) nearest-neighbor effects in dipeptides correlate with secondary structural preferences. Quantum mechanics (QM) calculations were used to define the electronic properties of individual amino acids and dipeptides. It was found that the hydrophilicity of an amino acid side-chain can be accurately represented as a function of the electron densities of its component atoms. In addition, the nature of an amino acid in the second position of a dipeptide affects the electronic properties (Mulliken populations and electron densities) of the main-chain atoms of the first residue. Certain electronic features of the dipeptides strongly correlated with propensity for secondary structure. Specifically, Mulliken population data at the Calpha atom and N atom predicted preference for alpha-helices versus coil and strand conformations, respectively. Analysis of dipeptides arrayed in either helical or extended structures revealed lengthening of main-chain bonds in the alpha-helical conformations. A thorough characterization of the electronic properties of amino acids and short peptide segments may provide a better understanding of the forces that determine secondary structure in proteins.     

Quantum-chemical investigation of the structure and the antioxidant properties of α-lipoic acid and its metabolites

Quantum-chemical computations were used to investigate the structure-antioxidant parameter relationships of α-lipoic acid and its natural metabolites bisnorlipoic acid and tetranorlipoic acid in their oxidized and reduced forms. The enantiomers of lipoic and dihydrolipoic acid were optimized using the B3LYP/6-311+G(3df,2p), B3LYP/aug-cc-pVDZ and MP2(full)/6-31+G(d,p) levels of theory as isolated molecules and in the presence of water. The geometries of the metabolites and the values of their antioxidant parameters (proton affinity, bond dissociation enthalpy, adiabatic ionization potential, spin density, and the highest occupied molecular orbital energy) were calculated at the B3LYP/6-311+G(3df,2p) level of theory. The results obtained reveal similarities between these structures: a pentatomic, nonaromatic ring is present in the oxidized forms, while an unbranched aliphatic chain (as found in saturated fatty acids) is present in both the oxidized and the reduced forms. Analysis of the spin density and the highest occupied molecular orbital energy revealed that the SH groups exhibited the greatest electron-donating activities. The values obtained for the proton affinity, bond dissociation enthalpy and adiabatic ionization potential indicate that the preferred antioxidant mechanisms for α-lipoic acid and its metabolites are sequential proton loss electron transfer in polar media and hydrogen atom transfer in vacuum.     

Theoretical study of borazine and its derivatives

The density functional theory is used to study the geometries, electronic structures, and aromaticity of borazine and its fused ring derivatives. Some new evidences for the ionic nature of B-N bond are found. Geometry studies show that the B-N bond lengths are equal. The lone pair VSCCs of the N atoms are found. As shown, the B-N bonds are of ionic nature based on their positive Laplacian. Magnatic shielding constants also are computed. The shielding and deshielding contributions are divided into Lewis and non-Lewis parts by the NCS-NBO method. It is demonstrated in the NICS studies that there are the ring current effects on borazine and its derivatives are very weak. The aromaticity of borazine is weakened with the fused ring number increasing.     

Conformational and stereoelectronic investigation of tryptamine. An AIM/NBO study

Due to the free radical scavenger properties of Tryptamine (TRA), as well as of others indole derivatives, it is in our interest to explore deeply the stereoelectronic aspects that would be relevant in their stabilization and antioxidant activity. In this work the conformational space of TRA was scanned using molecular dynamics complemented with functional density calculations at B3LYP/6-31 + G** level. Twenty one conformers of lowest energy were obtained, their electronic distributions were analyzed at a higher calculation level, thus improving the basis set (B3LYP/6-311++G**). A topological study based on Bader's theory (AIM: atoms in molecules) and natural bond orbital (NBO) framework was performed. The study was enriched by a deep analysis of maps of molecular electrostatic potential (MEP) through a coordinated NBO/AIM analysis. The conformational preferences were explained by hyperconjugative interactions, which were revealed by NBO data. Because radical scavenging by indolic compounds is strongly modulated by their functional residues our study was related to similar analysis done previously on Indole and 1H-indole-3-acetic acid (IAA). Therefore, the conformational space of TRA was studied from a new perspective focusing on a deep analysis of the geometric and electronic properties of TRA conformers. The changes of the electronic distribution introduced by the substituent and the conformational flexibility of the side chain were addressed. The results reported contribute to the understanding of the structure, stability and reactivity of TRA and others indole derivatives.