A new visualization scheme of chemical energy density and bonds in molecules

Covalent bond describes electron pairing in between a pair of atoms and molecules. The space is partitioned in mutually disjoint regions by using a new concept of the electronic drop region R_D, atmosphere region R_A, and the interface S (Tachibana in J Chem Phys 115:3497–3518, 2001). The covalent bond formation is then characterized by a new concept of the spindle structure. The spindle structure is a geometrical object of a region where principal electronic stress is positive along a line of principal axis of the electronic stress that connects a pair of the R_Ds of atoms and molecules. A new energy density partitioning scheme is obtained using the Rigged quantum electrodynamics (QED). The spindle structure of the stress tensor of chemical bond has been disclosed in the course of the covalent bond formation. The chemical energy density visualization scheme is applied to demonstrate the spindle structures of chemical bonds in H₂, C₂H₆, C₂H₄ and C₂H₂ systems.Figure Field theory of the energy density.     

Dynamic molecules: molecular dynamics for everyone. An internet-based access to molecular dynamic simulations: basic concepts

Molecular dynamics is a rapidly developing field of science and has become an established tool for studying the dynamic behavior of biomolecules. Although several high quality programs for performing molecular dynamic simulations are freely available, only well-trained scientists are currently able to make use of the broad scientific potential that molecular dynamic simulations offer to gain insight into structural questions at an atomic level. The "Dynamic Molecules" approach is the first internet portal that provides an interactive access to set up, perform and analyze molecular dynamic simulations. It is completely based on standard web technologies and uses only publicly available software. The aim is to open molecular dynamics techniques to a broader range of users including undergraduate students, teachers and scientists outside the bioinformatics field. The time-limiting factors are the availability of free capacity on the computing server to run the simulations and the time required to transport the history file through the internet for the animation mode. The interactive access mode of the portal is acceptable for animations of molecules having up to about 500 atoms.Figure Several main menus (see top) are provided to start "New Simulations", to "Display Simulations" and to "Analyze" statistical and geometrical properties of the molecule. Here the "Display Simulation" interface is shown. The Chime plugin is used to visualize molecular 3D structures and motions.     

Effect of Pt and Sn on the adsorption of n-heptane in γ-Al2O3 catalyst models

The Grand Canonical Monte-Carlo (GCMC) method has been used to carry out simulations of the adsorption of n-heptane in models of naphtha-reforming catalysts. Models used in the study differed in the number and distribution of metal atoms—Pt and Sn. The number of adsorbed n-heptane molecules grows linearly with increasing number of metal atoms. The effect of Pt content on the adsorption of n-heptane molecules is most distinct at approximately 100 kPa and within the lower range of the temperatures investigated. In the models of bimetallic catalysts, the effect of the two metals is additive.Figure Effect of Pt and Sn on number of n-heptane molecules adsorbed in Al₂O₃ catalyst in 773 K and 1000 kPa.      

The HF-SCF energy of HIV-1 MNgp120 V3 hairpin loop conformers

The purpose of this study is to analyze the structure of the V3 loop of the HIV-1 gp120 molecule at the atomic level. The total energy of each member of the antibody-complexed 16-mer V3 conformer data set of Sharon et al. (PDB 1NJ0) was determined by the Hartree–Fock-self-consistent field (HF-SCF) method and with the GROMOS96 force field. There was no correlation between the results of the classical GROMOS96 force field analysis and the ab initio HF-SCF quantum mechanical analysis of the energy of the V3-loop-peptide conformers. HF-SCF optimization (AM1) of conformer geometries yielded structures in which HIS315 is displaced from its original position in the combining site of human antibody fragment 447-52D, but with the hairpin turn intact. The hairpin shape of the V3 loop remained detectable, albeit distorted, even with perturbation by a lithium dicationic electrostatic force field and by substitution of the PRO320 at the crown of the V3 hairpin by a GLY. These data suggest that the hairpin conformation is at least partially stable to long-range electrostatic perturbations, either with or without PRO in the tip of the crown of the V3-hairpin loop. Figure Molecular geometry of HIV-1 V3 conformer model 5 and a GLY320 substituted model 5. Space-filling models were obtained with ViewMol3D [Sharon et al. (2002) PDB 1NJ0]). Red=oxygen, blue=nitrogen, black=carbon, white=hydrogen and purple=lithium. End-to-end distance (D) was obtained with ViewMol3D and is in Ångströms. Geometry optimized GLY320 Model 5, D=4.74 ÅThis revised version was published online in October 2004 with corrections to the Graphical Abstract.     

Molecular dynamics simulations of 14 HIV protease mutants in complexes with indinavir

The molecular mechanisms of HIV drug resistance were studied using molecular dynamics simulations of HIV-1 protease complexes with the clinical inhibitor indinavir. One nanosecond molecular dynamics simulations were run for solvated complexes of indinavir with wild type protease, a control variant and 12 drug resistant mutants. The quality of the simulations was assessed by comparison with crystallographic and inhibition data. Molecular mechanisms that contribute to drug resistance include structural stability and affinity for inhibitor. The mutants showed a range of structural variation from 70 to 140% of the wild type protease. The protease affinity for indinavir was estimated by calculating the averaged molecular mechanics interaction energy. A correlation coefficient of 0.96 was obtained with observed inhibition constants for wild type and four mutants. Based on this good agreement, the trends in binding were predicted for the other mutants and discussed in relation to the clinical data for indinavir resistance. Figure Poincare map representation for WT protease-indinavir complex. The side chain of Tyr 59 showing the positions of hydrogen atoms.This revised version was published online in October 2004 with corrections to the Graphical Abstract.     

Prototropic tautomerism of imidazolone in aqueous solution: a density functional approach using the combined discrete/self-consistent reaction field (SCRF) models

A systematic investigation of the proton transfer in the keto-amino/enol tautomerization of imidazolone was undertaken. Calculations in aqueous solution were performed using both combined discrete/self-consistent reaction field (SCRF) and SCRF methods. Complexes containing one to three water molecules around the hydrophilic site of imidazolone were used for the combined discrete/SCRF calculations. The DFT results predict that the barrier height for non-water-assisted intramolecular proton transfer is very high (214.8 kJmol^–1). Hydrogen bonding between imidazolone and the water molecule(s) will dramatically lower the barrier by a concerted multiple proton transfer mechanism. The proton transfer process through a eight-member ring formed by imidazolone and two water molecules is found to be more efficient and the calculated barrier height is ca. 61 kJmol^–1.Figure DFT calculations in aqueous solution predict the H-bonding between imidazolone(IZ) and the water molecule(s) will dramatically lower the tautomeric barrier by a concerted multiple proton transfer mechanism, in which an eight-member ring structure formed by IZ and 2H₂O is found to be more efficient and the barrier is 60.8 kJ mol^–1, much less than 214.8 kJ mol^–1 in the non-water-assisted mechanism.     

Theoretical study of the electron affinities of the alkaline-earth tetramers possessing T d symmetry: Be4 and Mg4

The electron affinities of beryllium and magnesium tetramers are calculated at the ROMP2 level employing the Dunning-type aug-cc-pVQZ basis set. The vertical electron detachment energy (VEDE) amounts to 1.685 eV for Be₄^– and 0.943 eV for Mg₄ . The decomposition of the VEDE into physical components and an atomic orbital population analysis are used to elucidate the nature of the outer electron binding in these anions.Figure The lowest unoccupied molecular orbitals in the ground state of Mg4 : a LUMO, symmetry A₁, b LUMO + 1, symmetry T₂; c the highest occupied molecular orbital (HOMO), symmetry A1 in the ground state of Mg₄^–.      

Molecular modeling of dipeptide and its analogous systems with water

Intermolecular hydrogen-bond interactions in the monohydrated complexes of formamide, N-methylacetamide and glycylglycine have been studied using ab initio and DFT methods. The geometries were optimized using second-order Møller–Plesset perturbation theory and the B3LYP DFT functional with the 6-311++G** basis set. It is observed that hydrogen-bond interactions at the carbonyl group of the peptide moiety are stronger than those at the amino group of the formamide and N-methylacetamide molecules. Because of the presence of cyclic hydrogen-bonding interactions in glycylglycine, the interaction at the amino group is higher than at the carbonyl. The ¹³C and ¹⁵N NMR shielding values were calculated for the non-hydrated and monohydrated complexes. Condensed Fukui functions have also been calculated for non-hydrated formamide, N-methylacetamide and glycylglycine molecules at the B3LYP/6-311++G** level of theory, and the results are discussed.Figure Structure of hydrated glycylglycine dipeptide     

An unusual feature of end-substituted model carbon (6,0) nanotubes

We have examined the effects of substituents on the computed electrostatic potentials V_S(r) and average local ionization energies on the surfaces of model carbon nanotubes of the types (5,5), (6,1) and (6,0). For the (5,5) and the (6,1), the effects upon both V_S(r) and of substituting a hydroxyl group at one end are primarily localized to that part of the system. For the (6,0) tube, however, a remarkable change is observed over its entire length, with V_S(r) showing a marked gradation from strongly positive at the substituted end to strongly negative at the other; correspondingly goes from higher to lower values. Replacing OH by another resonance- donor, NH₂, produces similar results in the (6,0) system, while the resonance withdrawing NO₂ does the opposite, but in equally striking fashion. We explain these observations by noting that the arrangement of the C–C bonds in the (6,0) tube facilitates charge delocalization over the full length and entire surface of the tube. Substituting NH₂ and NO₂ at opposite ends of the (6,0) tube greatly strengthens the gradations in both V_S(r) and The first hyperpolarizability of this system was found to be nine times that of para-nitroaniline, suggesting possible nonlinear optical applications.Figure HF/STO-5G electrostatic potential on outer surface of open (6,0) C₇₂H₁₀NH₂NO₂. The nitro group is at the right end of the tube, the amino group at the left. In eV: purple is less than 14, blue is between 14 and 15, green is between 15 and 16.5, yellow is between 16.5 and 17.5, and red is more than 17.5.      

Organization of rhodopsin molecules in native membranes of rod cells–an old theoretical model compared to new experimental data

It has been shown that rhodopsin forms an oligomer in the shape of long double rows of monomers. Because of the importance of rhodopsin as a template for all G protein-coupled receptors, its dimeric, tetrameric and higher-oligomeric structures also provide a useful pattern for similar structures in GPCRs. New experimental data published recently are discussed in the context of a proposed model of the rhodopsin oligomer 1N3M deposited in the protein data bank. The new rhodopsin structure at 2.2 Å resolution with all residues resolved as well as an electron cryomicroscopy structure from 2D crystals of rhodopsin are in agreement with the 1N3M model. Accommodation of movement of transmembrane helix VI, regarded as a major event during the activation of rhodopsin, in a steady structure of the oligomer is also discussed.Figure Superimposition of the 1U19 (red wire), 1GZM (purple wire) and 1N3M (blue wire) rhodopsin structures. Size of the wires is proportional to thermal factors of backbone C atoms, view parallel to the membrane.      

A multimeric model for murine anti-apoptotic protein Bcl-2 and structural insights for its regulation by post-translational modification.

A monomeric model for murine antiapoptotic protein Bcl-2 was constructed by comparative modeling with the software suite MPACK (EXDIS/DIAMOD/FANTOM) using human Bcl-xL as a template. The monomeric model shows that murine Bcl-2 is an all -helical protein with a central (helix 5) hydrophobic helix surrounded by amphipathic helices and an unstructured loop of 30 residues connecting helices 1 and 2. It has been previously shown that phosphorylation of Ser 70 located in this loop region regulates the anti-apoptotic activity of Bcl-2. Based on our current model, we propose that this phosphorylation may result in a conformational change that aids multimer formation. We constructed a model for the Bcl-2 homodimer based on the experimentally determined 3D structure of the Bcl-xL: Bad peptide complex. The model shows that it will require approximately a half turn in helix 2 to expose hydrophobic residues important for the formation of a multimer. Helices 5 and 6 of the monomeric subunit Bcl-2 have been proposed to form an ion-channel by associating with helices 5 and 6 of another monomeric subunit in the higher-order complex. In the multimeric model of Bcl-2, helices 5 and 6 of each subunit were placed distantly apart. From our model, we conclude that a global conformational change may be required to bring helices 5 and 6 together during ion-channel formation.Figure Hydrophobic interactions in the dimerization groove are shown. Helix 2' of monomer B interacts through residues V90, H91, L94, A97, G98, F101 and Y105 with the hydrophobic surface formed by residues in helices 3, 4, and 5 of the monomer A. Shown here is a lateral view of monomer A depicted in a surface model with hydrophobic regions in red. The backbone of the helix is shown using a neon representation in yellow and the interacting side chains are in blue.     

Molecular dynamics simulation of hepatitis C virus IRES IIId domain: structural behavior,electrostatic and energetic analysis

The dynamic behavior of the HCV IRES IIId domain is analyzed by means of a 2.6-ns molecular dynamics simulation, starting from an NMR structure. The simulation is carried out in explicit water with Na⁺ counterions, and particle-mesh Ewald summation is used for the electrostatic interactions. In this work, we analyze selected patterns of the helix that are crucial for IRES activity and that could be considered as targets for the intervention of inhibitors, such as the hexanucleotide terminal loop (more particularly its three consecutive guanines) and the loop-E motif. The simulation has allowed us to analyze the dynamics of the loop substructure and has revealed a behavior among the guanine bases that might explain the different role of the third guanine of the GGG triplet upon molecular recognition. The accessibility of the loop-E motif and the loop major and minor groove is also examined, as well as the effect of Na⁺ or Mg²⁺ counterion within the simulation. The electrostatic analysis reveals several ion pockets, not discussed in the experimental structure. The positions of these ions are useful for locating specific electrostatic recognition sites for potential inhibitor binding. Figure Superposition of 14 structures representative of the evolution of IRES IIId RNA along 2.6-ns MD simulation     

A study on the influence of molecular properties in the psychoactivity of cannabinoid compounds

Several molecular properties are calculated for a set of 26 cannabinoid compounds with the goal of connecting the psychoactivity of the compounds with an appropriate set of calculated properties. For this purpose we used quantum chemical (the AM1 semi-empirical method) and chemometric methods. The AM1 method was employed to calculate the set of quantum chemical molecular properties and the chemometric methods were employed with the aim of selecting the most relevant properties to be correlated with psychoactivity. The chemometric methods used were Principal Component Analysis (PCA), Hierarchical Cluster Analysis (HCA) and the K-Nearest Neighbor (KNN) method. The chemometric analysis showed that an electronic property (energy of LUMO), a hydrophobic property (log P), a steric property (volume of the substituent at the C4 position) and a topological property (Lovasz–Pelikan index) were the most important variables for the separation between the psychoactive and psychoinactive compounds. In order to validate our PCA, HCA and KNN results, eight new cannabinoid compounds (with known psychoactivity) were used in a prediction study and were classified correctly by the methods used in this work, indicating that our PCA, HCA and KNN models are able to predict reliable psychoactivity of cannabinoid compounds. Figure: ⁹-THC This revised version was published online in June 2005 with corrections to Table 1.     

Chiral recognition of aromatic compounds by β-cyclodextrin based on bimodal complexation

The chiral recognition of the selected aromatic chiral compounds by native -cyclodextrin (-CD) based on bimodal complexation was studied using a flexible docking algorithm FDOCK. A quantitative empirical free energy relationship model was employed to predict the complex stability constants and the preferred binding modes. The results showed that the calculated complex stability constants are in good agreement with the experimental data. Furthermore, the main force responsible for host-guest complexation is the van der Waals force and the chiral molecules are completely included into the -CD cavity. The chiral recognition for the selected aromatic chiral compounds is the result of the van der Waals force counterbalancing with the other effects, such as the electrostatic interaction and the hydrophobic effect.Figure The favorable structures for the inclusion complexes of (S)_phenylbutyric with -CD. View in the plane normal to the Z-axis. -CD is shown in surface and (S)_phenylbutyric in CPK representation.     

On the nature of bonding in HCOOH...Ar and HCOOH...Kr complexes

The chemical interaction in HCOOH...Ng (Ng=Ar, Kr) complex was analyzed by topological analysis of the electron density based on Atoms-In-Molecules theory. For all computationally stable equilibrium structures of 1:1 HCOOH...Ng complexes, an intermolecular bond path with a bond critical point was found and perturbation of formic acid (FA) atomic basins and electron density was observed. The intermolecular interaction between the two complex subunits can be classified, according to its nature, as a closed-shell van der Waals type of interaction. However, one of the computed structures (complex II), pictures a noble gas atom attached linearly to the acidic O–H tail of FA. In this particular case, the electron density at the intermolecular bond critical point was found to resemble a hydrogen-bonded system and thus, may be termed a hydrogen-bond-like interaction. This change in the nature of the interaction is also shown by large perturbations of the FA properties found for this complex structure. The structural and vibrational perturbations are larger than for the other three structures and they increase for the Kr complexes compared to the Ar complex.Figure Electron density analysis of HCOOH...Ng (Ng=Ar,Kr) complex.     

Determination of fuzzy logic membership functions using genetic algorithms: application to structure–odor modeling

Fuzzy logic has been used as a tool in structure–camphoraceous odor relationships. The data base studied included 99 molecules. The rules used to discriminate between camphor and non camphor molecules lead to 77% correct discrimination. Such rules account for the shape and the size of the molecule. Their adjustment by means of genetic algorithms led to 84% correct discrimination between camphor and non-camphor molecules.Figure Membership function for the chosen variables     

Analysis of the interactions of ribonuclease inhibitor with kanamycin

The interaction of ribonuclease inhibitor (RI) with kanamycin was studied by molecular modeling. The preliminary binding model was constructed using the Affinity module of the Insight II molecular modeling program and the key residues involved in the combination of RI binding to kanamycin were determined. Meanwhile, we determined relevant surface characteristics determining the interaction behavior. The modeling results indicated that electrostatic interactions and H-bond forces may work as major factors for the molecular interaction between kanamycin and RI. The above results are useful for elucidating the molecular principles upon which the selectivity of a kanamycin is based. The quartz-crystal microbalance (QCM) is a new method usually used to monitor the binding function of macromolecules with samples online in a flow-injection analysis (FIA) system. The experimental results demonstrate that kanamycin has an extraordinary affinity to the basic protein RI, and our result is consistent with the molecular modeling results. These principles can in turn be used to study the molecular recognition mechanism and design a mimic of kanamycin for the development of new RI binders.Figure Proposed binding model of kanamycin to RI obtained by the computer-aided docking method and optimized with molecular mechanics with the CVFF forcefield.     

Molecular aspects of the interaction between amphotericin B and a phospholipid bilayer: molecular dynamics studies

Amphotericin B (AmB) is a polyene macrolide antibiotic used to treat systemic fungal infections. The molecular mechanism of AmB action is still only partly characterized. AmB interacts with cell-membrane components and forms membrane channels that eventually lead to cell death. The interaction between AmB and the membrane surface can be regarded as the first (presumably crucial) step on the way to channel formation. In this study molecular dynamics simulations were performed for an AmB–lipid bilayer model in order to characterize the molecular aspects of AmB–membrane interactions. The system studied contained a box of 200 dimyristoylphosphatidylcholine (DMPC) molecules, a single AmB molecule placed on the surface of the lipid bilayer and 8,065 water molecules. Two molecular dynamics simulations (NVT ensemble), each lasting 1 ns, were performed for the model studied. Two different programs, CHARMM and NAMD2, were used in order to test simulation conditions. The analysis of MD trajectories brought interesting information concerning interactions between polar groups of AmB and both DMPC and water molecules. Our studies show that AmB preferentially took a vertical position, perpendicular to the membrane surface, with no propensity to enter the membrane. Our finding may suggest that a single AmB molecule entering the membrane is very unlikely.Figure The figure presents the whole structure of the system simulated—starting point. AmB is presented as a space-filling model, DMPC molecules—green sticks, water molecules—red sticks     

Quantitative structure-activity relationship by CoMFA for cyclic urea and nonpeptide-cyclic cyanoguanidine derivatives on wild type and mutant HIV-1 protease

Abstact 3D-QSAR studies using the Comparative Molecular Field Analysis (CoMFA) methodology were conducted to predict the inhibition constants, K_i, and the inhibitor concentrations, IC₉₀ of 127 symmetrical and unsymmetrical cyclic urea and cyclic cyanoguanidine derivatives containing different substituent groups such as: benzyl, isopropyl, 4-hydroxybenzyl, ketone, oxime, pyrazole, imidazole, triazole and having anti-HIV-1 protease activities. A significant cross-validated correlation coefficient (q²) of 0.63 and a fitted correlation coefficient r² of 0.70 were obtained, indicating that the models can predict the anti-protease activity from poorly to highly active compounds reliably. The best predictions were obtained for: XV643 (predicted log 1/K_i=9.86), a 3,5-dimethoxy-benzyl cyclic urea derivate (molec60, predicted log 1/K_i=8.57) and a benzyl cyclic urea derivate (molec 61, predicted log 1/IC₉₀=6.87). Using the CoMFA method, we also predicted the biological activity of 14 cyclic urea derivatives that inhibit the HIV-1 protease mutants V82A, V82I and V82F. The predicted biological activities of the: (i) XNO63 (inhibitory activity on the mutant HIV-1 PR V82A), (ii) SB570 (inhibiting the mutant HIV-1 PR V82I) and also (iii) XV652 (during the interaction with the mutant HIV-1 PR V82F) were in good agreement with the experimental values.Figure Stereoview of the contour plots of the CoMFA steric and electrostatic fields. The favorable (indicated by blue polyhedra) and unfavorable (represented by red polyhedra) electrostatic areas and also the favorable (shown by green polyhedra) and unfavorable (shown by yellow polyhedra) steric areas formed around the most active molecule, 6a.     

Structure of glutathione S-transferase of the filarial parasite Wuchereria bancrofti: a target for drug development against adult worm

A three dimensional structural model of Glutathione-S-transferase (GST) of the lymphatic filarial parasite Wuchereria bancrofti (wb) was constructed by homology modeling. The three dimensional X-ray crystal structure of porcine -class GST with PDB ID: 2gsr-A chain protein with 42% sequential and functional homology was used as the template. The model of wbGST built by MODELLER6v2 was analyzed by the PROCHECK programs. Ramachandran plot analysis showed that 93.5% of the residues are in the core region followed by 5.4 and 1.1% residues in the allowed and generously allowed regions, respectively. None of the non-glycine residues is in disallowed regions. The PROSA II z-score and the energy graph for the final model further confirmed the quality of the modeled structure. The computationally modeled three-dimensional (3D) structure of wbGST has been submitted to the Protein Data Bank (PDB) (PDB ID: 1SFM and RCSB ID: RCSB021668). 1SFM was used for docking with GST inhibitors by Hex4.2 macromolecular docking using spherical polar Fourier correlations.Figure: A three-dimensional (3D) structure of Glutathione-S-transferase (GST) of the lymphatic filarial parasite Wuchereria bancrofti (wb) was constructed by homology modeling. This modeled 3D structure of wbGST has been submitted to the Protein Data Bank (PDB) (PDB ID: 1SFM and RCSB ID: RCSB021668).