Probing the effect of electric field in 9,10-dimethoxy-2,6-bis(2-p-tolylethynyl)anthracene molecular nanowire using quantum chemical and charge density analysis

The effect of electric field (EF) in a newly designed molecular nanowire 9,10-dimethoxy-2,6-bis(2-p-tolylethynyl)anthracene has been analysed theoretically from the structural and electronic charge transport properties using quantum chemical and charge density calculations. The applied EF (0–0.36 VÅ^? 1) alters the molecular conformation, charge density distribution, electrostatic properties and the electronic energy levels of the molecule. Furthermore, the applied EF decreases the highest occupied molecular orbital–lowest unoccupied molecular orbital gap significantly from 1.775 to 0.258 eV and it also induces polarisation in the molecule, which leads to increase the dipole moment of the molecule. The electrostatic potential for various levels of applied EF reveals the charge-accumulated regions of the molecule. The I–V characteristics of the molecule have been studied against various applied fields using Landauer formalism.     

Model for the conformational analysis of hydrated peptides. Effect of hydration on the conformational stability of the terminally blocked residues of the 20 naturally occurring amino acids

The effects of hydration are included in empirical conformational energy computations on oligopeptides by means of a modified hydration-shell model. Free energy terms are introduced to account for “specific hydration” due to water–solute hydrogen bonding and for “nonspecific hydration” describing the interaction of the solute with water molecules in a first-neighbor shell. The dielectric constant has been doubled (over the value used for calculations in the absence of water) to take into account the presence of solvent. Computations were carried out for the N-acetyl-N′-methylamides of the 20 naturally occurring amino acids. Conformational energy maps are compared with similar maps calculated in the absence of hydration. Minimum-energy conformations are located and compared with the corresponding minima for unhydrated peptides in terms of ordering with respect to potential energy, the dihedral angles at the minima, and the presence of intramolecular hydrogen bonds. The Boltzmann factors for various conformational regions are altered significantly on hydration in some cases. These changes can be explained in terms of differences in the hydration free energy terms for various conformations.     

The conformational and vibrational behavior of the inhibitory neuropeptide derived from beta-endorphin

In this study, conformational behavior, structural, and vibrational characterization of the carboxy terminal dipeptide of β-endorphin (glycy-l-glutamine, glycyl-glutamine, beta-endorphin_30-31), which is an inhibitory neuropeptide synthesized from beta-endorphin_1-31 in brain stem regions, has been investigated. The theoretically possible stable conformers were searched by means of molecular mechanics method to determine their energetically preferred conformations. The 360 different conformations were calculated with the φ, Ψ, χ dihedral angles using the Ramachandran maps. The most stable conformation of the title molecule is characterized by the extended backbone shape (e) in the BR conformational range with ?.78 kcal/mol energy. The cis- and trans-dimeric forms of the dipeptide were also formed and energetically preferred conformations of dimers were investigated. The experimental methods (FT-IR, micro-Raman spectroscopies) coupled with quantum chemical calculations based on density functional theory (DFT) have been used to identify the geometrical, energetic, and vibrational characteristics of the dipeptide. The assignment of the vibrational spectra was performed based on the potential energy distribution of the vibrational modes. To investigate the electronic properties, such as nonlinear optical properties, the electric dipole moment, the mean polarizability, the mean first hyperpolarizability, and HOMO–LUMO energy gaps were computed using the DFT with the B3LYP/6-31++G(d,p) basis set combination. The second-order interaction energies were derived from natural bonding orbital analysis. The focus of this study is to determine possible stable conformation on inhibitory neuropeptide and to investigate molecular geometry, molecular vibrations of monomeric and dimeric forms, and hydrogen bonding interactions of glycy-l-glutamine dipeptide.     

Electronic structure and conformational properties of <Emphasis Type="Italic">1H</Emphasis>-indole-3-acetic acid

The conformational space of 1H-Indole-3-Acetic Acid (IAA) was scanned using molecular dynamics at semiempirical level, and complemented with functional density calculations at B3LYP/6-31G** level, 14 conformers of lowest energy were obtained. Electronic distributions were analyzed at a higher calculation level, thus improving the basis set (B3LYP/6-311++G**). A topological study based on Bader’s theory (AIM: atoms in molecules) and natural bond orbital (NBO) framework performed with the aim to analyze the stability and reactivity of the conformers allowed the understanding of electronic aspects relevant in the study of the antioxidant properties of IAA. Intramolecular hydrogen bonds were found and were characterized as blue-shifting hydrogen bonding interactions. Furthermore, molecular electrostatic potential maps (MEPs) were obtained and analyzed in the light of AIM and NBO results, thus showing subtle but essential features related not only to reactivity but also with intramolecular weak interactions, charge delocalization and structure stabilization.     

Molecular structure and spectroscopic (FT-IR,FT-Raman, 13C, 1H NMR and UV) studies of 3,4-dihydroxy-l-phenylalanine using density functional theory

Fourier transform infrared (FT-IR) and Fourier transform Raman (FT-Raman) spectra of 3,4-dihydroxy-l-phenylalanine (3,4-DPA) in solid phase were recorded and analysed in this research. Along with this, the IR spectra in CHCl₃ and the use of acetone as solvents of 3,4-DPA were also recorded. The equilibrium geometry, bonding features and harmonic vibrational frequencies were investigated with the help of density functional theory (DFT) method. The ¹H and ¹³C nuclear magnetic resonance chemical shifts of the molecule were calculated using the gauge including atomic orbital method and compared with experimental results. Stability of the molecule arising from hyperconjugative interactions and charge delocalisation was analysed using natural bond orbital analysis. The results show that charge in electron density (E _D) in the σ* and π* antibonding orbitals and second-order delocalisation energies E(2) confirms the occurrence of intramolecular charge transfer within the molecule. UV–vis spectrum of the compound was recorded and the electronic properties, such as HOMO and LUMO energies, were analysed using the time-dependent (TD)-DFT approach. Finally, the calculation results were applied to simulate infrared and Raman spectra of the title compound, which showed good agreement with the observed spectra.     

Computer simulation of the conformational properties of retro–inverso peptides. I. Empirical force field calculations of rigid and flexible geometries of N-acetylglycine-N′- methylamide,bis(acetamido) methane,and N,N′- dimethylmalonamide and their corresponding Cα-methylated analogs

Rigid and flexible geometry calculations are described for N-acetylglycine-N′-methylamide, N-acetylalanine-N′-methylamide, and their retro-inverso analogs, bis(acetamido) methane, 1,1-bis(acetamido) ethane, N,N′-dimethylmalonamide, and N,N′-dimethyl-2-methyl-malonamide. The significance of relaxing all degrees of freedom, especially angular flexibility is demonstrated. The flexible geometry approach yields energy maps similar to those from rigid geometry, but the energy barriers between minima are substantially reduced, leading in general, to more probable transitions and a higher volume of accessible conformational space. Whereas the glycine and alanine derivatives exhibit their lowest energy minima in the C region, the gem-diaminoalkyl and malonyl residues show their lowest minima in the “α-helical” regions. With respect to the effect of side chains (H versus CH₃), the greatest conformational influence appears with the gem-diaminoalkyl residues. These results indicate significantly different conformational behavior of retro peptides and the implications of these pairwise incorporations of retro-inverso residues in peptide chains, are discussed.     

Structural analysis,spectroscopic characterization and molecular docking studies of the cyclic heptapeptide

The theoretically possible stable conformer of the cyclic heptapeptide, that has significant anti-metastatic activity, was examined by conformational analysis followed by DFT calculations. Experimental infrared and Raman spectroscopy, together with theoretical DFT (6-31G (d,p) basis set)-based quantum chemical calculations, have been used to understand the structural and spectral characteristics of cyclo(Gly-Arg-Gly-Asp-Ser-Pro-Ala) {cyclo(GRGDSPA)}. A complete analysis of the vibrational spectrum has been reported on the basis of potential energy distribution (PED%) data of the vibrational modes. Finally, the calculation results were applied to simulate infrared and Raman spectra of the title compound. The simulated spectra satisfactorily coincide with the experimental spectra. In addition, molecular electrostatic potential and frontier molecular orbital analysis were investigated using theoretical calculations. The stability of the molecule, arising from hyperconjugative interaction and charge delocalization, has been analyzed using natural bond orbital analysis and a high E(2) value reveals the presence of strong interaction between donors and acceptors. Molecular docking studies with fibronectin were performed on cyclo(GRGDSPA) in order to understand its inhibitory nature. The results indicate that the docked ligand {cyclo(GRGDSPA)} forms a stable complex with human fibronectin and gives a binding affinity value of ?7.7 kcal/mol, which points out that cyclo(GRGDSPA) might exhibit inhibitory activity against the attachment of melanoma cells to human fibronectin.     

Structural and vibrational spectroscopic elucidation of sulpiride in solid state

The study on the conformational and vibrational behaviors of sulpiride molecule which is known as a neuroleptic or antipsychotic drug that is widely used clinically in the treatment of schizophrenic or depressive disorders is an important scientific and practical task. In here, a careful enough study of monomer and dimeric forms of sulpiridine {5-(aminosulfonyl)-N-[(1-ethyl-2-pyrrolidinyl) ethyl]-2-methoxy-benzamide (C₁₅H₂₃N₃O₄S)} is undertaken by density functional theory (DFTB3LYP) method with the B3LYP/6-31G(d,p) basis set. The conformations of free molecule were searched by means of torsion potential energy surfaces scan studies through dihedral angles D₁ (8?N, 18C, 20C, 23?N), D₂ (18C, 20C, 23?N, 25C) and D₃ (28C, 30C, 41S, 44?N) in electronically ground state, employing 6-31G basic set. The final geometrical parameters for the obtained stable conformers were determined by means of geometry optimization, carried out at DFT/B3LYP/6-31G(d,p) theory level. Afterwards, the possible dimer forms of the molecule were formed and their energetically preferred conformations were investigated. Moreover, the effect of basis set superposition error on the structure and energy of the three energetically favourable sulpiride dimers has been determined. The optimized structural parameters of the most stable monomer and three low energy dimer forms were used in the vibrational wavenumber calculations. Raman and IR (4000–400?cm^?1) spectra of sulpiride have been recorded in the solid state. The assignment of the bands was performed based on the potential energy distribution data. The natural bond orbital analysis has been performed on both monomer and dimer geometries in order to elucidate delocalization of electron density within the molecule. The predicted frontier molecular orbital energies at DFT/B3LYP/6-31G(d,p) theory level show that charge transfer occurs within the molecule. The first-order hyperpolarizability (β₀) and related properties (μ and α) of the title molecule were also calculated.     

Theoretical Study of Antagonists and Inhibitors of Mammalian Adenosine Deaminase: I. Adenosine and Its Aza- and Deazaanalogues

Aza- and deazaanalogues of adenosine, including their 1-protonated forms (except for that of 1-deazaadenosine), were studied by computer computation to find a relationship between their molecular structures and substrate properties for the mammalian adenosine deaminase. The atomic charge distribution and maps of the electrostatic potential around their van der Waals molecular surface were calculated for these compounds using the ab initio STO-3G method. The conformational studies are carried out by the MM⁺ method of molecular mechanics. The mechanism that determines the substrate selectivity of mammalian adenosine deaminase is discussed.     

Understanding the conformational flexibility and electrostatic properties of curcumin in the active site of rhAChE via molecular docking,molecular dynamics,and charge density analysis

Acetylcholinesterase (AChE) is an important enzyme responsible for Alzheimer’s disease, as per report, keto-enol form of curcumin inhibits this enzyme. The present study aims to understand the binding mechanism of keto-enol curcumin with the recombinant human Acetylcholinesterase (rhAChE) from its conformational flexibility, intermolecular interactions, charge density distribution, and the electrostatic properties at the active site of rhAChE. To accomplish this, a molecular docking analysis of curcumin with the rhAChE was performed, which gives the structure and conformation of curcumin in the active site of rhAChE. Further, the charge density distribution and the electrostatic properties of curcumin molecule (lifted from the active site of rhAChE) were determined from the high level density functional theory (DFT) calculations coupled with the charge density analysis. On the other hand, the curcumin molecule was optimized (gas phase) using DFT method and further, the structure and charge density analysis were also carried out. On comparing the conformation, charge density distribution and the electrostatic potential of the active site form of curcumin with the corresponding gas phase form reveals that the above said properties are significantly altered when curcumin is present in the active site of rhAChE. The conformational stability and the interaction of curcumin in the active site are also studied using molecular dynamics simulation, which shows a large variation in the conformational geometry of curcumin as well as the intermolecular interactions.     

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.     

Conformational flexibility of peptides containing α,β-unsaturated amino acid residues. I. Conformational analysis of N-acetyl-N′-methylamides of dehydroalanine and N-methyldehydroalanine

Conformational energy calculations on the N-acetyl-N′-methylamides of dehydroalanine and N-methyldehydroalanine indicate that their conformational behavior is very different from that of the corresponding saturated compounds. The conformational data in the literature from x-ray and nmr investigations on peptides containing α,β-unsaturated residues are discussed on the basis of these theoretical results.     

GROmaρs: A GROMACS-Based Toolset to Analyze Density Maps Derived from Molecular Dynamics Simulations

We introduce a computational toolset, named GROmaρs, to obtain and compare time-averaged density maps from molecular dynamics simulations. GROmaρs efficiently computes density maps by fast multi-Gaussian spreading of atomic densities onto a three-dimensional grid. It complements existing map-based tools by enabling spatial inspection of atomic average localization during the simulations. Most importantly, it allows the comparison between computed and reference maps (e.g., experimental) through calculation of difference maps and local and time-resolved global correlation. These comparison operations proved useful to quantitatively contrast perturbed and control simulation data sets and to examine how much biomolecular systems resemble both synthetic and experimental density maps. This was especially advantageous for multimolecule systems in which standard comparisons like RMSDs are difficult to compute. In addition, GROmaρs incorporates absolute and relative spatial free-energy estimates to provide an energetic picture of atomistic localization. This is an open-source GROMACS-based toolset, thus allowing for static or dynamic selection of atoms or even coarse-grained beads for the density calculation. Furthermore, masking of regions was implemented to speed up calculations and to facilitate the comparison with experimental maps. Beyond map comparison, GROmaρs provides a straightforward method to detect solvent cavities and average charge distribution in biomolecular systems. We employed all these functionalities to inspect the localization of lipid and water molecules in aquaporin systems, the binding of cholesterol to the G protein coupled chemokine receptor type 4, and the identification of permeation pathways through the dermicidin antimicrobial channel. Based on these examples, we anticipate a high applicability of GROmaρs for the analysis of molecular dynamics simulations and their comparison with experimentally determined densities.     

Electrostatic components of drug-receptor recognition. II. The DNA-binding antibiotic actinomycin

The semiempirical molecular orbital method PCILO has been used to study the electronic charge distribution and conformation of the antibiotic actinomycin. Molecular electrostatic potential fields have been constructed in three-dimensional space round the molecule and displayed stereoscopically together with contour maps in the orthogonal planes relative to that of the chromophore. It was found that the drug has a considerable dipole moment and that the electrostatic fields can be separated into two large regions; an area of positive potential surrounds the chromophore and one of negative potential is directed away from the peptide rings. In this and the preceding paper these observations are discussed with respect to drug-receptor recognition phenomena. Pattern-matching of complementary electrostatic fields between the drug and polynucleotide receptor can be discerned.     

Optical and conformational properties of the 7-(beta-D-ribofuranosyl) purines

The optical and conformational properties of certain 7-ribosyl purines have been studied using several experimental and theoretical methods. The CD, MCD, and absorption spectra of the 7-ribosyl purines and their cations have been determined experimentally with some solvent effects being noted. Theoretically we have used the SCF-CI and CNDO molecular orbital calculations to determine the optical properties of the bases; Simpson's bond exciton theory to determine the optical properties of the bases; Simpson's bond exciton theory to determine the optical properties of the ribose moiety; the optical activity was determined by the dipole coupling and electricmagnetic coupling equations; and the conformational energy calculations are included as a basis for excluding highly improbable conformations. On the basis of these calculations, we concluded that the oxy derivatives and the amino derivatives are predominately in the “extreme anti” conformation and “standard anti” conformation, respectively. However, the conclusions are tentative since at the present time the band assignments are not unequivocal. At least one n-π* transition has been identified in the spectra which is not notably sensitive to pH and solvent effects. Several implications of this observation are discussed.     

Comparison of non-relativistic and quasi-relativistic SCF-Xα scattered-wave calculations of uranocene,thorocene and cerocene

The results of non-relativistic and quasi-relativistic SCF-Xα-SW calculations of uranocene, thorocene and cerocene are analyzed on the basis of selected orbital energies and charge distributions. The character of the ring–metal bonding is discussed. The f-orbital expansion found in previous non-relativistic molecular calculations is interpreted as a consequence of the Xα approximation by comparison of Xα model ground states of the corresponding metal atoms.     

Computer simulation of the conformational properties of retro-inverso peptides. II. Ab initio study, spatial electron distribution, and population analysis of N-formylglycine methylamide, N-formyl N'-acetyldiaminomethane, and N-methylmalonamide

Ab initio minimal and split-valence basis set calculations have been performed on compounds that are involved in retro–inverso modifications, i.e., gem-diaminoalkyl and malonyl structures. These calculations are compared with empirical force field calculations and the minor differences discussed. All calculations agree that the preferred helical conformation of the isolated gem-diaminoalkyl and malonyl derivatives of residues found in the retro-inverso modified peptides is 5–8 kcal/mol lower than the C conformation preferred by the isolated peptide residues. Population analysis and contour plots of the charge distribution are used to help explain the differences between the model compounds.     

Experimental and theoretical investigations on the folding modes of depsipeptide molecules

The coupling of N-acyl-α-amino-acids with α-hydroxyacid-methyl amides results in depsipeptide molecules containing two chiral centers and one ester function inserted between two amide functions. Their conformational features have been investigated by IR spectroscopy, proton magnetic resonance, X-ray diffraction, and theoretical P.C.I.L.O. calculations. It is shown that most of these molecules are folded by an intramolecular 4 → 1 hydrogen bonding. Two folded conformations, similar to the well known β turn in peptides, are described, the stability of which depends on the configurational sequence in the investigated molecule. LL and LD species are folded in two different ways whereas LG sequences containing an achiral hydroxy-acid residue accommodate both of them. The presence of a N-terminal achiral amino acid noticeably decreases the folding ratio. The above conclusions are then compared with the conformational features of homologous tripeptide molecules.     

Quantitative structure activity relationships of some pyridine derivatives as corrosion inhibitors of steel in acidic medium

Quantum chemical calculations using the density functional theory (B3LYP/6-31G* DFT) and semi-empirical AM1 methods were performed on ten pyridine derivatives used as corrosion inhibitors for mild steel in acidic medium to determine the relationship between molecular structure and their inhibition efficiencies. Quantum chemical parameters such as total negative charge (TNC) on the molecule, energy of highest occupied molecular orbital (E _HOMO), energy of lowest unoccupied molecular orbital (E _LUMO) and dipole moment (μ) as well as linear solvation energy terms, molecular volume (Vi) and dipolar-polarization (π*) were correlated to corrosion inhibition efficiency of ten pyridine derivatives. A possible correlation between corrosion inhibition efficiencies and structural properties was searched to reduce the number of compounds to be selected for testing from a library of compounds. It was found that theoretical data support the experimental results. The results were used to predict the corrosion inhibition of 24 related pyridine derivatives.     

Optical rotatory properties of L-cystine conformational isomers. Theoretical analysis

The near ultraviolet chiroptical properties of L -cystine conformational isomers are examined on a static, “one-electron” model in which the disulfide moiety is the chromophoric group and the atoms of the L -alanyl fragments are treated as perturbers. The zeroth order chromophoric wave functions are calculated on a semiempirical molecular orbital model in which excited states are constructed in the virtual orbital-configuration interaction approximation. The perturbing environment is represented by point charges located at the atomic centers of the L -alanyl fragments. Chromophore–perturber interactions are expressed as charge–multipole moments with only the charge–dipole and charged–quadrupole terms being retained in the calculations. Vicinal contributions to the rotatory strengths of the five lowest energy disulfide transitions are computed for 30 conformational isomers of the L -cystine dizwitterion. The results provide support for the view that vicinal or peripheral effects can account entirely for the observed near ultraviolet (λ > 230 nm) chiroptical properties of L -cystine and its derivatives and that these properties are diagnostic of conformational features external to the disulfide moiety.