Structure-Activity Relationship for Some 2′,3′-Dideoxynucleoside Anti-HIV Drugs Using Molecular Electrostatic Potential Mapping

Molecular electrostatic potential(MEP) maps of azido thymidine (AZT), some of its analogs and derivatives and certain other 2′,3′-dideoxy nucleosides having different anti-HIV activities have been studied. The optimised hybridization displacement charges (HDC) combined with MNDO Löwdin charges, continuosly distributed in three dimension spherically symmetrically as a Slater cloud at each site were used to compute the MEP maps. The negative MEP region near the O5′ sites of these molecules appears to be of primary importance from the point of view of their anti-HIV activity. The roles of the azido group in AZT and fluorine atoms substituted at different positions in the sugar moiety have been evaluated. The azido group in AZT behaves as a strongly electronegative group.     

Hybridization-displaced charges for amino-acids: a new model using two point charges per atom along with bond-center charges

A new charge distribution is proposed for the amino acids where each atom is associated with two point charges while each bond center is associated with one point charge. Centroids of charges arising due to atomic orbital hybridization called hybridization-displaced charges (HDC) and those located at the atomic sites and bond centers obtained by a modified form of the Mulliken scheme were combined. The density matrix calculations required for this analysis were performed at the B3LYP/6-31G** level of density functional theory. The combination of HDC centroids with the modified Mulliken charges was found to yield dipole moments and surface molecular electrostatic potentials (MEP) of the amino acids in good agreement with those obtained by rigorous DFT calculations or those obtained using the MEP-fitted CHelpG charges. This study shows that the combination of HDC centroids with the modified Mulliken charges is significantly superior to the conventional Mulliken charges.     

A theoretical study of some new analogues of the anti-cancer drug camptothecin

The enzyme topoisomerase I (topo I), which is essential for cell replication, transiently causes a DNA single strand break and makes a complex with it. The anti-cancer agent camptothecin (CPT) binds to the topo I–DNA complex and stabilizes it, preventing resealing of the broken DNA strand and cell growth. Considering the structural factors of CPT that are believed to be involved in stabilizing the topo I–DNA complex via hydrogen bonding and stacking interactions, designs of two new analogues of CPT (topo I inhibitors) have been suggested. The molecular geometries of CPT, two of its analogues and certain other related molecules included in the study were fully optimized in both gas phase and aqueous media at the B3LYP/6-311++G(d,p) level of density functional theory. Solvation effects of aqueous media were treated using the polarizable continuum model (PCM). Net CHelpG charges and surface molecular electrostatic potentials (MEP) near the atomic sites of the molecules were studied. Structural analogy and surface MEP values suggests that the two new CPT analogues studied here would be potent topoisomerase I inhibitors. Figure Optimized structures of CPT and two of its new analogues, 10 and 11.     

Improved electrostatic properties using combined Mulliken and hybridization-displaced charges for radicals

Effects of explicit consideration of charges displaced from atomic sites due to atomic orbital hybridization called hybridization-displaced charges (HDC) on dipole moments and surface molecular electrostatic potentials of certain radicals and their complexes with closed-shell molecules have been studied. HDC were computed for several radicals and their complexes at the B3LYP/6–31G** level of theory. At this level, HDC consist of three point charges associated with hydrogen atoms and seven point charges associated with heavy atoms belonging to the second row of the periodic table. HDC are so calculated that the contribution of each atom to the component of molecular dipole moment arising due to atomic orbital hybridization is preserved. It is found that dipole moments and electrostatic potentials of the systems studied here can be obtained with a significantly improved accuracy using a combination of Mulliken charges and HDC over that obtained by Mulliken charges only. Figure Surface MEP map of H₂O-HO· radical complex obtained using Mulliken charges combined with HDC     

A Molecular Electrostatic Potential Mapping Study of Some Fluoroquinolone Anti-Bacterial Agents

Molecular electrostatic potential (MEP) maps of some fluoroquinolones having varying degrees of activity against the bacterium Staphylococcus Aureus have been studied using the optimized hybridization displacement charges (HDC) combined with Löwdin charges obtained by the AM1 method. The roles of different substitutions at the N₁-position in the parent quinolone ring have been studied. The conformation of the carboxylic group attached to the quinolone ring was shown to be such that there is an intramolecular hydrogen bonding between the hydrogen atom of this group and the oxygen atom of the carbonyl group of the quinolone moiety. The carbonyl oxygen atom of the quinolone moiety, hydroxyl oxygen atom of the carboxylic group and the terminal nitrogen atom of the piperazin ring attached to the quinolone ring appear to be involved in the action of the drugs through electrostatic interactions while the N₁-alkyl substituents seem to be involved in the same through hydrophobic interactions.     

Optimization of binding electrostatics: charge complementarity in the barnase-barstar protein complex

Theoretical and experimental studies have shown that the large desolvation penalty required for polar and charged groups frequently precludes their involvement in electrostatic interactions that contribute strongly to net stability in the folding or binding of proteins in aqueous solution near room temperature. We have previously developed a theoretical framework for computing optimized electrostatic interactions and illustrated use of the algorithm with simplified geometries. Given a receptor and model assumptions, the method computes the ligand-charge distribution that provides the most favorable balance of desolvation and interaction effects on binding. In this paper the method has been extended to treat complexes using actual molecular shapes. The barnase-barstar protein complex was investigated with barnase treated as a target receptor. The atomic point charges of barstar were varied to optimize the electrostatic binding free energy. Barnase and natural barstar form a tight complex (K(d) approximately 10(-14) M) with many charged and polar groups near the interface that make this a particularly relevant system for investigating the role of electrostatic effects on binding. The results show that sets of barstar charges (resulting from optimization with different constraints) can be found that give rise to relatively large predicted improvements in electrostatic binding free energy. Principles for enhancing the effect of electrostatic interactions in molecular binding in aqueous environments are discussed in light of the optima. Our findings suggest that, in general, the enhancements in electrostatic binding free energy resulting from modification of polar and charged groups can be substantial. Moreover, a recently proposed definition of electrostatic complementarity is shown to be a useful tool for examining binding interfaces. Finally, calculational results suggest that wild-type barstar is closer to being affinity optimized than is barnase for their mutual binding, consistent with the known roles of these proteins.     

Beyond CHELP: improved potential derived charges for sugars

The partitioning of the overall molecular charge distribution into atom centered monopole charges, while quantum mechanically ill-defined, is nevertheless a technique which finds applications in several broad classes of chemical problems. Charges derived from fits to electrostatic potentials have an intuitive appeal since, in principle, these could be derived from either theoretical or experimental data. It has been noted, however, that such potential derived charges can be conformationally dependent in ways that do not appear to reflect the changes in the molecular wavefunction. Both the algorithm used for selecting points at which the molecular electrostatic potential will be fit and the density of points used in the fit have been suggested to influence the resultant charges. Recently [Stouch TR, Williams DE (1992) J Comp Chem 13: 622–32; Stouch TR, Williams DE (1993) J Comp Chem 14: 858–66] it has been noted that numerical difficulties may make it impossible to fit all the atomic charges in a molecule. Singular value decomposition (SVD) of the linear least squares matrices used in fitting atom based monopoles to molecular electrostatic potentials provides a tool for evaluating the integrity of the calculated charges. Based on the SVD analysis for a selected group of molecules we have noted particularly that increasing the molecular size reduces the fraction of charges which can be validly assigned. Users of PD derived charges, especially those who are using those charges for tasks other than reproduction of the MEP, should be aware that there is a high probability that a significant portion of those charges are statistically unreliable. Therefore, charges in many biological molecules, such as sugars, prove to be difficult to obtain by potential derived (PD) methods such as CHELP or CHELPG. Results from the SVD can be used to both assess PD charges and to generate an improved, albeit incomplete, set. Improved PD fits are presented for a series of simple saccharides. Abbreviations: HF, Hartree-Fock; LLS, linear least squares; MEP, molecular electrostatic potential; PD, potential derived; SVD, singular value decomposition This revised version was published online in November 2006 with corrections to the Cover Date.     

Enzyme selectivity of new cyclooxygenase-2/5 lipoxygenase inhibitors using molecular modeling approach

We have studied the conformational flexibility of three 5-keto-substituted 7-tert-butyl-2,3-dihydro-3,3-dimethylbenzofurans (DHDMBFs) which show dual cyclooxygenase (COX) and 5-lipoxygenase (LOX) inhibition and are potential candidates as antiinflammatory agents and analgesics. The conformations were studied by systematic search, molecular mechanics (MM) and simulated annealing molecular dynamics (SAMD) techniques. We also studied several structure based parameters and distribution of molecular electrostatic potential (MEP) around these molecules. All the three compounds were docked in the active cavity of cyclooxygenase-2 (COX-2) using graphical and energy grid search techniques. The complex geometries were optimized by MM. The results on conformational flexibility, inter-atomic distances and angles, MEP distribution and points of contacts with peptide side chains in active cavity have been used to understand the mechanistic cause of differential action of these molecules.     

Multipole electrostatic potential derived atomic charges in NDDO-methods with <Emphasis Type="Italic">spd</Emphasis>-basis sets

The recently introduced multipole approach for computing the molecular electrostatic potential (MEP) within the semiempirical neglect of diatomic differential overlap (NDDO) framework [Horn AHC, Lin Jr-H., Clark T (2005) Theor Chem Acc 114:159–168] has been used to obtain atomic charges of nearly ab initio quality by scaling the semiempirical MEP. The parameterization set comprised a total of 797 compounds and included not only the newly parameterized AM1* elements Al, Si, P, S, Cl, Ti, Zr, and Mo but also the standard AM1 elements H, C, N, O and F. For comparison, the ZDO-approximated MEP was also calculated analytically in the spd-basis. For the AM1*-optimized structures, single-point calculations at the B3LYP, HF and MP2 levels with the 6-31G(d) and LanL2DZP basis sets were performed to obtain the MEP. The regression analysis of all 12 combinations of semiempirical and ab initio MEP data yielded correlation coefficients of at least 0.99 in all cases. Scaling the analytical and multipole-derived semiempirical MEP by the regression coefficients yielded mean unsigned errors below 2.6 and 1.9 kcal mol⁻¹, respectively. Subsequently, for 22 drug molecules from the World Drug Index, atomic charges were computed according to the RESP procedure using XX/6-31G(d) (XX=B3LYP, HF, MP2) and scaled AM1* multipole MEP; the correlation coefficients obtained are 0.83, 0.85 and 0.83, respectively. Figure: Schematic representation of the atomic charge generation: The molecular electrostatic potential (MEP) is calculated using the AM1* Hamiltonian; then the semiempirical MEP is scaled to DFT or ab initio level, and atomic charges are generated subsequently by the restraint electrostatic potential (RESP) fit method. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible to authorized users. Proceedings of “Modeling Interactions in Biomolecules II”, Prague, September 5th–9th, 2005.     

Molecular similarity analysis between insect juvenile hormone and N, N-diethyl-m-toluamide (DEET) analogs may aid design of novel insect repellents

Molecular similarity analysis of stereoelectronic properties between natural insect juvenile hormone (JH), -a synthetic insect juvenile hormone mimic (JH-mimic, undecen-2-yl carbamate), and N, N-diethyl-m-toluamide (DEET) and its analogs reveals similarities that may aid the design of more efficacious insect repellents and give a better insight into the mechanism of repellent action. The study involves quantum chemical calculations using the AM1 semi-empirical computational method enabling a conformational search for the lowest and most abundant energy conformers of JH, JH-mimic, and 15 DEET compounds, followed by complete geometry optimization of the conformers. Similarity analyses of stereoelectronic properties such as structural parameters, atomic charges, dipole moments, molecular electrostatic potentials, and highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energies were performed on JH, JH-mimic and the DEET compounds. The similarity of stereoelectronic attributes of the amide/ester moiety, the negative electrostatic potential regions beyond the van der Waals surface, and the large distribution of hydrophobic regions in the compounds appear to be the three important factors leading to a similar interaction with the JH receptor. The similarity of electrostatic profiles beyond the van der Waals surface is likely to play a crucial role in molecular recognition interaction with the JH receptor from a distance. This also suggests electrostatic bioisosterism of the amide group of the DEET compounds and JH-mimic and, thus, a model for molecular recognition at the JH receptor. The insect repellent property of the DEET analogs may thus be attributed to a conflict of complementarity for the JH receptor binding sites.     

An experimental and theoretical approach to 5alpha-cholestan-6-spiro-1',2',4'-triazolidine-3'-one

The 5alpha-cholestan-6-one semicarbazone (1) on reaction with hydrogen peroxide at 0 degrees C affords selectively 5alpha-cholestan-6-spiro-1',2',4'-triazolidine-3'-one. (2) The structural assignment of the product was confirmed on the basis of its elemental, analytical and spectral analysis. The Hartree-Fock method using 6-31G* basis set was employed in order to explore the reaction mechanism. The results of the computational study show that the reaction proceeds through two radical intermediates formation. The different characteristics involved during the reaction were explained, firstly, the lower energy conformation of each molecule using total energy, hardness and dipole moment, and secondly, the explanation of the free radical mechanism, using frontier molecular orbital (FMO) theory, encoded electrostatic potential, spin electronic density and atomic charges. The localization of highest occupied molecular orbital (HOMO) or alpha-HOMO, lowest unoccupied molecular orbital (LUMO) or alpha-LUMO and the flow of atomic charges are in good agreement to support the present mechanism of the reaction. Stability and feasibility of all the optimized structures were supported by their respective fundamental frequencies and energy minima.     

On contribution of known atomic partial charges of protein backbone in electrostatic potential density maps

Partial charges of atoms in a molecule and electrostatic potential (ESP) density for that molecule are known to bear a strong correlation. In order to generate a set of point‐field force field parameters for molecular dynamics, Kollman and coworkers have extracted atomic partial charges for each of all 20 amino acids using restrained partial charge‐fitting procedures from theoretical ESP density obtained from condensed‐state quantum mechanics. The magnitude of atomic partial charges for neutral peptide backbone they have obtained is similar to that of partial atomic charges for ionized carboxylate side chain atoms. In this study, the effect of these known atomic partial charges on ESP is examined using computer simulations and compared with the experimental ESP density recently obtained for proteins using electron microscopy. It is found that the observed ESP density maps are most consistent with the simulations that include atomic partial charges of protein backbone. Therefore, atomic partial charges are integral part of atomic properties in protein molecules and should be included in model refinement.     

Molecular surface electrostatic potentials as guides to Si-O-N angle contraction: tunable σ-holes

We have demonstrated that the variation in the experimentally-determined Si-O-N angles in XYZSi-O-N(CH₃)₂ molecules, which depends upon the positions and natures of the substituents X, Y and Z, can be explained in terms of computed electrostatic potentials on the molecular surfaces of the corresponding XYZSi-H molecules. The latter framework has been used as a model for what the nitrogen lone pair in the XYZSi-O-N(CH₃)₂ molecules sees. Both optimized geometries and electrostatic potentials of our model XYZSi-H systems have been obtained at the B3PW91/6-31G(d,p) level. We propose that the driving force for the observed Si-O-N angle contraction in XYZSi-O-N(CH₃)₂ molecules is largely the electrostatic attraction between a positive σ-hole on the silicon and the lone pair of the nitrogen. Negative regions that may be near the silicon σ-hole, arising from substituents with negative potentials, also play an important role, as they impede the approach of the nitrogen lone pair. These two factors work in synergy and attest to the electrostatically-driven nature of the Si---N intramolecular interactions, highlighting their tunability.     

Experimental and DFT computational studies on 5-benzyl-4-(3,4-dimethoxyphenethyl)-2H-1,2,4-triazol-3(4H)-one

The triazole compound, 5-benzyl-4-(3,4-dimethoxyphenethyl)-2H-1,2,4-triazol-3(4H)-one, has been synthesized and characterized by ¹H-NMR, ¹³C-NMR, IR, and X-ray single-crystal determination. The compound crystallizes in the monoclinic space group P2₁ with a?=?11.8844(3) Å, b?=?17.5087(4) Å, c?=?17.3648(6) Å, β?=?99.990(2)? and Z?=?8. In addition to the molecular geometry from X-ray experiment, the molecular geometry, vibrational frequencies and gauge including atomic orbital (GIAO) ¹H- and ¹³C-NMR chemical shift values of the title compound in the ground state have been calculated using the density functional method (B3LYP) with 6-31G(d,p) basis set. The calculated results show that the optimized geometries can well reproduce the crystal structure and the theoretical vibrational frequencies and chemical shift values show good agreement with experimental ones. Besides, molecular electrostatic potential (MEP), natural bond orbital (NBO), and frontier molecular orbitals (FMO) analysis of the title compound were performed by the B3LYP/6-31G(d,p) method.     

Molecular polarization maps as a tool for studies of intermolecular interactions and chemical reactivity

Maps for the interaction energy of acetone, pyrrole, furan, and pyridine with a positive unitary charge were computed using ab initio techniques, together with their molecular electrostatic potentials at the same points. The difference between the interaction and electrostatic potential maps yielded polarization maps for the molecules. Finally, maps for the interaction with a negative charge were obtained as the difference between the polarization and electrostatic potential maps.The calculations were carried out for three planes, 2 Bohr radii, 4 Bohr radii, and 8 Bohr radii from the plane containing the heavy atoms for all the molecules. At larger distances, the interaction and electrostatic maps resemble each other qualitatively; however, at shorter distances, where the polarization effects are more significant, the differences between the maps are notable.Interaction and polarization maps can be routinely evaluated for medium-sized molecules, and are likely to become an important tool in drug design and chemical reactivity.     

New Multicentre Point Charge Models for Molecular Electrostatic Potentials from Semiempirical M0-Calculations

Two quasi-multipole electrostatic models for molecular charge distributions are presented. They assign arrays of point charges to nonhydrogen atoms on the basis of hybrid orbitals or localised molecular orbitals. When used with common semiempirical MO-techniques, they reproduce natural atomic orbital derived point charge (NAO-PC) and ab initio molecular potentials well. The localised orbital technique (LMO-PC) is intuitively more attractive than the hybrid orbital-point charge (HO-PC) method, although the former is more CPU-intensive.Electronic Supplementary Material available.