Diels–Alder reactivity of benzannulated isobenzofurans as assessed by density functional theory

Density functional theory (DFT) calculations at the B3LYP/6-31G* level for isobenzofuran 1 and eleven benzannulated derivatives of types 2 and 3 have been performed in order to compare their relative reactivities as dienes in Diels–Alder reactions. The transition state (TS) energies for their reactions with ethylene have been determined and shown to form a linear correlation between activation energies and structure count (SC) ratios. TS energies as a method for comparison of diene reactivities can be applied to IBFs bearing substituents on the ring as well as those containing heteroatoms, for which the SC ratio method failed. Different measures of aromaticity of benzannulated IBFs indicated a decrease in aromaticity going from 4 to 14, which is also reflected in their reactivity as a dienes in Diels–Alder reaction. Figure Isobenzofuran 1 and benzannulated isobenzofurans 2 and 3      

Analysis of the effect of electrostatic energy truncation in molecular dynamics simulations of immunoglobulin G light chain dimer

Molecular dynamics (MD) simulations of immunoglobulin G (IgG) light chain dimer using particle mesh Ewald (PME) and cutoff methods of treating electrostatic interactions were performed. The results indicate that structural parameters (RMSD, radius of gyration, solvent accessible surface) are very similar for both schemes; however, PME simulation shows increased mobility of side chains. This leads to larger fluctuations in the distance between the monomers in the dimer molecule, and, as a consequence, results in decreased number of interactions across the dimer interface. The wall clock time of the simulations was also compared. It was shown that the PME method is approximately 30% faster than the cutoff method for the system studied on a single processor.Figure Backbone order parameters for PME (red) and cutoff (green) calculations. Thick, horizontal lines show stable secondary structures     

Molecular modeling of B-DNA site recognition by Ru intercalators: molecular shape selection

In this work, molecular modeling methods have been applied to the interaction characterization of polypyridyl transitional-metal complexes with the oligonucleotide (B-DNA fragment). In order to explore the factors governing the groove recognition and intercalative depth, we establish a simple and practical docking method (step-by-step docking operation) to obtain potential curves while making complexes inset into B-DNA along an assigned path. Energy values in the potential curve are obtained from energy minimization of binding geometries. Modeling results clearly show that the optimum binding conformation corresponding to the global minimum in the potential curve for each complex is found to correlate well with the experimental results. Our results also confirm that minor changes of the ligand structure can lead to profound influences on binding geometries, so the molecular shape of the complexes is a predominant factor in governing the binding mode. Moreover, we find that the vdW force and water molecular effect are strongly associated with molecular-shape selection in our model. These results complement and extend the knowledge of the nature of these complexes binding to B-DNA.Figure Schematic illustration of metal complexes bound to B-DNA. The complexes are intercalated into the A₅T₆/T₆A₅ base step via a head-on fashion     

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.     

Characterization of molecular structures and properties of polyurethanes using molecular dynamics simulations

Thermotropic polyurethanes with mesogenic groups in side chains were prepared from two diisocyanates and four diols with stoichiometric ratios of reactive isocyanate (NCO) and hydroxy (OH) groups. Their thermal behavior was determined by differential scanning calorimetry. The effect of structure modifications of the diisocyanates and diols, in particular changes in the mesogen, were investigated. Introduction of mesogenic segments into the polymers suppresses the ordering. Stiff end substituents (phenyl and alkoxy groups) of the mesogens stabilize the mesophases to such an extent that the negative influence of long polymer chains is compensated and the liquid-crystalline properties are recovered. All-atom molecular dynamics simulations in the Cerius² modeling environment were carried out to characterize the structures of the polymers. Analysis of the dynamic trajectories at 20, 100, 120 and 170 °C revealed changes in conformation of macromolecules, which correlate with DSC measurements.Figure Example of structure relaxation of D4/TDI molecule at indicated simulation times (temperature 20 °C): a complete structure; b backbone structure; c top view of molecule     

Active site modeling in copper azurin molecular dynamics simulations

Active site modeling in molecular dynamics simulations is investigated for the reduced state of copper azurin. Five simulation runs (5 ns each) were performed at room temperature to study the consequences of a mixed electrostatic/constrained modeling for the coordination between the metal and the polypeptide chain, using for the ligand residues a set of charges that is modified with respect to the apo form of the protein by the presence of the copper ion.The results show that the different charge values do not lead to relevant effects on the geometry of the active site of the protein, as long as bond distance constraints are used for all the five ligand atoms. The distance constraint on the O atom of Gly45 can be removed without altering the active site geometry. The coordination between Cu and the other axial ligand Met121 is outlined as being flexible. Differences are found between the bonds of the copper ion with the two apparently equivalent N¹ atoms of His46 and His117.The overall findings are discussed in connection with the issue of determining a model for the active site of azurin suitable to be used in molecular dynamics simulations under unfolding conditions. Figure Model of azurin active site. Copper ligand residues are cut off at C position except Gly45, for which the portion of backbone connecting it to His46 is shown. Only polar H atoms are shown. All atoms are in standard colors (Cu in violet), and the five ligands are labeled     

Local molecular properties and their use in predicting reactivity

Expressions for the local electron affinity, electronegativity and hardness are derived in analogy to the local ionization energy introduced by Sjoberg, Murray and Politzer. The local polarizability is also defined based on an additive atomic orbital polarizability model that uses Rivail's variational technique. The characteristics of these local properties at molecular surfaces and their relevance to electrophilic aromatic substitution, to S_N2 reactivity and to the nucleophilicity of enolate ions are discussed.Figure The local ionization energy at the SES surfaces of methyl benzoate. The color scale ranges from 375 (blue) to 550 kcal mol^–1 (red). The blue areas are those for which interaction with an acceptor is most favorable.     

Reinvestigation of molecular structure and barrier to internal rotation of pyridinium N-phenolate betaine dye

In the present paper, the results of a systematic theoretical study of the molecular structure of 4-(1-pyridinium-1-yl)phenolate betaine are reported. The ground-state molecular structure and the barrier to internal rotation of the betaine dye molecule were calculated ab inito (with Hartree–Fock theory and the second-order of Möller–Plesset method) and with density functional theory (DFT). In order to estimate the complete basis set limit, the calculations of barriers to internal rotations were performed using correlation–consistent basis sets with a maximal cardinal number of four. It was determined that electron correlation is crucial in order to obtain reliable geometries and rotational barriers of the molecule investigated. For the sake of comparison, the results of calculations using the AM1 Hamiltonian are also presented.Figure Investigated betaine dye.     

Binding affinity of hydroxamate inhibitors of matrix metalloproteinase-2

Here we report molecular dynamics (MD) and free energy perturbation (FEP) simulations applied to hydroxamate-matrix metalloproteinase-2 (MMP-2) complex systems. We have developed some new force field parameters for the hydroxamate functional group that were not included in the AMBER94 force field but were necessary in our simulations. For the representation of the active zinc center, a bonded model was adopted in which restrained electrostatic potential fitting (RESP) charges were used as the electrostatic representation of this model. Using the resulted bonded model, FEP simulations predict the relative binding free energy in good agreement with the experimental value. By analyzing the molecular dynamics (MD) trajectories of the two complex systems, we can provide an explanation of why one of the two inhibitors is favored over the other. The results provide a chemical insight into the interactions between inhibitor and enzyme, and can indicate changes in the inhibitor that would enhance inhibitor–enzyme interactions.Figure The scheme of the binding site     

Molecular dynamics exposes α-helices in myelin basic protein

Molecular dynamics simulations of models of unmodified and deiminated MBP (myelin basic protein) have been performed on solvated structures with added counterions, for 10 ns using AMBER (assisted model building with energy refinement). The protein structures became extended, and a considerable number of -helical segments formed spontaneously. The degree of molecular extension was greater in the deiminated species, and the -helices were more transient. These structural disruptions may be operative in vivo during multiple sclerosis.Figure A model of the C1 isoform of myelin basic protein showing major -helical segments that were stable over the last 1 ns of a molecular dynamics simulation. The -helices are colored red, the -strands are colored yellow, the -turns are colored blue, and random coils are colored green     

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     

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.     

A computational approach to the synthesis of dirithromycin

Dirithromycin is a macrolide antibiotic derived from erythromycin A. Dirithromycin is synthesized by the condensation of 9(S)-erythromycylamine with 2-(2-methoxyethoxy)-acetaldehyde. To gain insight into the synthesis, the condensation mechanism has been analyzed computationally by the AM1 method in the gas phase. First, the formation of the Schiff bases of dirithromycin and epidirithromycin from 9(S)-erythromycylamine and 2-(2-methoxyethoxy)-acetaldehyde were modeled. Then, the tautomerization of the Schiff bases to dirithromycin and epidirithromycin were considered. Finally, the epimerization of the Schiff base of epidirithromycin to the Schiff base of dirithromycin was investigated. Our results show that, even though carbinolamine forms faster for epidirithromycin than the corresponding structure for dirithromycin, dirithromycin is the major product of the synthesis. Figure Synthesis of dirithromycin     

3D QSAR studies on peroxisome proliferator-activated receptor γ agonists using CoMFA and CoMSIA

The peroxisome proliferator-activated receptors (PPARs) have increasingly become attractive targets for developing novel anti-type 2 diabetic drugs. We employed comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) to study three-dimensional quantitative structure–activity relationship (3D QSAR) based on existing agonists of PPAR (including five thiazolidinediones and 74 tyrosine-based compounds). Predictive 3D QSAR models with conventional r² and cross-validated coefficient (q²) values up to 0.974 and 0.642 for CoMFA and 0.979 and 0.686 for COMSIA were established using the SYBYL package. These models were validated by a test set containing 18 compounds. The CoMFA and CoMSIA field distributions are in general agreement with the structural characteristics of the binding pockets of PPAR, which demonstrates that the 3D QSAR models built here are very useful in predicting activities of novel compounds for activating PPAR.      

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.     

Modeling of the phase equilibria of polystyrene in methylcyclohexane with semi-empirical quantum mechanical methods I

A method for calculating interaction parameters traditionally used in phase-equilibrium computations in low-molecular systems has been extended for the prediction of solvent activities of aromatic polymer solutions (polystyrene+methylcyclohexane). Using ethylbenzene as a model compound for the repeating unit of the polymer, the intermolecular interaction energies between the solvent molecule and the polymer were simulated. The semiempirical quantum chemical method AM1, and a method for sampling relevant internal orientations for a pair of molecules developed previously were used. Interaction energies are determined for three molecular pairs, the solvent and the model molecule, two solvent molecules and two model molecules, and used to calculated UNIQUAC interaction parameters, a _ij and a _ji . Using these parameters, the solvent activities of the polystyrene 90,000 amu+methylcyclohexane system, and the total vapor pressures of the methylcyclohexane+ethylbenzene system were calculated. The latter system was compared to experimental data, giving qualitative agreement. Figure Solvent activities for the methylcylcohexane(1)+polystyrene(2) system at 316 K. Parameters a _ij (blue line) obtained with the AM1 method; parameters a _ij (pink line) from VLE data for the ethylbenzene+methylcyclohexane system. The abscissa is the polymer weight fraction defined as ₂(x ₁)=(1–x ₁)M ₂/[x ₁ M ₁+(1–x ₁)M ₂], where x ₁ is the solvent mole fraction and M _i are the molecular weights of the components.An erratum to this article can be found at     

Properties of star-branched and linear chains in confined space. A Monte-Carlo study

We have studied the properties of simple models of linear and star-branched polymer chains confined in a slit formed by two parallel impenetrable walls. The polymer chains consisted of identical united atoms (homopolymers) and were restricted to a simple cubic lattice. Two macromolecular architectures of the chain: linear and regular stars with three branches of equal length, were studied. The excluded volume was the only potential introduced into the model and thus the system was athermal. Monte-Carlo simulations with the sampling algorithm based on the chains local changes of conformation were carried out for chains with different lengths as well as for different distances between the confining surfaces. We found that the properties of model chains differ for both macromolecular architectures but a universal behavior for both kinds of chains was also found. Investigation of the frequency of chain-wall contacts shows that the ends of the chains are much more mobile than the rest of the chain, especially in the vicinity of the branching point in star polymers.Figure The scheme of a star-branched (left) and a linear (right) chain located between two parallel impenetrable surfaces.     

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     

Predicting anti-HIV activity of 2,3-diaryl-1,3-thiazolidin-4-ones: computational approach using reformed eccentric connectivity index

The eccentric connectivity index, which has recently been employed successfully for the development of numerous mathematical models for the prediction of biological activities of diverse nature, has been reformed to overcome its limitations caused by degeneracy and insensitivity towards heteroatoms. The reformed eccentric connectivity index, termed the eccentric connectivity topochemical index, overcomes the limitations of the eccentric connectivity index by exhibiting very low degeneracy and displaying sensitivity to both the presence and relative position of heteroatoms without compromizing the discriminating power of the eccentric connectivity index. The relationship of the eccentric connectivity topochemical index, eccentric connectivity index and Wieners index with regard to the anti-HIV activity of 2, 3-diaryl-1, 3-thiazolidin-4-one derivatives was subsequently investigated. The values of the eccentric connectivity topochemical index, the eccentric connectivity index and Wieners index of each of 31 analogues comprizing the data set were computed using in-house computer program. Resultant data was analyzed and suitable models developed after identification of active ranges. Subsequently, each derivative was assigned a biological activity using these models, which was then compared with the reported anti-HIV activity. The accuracy of prediction using these models was found to vary from 81 to 90%. The proposed index offers a vast potential for virtual screening of combinatorial libraries, structure property/activity studies and drug design.Figure Basic structure of 2,3-diaryl-1, 3-thiazoidin-4-ones.     

Density functional and docking studies of retinoids for cancer treatment

The retinoic acid receptor (RAR) and retinoid X receptor (RXR) are members of the nuclear receptor superfamily. The ligand-binding domain contains the ligand-dependent activation function. The isotypes RAR, and are distinct pharmacological targets for retinoids involved in the treatment of various cancers and skin diseases. There is thus considerable interest in synthetic retinoids with isotype selectivity and reduced side effects. In this work we have focused on the retinoid acid receptor and three of its panagonists. We have carried out density functional geometry optimizations at the B3LYP/6-31G* level, computed two types of atomic charges and also electrostatic potentials. A docking program was used to investigate the interactions between the receptor and the three ligands. A theoretically more potent inhibitor, which was obtained by modifying one of the retinoic acids investigated, is proposed. Figure Superposition of the crystal structure (Å) of the 1FCX ligand with the proposed new inhibitor