Molecular Dynamics Simulations of a β-Hairpin Fragment of Protein G by Means of Atom-Bond Electronegativity Equalization Method Fused into Molecular Mechanics (ABEEMσπ/MM)

There are some controversial opinions about the origin of folding β‐hairpin stability in aqueous solution. In this study, the structural and dynamic behavior of a 16‐residue β‐hairpin from B1 domain of protein G has been investigated at 280, 300, 350 and 450 K using molecular dynamics (MD) simulations by means of Atom‐Bond Electronegativity Equalization Method Fused into Molecular Mechanics i.e., ABEEMδπ/MM and the explicit ABEEM‐7P water solvent model. In addition, a 300 K simulation of one mutant having the aromatic residues substituted with alanines has been performed. The hydrophobic surface area, hydrophilic surface area and some structural properties have been used to measure the role of the hydrophobic interactions. It is found that the aromatic residues substituted with alanines have shown an evident destabilization of the structure and unfolding started after 1.5 ns. It is also found that the number of the main chain hydrogen bonds have different distributions through three different simulations. All above demonstrate that the hydrophobic interactions and the main chain hydrogen bonds play an important role in the stability of the folding structure of β‐hairpin in solution. Furthermore, through the structural analyses of the β‐hairpin structures from four temperature simulations and the comparison with other MD simulations of β‐hairpin peptides, the new ABEEMδπ force field can reproduce the structural data in good agreement with the experimental data.     

Performance assessment of density-functional methods for study of charge-transfer complexes

Various density functionals are applied to a number of weakly bound intermolecular pi-pi charge-transfer (CT) complexes. Most functionals, including the recently developed mPWPW91 and mPW1PW91, grossly underestimate experimental excitation energies; good agreement is obtained only with the half-and-half hybrid BH&HLYP functional. PW91PW91 provides the best agreement with intermolecular distances measured in crystal, while the BH&HLYP values are about 0.1 A too long. Various hybrid functionals with nonlocal exchange correction provide binding energies that compare favorably with the experimental heats of formation measured in solution.     

具有分子内折叠结构锌配合物的合成，与DNA键合及其对DNA切割作用研究

本文报道了一个锌配合物[Zn(bipy)(pmal)(H2O)]·2H2O（其中bipy = 2,2’-联吡啶, pmal = 苯基丙二酸）的合成,晶体结构及其光谱学研究。并通过单晶X射线衍射,元素分析,红外光谱等手段对它进行了表征。同时,利用紫外光谱和荧光光谱方法考察了该配合物与小牛胸腺DNA的键和作用。琼脂糖凝胶电泳实验的结果说明该配合物的平面结构对pBR 322DNA切割作用显著。     

Synthesis, Characterization and Weak Intramolecular Interactions of Porphyrins Bearing Nucleobases

5,10, 15-Triphenyl-20-{2- [α- (adenine-9 ) acetylamino]} phenyl porphyrin ( 1 ), 5,10, 15-triphenyl-20-{2-[α-(cytosine-1)acetylamino]} phenyl porphyrin (2), 5, 10, 15-triphenyl-20-{4-[α-(cytosine-1)ethoxy]} phenyl porphyrin (3) and their zinc complexes Zn-1, Zn-2 and Zn-3 have been prepared and characterized by ^1H NMR spectra, elemental analyses, electronic absorption spectra and mass spectra (FAB). Intramolecular π-π interactions and intramolecular metal-~ interaction for 1, 2, Zn-1,and Zn-2 have been investigated by several methods. ^1H NMR studies demonstrate that the porphyrin π-system in 1 and 2 is parallel to the adenine and the cytosine aromatic ring, respectively. The electronic absorption spectral properties of free porphyrin derivatives and their zinc complexes have been compared with those of H2TPP and ZnTPP. The results show that the UV-vis spectra of 1 and 2 are the same as that of H2TPP,whereas the spectra of their zinc complexes show 7 nm red shifts of the Soret bands compared to that of ZnTPP. The emission spectra of Zn-1 and Zn-2 are independent of excitation wavelength. From combination of the evidence of absorption and emission spectra it is suggested the existence of intramolecular metal-π interaction in Zn-1 and Zn-2. The results of conformational analysis agreed quite nicely with that of experiments, thus it was further to validate the experimental conclusions.     

Effect of octabromination of a tetrakis(4-carboxyphenyl)porphine derivative bound to silica gels on HPLC retention behaviors of polyaromatic hydrocarbons

The effect of bromination of Cu-porphyrin-derivative-immobilized silica gels (Cu-TCPP_D) was examined by comparing the retention behaviors of polyaromatic hydrocarbons (PAHs) on Cu-TCPP_D and Cu-octabromotetrakis(4-carboxyphenyl)porphine-derivative-immobilized silica gels (Cu-OBTCPP_D) columns. It was revealed that bromination affects strongly the π–π electron interactions caused from hydration energy in a polar eluent (80% methanol) possibly as a result of a destruction of planar structure of porphine-ring by bromination. It was also revealed that bromination enhances π–d interactions as well as the π–π electron interactions in a broad sense (e.g., dispersion forces) in a non-polar eluent (n-hexane). However, the bromination did not exert much influence on electrostatic interactions caused from polarization of mono-halogenated benzenes.     

Violet-to-blue tunable emission of aryl-substituted dispirofluorene-indenofluorene isomers by conformationally-controllable intramolecular excimer formation

Two series of DiSpiroFluorene-IndenoFluorene (DSF-IF) positional isomers, namely dispiro[2,7-diarylfluorene-9',6,9',12-indeno[1,2-b]fluorenes], (1,2-b)-DSF-IFs 1 and dispiro[2,7-diarylfluorene-9',6,9',12-indeno[2,1-a]fluorenes], (2,1-a)-DSF-IFs 2 have been synthesized. These violet-to-blue fluorescent emitters possess a 3π-2spiro architecture, which combines via two spiro links two different indenofluorene cores, that is, (1,2-b)-IF or (2,1-a)-IF and 2,7-substituted-diaryl-fluorene units. Due to their different geometric profiles, the two families of positional isomers present drastically different properties. The marked difference observed between the properties of (1,2-b)-DSF-IF (1) and (2,1-a)-DSF-IF (2) is discussed in terms of intramolecular π-π interactions occurring in (2,1-a)-DSF-IF (2) leading to conformationally-controllable intramolecular excimer formation. Indeed, the original geometry of the (2,1-a)-DSF-IF (2) family, with face-to-face "aryl-fluorene-aryl" moieties, leads to remarkable excimer emission through intramolecular π-π interactions in the excited state. Furthermore, the emission wavelengths can be gradually modulated by the control of the steric hindrance between the adjacent substituted phenyl rings. Thus, through a comparative and detailed study of the (1)H NMR, electrochemical and photophysical properties of DSF-IFs 1 and 2, we have evidenced the intramolecular π-π interactions occurring between the two "aryl-fluorene-aryl" moieties in the ground state and in the excited state. These properties have been finally correlated to the spectacular conformational change modeled by density functional theory (DFT) calculations. Indeed, the two "aryl-fluorene-aryl" moieties switch from a staggered conformation in the ground state to an eclipsed conformation in the first excited state.     

Competition between π–π stacking and hydrogen bonding in (1:2) picrates of 1,5-diamino-3-oxapentane, 1,8-diamino-3,6-dioxaoctane and 1,5-diamino-3-azapentane—solid state studies

Crystallographic studies of (2:1) salts of picric acid with 1,5-diamino-3-oxapentane (1OPICR), 1,8-diamino-3,6-dioxaoctane (2OPICR) and 1,5-diamino-3-azapentane (1NPICR) showed significant conformational change of the picrate ion due to numerous electrostatic, H-bonding and π–π stacking interactions present in the crystal lattice. In particular, intermolecular N–HO H-bonds were found to cause significant twisting of the o-NO₂ groups from the plane of the benzene ring, whereas overlapping of the picrate ions due to electrostatic interactions and π–π stacking caused flattening of the molecule. Analysis of the geometry of 74 picrate ions found in the Cambridge Crystallographic Database, in their various crystallochemical environments, showed that competition between essentially weak but numerous intermolecular interactions of different types led to systematic changes in geometric parameters within the picrate ion. In particular, relations found between the C1–C2–N–O (C1–C6–N–O) torsion angle and the endocyclic C1–C2–C3 (C1–C6–C5) valence angle can be explained on the basis of competition between resonance effects of the o-NO₂ group and π–π stacking.     

蒽与苯甲酸及其衍生物之间的作用机理研究

用紫外差光谱、红外光谱法和荧光猝灭法研究了典型多环芳烃蒽与苯甲酸及其羟基取代衍生物邻羟基苯甲酸和对羟基苯甲酸之间的作用机理。实验结果表明,在此芳香羧酸与蒽之间存在定向的、特异性的作用,其作用方式受到反应物结构和环境酸度的影响。对苯甲酸和邻羟基苯甲酸而言,当pH<pKa时,二者之间以蒽离域大π电子与羧基质子之间的π-H氢键作用为主,pH>pKa时,π-π电子给体受体作用逐渐成为主要结合方式。邻位羟基的存在使得苯甲酸与蒽的作用强度明显降低。对羟基苯甲酸特殊的D-π-A型分子结构使得它在溶液中形成平面多分子聚集体,这种多分子聚集体的生成使得对羟基苯甲酸与蒽的结合方式不随酸度变化,在pH2.0~10.0的范围内均以π-π电子给体受体作用相结合,且结合强度大于苯甲酸和邻羟基苯甲酸。     

The Bright Future of Unconventional σ/π‐Hole Interactions

Non‐covalent interactions play a crucial role in (supramolecular) chemistry and much of biology. Supramolecular forces can indeed determine the structure and function of a host–guest system. Many sensors, for example, rely on reversible bonding with the analyte. Natural machineries also often have a significant non‐covalent component (e.g. protein folding, recognition) and rational interference in such ‘living’ devices can have pharmacological implications. For the rational design/tweaking of supramolecular systems it is helpful to know what supramolecular synthons are available and to understand the forces that make these synthons stick to one another. In this review we focus on σ‐hole and π‐hole interactions. A σ‐ or π‐hole can be seen as positive electrostatic potential on unpopulated σ* or π⁽*⁾ orbitals, which are thus capable of interacting with some electron dense region. A σ‐hole is typically located along the vector of a covalent bond such as X?H or X?Hlg (X=any atom, Hlg=halogen), which are respectively known as hydrogen and halogen bond donors. Only recently it has become clear that σ‐holes can also be found along a covalent bond with chalcogen (X?Ch), pnictogen (X?Pn) and tetrel (X?Tr) atoms. Interactions with these synthons are named chalcogen, pnigtogen and tetrel interactions. A π‐hole is typically located perpendicular to the molecular framework of diatomic π‐systems such as carbonyls, or conjugated π‐systems such as hexafluorobenzene. Anion–π and lone‐pair–π interactions are examples of named π‐hole interactions between conjugated π‐systems and anions or lone‐pair electrons respectively. While the above nomenclature indicates the distinct chemical identity of the supramolecular synthon acting as Lewis acid, it is worth stressing that the underlying physics is very similar. This implies that interactions that are now not so well‐established might turn out to be equally useful as conventional hydrogen and halogen bonds. In summary, we describe the physical nature of σ‐ and π‐hole interactions, present a selection of inquiries that utilise σ‐ and π‐holes, and give an overview of analyses of structural databases (CSD/PDB) that demonstrate how prevalent these interactions already are in solid‐state structures.     

Asymmetric Epoxidation of Terminal Olefins with Binaphthyl Strapped Porphyrin Catalysts: (‐( Stacking Interaction and Steric Effects on the Enantioselectivities

将两种具有相似手性联萘空腔的手性卟啉1b和2b作为催化剂,对苯乙烯衍生物及非芳香烯烃底物进行了催化不对称环氧化反应。应用理论计算方法研究了底物同两种催化剂之间的作用机理;应用1H NMR探讨了催化剂的手性空腔结构;结果表明芳香烯烃同催化剂的手性联萘基团之间的p-p相互作用是决定该类催化剂对映选择性的重要因素,同时催化剂手性结构的立体位阻效应使该类催化剂具有良好的ee值。     

Typical aromatic noncovalent interactions in proteins: A theoretical study using phenylalanine

A systematic study of CH ··· π, OH ··· π, NH ··· π, and cation ··· π interactions has been done using complexes of phenylalanine in its cationic, anionic, neutral, and zwitterionic forms with CH₄, H₂O, NH₃, and NH at B3LYP, MP2, MPWB1K, and M06‐2X levels of theory. All noncovalent interactions are identified by the presence of bond critical points (bcps) of electron density (ρ( r )) and the values of ρ( r ) showed linear relationship to the binding energies (E_total). The estimated E_total from supermolecule, fragmentation, and ρ( r ) approaches suggest that cation ··· π interactions are in the range of 36 to 46 kcal/mol, whereas OH ··· π, and NH ··· π interactions have comparable strengths of 6 to 27 kcal/mol and CH ··· π interactions are the weakest (0.62–2.55 kcal/mol). Among different forms of phenylalanine, cationic form generally showed the highest noncovalent interactions at all levels of theory. Cooperativity of multiple interactions is analyzed on the basis of ρ( r ) at bcps which suggests that OH ··· π and NH ··· π interactions show positive, whereas CH ··· π and cation ··· π interactions exhibit negative cooperativity with respect to the side chain hydrogen bond interactions. In general, side chain interactions are strengthened as a result of aromatic interaction. Solvation has no significant effect on the overall geometry of the complex though slight weakening of noncovalent interactions by 1–2 kcal/mol is observed. An assessment of the four levels of theory studied herein suggests that both MPWB1K and M06‐2X give better performance for noncovalent interactions. The results also support the fact that B3LYP is inadequate for the study of weak interactions. © 2008 Wiley Periodicals, Inc. J Comput Chem 2009     

两种吩噻嗪衍生物的合成，光学性质及其电化学性质

The synthesis, luminescence and electrochemical properties of two novel phenothiazine derivatives were presented, which exhibit good one-photon fluorescence emission. The quantum yields and solvent effects of them were studied in detail. Compound 3,7-bis（2-4＇-imidazolylbenzylidenehydrazonoethyl）-（0-ethylphenothiazine） （5） was used as an effective initiator, of which two-photon fluorescence spectra were investigated under 800 nm fs laser pulse and the measured two-photon absorption cross-section was 18 × 10^-50 cm^4·s per photon. Two-photon initiated polymerization microfabrication experiments were carried out. X-ray diffraction analyses revealed that the shape of the compound 3,7-bis（2-pyrid-4＇-ylmethylidenehydrazonoethyl）-（0-ethylphenothiazine） （4） looks like a butterfly with nearly planar wings. The dihedral angle of the two benzene rings is 37.6° and there is an obvious π-π stacking interaction between the molecules in the crystal.     

Study on the Cation-π Interactions between Ammonium Ion and Aromatic π Systems

The nature and strength of the cation-π interactions between NH4^＋ and toluene, p-cresol, or Me-indole were studied in terms of the topological properties of molecular charge density and binding energy decomposition. The results display that the diversity in the distribution pattern of bond and cage critical points reflects the profound influence of the number and nature of substituent on the electron density of the aromatic rings. On the other hand, the energy decomposition shows that dispersion and repulsive exchange forces play an important role in the organic cation （NH4^＋）-π interaction, although the electrostatic and induction forces dominate the interaction. In addition, it is intriguing that there is an excellent correlation between the electrostatic energy and ellipticity at the bond critical point of the aromatic π systems, which would be helpful to further understand the electrostatic interaction in the cation-π complexes.     

A generalized Langevin dynamics approach to model solvent dynamics effects on proteins via a solvent-accessible surface. The carboxypeptidase A inhibitor protein as a model

A generalized Langevin dynamics (GLD) scheme is derived for (bio)macromolecules having internal structure, arbitrary shapes and a size larger than solvent molecules (i.e. proteins). The concept of solvent-accessible surface area (SASA) is used to incorporate solvent effects via external forces thereby avoiding its explicit molecular representation. A simulation algorithm is implemented in the GROMOS molecular dynamics (MD) program including random forces and memory effects, while solvation effects enter via derivatives of the surface area. The potato carboxypeptidase inhibitor (PCI), a small protein, is used to numerically test the approach. This molecule has N- and C-terminal tails whose structure and fluctuations are solvent dependent. A 1-ns MD trajectory was analyzed in depth. X-ray and NMR structures are used in conjunction with MD simulations with and without explicit solvent to gauge the quality of the results. All the analyses showed that the GLD simulation approached the results obtained for the MD simulation with explicit simple-point-charge-model water molecules. The SASAs of the polar atoms show a natural exposure towards the solvent direction. A FLS solvent simulation was completed in order to sense memory effects. The approach and results presented here could be of great value for developing alternatives to the use of explicit solvent molecules in the MD simulation of proteins, expanding its use and the time-scale explored. Received: 2 February 2000 / Revised: 12 March 2000 / Accepted: 26 May 2000 / Published online: 2 November 2000     

Synthesis, crystal structure and luminescence properties of europium complexes with a new terpyridine-like ligand

New 1:1, 1:2 and 1:3 europium (Ⅲ) complexes with 2,4,6-tris(3,5-dimethylpyrazol-1-yl)-1, 3,5-triazine have been prepared, and their luminescence has been demonstrated to be sensitive to the π-π stacking and H-bonding interactions.     

Molecular docking approach on the molecular interactions involving beta-lactoglobulin (βLG)-4-Dicyanomethylene2,6-Dimethyl-4-Hpyran (DDP) dye in the presence of an antibiotic,norfloxacin

Molecular docking (MD) techniques were employed in the determination of binding stability of Beta-lactoglobulin (βLG) with 4-Dicyanomethylene2,6-Dimethyl-4-Hpyran (DDP) dye in the presence of quinolone-based antibiotic, norfloxacin (NOR). βLG acts as the host and drug/dye is employed as the guest molecule. The energetics and molecular interaction parameters of βLG-DDP were found to be stable than βLG-NOR. DDP resides in three distinguishably different binding sites of βLG, however the binding energies does not differ considerably. Interestingly, the energetics of βLG-NOR conformers differs significantly and reveals the existence of two distinguishable conformers. NOR acts as a hydrogen-bonding (HB) acceptor and the amino acid (AA) residues of βLG as the donor. The extent of HB interactions predominates over hydrophobic interactions in βLG-NOR than βLG-DDP complex. MD studies authenticates no direct binding of dye-drug inside the domains of βLG, when docked simultaneously. On the contrary, docking of dye to βLG-NOR complex enhances the binding stability. Several molecular interactions govern the stability of the complex and play a major role on the binding affinity of the host-guest complex in the presence of two competitive ligands. MD techniques authenticate that DDP/NOR is bound to several AA residues through conventional and non-conventional HB interactions accompanied with hydrophobic, pi-pi, pi-alkyl and van der Waals forces of interactions. Further, MD methods is employed as an efficient tool and a non-evasive technique in establishing the stability of 1:1complex of drug/dye with βLG in the present study.     

Interplay between anion-pi and hydrogen bonding interactions

The interplay between two important noncovalent interactions involving aromatic rings is studied by means of high level ab initio calculations. They demonstrate that synergistic effects are present in complexes where anion-pi and hydrogen bonding interactions coexist. These synergistic effects have been studied using the "atoms-in-molecules" theory and the Molecular Interaction Potential with polarization partition scheme. The present study examines how these two interactions mutually influence each other.     

Complete basis set limit of Ab initio binding energies and geometrical parameters for various typical types of complexes

Using basis‐set extrapolation schemes for a given data set, we evaluated the binding energies and geometries at the complete basis set (CBS) limit at the levels of the second order Møller–Plesset perturbation theory (MP2) and the coupled cluster theory with singles, doubles, and perturbative triples excitations [CCSD(T)]. The systems include the hydrogen bonding (water dimer), aromatic interaction (benzene dimer), π–H interaction (benzene–water), cation–water, anion–water, π–cation interaction (cation–benzene), and π–anion interaction (anion–triazine). One extrapolation method is to exploit both BSSE‐corrected and BSSE‐uncorrected binding energies for the aug‐cc‐pVNZ (N = 2, 3, 4, …) basis set in consideration that both binding energies give the same CBS limit (CBS^B). Another CBS limit (CBS^C) is to use the commonly known extrapolation approach to exploit that the electron correlation energy is proportional to N^?3. Since both methods are complementary, they are useful for estimating the errors and trend of the asymptotic values. There is no significant difference between both methods. Overall, the values of CBS^C are found to be robust because of their consistency. However, for small N (in particular, for N = 2, 3), CBS is found to be slightly better for water–water interactions and cation–water and cation–benzene interactions, whereas CBS is found to be more reliable for bezene–water and anion–water interactions. We also note that the MP2 CBS limit value based on N = 2 and 3 combined with the difference between CCSD(T) and MP2 at N = 2 would be exploited to obtain a CCSD(T)/CBS value for aromatic–aromatic interactions and anion–π interactions, but not for cationic complexes. © 2007 Wiley Periodicals, Inc. J Comput Chem, 2008     

Enhanced sampling method in molecular simulations using genetic algorithm for biomolecular systems

We propose a molecular simulation method using genetic algorithm (GA) for biomolecular systems to obtain ensemble averages efficiently. In this method, we incorporate the genetic crossover, which is one of the operations of GA, to any simulation method such as conventional molecular dynamics (MD), Monte Carlo, and other simulation methods. The genetic crossover proposes candidate conformations by exchanging parts of conformations of a target molecule between a pair of conformations during the simulation. If the candidate conformations are accepted, the simulation resumes from the accepted ones. While conventional simulations are based on local update of conformations, the genetic crossover introduces global update of conformations. As an example of the present approach, we incorporated genetic crossover to MD simulations. We tested the validity of the method by calculating ensemble averages and the sampling efficiency by using two kinds of peptides, ALA3 and (AAQAA)₃. The results show that for ALA3 system, the distribution probabilities of backbone dihedral angles are in good agreement with those of the conventional MD and replica-exchange MD simulations. In the case of (AAQAA)₃ system, our method showed lower structural correlation of α-helix structures than the other two methods and more flexibility in the backbone ψ angles than the conventional MD simulation. These results suggest that our method gives more efficient conformational sampling than conventional simulation methods based on local update of conformations. © 2018 Wiley Periodicals, Inc.