Halogen bonding, a specific intermolecular noncovalent interaction, plays crucial roles in fields as diverse as molecular recognition, crystal engineering, and biological systems. This paper presents an ab initio investigation of a series of dimeric complexes formed between bromobenzene and several electron donors. Such small model systems are selected to mimic halogen bonding interactions found within crystal structures as well as within biological molecules. In all cases, the intermolecular distances are shown to be equal to or below sums of van der Waals radii of the atoms involved. Halogen bonding energies, calculated at the MP2/aug-cc-pVDZ level, span over a wide range, from -1.52 to -15.53 kcal/mol. The interactions become comparable to, or even prevail over, classical hydrogen bonding. For charge-assisted halogen bonds, calculations have shown that the strength decreases in the order OH- > F- > HCO2- > Cl- > Br-, while for neutral systems, their relative strengths attenuate in the order H2CS > H2CO > NH3 > H2S > H2O. These results agree with those of the quantum theory of atoms in molecules (QTAIM) since bond critical points (BCPs) are identified for these halogen bonds. The QTAIM analysis also suggests that strong halogen bonds are more covalent in nature, while weak ones are mostly electrostatic interactions. The electron densities at the BCPs are recommended as a good measure of the halogen bond strength. Finally, natural bond orbital (NBO) analysis has been applied to gain more insights into the origin of halogen bonding interactions.

Ab initio calculations have been performed on a series of complexes formed between halogen-containing molecules and ammonia to gain a deeper insight into the nature of halogen bonding. It appears that the dihalogen molecules form the strongest halogen-bonded complexes with ammonia, followed by HOX; the charge-transfer-type contribution has been demonstrated to dominate the halogen bonding in these complexes. For the complexes involving carbon-bound halogen molecules, our calculations clearly indicate that electrostatic interactions are mainly responsible for their binding energies. Whereas the halogen-bond strength is significantly enhanced by progressive fluorine substitution, the substitution of a hydrogen atom by a methyl group in the CH(3)X...NH(3) complex weakened the halogen bonding. Moreover, remote substituent effects have also been noted in the complexes of halobenzenes with different para substituents. The influence of the hybridization state of the carbon atom bonded to the halogen atom has also been examined and the results reveal that halogen-bond strengths decrease in the order HC triple bond CX > H(2)C=CHX approximately O=CHX approximately C(6)H(5)X > CH(3)X. In addition, several excellent linear correlations have been established between the interaction energies and both the amount of charge transfer and the electrostatic potentials corresponding to an electron density of 0.002 au along the R-X axis; these correlations provide good models with which to evaluate the electron-accepting abilities of the covalently bonded halogen atoms. Finally, some positively charged halogen-bonded systems have been investigated and the effect of the charge has been discussed.

A systematic theoretical investigation on a series of dimeric complexes formed between some halocarbon molecules and electron donors has been carried out by employing both ab initio and density functional methods. Full geometry optimizations are performed at the Moller-Plesset second-order perturbation (MP2) level of theory with the Dunning's correlation-consistent basis set, aug-cc-pVDZ. Binding energies are extrapolated to the complete basis set (CBS) limit by means of two most commonly used extrapolation methods and the aug-cc-pVXZ (X = D, T, Q) basis sets series. The coupled cluster with single, double, and noniterative triple excitations [CCSD(T)] correction term, determined as a difference between CCSD(T) and MP2 binding energies, is estimated with the aug-cc-pVDZ basis set. In general, the inclusion of higher-order electron correlation effects leads to a repulsive correction with respect to those predicted at the MP2 level. The calculations described herein have shown that the CCSD(T) CBS limits yield binding energies with a range of -0.89 to -4.38 kcal/mol for the halogen-bonded complexes under study. The performance of several density functional theory (DFT) methods has been evaluated comparing the results with those obtained from MP2 and CCSD(T). It is shown that PBEKCIS, B97-1, and MPWLYP functionals provide accuracies close to the computationally very expensive ab initio methods.

Glycogen synthase kinase-3 (GSK-3) is a kind of serine-threonine protein kinase. It places important roles in several signaling pathways and it is a key therapeutic target for a number of diseases, such as diabetes, cancer, Alzheimer's disease and chronic inflammation. Mg(2+) ions which interact with ATP are conserved in GSK. They are important in phosphoryl transfer. Li(+) is an inhibitor for GSK-3. It is used to treat bipolar mood disorder. This paper illustrates the effect of Li(+) on GSK-3. When Mg(I)(2+) is replaced by Li(+), the atom fluctuation of GSK-3 will rise, and the in-line phosphoryl transfer mechanism is probably demolished and the binding of pre-phosphorylated substrates may be disturbed. All the results we obtained clearly suggest that inhibition to GSK-3 is caused by the Mg(I)(2+) replacement with Li(+).

Triangular halogen trimers (RX)3, where X = Br, I and R represents H, H3C, H2FC, HF2C, F3C, CH2=CH, CH[triple bond]C, and Ph, have been investigated using the density functional theory in the Perdew, Burke, and Ernzerhof method. We report herein the optimized geometries of the stable structures, their vibrational frequencies, and binding energies with the two- and three-body terms. All trimer structures possess a cyclic array of halogen atoms in the type II approach by virtue of the nonspherical atomic charge distribution around the halogens. The Br...Br interactions in trimers are very weak, whereas the I...I interactions in trimers are relatively stronger. While all bromine trimers and most of iodine trimers are predicted to be noncooperative, three of iodine trimers show weak cooperativity. The analysis of vibration modes reveals that all halogen trimers exhibit no especially remarkable frequency shifts. It is also shown that the electrostatic contribution plays a major role in the halogen...halogen interactions in halogen trimers. In contrast to bromine trimers, the relative contribution of charge-transfer component to the halogen...halogen interactions becomes more important for iodine trimers.

In this study, structural topology scale (ST-scale) was recruited as a novel structural topological descriptor derived from principal component analysis on 827 structural variables of 167 amino acids. By using partial least squares (PLS), we applied ST-scale for the study of quantitative sequence-activity models (QSAMs) on three peptide datasets (58 angiotensin-converting enzyme (ACE) inhibitors, 34 antimicrobial peptides (AMPs) and 89 elastase substrates (ES)). The results of QSAMs were superior to that of the earlier studies, with determination coefficient (r(2)) and cross-validated (q(2)) equal to 0.855, 0.774; 0.79, 0.371 (OSC-PLS: 0.995, 0.848) and 0.846, 0.747, respectively. Therefore, ST-scale descriptors were considered to be competent to extract information from 827 structural variables and relate with their bioactivities.

No abstract

The phenomenon that SARS coronavirus main protease (SARS M(pro)) dimer is the main functional form has been confirmed by experiment. However, because of the absence of structural information of the monomer, the reasons for this remain unknown. To investigate it, two molecular dynamics (MD) simulations in water for dimer and monomer models have been carried out, using the crystal structure of protomer A of the dimer as the starting structure for the monomer. During the MD simulation of dimer, three interest phenomena of protomer A have been observed: (i) the distance between NE2 of His41 and SG of Cys145 averages 3.72 A, which agrees well with the experimental observations made by X-ray crystallography; (ii) His163 and Glu166 form the "tooth" conformational properties, resulting in the specificity for glutamine at substrate P1 site; and (iii) the substrate-binding pocket formed by loop 140-146 and loop 184-197 is large enough to accommodate the substrate analog. However, during the MD simulation of the monomer complex, the three structural characteristics are all absent, which results directly in the inactivation of the monomer. Throughout the MD simulation of the dimer, the N-terminus of protomer B forms stable hydrogen bonds with Phe140 and Glu166, through which His163, Glu166, and loop 140-146 are kept active form. Furthermore, a water-bridge has been found between the N-terminus of protomer B and Gly170, which stabilizes His172 and avoids it moving toward Tyr161 to disrupt the H-bond between Tyr161 and His163, stabilizing the conformation of His163. The interactions between the N-terminus and another monomer maintain the activity of dimer.