Analysis of the binding selectivity and inhibiting mechanism of chlorogenic acid isomers and their interaction with grass carp endogenous lipase using multi-spectroscopic, inhibition kinetics and modeling methods

Polyphenols are inhibitors for lipase, but the binding selectivity and mechanism of polyphenol isomers and how they interact with lipase are not clear. Here, chlorogenic acid (CGA) isomers, neochlorogenic acid (NCGA) and cryptochlorogenic acid (CCGA) were used to explore the binding selectivity and mechanism of lipase. An inhibition assay indicated that both CGA isomers had dose-dependent inhibitory effects on lipase; however, the inhibitory effect of NCGA was better (IC₅₀: 0.647 mg/mL) than that of CCGA (IC₅₀: 0.677 mg/mL). NCGA and CCGA formed complexes with lipase at a molar ratio of 1:1, and the electrostatic interaction force plays a major role in the lipase-CCGA system. Molecular dynamics studies demonstrated that NCGA had a greater impact on the structure of lipase. The multi-spectroscopic and modeling results explained the effects of micro-structural changes on the binding site, the interaction force and the inhibition rate of the isomers when they combined with lipase.

In silico decryption of serotonin-receptor binding: local non-covalent interactions and long-range conformational changes

Serotonin-receptor binding is the key step in the process behind serotonin functionality, including several psychological and physiological behaviours. This study is focused on identifying the main non-covalent interactions controlling the stability of serotonin-receptor complexes as well as the main conformational changes in the receptor due to serotonin-receptor binding using classical molecular dynamics simulations and quantum chemical calculations. A qualitative analysis based on two order parameters ((i) the centre of mass distance and (ii) the angle between the surface normals of each aromatic residue and serotonin in the binding site) on the serotonin-receptor complex trajectory suggests that the T-type stacking interaction is predominant in the binding site. Quantum chemical calculations of the stacking interaction energy provide the quantitative contributions of important aromatic residues to the stabilization of the complex. Furthermore, a three body stacking interaction (named 'L'-type) was observed and likely contributes to the stability of the complex. Direct and water-mediated hydrogen bonding between the residues in the binding site and serotonin contributes to the complex stability. Principal component analysis of the molecular dynamics simulation trajectory of the serotonin-receptor complex and the apo-receptor in water indicates that the whole receptor is significantly stabilized due to serotonin binding. An analysis based on the dynamic cross correlation function reflects the strong correlation between trans-membrane (TM)3, TM5, TM6 (containing residues responsible for the stacking interaction and hydrogen bonding) and mini-G0 which may participate in signal transduction leading to the functionality of serotonin.

Adenine nucleobase directed supramolecular architectures based on ferrimagnetic heptanuclear copper(II) entities and benzenecarboxylate anions

Two planar organic anions, benzoate and benzene-1,4-dicarboxylate (terephthalate), have been selected as potential π-stacking intercalators among ferrimagnetic [Cu₇(μ-adeninato)₆(μ₃-OH)₆(μ-H₂O)₆]²⁺ heptameric discrete entities. The resulting supramolecular architecture is highly dependent on the negative charge density distribution, mainly located in the carboxylate groups of the organic anions. In this sense, the benzoate anion, with just one carboxylate group, does not allow its intercalation between the adeninato ligands as it would imply a high steric hindrance among the heptameric entities. As a consequence, these benzoate anions are located inside the voids of the crystal structure reducing the accessible volume of compound [Cu₇(μ-adeninato)₆(μ₃-OH)₆(μ-H₂O)₆](benzoate)₂·~17H₂O (1). On the contrary, the terephthalate anion, containing two carboxylate groups at opposite sites, adopts a π-stacking sandwich arrangement between two adeninato ligands that affords the porous open structure of formula [Cu₇(μ-adeninato)₆(μ₃-OH)₆(μ-H₂O)₆](terephthalate)·nH₂O (2a, 2b; n: 12 and 24, respectively). In addition to that, the less directional nature of the π-stacking interactions in comparison to the complementary hydrogen bonding based supramolecular metal-organic frameworks (SMOFs), suits them with a flexible architecture able to reversibly adsorb/desorb water (up to a 25-30% at 20 °C) altogether with the expansion/shrinkage of the crystal structure. The bridging adeninato and hydroxido ligands are effective magnetic exchange mediators to provide a S_T = 5/2 ferrimagnetic state for the heptanuclear entity.

DL-3-Aminoisobutyric acid: vibrational, NBO and AIM analysis of N-H⋯O bonded-zwitterionic dimer model

A zwitterionic dimer model constructed of inter-molecular -N-H⋯O bonding has been proposed for the solid sample of DL-3-Aminoisobutyric acid consistent with IR absorption and Raman spectral features measured in the 3500-400/50 cm^-1. This zwitterionic dimer model in water as solvent has been computed at B3LYP/6-311++G(d,p) and B3LYP-D3/6-311++G(d,p) levels including Grimme's dispersion correction associated with the -N-H⋯O interaction and SCRF-SMD method. Of the several possible monomer and dimer conformational structures, the most stable dimer constructed of two zwitterion monomer units has produced vibrational modes due to the -NH₃ ⁺ cation and -CO₂‾ anion involved in the -N-H⋯O bonding in fair agreement with the observed broad but composite IR modal features near the 3500-2000 cm^-1. Except for the frequency of asymmetric stretching mode of the -NH₃ ⁺ cation, its symmetric and bending modes agree with the observed values. As for the -CO₂‾ anion, the frequencies of all of its modes are in good agreement with the experiment. Natural bond orbital (NBO), molecular electrostatic potential (MEP), atoms-in-molecules (AIM) and non-covalent interaction (NCI) analyses have been used to understand electronic characterization of the -N-H⋯O bonding.

Biointeractions of Herbicide Atrazine with Human Serum Albumin: UV-Vis, Fluorescence and Circular Dichroism Approaches

The herbicide atrazine is widely used across the globe, which is a great concern. To investigate its potential toxicity in the human body, human serum albumin (HSA) was selected as a model protein. The interaction between atrazine and HSA was investigated using steady-state fluorescence spectroscopy, synchronous fluorescence spectroscopy, UV-Vis spectroscopy, three-dimensional (3D) fluorescence spectroscopy and circular dichroism (CD) spectroscopy. The intrinsic fluorescence of HSA was quenched by the atrazine through a static quenching mechanism. Fluorescence spectra at two excitation wavelengths (280 and 295 nm) showed that the fluorescence quenched in HSA was mainly contributed to by tryptophan residues. In addition, the atrazine bound to HSA, which induced changes in the conformation and secondary structure of HSA and caused an energy transfer. Thermodynamic parameters revealed that this binding is spontaneous. Moreover, electrostatic interactions play a major role in the combination of atrazine and HSA. One atrazine molecule can only bind to one HSA molecule to form a complex, and the atrazine molecule is bound at site II (subdomain IIIA) of HSA. This study furthers the understanding of the potential effects posed by atrazine on humans at the molecular level.

Using Density Functional Theory To Study Neutral and Ionized Stacked Thymine Dimers

Stacking interactions in thymine dimers are studied with density functional theory. According to our calculations, six dimers of comparable stability can be prepared at low temperature, but dimerization is impossible at room temperature due to the large entropy contribution that accompanies it. Analysis of vibrational anharmonic coupling terms shows that each of the dimers exhibits distinct vibrational dynamics. Properties of electron density in the intermolecular region are used to analyze neutral stacked species and their ionized forms. Bond paths and critical points in the intermolecular region are identified, but a simple relationship between binding energy and total electron density in the intermolecular critical points could not be found due to an uneven electron distribution in the binding region. The reduced density gradient was confirmed to be a useful tool for analysis of weak stacking interactions. Those interactions also affect vertical and adiabatic ionization energies, which are computed to be slightly lower for the dimers compared to the monomer.

Non-covalent interaction of benzene with methanol and diethyl ether solid surfaces

We have investigated the interactions involved at the interface of binary, layered ices (benzene on methanol and on diethyl ether) by means of laboratory experiments and ab initio calculations on model clusters.

Differential Editosome Protein Function between Life Cycle Stages of Trypanosoma brucei

Uridine insertion and deletion RNA editing generates functional mitochondrial mRNAs in Trypanosoma brucei. The mRNAs are differentially edited in bloodstream form (BF) and procyclic form (PF) life cycle stages, and this correlates with the differential utilization of glycolysis and oxidative phosphorylation between the stages. The mechanism that controls this differential editing is unknown. Editing is catalyzed by multiprotein ∼20S editosomes that contain endonuclease, 3'-terminal uridylyltransferase, exonuclease, and ligase activities. These editosomes also contain KREPB5 and KREPA3 proteins, which have no functional catalytic motifs, but they are essential for parasite viability, editing, and editosome integrity in BF cells. We show here that repression of KREPB5 or KREPA3 is also lethal in PF, but the effects on editosome structure differ from those in BF. In addition, we found that point mutations in KREPB5 or KREPA3 differentially affect cell growth, editosome integrity, and RNA editing between BF and PF stages. These results indicate that the functions of KREPB5 and KREPA3 editosome proteins are adjusted between the life cycle stages. This implies that these proteins are involved in the processes that control differential editing and that the 20S editosomes differ between the life cycle stages.

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.

Energy interactions in amyloid-like fibrils from NNQQNY

We use large-scale MP2 calculations to analyze the interactions appearing in amyloid fibers, which are difficult to determine experimentally. To this end, dimers and trimers of the hexapeptide NNQQNY from the yeast prion-like protein Sup35 were considered as model systems. We studied the energy interactions present in the three levels of organization in which the formation of amyloid fibrils is structured. The structural changes in the hydrogen bonds were studied too. It was found that the most energetic process is the formation of the β-sheet, which is equally due to both hydrogen bonds and van der Waals interactions. The aromatic rings help stabilize these aggregates through stacking of the aromatic rings of tyrosine, the stability produced by the aromatics residues increasing with their aromaticity. The formation of the basic unit of the assembled proto-fiber, the steric zipper, is less energetic and is associated to both dispersion forces and hydrogen bonds. The interactions between pair of β-sheets across the peptide-to-peptide contact through the tyrosine rings are cooperative and due to dispersion effects. Moreover, the strength of this interaction can rationalize the variation of mobility of the aromatic ring in the tyrosine units found in solid NMR experiments.

DNA-protein π-interactions in nature: abundance, structure, composition and strength of contacts between aromatic amino acids and DNA nucleobases or deoxyribose sugar

Four hundred twenty-eight high-resolution DNA-protein complexes were chosen for a bioinformatics study. Although 164 crystal structures (38% of those searched) contained no interactions, 574 discrete π-contacts between the aromatic amino acids and the DNA nucleobases or deoxyribose were identified using strict criteria, including visual inspection. The abundance and structure of the interactions were determined by unequivocally classifying the contacts as either π-π stacking, π-π T-shaped or sugar-π contacts. Three hundred forty-four nucleobase-amino acid π-π contacts (60% of all interactions identified) were identified in 175 of the crystal structures searched. Unprecedented in the literature, 230 DNA-protein sugar-π contacts (40% of all interactions identified) were identified in 137 crystal structures, which involve C-H···π and/or lone-pair···π interactions, contain any amino acid and can be classified according to sugar atoms involved. Both π-π and sugar-π interactions display a range of relative monomer orientations and therefore interaction energies (up to -50 (-70) kJ mol(-1) for neutral (charged) interactions as determined using quantum chemical calculations). In general, DNA-protein π-interactions are more prevalent than perhaps currently accepted and the role of such interactions in many biological processes may yet to be uncovered.

On the non-additivity of the substituent effect in ortho-, meta- and para-homo-disubstituted benzenes

The non-additivity of the substituent effect in para-, meta-, and ortho- homo-disubstituted benzenes on π-valence orbitals is smaller than that on σ-ones. The former increases while the latter decreases with π-electron-donating character of the substituent which demonstrates the role of hyperconjugation in the substituent effect.

Examination of tyrosine/adenine stacking interactions in protein complexes

The π-stacking interactions between tyrosine amino acid side chains and adenine-bearing ligands are examined. Crystalline protein structures from the protein data bank (PDB) exhibiting face-to-face tyrosine/adenine arrangements were used to construct 20 unique 4-methylphenol/N9-methyladenine (p-cresol/9MeA) model systems. Full geometry optimization of the 20 crystal structures with the M06-2X density functional theory method identified 11 unique low-energy conformations. CCSD(T) complete basis set (CBS) limit interaction energies were estimated for all of the structures to determine the magnitude of the interaction between the two ring systems. CCSD(T) computations with double-ζ basis sets (e.g., 6-31G*(0.25) and aug-cc-pVDZ) indicate that the MP2 method overbinds by as much as 3.07 kcal mol(-1) for the crystal structures and 3.90 kcal mol(-1) for the optimized structures. In the 20 crystal structures, the estimated CCSD(T) CBS limit interaction energy ranges from -4.00 to -6.83 kcal mol(-1), with an average interaction energy of -5.47 kcal mol(-1), values remarkably similar to the corresponding data for phenylalanine/adenine stacking interactions. Geometry optimization significantly increases the interaction energies of the p-cresol/9MeA model systems. The average estimated CCSD(T) CBS limit interaction energy of the 11 optimized structures is 3.23 kcal mol(-1) larger than that for the 20 crystal structures.

Significant strength of charged DNA-protein π-π interactions: a preliminary study of cytosine

The present work characterized the preferred gas-phase structure and optimum interaction energy of both parallel stacked and perpendicular T-shaped dimers between cytosine (C), as a representative nucleobase, and aspartic/glutamic acid (DE), aspartate/glutamate (DE(-)) or arginine (R(+)), using detailed M06-2X/6-31+G(d,p) potential energy surface scans as a function of the relative monomer orientation. Through comparison to previous literature on the π-π interactions between the DNA nucleobases and the aromatic amino acid residues, this work will allow for comparisons between DNA-protein interactions involving aromatic and acyclic R-side chains, as well as comparisons of the relative geometric dependence and magnitude of π-π (C:DE), πcation-π (C:R(+)), and πanion-π (C:DE(-)) interactions. Our results show that the preferred relative monomer orientation is highly dependent on the monomer composition and charge, and is dictated by electrostatic-driven interactions. More importantly, for the first time, we report that the π-π interactions between cytosine and (neutral) aspartic/glutamic acid are up to approximately -40 kJ mol(-1), while the πcation-π or πanion-π interactions between cytosine and arginine or aspartate/glutamate are up to approximately -90 and -99 kJ mol(-1), respectively. An extensive investigation of the effects of the computational methodology implemented, including comparisons to detailed CCSD(T)/CBS potential energy surfaces and interaction energies, supports the use of M06-2X, as well as ωB97X-D, to study DNA-protein π-π interactions of varying composition and charge. Most importantly, the CCSD(T)/CBS results verify the strong nature of these DNA-protein π-π interactions, as well as the unique nature of the πcation-π and πanion-π counterparts. Therefore, our results emphasize that a wide variety of different types of noncovalent interactions between both cyclic and acyclic π-containing components can significantly contribute to the stability of DNA-protein complexes and likely play a larger role in biology than currently accepted.

Evaluating how discrete water molecules affect protein-DNA π-π and π(+)-π stacking and T-shaped interactions: the case of histidine-adenine dimers

Changes in the magnitude of (M06-2X/6-31+G(d,p)) π-π stacking and T-shaped (nucleobase-edge and amino acid-edge) interactions between (neutral or protonated) histidine (His) and adenine (A) dimers upon microsolvation with up to four discrete water molecules were determined. A variety of histidine-water interactions were considered including conventional (N-H···O, N···H-O, C-H···O) hydrogen bonding and nonconventional (X-H···π (neutral His) or lone-pair···π (protonated His)) contacts. Overall, the effects of discrete His-H(2)O interactions on the neutral histidine-adenine π-π contacts are negligible (<3 kJ mol(-1) or 15%) regardless of the type of water binding, the number of water molecules bound, or the His-A dimer (stacked or (amino acid- or nucleobase-edge) T-shaped) configuration. This suggests that previously reported gas-phase binding strengths for a variety of neutral amino acid-nucleobase dimers are likely relevant for a wide variety of (microsolvated) environments. In contrast, the presence of water decreases the histidine-adenine π(+)-π interaction by up to 15 kJ mol(-1) (or 30%) for all water binding modes and orientations, as well as different stacked and T-shaped His(+)-A dimers. Regardless of the larger effect of discrete histidine-water interactions on the magnitude of the π(+)-π compared with π-π interactions, the π(+)-π binding strengths remain substantially larger than the corresponding π-π contacts. These findings emphasize the distinct nature of π(+)-π and π-π interactions and suggest that π(+)-π contacts can provide significant stabilization in biological systems relative to π-π contacts under many different environmental conditions.

Effects of the biological backbone on stacking interactions at DNA-protein interfaces: the interplay between the backbone···π and π···π components

The (gas-phase) MP2/6-31G*(0.25) π···π stacking interactions between the five natural bases and the aromatic amino acids calculated using (truncated) monomers composed of conjugated rings and/or (extended) monomers containing the biological backbone (either the protein backbone or deoxyribose sugar) were previously compared. Although preliminary energetic results indicated that the protein backbone strengthens, while the deoxyribose sugar either strengthens or weakens, the interaction calculated using truncated models, the reasons for these effects were unknown. The present work explains these observations by dissecting the interaction energy of the extended complexes into individual backbone···π and π···π components. Our calculations reveal that the total interaction energy of the extended complex can be predicted as a sum of the backbone···π and π···π components, which indicates that the biological backbone does not significantly affect the ring system through π-polarization. Instead, we find that the backbone can indirectly affect the magnitude of the π···π contribution by changing the relative ring orientations in extended dimers compared with truncated dimers. Furthermore, the strengths of the individual backbone···π contributions are determined to be significant (up to 18 kJ mol(-1)). Therefore, the origin of the energetic change upon model extension is found to result from a balance between an additional (attractive) backbone···π component and differences in the strength of the π···π interaction. In addition, to understand the effects of the biological backbone on the stacking interactions at DNA-protein interfaces in nature, we analyzed the stacking interactions found in select DNA-protein crystal structures, and verified that an additive approach can be used to examine the strength of these interactions in biological complexes. Interestingly, although the presence of attractive backbone···π contacts is qualitatively confirmed using the quantum theory of atoms in molecules (QTAIM), QTAIM electron density analysis is unable to quantitatively predict the additive relationship of these interactions. Most importantly, this work reveals that both the backbone···π and π···π components must be carefully considered to accurately determine the overall stability of DNA-protein assemblies.

A preliminary investigation of the additivity of pi-pi or pi+-pi stacking and T-shaped interactions between natural or damaged DNA nucleobases and histidine

Previous computational studies have examined pi-pi and pi(+)-pi stacking and T-shaped interactions in nucleobase-amino acid dimers, yet it is important to investigate how additional amino acids affect these interactions since simultaneous contacts often appear in nature. Therefore, this paper investigates the geometries and binding strengths of amino acid-nucleobase-amino acid trimers, which are compared to the corresponding nucleobase-amino acid dimer interactions. We concentrate on systems containing the natural nucleobase adenine or its (cationic) damaged counterpart, 3-methyladenine, and the aromatic amino acid histidine, in both the neutral and protonated forms. This choice of molecules provides information about pi-pi and pi(+)-pi stacking and T-shaped interactions in asymmetric, biologically relevant systems. We determined that both stacked and T-shaped interactions, as well as both pi-pi and pi(+)-pi interactions, exhibit geometric additivity. To investigate the energetic additivity in our trimers, the synergy (E(syn)) and the additivity (E(add)) energy were examined. E(add) reveals that it is important to consider the interaction between the two amino acids when examining the additivity of nucleobase-amino acid interactions. Additionally, E(syn) and E(add) indicate that pi(+)-pi interactions are quite different from pi-pi interactions. The magnitude of E(add) is generally less than 2 kJ mol(-1), which suggests that these interactions are additive. However, the interaction energy analysis does not provide information about the individual interactions in the trimers. Therefore, the quantum theory of atoms in molecules (QTAIM) was implemented. We find inconsistent conclusions from our QTAIM analysis and interaction energy evaluation. However, the magnitudes of the differences between the dimer and trimer critical point properties are extremely small and therefore may not be able to yield conclusive descriptions of differences (if any) between the dimer and trimer interactions. We hypothesize that, due to the limited number of investigations of this type, it is currently unclear how QTAIM can improve our understanding of pi-pi and pi(+)-pi dimers and trimers. Therefore, future work must systematically alter the pi-pi or pi(+)-pi system to definitively determine how the geometry, symmetry, and system size alter the QTAIM analysis, which can then be used to understand biologically relevant complexes.

Effects of the biological backbone on DNA-protein stacking interactions

The pi-pi stacking (face-to-face) interactions between the five natural DNA or RNA nucleobases and the four aromatic amino acids were compared using three different types of dimers: (1) a truncated nucleoside (nucleobase) stacked with a truncated amino acid; (2) a truncated nucleoside (nucleobase) stacked with an extended amino acid; and (3) a nucleoside (extended nucleobase) stacked with a truncated amino acid. Systematic (MP2/6-31G*(0.25)) potential energy surface scans reveal important information about the effects of the deoxyribose sugar and protein backbone on the structure and binding energy between truncated nucleobase and amino acid models that are typically implemented in the literature. Most notably, electrostatic and steric interactions arising from the bulkiness of the biological backbones can change the preferred relative orientations of DNA and protein pi-systems. More importantly, the protein backbone can strengthen the stacking energy (by up to 10 kJ mol(-1)), while the deoxyribose moiety can strengthen or weaken the stacking interaction depending on the positioning of the amino acid relative to the sugar residue. These effects are likely due to additional interactions between the amino acid or nucleobase ring and the backbone in the extended monomer rather than significant changes in the properties of the biological pi-systems upon model extension. Since the present work reveals that all calculated DNA-protein stacking interactions are significant and approach the strength of other noncovalent interactions between biomolecules, both pi-pi and backbone-pi interactions must be considered when attempting to gain a complete picture of DNA-protein binding.