Environment influences on the aromatic character of nucleobases and amino acids

Mechanistic insights into the conformational switch in profilin-1 subject to collective effects of mutation and histidine tautomerism

Mutations and histidine (His) tautomerism in profilin-1 (PFN1) are associated with amyotrophic lateral sclerosis (ALS). The conformational changes in PFN1 caused by the collective effects of G117V mutation and His tautomeric isomers εε, εδ, δε, and δδ were clarified using molecular dynamics (MD) simulations. The predominant structural variations were seen in α-helices, β-sheets, turns, and coils and the His tautomer's unique degree of disruption was seen in these conformations. The content of α-helices was 23.2 % in the εε and δδ isomers, but the observed α-helices in the isomers εδ and δε were 20.3 % and 21.7 % respectively. The percentage of β-sheet was found to be higher (34.1) in the εε isomer than in the εδ, δε, and δδ isomers, and the values were 30.4, 29.7, and 31.9, respectively. Intermolecular water dynamics analysis discloses that His 133 can form an intramolecular H-bond interaction (N^α-H---N^δ), confirming the experimental observations in the simulations of εε, δε, and δδ isomers of G117V PFN1 mutant. It was concluded that these solvent molecules are crucial for aggregation and must be considered in future research on the PFN1 associated with ALS. Overall, the study offers a thorough microscopic understanding of the pathogenic mechanisms behind conformational changes that cause aggregation illnesses like ALS.

Histidine tautomerism dependent conformational transitions driven aggregation of profilin-1: Implications in amyotrophic lateral sclerosis

Aggregation of profilin-1 (PFN1) causes a fatal neurodegenerative disease, familial amyotrophic lateral sclerosis (fALS). Histidine (His) tautomerism has been linked to the formation of fibril aggregation causing neurodegenerative disease. Characterization of intermediate species that form during aggregation is crucial, however, this has proven very challenging for experimentalists due to their transient nature. Hence, molecular dynamics (MD) simulations have been performed on the His tautomeric isomers εε, εδ, δε, and δδ of PFN1 to explain the structural changes and to correlate them with its aggregation propensity. MD simulations show that His133 presumably plays a major role in the aggregation of PFN1 upon His tautomerism compared to His119. Further, the formation of a new 3₁₀-helix is observed in εε and δε but 3₁₀-helix is not observed in δδ and εδ isomers. In addition, our findings unveil that β-sheet dominating conformations are observed in His119(δ)-His133(δ) δδ isomer of PFN1 with significant antiparallel β-sheets between residues T15-G23, S29-A33, L63-L65, Q68-S76, F83-T89, T97-T105, and K107-K115, suggesting a novel aggregation mechanism possibly occur for the formation of PFN1 aggregates. Overall, these results propose that MD simulations of PFN1 His tautomers can provide a detailed microscopic understanding of the aggregation mechanisms which are hard to probe through experiments.

Effect of Alkali Metal Cations on Length and Strength of Hydrogen Bonds in DNA Base Pairs

For many years, non-covalently bonded complexes of nucleobases have attracted considerable interest. However, there is a lack of information about the nature of hydrogen bonding between nucleobases when the bonding is affected by metal coordination to one of the nucleobases, and how the individual hydrogen bonds and aromaticity of nucleobases respond to the presence of the metal cation. Here we report a DFT computational study of nucleobase pairs interacting with alkali metal cations. The metal cations contribute to the stabilization of the base pairs to varying degrees depending on their position. The energy decomposition analysis revealed that the nature of bonding between nucleobases does not change much upon metal coordination. The effect of the cations on individual hydrogen bonds were described by changes in VDD charges on frontier atoms, H-bond length, bond energy from NBO analysis, and the delocalization index from QTAIM calculations. The aromaticity changes were determined by a HOMA index.

Microsolvation of Histidine—A Theoretical Study of Intermolecular Interactions Based on AIM and SAPT Approaches

Histidine is unique among amino acids because of its rich tautomeric properties. It participates in essential enzymatic centers, such as catalytic triads. The main aim of the study is the modeling of the change of molecular properties between the gas phase and solution using microsolvation models. We investigate histidine in its three protonation states, microsolvated with 1:6 water molecules. These clusters are studied computationally, in the gas phase and with water as a solvent (Polarizable Continuum Model, PCM) within the Density Functional Theory (DFT) framework. The structural analysis reveals the presence of intra- and intermolecular hydrogen bonds. The Atoms-in-Molecules (AIM) theory is employed to determine the impact of solvation on the charge flow within the histidine, with emphasis on the similarity of the two imidazole nitrogen atoms—topologically not equivalent, they are revealed as electronically similar due to the heterocyclic ring aromaticity. Finally, the Symmetry-Adapted Perturbation Theory (SAPT) is used to examine the stability of the microsolvation clusters. While electrostatic and exchange terms dominate in magnitude over polarization and dispersion, the sum of electrostatic and exchange term is close to zero. This makes polarization the factor governing the actual interaction energy. The most important finding of this study is that even with microsolvation, the polarization induced by the presence of implicit solvent is still significant. Therefore, we recommend combined approaches, mixing explicit water molecules with implicit models.

Nucleotide conjugated (ZnO)3 cluster: Interaction and optical characteristics using TDDFT

Binding of four DNA nucleotide units with (ZnO)₃ cluster in an aqueous phase has been investigated using density functional theory (DFT) and time dependent-density functional theory (TDDFT) method and the stability order for (ZnO)₃-nucleobases/sugar/phosphate systems is predicted as phosphate > C > A > S > T ∼ G. The order of binding energy for (ZnO)₃-nucleotide hybrid systems is observed to be (ZnO)₃ + nuc-C ˃ (ZnO)₃ + nuc-A ˃ (ZnO)₃ + nuc-G ˃ (ZnO)₃ + nuc-T. The binding of nucleotide units with the cluster has been explained on the basis of molecular electrostatic potential (MEP) plots, hydrogen bonding, glycosidic torsion angles, density of state (DOS) plots. The photophysical properties of (ZnO)₃-nucleotide complexes have been studied using TDDFT approach. Among all (ZnO)₃-nucleotide complexes, the absorption spectra of (ZnO)₃ + nuc-A and (ZnO)₃ + nuc-C complexes are seen to undergo red shift with respect to their bare nucleotide units that would be useful in the optical sensing of the respective nucleotides of DNA. It is interesting to note that binding of the nucleotide unit with the cluster makes it fluorescent, the study reports the fluorescence activity of (ZnO)₃ + nuc-T complex.

Interaction of nucleobases with silicon doped and defective silicon doped graphene and optical properties

The interaction of nucleobases (NBs) with the surface of silicon doped graphene (SiGr) and defective silicon doped graphene (dSiGr) has been studied using electronic structure methods.

π-Cooperativity effect on the base stacking interactions in DNA: is there a novel stabilization factor coupled with base pairing H-bonds?

The results from absolutely localized molecular orbital (ALMO)-energy decomposition analysis (EDA) and ALMO-charge transfer analysis (CTA) at M06-2X/cc-pVTZ level reveal that double-proton transfer (DPT) reactions through base pairing H-bonds have nonignorable effects on the stacking energies of dinucleotide steps, which introduces us to a novel stabilization (or destabilization) factor in the DNA duplex. Thus, intra- and inter-strand base stacking interactions are coalesced with each other mediated by H-bridged quasirings between base pairs. Changes in stacking energies of dinucleotide steps depending on the positions of H atoms are due to variations in local aromaticities of individual nucleobases, manifesting π-cooperativity effects. CT analyses show that dispersion forces in dinucleotide steps can lead to radical changes in the redox properties of nucleobases, in particular those of adenine and guanine stacked dimers in a strand. Besides Watson-Crick rules, novel base pairing rules were propounded by considering CT results. According to these, additional base pairing through π-stacks of nucleobases in dinucleotide steps does not cause any intrinsic oxidative damage to the associated nucleobases throughout DPT.

Structural, electronic and energetic consequences of epigenetic cytosine modifications

The hydrogen bonding patterns of cytosine and its seven C5-modifed analogues paired with canonical guanine were studied using the first principle approach. Both global minima and biologically relevant conformations were studied. The former resulted from full gradient geometry optimizations of hydrogen bonded pairs, while the latter were obtained based on 125 d(GpC) dinucleotides found in the PDB database. The obtained energetic, electronic and structural data lead to the conclusion that the epigenetically relevant modification of cytosine may have serious consequences on hydrogen bonding with guanine. First of all, the significant substituent effects were observed for such trends as charges on sites involved in hydrogen bonding, the total intermolecular interaction energy or electron densities at bond critical points. Moreover, the molecular orbital polarization contribution resulting from energy decomposition expressed in terms of absolutely localized molecular orbitals exhibited an inverse linear correlation with frozen density contributions. A substituent effect on the amount of charge transfer from pyrimidine toward guanine was also observed. The increase of intermolecular interactions of guanine with modified cytosine is associated with the increase of the electro-donating character of the C5-substituent. However, only pairs involving 5-methylcytosine are more stable than those formed by canonical cytosine. Furthermore, the energy differences observed for global minima also remain important for a broad range of displacement and angular parameters defining pair conformations in model d(GpC) dinucleotides. Due to the sensitivities of intermolecular interactions to mutual arrangements of monomers the modification of cytosine at the C5 site can significantly alter the actual energy profiles. Consequently, it may be anticipated that the modified dinucleotides will adopt different conformations than a standard G-C pair in a B-DNA double helix.

Pressure-imposed changes of benzoic acid crystals

Structural and energetic properties of benzoic acid crystals at pressure elevated from ambient condition up to 2.21 GPa were characterized. The directly observed variations of cell parameters and consequently cell volume are associated with many other changes including energetic, geometric, and electronic characteristics. First of all the non-monotonous change of lattice energy are noticed with the rise of pressure since the increase of stabilization up to 1GPa is followed by systematic decrease of lattice energies after extending the hydrostatic compression. There is also an observed increase of C₂²(8) synthon stabilization interaction with increase of pressure. The lattice response rather than interaction within synthons are source of observed pressure-related trend of lattice energy changes. The energy decomposition analysis revealed that the total steric interactions determine the overall trend of lattice energy change with the rise of pressure. Besides geometric aromaticity index was used as a measure of geometric changes. Serious discrepancies were noticed between HOMA values computed with the use of experimental and optimized geometries of the ring. Even inclusion of uncertainties of experimental geometries related to limited precision of X-ray diffraction measurements does not cancel mentioned discrepancies. Although HOMA exhibit similar trends at modest pressures the diversity became surprisingly high at more extreme conditions. This might suggest limitations of periodic DFT computations at elevated pressures and the experimentally observed breaking of molecules at very high pressures will probably not be accounted properly in this approach. Also limitation of direct use of experimental geometries were highlighted.

Geometric and energetic consequences of prototropy for adenine and its structural models – a review

Prototropy for adenine and its convenient models causes parallel changes of geometric (HOMED) and energetic (ΔE) parameters for neutral tautomers.

Geometric consequences of electron delocalization for adenine tautomers in aqueous solution

Geometric consequences of electron delocalization were studied for all possible adenine tautomers in aqueous solution by means of ab initio methods {PCM(water)//DFT(B3LYP)/6-311+G(d,p)} and compared to those in the gas phase {DFT(B3LYP)/6-311+G(d,p)}. To measure the consequences of any type of resonance conjugation (π-π, n-π, and σ-π), the geometry-based harmonic oscillator model of electron delocalization (HOMED) index, recently extended to the isolated (DFT) and hydrated (PCM//DFT) molecules, was applied to the molecular fragments (imidazole, pyrimidine, 4-aminopyrimidine, and purine) and also to the whole tautomeric system. For individual tautomers, the resonance conjugations and consequently the bond lengths strongly depend on the position of the labile protons. The HOMED indices are larger for tautomers (or their fragments) possessing the labile proton(s) at the N rather than C atom. Solvent interactions with adenine tautomers slightly increase the resonance conjugations. Consequently, they slightly shorten the single bonds and lengthen the double bonds. When going from the gas phase to water solution, the HOMED indices increase (by less than 0.15 units). There is a good relation between the HOMED indices estimated in water solution and those in the gas phase for the neutral and ionized forms of adenine. Subtle effects, being a consequence of intramolecular interactions between the neighboring groups, are so strongly reduced by solvent that the relation between the HOMED indices and the relative energies for the neutral adenine tautomers seems to be better in water solution than in the gas phase.

Effect of the H-bonding on aromaticity of purine tautomers

Four tautomers of purine (1-H, 3-H, 7-H, and 9-H) and their equilibrium H-bonded complexes with F(-) and HF for acidic and basic centers, respectively, were optimized by means of the B3LYP/6-311++G(d,p) level of theory. Purine tautomer stability increases in the following series: 1-H < 3-H < 7-H < 9-H, consistent with increasing aromaticity. Furthermore, the presence of a hydrogen bond with HF does not change this order. For neutral H-bonded complexes, the strongest and the weakest intermolecular interactions occur (-14.12 and -10.49 kcal/mol) for less stable purine tautomers when the proton acceptor is located in the five- and six-membered rings, respectively. For 9-H and 7-H tautomers the order is reversed. The H-bond energy for the imidazole complex with HF amounts to -14.03 kcal/mol; hence, in the latter case, the fusion of imidazole to pyrimidine decreases its basicity. The ionic H-bonds of N(-)···HF type are stronger by ~10 kcal/mol than the neutral N···HF intermolecular interactions. The hydrogen bond N(-)···HF energies in pyrrole and imidazole are -32.28 and -30.03 kcal/mol, respectively, and are substantially stronger than those observed in purine complexes. The aromaticity of each individual ring and of the whole molecule for all tautomers in ionic complexes is very similar to that observed for the anion of purine. This is not the case for neutral complexes and purine as a reference. The N···HF bonds perturb much more the π-electron structure of five-membered rings than that of the six-membered ones. The H-bonding complexes for 7-H and 9-H tautomers are characterized by higher aromaticity and a much lower range of HOMA variability.