Density functional theory study on (LiNH<Subscript>2</Subscript>)<Subscript>n</Subscript> (<Emphasis Type="Italic">n</Emphasis> = 1–5) clusters

Geometrical structures and relative stabilities of (LiNH₂)_n (n = 1–5) clusters were studied using density functional theory (DFT) at the B3LYP/6-31G* and B3LYP/6-31++G* levels. The electronic structures, vibrational properties, N–H bond dissociation energies (BDE), thermodynamic properties, bond properties and ionization potentials were analyzed for the most stable isomers. The calculated results show that the Li–N and Li–Li bonds can be formed more easily than those of the Li–H or N–H bonds in the clusters, in which NH₂ is bound to the framework of Li atomic clusters with fused rings. The average binding energies for each LiNH₂ unit increase gradually from 142 kJ mol⁻¹ up to about 180 kJ mol⁻¹ with increasing n. Natural bond orbital (NBO) analysis suggests that the bonds between Li and NH₂ are of strong ionicity. Three-center–two-electron Li–N–Li bonding exists in the (LiNH₂)₂ dimer. The N–H BDE values indicate that the change in N–H BDE values from the monomer a1 to the singlet-state clusters is small. The N–H bonds in singlet state clusters are stable, while the N–H bonds in triplet clusters dissociate easily. A study of their thermodynamic properties suggests that monomer a1 forms clusters (b1, c1, d2 and e1) easily at low temperature, and clusters with fewer numbers of rings tend to transfer to ones with more rings at low temperature. E _g, E _HOMO and E _av decrease gradually, and become constant. Ring-like (LiNH₂)_3,4 clusters possess higher ionization energy (VIE) and E _g, but lower values of E _HOMO. Ring-like (LiNH₂)_3,4 clusters are more stable than other types. A comparison of structures and spectra between clusters and crystal showed that the NH₂ moiety in clusters has a structure and spectral features similar to those of the crystal.     

A DFT study on the sensing behavior of a BC2N nanotube toward formaldehyde

We investigated the viability of using a BC₂N nanotube to detect formaldehyde (H₂CO) molecule by means of B3LYP and M06 density functionals. The results indicate that the molecule is weakly adsorbed on the intrinsic BC₂N nanotube releasing energy of 0.8 kcal mol^-1 (at B3LYP/6-31G(d)) without significant effect on the HOMO-LUMO energy gap and electrical conductivity of the tube. Thus, H₂CO cannot be detected using this intrinsic nanotube. To overcome this problem, a carbon atom of the tube wall was substituted by a Si atom. It was demonstrated that the Si-doped tube cannot only strongly adsorb the H₂CO molecule, but also may effectively detect its presence because of the increase in the electric conductivity of the tube.     

Single-Chain Fv of Anti-idiotype 11-1G10 Antibody Interacts with Antibody NC41 Single-Chain Fv with a Higher Affinity than the Affinity for the Interaction of the Parent Fab Fragments

A single-chain Fv (scFv) fragment of anti-idiotype antibody 11-1G10, which recognizes an idiotope of anti-neuraminidase antibody NC41, was constructed by joining V_H and V_L domains with a (Gly₄Ser)₃ linker, with a pelB leader sequence, and two C-terminal FLAG tag sequences, and expressed in E. coli (10 mg/L). The 11-1G10 scFv was isolated by affinity chromatography on an anti-FLAG M2 antibody column as a 2:1 mixture of monomer and dimer forms which were separated by Superdex 75 chromatography; monomer (at 100 g/ml) was stable for 7 days at 21°C and 30 days at 4°C, whereas the dimer slowly dissociated to monomer to yield a 2:1 monomer–dimer equilibrium mixture after 30 days at 4°C. The dimer was bivalent, with each combining site binding an NC41 Fab to yield a stable complex of M _r 156,000. Binding affinities, determined in solution using a BIAcore biosensor, showed that the affinity for the interaction of 11-1G10 scFv monomer with NC41 scFv monomer was five- to six-fold higher than the interaction of the parent Fab pair. This is the first example of an scFv derived from a monoclonal antibody with a higher affinity than its parent Fab.     

Theoretical potential scans and vibrational spectra of vinyl selenonyl halides CH2=CH–SeO2X (X is F, Cl and Br)

The structure and conformational stability of vinyl selenonyl fluoride, chloride and bromide CH₂=CH–SeO₂X (X is F, Cl and Br) were investigated using density functional B3LYP/6-311+G** and ab initio MP2/6-311+G** calculations. From the calculations the molecules were predicted to exist only in the non-planar gauche conformation with the vinyl C=C group almost eclipsing one of the selenonyl Se=O bonds as a result of conjugation between the two moieties. Single-minimum potential scans were calculated at the DFT level for the molecules. The vibrational frequencies were computed using B3LYP/6-311+G**. Normal coordinate calculations were then carried out and potential energy distributions were calculated for the three molecules in the gauche conformation.Figure Potential function for the asymmetric torsion in vinyl selenonyl fluoride (dotted line), chloride (dashed line) and bromide (solid line) as determined at the DFT-B3LYP/6-311+G** level     

Ab initio study of phenyl benzoate: structure,conformational analysis,dipole moment,IR and Raman vibrational spectra

The molecular structure (bond distances and angles), conformational properties, dipole moment and vibrational spectroscopic data (vibrational frequencies, IR and Raman intensities) of phenyl benzoate were calculated using Hartree–Fock (HF), density functional (DFT), and second order Møller–Plesset perturbation theory (MP2) with basis sets ranging from 6-31G* to 6-311++G**. The theoretical results are discussed mainly in terms of comparisons with available experimental data. For geometric data, good agreement between theory and experiment is obtained for the MP2, B3LYP and B3PW91 levels with basis sets including diffuse functions. The B3LYP/6-31+G* theory level estimates the shape of the experimental functions for phenyl torsion around the Ph–O and Ph–C bonds well, but reproduces the height of the rotational barriers poorly. The B3LYP/6-31+G* harmonic force constants were scaled by applying the scaled quantum mechanical force field (SQM) technique. The calculated vibrational spectra were interpreted and band assignments were reported. They are in excellent agreement with experimental IR and Raman spectra.Figure Calculated and experimental (GED) potential energy functions for torsional motion of phenyl benzoate relative to the minimum value. a The potential function for torsion about the O₃–C₄ bond. b The potential function for torsion about the C₂–C₁₀ bond.     

Normal coordinate analyses of 3,5-dichlorophenylcyanamide

The structure of 3,5-dichlorophenylcyanamide c-C₆H₃Cl₂–NHCN was investigated by DFT-B3LYP and ab initio MP2 calculations with the 6-311+G** basis set. The planar to perpendicular rotational barrier was calculated to be of about 5 kcal mol^–1 at both levels of calculation. The stability of the planar structure of the molecule was explained on the basis of conjugation effects between the cyanamide–NHCN moiety and the phenyl c-C₆H₅ ring in agreement with earlier NMR results. The CNC and the HNC bond angles were calculated to be about 120° especially by MP2 calculation, which is consistent with sp² (planar –NH–CN group) and not sp³ (pyramidal –NH–CN group) structure. The vibrational frequencies of the d₀, d₁ and d₃ species of 3,5-dichlorophenylcyanamide and the potential energy distributions among symmetry coordinates of the normal modes of the parent molecule were computed at the DFT-B3LYP level. The calculated infrared and Raman spectra of the molecule were plotted. Complete vibrational assignments were made on the basis of isotopic substitution and normal coordinate calculations.Figure Potential curves for the internal rotation in 3,5-dichlorophenylcyanamide as determined by DFT-B3LYP/6-311+G** (solid) and MP2/6-311+G** (dotted) calculations     

Molecular structure,vibrational spectra and HOMO,LUMO analysis of 4-piperidone by density functional theory and ab initio Hartree–Fock calculations

Vibrational frequencies and geometrical parameters of 4-piperidone (4-PID) in the ground state have been calculated by using the Hartree–Fock (HF) and density functional methods (B3LYP) with 6-311++G(d,p) and 6-311+G(3df,2p) basis sets. These methods are proposed as a tool to be applied in the structural characterisation of 4-PID (C₅H₉NO). The title molecule has C _s point group symmetry, thus providing useful support in the interpretation of experimental IR and Raman data. The DFT-B3LYP/6-311+G(3df,2p) calculations have been found more reliable than the ab initio HF/6-311++G(d,p) calculations for the vibrational study of 4-PID. The calculated highest occupied molecular orbital and lowest unoccupied molecular orbital energies show that charge transfer occurs within the molecule. The theoretical spectrograms for FT-IR and FT-Raman spectra of the title molecule have been constructed.     

Studies on molecular structure and tautomerism of a vitamin B<Subscript>6</Subscript> analog with density functional theory

This work presents a computational study on the molecular structure and tautomeric equilibria of a novel Schiff base L derived from pyridoxal (PL) and o-phenylenediamine by using the density functional method B3LYP with basis sets 6-31 G(d,p), 6-31++G(d,p), 6-311 G(d,p) and 6-311++G(d,p). The optimized geometrical parameters obtained by B3LYP/6-31 G(d,p) method showed the best agreement with the experimental values. Tautomeric stability study of L inferred that the enolimine form is more stable than its ketoenamine form in both gas phase and solution. However, protonation of the pyridoxal nitrogen atom (LH) have accelerated the formation of ketoenamine form, and therefore, both ketoenamine and enolimine forms could be present in acidic media.     

Measurement of diffusion coefficients of meningococcal polysaccharide by optical mixing spectroscopy. I. A preliminary characterization on the aggregation of the group C polysaccharide

The angular distribution of scattered intensity and decay times of concentration fluctuations have been measured by means of digital photon counting and single-clipped photon correlation for solutions of Group C meningococcal polysaccharides at 31°C. The z-average diffusion coefficient 〈D〉_z and its second moment 〈D²〉_z have been determined from the time-dependent correlation function using the cumulant expansion technique. Very low observed values of 〈D〉_z and the tremendous width of the polydispersity index, which is the z-average normalized variance, suggest a higher degree of aggregation than the monomer–dimer type self-association at finite concentrations.     

ONIOM DFT/PM3 calculation on the interaction between STI-571 and abelson tyrosine kinase

The ONIOM2 (B3LYP/6–31G (d, p): PM3) and B3LYP/6–31G (d, p) methods were applied to investigate the interaction between STI-571 and abelson tyrosine kinase binding site. The complex of N-[4-methyl-3-(4-pyridin-3-yl-pyrimidin-2-ylamino)- phenyl]-benzamide (part of STI-571) and related 16 amino acid residues were found at B3LYP/6–31G (d, p) level to have hydrogen bonds and π....π stacking interaction, their binding energy via HAF optimization was −20.4 kcal mol⁻¹. The results derived from this study agreed well with the reported observation. Figure Optimized structure of STI-571 and Thr315 in abelson tyrosine kinase based on ONIOM2 method     

Computational Study of the Sulfonylated Amino Acid Hydroxamates Binding to the Zinc Ion within the Active Site of Carbonic Anhydrase

Abstract

The B3LYP/6–311+G(d,p) method and three ONIOM extrapolation methods ONI-OM (B3LYP/6–311+G(d,p): AM1); ONIOM(B3LYP/6–311+G(d,p): MNDO); ONIOM (B3LYP/6–311+G(d,p): HF/3-21G(d)) were used to characterize the complexes of Zn²⁺ cation with anionic sulfonylated amino acid hydroxamates (RSO₂NH-AA-CON(-)OH), possessing an unsubstituted RSO₂NH—amino acyl moiety. According to the R moiety we distinguish between pentafluorophenyl and 4-methoxyphenyl derivates. The amino acid hydroxamates included in the study were the Gly, Ala, and Leu derivates. Of the inhibitors investigated, the weakest zinc affinity exhibits the pentafluorophenyl derivate with Gly amino acid and the strongest affinity the 4-methoxyphenyl derivate with Leu amino acid. The inhibitors form bidentate coordination bonds with the zinc cation by means of the sulfonyl oxygen and the ionized hydroxamate nitrogen atoms, respectively. The zinc affinities computed using the B3LYP/6–311 +G(d,p)//HF/6–31 +G(d,p) method are in very good agreement with the full density functional theory (DFT) B3LYP/6–311+G(d,p)//B3LYP/6- 311+G(d,p) method and this method can be adopted to model larger complexes of inhibitors with the active site of carbonic anhydrase.     

DFT studies on the intrinsic conformational properties of non-ionic pyrrolysine in gas phase

B3LYP/6-31G(d,p) level of theory is used to carry out a detailed gas phase conformational analysis of non-ionized (neutral) pyrrolysine molecule about its nine internal back-bone torsional angles. A total of 13 minima are detected from potential energy surface exploration corresponding to the nine internal back-bone torsional angles. These minima are then subjected to full geometry optimization and vibrational frequency calculations at B3LYP/6-31++G(d,p) level. Characteristic intramolecular hydrogen bonds present in each conformer, their relative energies, theoretically predicted vibrational spectra, rotational constants and dipole moments are systematically reported. Single point calculations are carried out at B3LYP/6-311++G(d,p) and MP2/6-31++G(d,p) levels. Six types of intramolecular H-bonds, viz. O…H–O, N…H-O, O…H–N, N…H–N, O…H–C and N…H–C, are found to exist in the pyrrolysine conformers; all of which contribute to the stability of the conformers. The vibrational frequencies are found to shift invariably toward the lower side of frequency scale corresponding to the presence of intramolecular H-bond interactions in the conformers.     

Density functional theory study of the potassium complexation of an unsymmetrical 1,3-alternate calix[4]-crown-5-N-azacrown-5 bearing two different crown rings

Theoretical studies of an unsymmetrical calix[4]-crown-5-N-azacrown-5 (1) in a fixed 1,3-alternate conformation and the complexes 1·K⁺(a), 1·K⁺(b), 1·K⁺(c) and 1·K⁺K⁺ were performed using density functional theory (DFT) at the B3LYP/6-31G* level. The fully optimized geometric structures of the free macroligand and its 1:1 and 1:2 complexes, as obtained from DFT calculations, were used to perform natural bond orbital (NBO) analysis. The two main types of driving force metal–ligand and cation–π interactions were investigated. NBO analysis indicated that the stabilization interaction energies (E ₂) for O…K⁺ and N…K⁺ are larger than the other intermolecular interactions in each complex. The significant increase in electron density in the RY* or LP* orbitals of K⁺ results in strong host–guest interactions. In addition, the intermolecular interaction thermal energies (ΔE, ΔH, ΔG) were calculated by frequency analysis at the B3LYP/6-31G* level. For all structures, the most pronounced changes in the geometric parameters upon interaction are observed in the calix[4]arene molecule. The results indicate that both the intermolecular electrostatic interactions and the cation–π interactions between the metal ion and π orbitals of the two pairs that face the inverted benzene rings play a significant role.     

Ionic interactions for substituted MCH1R inhibitors studied by pKa values

MCH1R inhibitors with the quinoline moiety having the aromatic amine and aliphatic amine chain were selected, and then the effect of substituents of the quinoline ring on the ionic interaction were studied by calculating pK_a values for these amines at the B3LYP/6-311++G(d,p)//B3LYP/6-31+G(d) level in the gas phase and in water. For substituent with C, N, O, and S atoms next to the quinoline ring, respectively, the pK_a values of aromatic amines are estimated to be 8.98, 12.19, 4.64, and 4.33 and those of the aliphatic amines are 12.65, 10.82, 9.94, and 11.55, respectively.     

Biological exploration of a novel 1,2,4-triazole-indole hybrid molecule as antifungal agent

Abstract

(2-(2,4-Dichlorophenyl)-3-(1H-indol-1-yl)-1-(1,2,4-1H-triazol-1-yl)propan-2-ol (8?g), a new 1,2,4-triazole-indole hybrid molecule, showed a broad-spectrum activity against Candida, particularly against low fluconazole-susceptible species. Its activity was higher than fluconazole and similar to voriconazole on C. glabrata (MIC₉₀ = 0.25, 64 and 1?µg/mL, respectively), C. krusei (MIC₉₀ = 0.125, 64 and 0.125?µg/mL, respectively) and C. albicans (MIC₉₀ = 0.5, 8 and 0.25?µg/mL, respectively). The action mechanisms of 8?g were also identified as inhibition of ergosterol biosynthesis and phospholipase A2-like activity. At concentration as low as 4 ng/mL, 8g inhibited ergosterol production by 82% and induced production of 14a-methyl sterols, that is comparable to the results obtained with fluconazole at higher concentration. 8?g demonstrated moderate inhibitory effect on phospholipase A2-like activity being a putative virulence factor. Due to a low MRC5 cytotoxicity, this compound presents a high therapeutic index. These results pointed out that 8?g is a new lead antifungal candidate with potent ergosterol biosynthesis inhibition.     

Ab initio and DFT study on the electrophilic addition reaction of bromine to tetracyclo[5.3.0.02,6.03,10]deca–4,8–diene

The electronic and geometric structures of tetracyclo[5.3.0.0^2,6.0^3,10]deca-4,8-diene (hypostrophene) have been investigated by ab initio and DFT/B3LYP methods using the 6-31G* and 6-311G* basis sets. The double bonds of hypostrophene are endo-pyramidalized. The cationic intermediates and products formed in the addition reaction have been investigated using the HF/6-311G*, HF/6-311G**, and B3LYP/6-311G* methods. The bridged bromonium cation was more stable than the U-type cation. Considering that the bridged cation does not isomerize to the less stable U-type cation, it is not possible for the U-type product to be obtained in the reaction. The bridged bromonium cation transformed into the more stable N-type cation and the N-type product was obtained via this cation. The thermodynamic stability of the exo, exo and exo, endo isomers of the N-type dibromide molecule were almost identical. The N-type product was 16.6 kcal mol⁻¹ more stable than the U-type product. Figure General energy diagram of the hypostrophene–bromine (HS–Br₂) system (kcal mol⁻¹) (MP2/6-311G*//HF/6-311G*)     

An Investigation of the Structural Diversities of Lithiated HMPA Complexes of o-Mercaptopyridine and Trithiocyanuric Acid: Syntheses, Crystal structures and Model Molecular Orbital Calculations

Reaction of o-mercaptopyridine (o-MPH) and trithiocyanuric acid (TTCyH₃) with one equivalent of BuⁿLi in the presence of HMPA yields the mono-lithiated salts MPLi.HMPA (1) and TTCyH₂Li.2HMPA (2) respectively, which have been characterised by NMR spectroscopy and X-ray crystallography. Reaction of three equivalents of BuⁿLi with anhydrous TTCyH₃ in THF yields the tri-lithiated species TTCyLi₃.4THF (3). In all three compounds the lithium centres have N,S-bridged coordination modes. Whereas 1 is dimeric in the solid state, 2 has an unusual monomeric structure and 3, which is a very rare example of a structurally characterised tri-lithiated compound, has an unprecedented polymeric structure incorporating (NCSLi) _n (n = 1, 2) rings. The structural diversities displayed by 1 and 2 have been probed, and thereby in part rationalised, by ab initio (6-31G*/RHF, 6-31G**/RHF and 6-31G*/MP2 levels) MO calculations on both their thio-keto and thiol isomers and on their uncomplexed and complexed lithiated derivatives. In particular, the optimised structures predict and reproduce the N,S-bridging coordination modes found for lithium and explain why structure 1 is dimeric whereas 2 is monomeric.Electronic Supplementary Material available.     

Subunit interaction enhances enzyme activity and stability of sweet potato cytosolic Cu/Zn-superoxide dismutase purified by a His-tagged recombinant protein method

The coding region of copper/zinc-superoxide dismutase (Cu/Zn-SOD) cDNA from sweet potato, Ipomoea batatas (L.) Lam. cv. Tainong 57, was introduced into an expression vector, pET-20b(+). The Cu/Zn-SOD purified by His-tagged technique showed two active forms (dimer and monomer). The amount of proteins of dimer and monomer appeared to be equal, but the activity of dimeric form was seven times higher than that of monomeric form. The enzyme was dissociated into monomer by imidazole buffer above 1.0 M, acidic pH (below 3.0), or SDS (above 1%). The enzyme is quite stable. The enzyme activity is not affected at 85 °C for 20 min, in alkali pH 11.2, or in 0.1 M EDTA and also quite resistant to proteolytic attack. Dimer is more stable than monomer. The thermal inactivation rate constant k _dcalculated for the monomer at 85 °C was 0.029 min^-1 and the half-life for inactivation was about 28 min. In contrast, there is no significant change of dimer activity after 40 min at 85 °C. The enzyme dimer and monomer retained 83% and 58% of original activity, respectively, after 3 h incubation with trypsin at 37 °C, while those retained 100% and 31% of original activity with chymotrypsin under the same condition. These results suggest subunit interaction might change the enzyme conformation and greatly improve the catalytic activity and stability of the enzyme. It is also possible that the intersubunit contacts stabilize a particular optimal conformation of the protein or the dimeric structure enhances catalytic activity by increasing the electrostatic steering of substrate into the active site.     

Synthesis of a Carboxamide Linked T∗T Dimer with an Acyclic Nucleoside Unit and Its Incorporation in Oligodeoxynucleotides

Abstract

A T?T dimer with ? representing a 2′-OCH₂CH₂NHC(O)-4′ linkage connecting two nucleoside units was prepared by condensation of (S)-1-[2-(2-aminoethoxy)-3-(4,4′-dimethoxytrityloxy)propyl]thymine with 1,2-dideoxy-1-thyminyl-β-D-erythro-pento-furanuronic acid. The T?T dimer was incorporated in oligodeoxynucleotides and investigated for hybridization to DNA.