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     

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     

AIPAR: ab initio parametrization of intermolecular potentials for computer simulations

An unambiguous, fully ab initio and automated technique denoted AIPAR (ab initio parametrization) implemented in the SJBR program has been proposed to yield intermolecular interaction potentials between polar molecules and water. The AIPAR procedure has been applied to several organic molecules covering a wide range of structure and functional groups, namely methanol, acetone (propanone), methanethiol (methyl mercaptan), imidazole (1,3-diazole), oxazole and furan. The AIPAR-derived sets of parameters compare well with the empirical OPLS ones, mainly when the all-atoms model is employed in the OPLS procedure. Monte Carlo simulations were performed for an aqueous solution of methanol and for an equimolar binary mixture methanol–water using the AIPAR and OPLS parameters. The thermodynamic and geometric results obtained with the parameters obtained with the AIPAR procedure compare favorably with the OPLS simulations, even for the binary mixture, demonstrating the precision, robustness and transferability of the parameters obtained with the AIPAR procedure.Figure Superimposed configurations of water (without the hydrogen atoms) around the methanol molecule obtained with the AGOA procedure.     

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     

Quantitative structure–activity relationship study on some benzodiazepine derivatives as anti-Alzheimer agents

A QSAR study was performed in an attempt to explore the pharmacophore of some benzodiazepine derivatives as anti-Alzheimer agents for the inhibition of -secretase. The study, which used the electrotopological state atom (ETSA) index, which encodes electronic and topological information, reveals the importance of two phenyl rings—one substituted and another unsubstituted, for the inhibition of the enzyme. Fluorine substitution on the substituted phenyl ring has an important contribution to the activity. R-configurations of the aliphatic chain substituents provide the exact conformation of the molecules to enter into the binding pockets of the receptor(s).Figure General structure of benzodiazepine containing -secretase inhibitors     

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.     

Modeling of adsorption in nanopores

Adsorption in nonporous materials has been studied using Grand Canonical Monte Carlo simulations. We discuss three types of materials: (a) a model of cylindrical pores with smooth walls, representing MCM-41 like materials, (b) a model of cylindrical pores with regular structured walls (model of carbon nanotubes) and (c) a material with crystalline wall structure (zeolites). Typical problems related to the stability of adsorbed layers have been analyzed. We have shown that the mechanism of adsorption is strongly dependent on the structure of the pore walls. In the case of amorphous walls it may lead to metastable configurations. In nanotubes, the ordered corrugation structure of walls determines the low temperature structure of the adsorbed system. In 3D ordered porous system, such as zeolites, the mechanism of adsorption is mostly determined by characteristic sites of adsorption.Figure Adsorbed atoms and energy fluctuations at the pressure of the first layer formation of krypton atoms: (a) instantaneous numbers of adsorbed atoms (per nm² of the pore wall) as a function of the time of simulation (Monte Carlo steps) observed in a relatively long run, (b) the bimodal distribution of the energy fluctuations is a consequence of the behavior of the systems as shown in (a).     

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     

Theoretical studies of the conformational behavior of chain molecules containing polar groups: simulations of a poly(vinylidene fluoride) model

Atomistic Monte Carlo simulations have been conducted to elucidate the conformational behavior of a single chain molecule containing polar functional groups. Here, we resort to an atomistic poly(vinylidene fluoride) (PVDF) chain model as a representative example. The model is modified in such a way that bond lengths and bond angles are fixed, aiming to manifest the role of dipolar interactions. For a given chain length, chain conformation is sensitive to two environmental parameters, temperature and dielectric constant. The mean chain size increases when temperature and/or dielectric constant are increased. The conformational behavior is further characterized by chain size distribution function, and our findings show that temperature induced conformational transition for a chain molecule can be discrete or continuous, depending on its chain length. Also, the dipolar interactions in PVDF are effectively attractive, and enhance chain contraction. As a result, when the strength of dipolar interactions is increased, the discrete conformational transition shifts toward longer chains; and for a given chain length, such a transition occurs at higher temperatures.Figure Variation of R² with temperature for different dielectric constants =1 and 8, denoted by dotted and solid lines, respectively, and for different chain lengths M=8 and 12, as marked. Lines are meant for eye guidance     

Reasoning of spike glycoproteins being more vulnerable to mutations among 158 coronavirus proteins from different species

In this study, we used the probabilistic models developed by us over the last several years to analyze 158 proteins from coronaviruses in order to determine which protein is more vulnerable to mutations. The results provide three lines of evidence suggesting that the spike glycoprotein is different from the other coronavirus proteins: (1) the spike glycoprotein is more sensitive to mutations, this is the current state of the spike glycoprotein, (2) the spike glycoprotein has undergone more mutations in the past, this is the history of spike glycoprotein, and (3) the spike glycoprotein has a bigger potential towards future mutations, this is the future of spike glycoprotein. Furthermore, this study gives a clue on the species susceptibility regarding different proteins.Figure Predictable and unpredictable portions in coronavirus proteins. The data are presented as median with interquartile range. * the predictable and unpredictable portions in spike glycoprotein group are statistically different from any other protein groups at p<0.05 level, except for hemagglutinin-esterase precursor group. # the predictable and unpredictable portions in spike glycoprotein group are statistically different from hemagglutinin-esterase precursor, membrane protein and nucleocapsid protein groups at p<0.05 level. the predictable and unpredictable portions in spike glycoprotein group are statistically different from hemagglutinin-esterase precursor, and membrane protein groups at p<0.05 level.Electronic Supplementary Material is available for this article if you access the article at .     

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     

An ab initio study of di- and trifluorobenzene–benzene complexes as relevant to carbonic anhydrase II–drug interactions

Ab initio calculations at the MP2/6-31G* level have shown that variously substituted di- and trifluorobenzenes form non-covalent complexes with benzene that adopt either aromatic–aromatic or H---F binding, the choice being determined by the pattern of fluorination. The binding energies of these structures are from 3.4 to 4.5 kcal mol^–1. This range is large enough to account for observed variations in the binding affinity of a library of fluoroaromatic inhibitors of carbonic anhydrase. This enzyme has an aromatic amino acid at a central position in the active site. The diverse modes of binding of the dimers also suggest that aggregates of fluorobenzenes might adopt specified 3-dimensional shapes in the solid state. Figure Observed structure for 1,2-difluorobenzene     

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.     

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.     

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     

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.     

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.     

Theoretical mechanism of the formation of cholesteryl chloride from cholesterol and thionyl chloride

This work describes a theoretical approach to the substitution reaction mechanism involving the conversion of cholesterol to cholesteryl chloride. Two chlorosulfite ester molecules were formed as intermediates. An iso-steroid was found as the transition state. The final product was cholesteryl chloride and the side products were HCl and SO₂. Calculations were carried out at high level Hartree–Fock theory, using the 6–31G* basis set. From the electronic structure of the reactants, the most important physicochemical properties involved in the reaction pathway were used. Thus, to determine the participation of each molecule and to explain the mechanism of reaction; the total energy, HOMO and LUMO, atomic orbital contribution to frontier orbitals formation, electrostatic potentials, atomic charges, hardness and dipole moment were used. Characterization of intermediates and transition state was supported by their respective energy minima, fundamental frequencies and equilibrium geometry.Figure Synopsis of the reaction pathway. The reaction starts when the lone pair of the Ch oxygen interacts with the sulfur atom, releasing a chloride ion. As a result, the first intermediate is formed. Next, in the first intermediate the nucleophilic chloride ion bonds the electrophilic hydrogen atom, releasing HCl and yielding the second intermediate. In the second intermediate, the electrophilic H-atom from HCl bonds with the lone pair of the Cl atom adjacent to the sulfur atom, restoring HCl. Concurrently, SO₂ is liberated and causes the formation of the C3-C5 partial bond and breaking of the C5-C6 -bond leading to the transition state. In the transition state, the electrophilic H from HCl bonds with the Cl lone pair at C6-Cl, forming HCl again and leaving the C6 atom electron-deficient, which restores the C5-C6 -bond and breaks the C3-C5 partial bond. Finally, the electrophilic C3 atom and the nucleophilic Cl atom form a bond, yielding cholesteryl chloride. HCl and SO₂ are also formed as side products. The arrows show the rearrangement of electrons.     

Homology modeling and PAPS ligand (cofactor) binding study of bovine phenol sulfotransferase

In order to understand the mechanisms of ligand binding and the interaction between the ligand and the bovine phenol sulfotransferase, (bSULT1A1, EC 2.8.2.1) a three-dimensional (3D) model of the bSULT1A1 is generated based on the crystal structure of the estrogen sulfotransferase (PDB code 1AQU) by using the InsightII/Homology module. With the aid of the molecular mechanics and molecular dynamics methods, the final refined model is obtained and is further assessed by Profile-3D and ProStat, which show that the refined model is reliable. With this model, a flexible docking study is performed and the results indicate that 3-phosphoadenosine-5- phosphosulfate (PAPS) is a more preferred ligand than coenzyme A (CoA), and that His108 forms hydrogen bond with PAPS, which is in good agreement with the experimental results. From these docking studies, we also suggest that Phe255, Phe24 and Tyr169 in bSULT1A1 are three important determinant residues in binding as they have strong van-der-Waals contacts with the ligand. The hydrogen–bonding interactions also play an important role for the stability of the complex. Our results may be helpful for further experimental investigations.Figure The final 3D-structure of bSULT1A1. The structure is obtained by energy minimizing an average conformation over the last 100 ps of MD simulation. The -helix is represented in red and the -sheet in yellow.     

Normal coordinate analyses of 3,5-dichlorophenylcyanamide

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