Generation and Analysis of Interaction Energy Maps of p-Substituted Benzoic Acid [(5-Nitro-thiophen-2-yl)-methylene]-Hydrazides Active against MRSA

In this preliminary study eighteen p-substituted benzoic acid [(5-nitro-thiophen-2-yl)-methylene]-hydrazides with antimicrobial activity were evaluated against multidrug-resistant Staphylococcus aureus, correlating the three-dimensional characteristics of the ligands with their respective bioactivities. The computer programs Sybyl and CORINA were used, respectively, for the design and three-dimensional conversion of the ligands. Molecular interaction fields were calculated using GRID program. Calculations using Volsurf resulted in a statistically consistent model with 48 structural descriptors showing that hydrophobicity is a fundamental property in the analyzed biological response.

Casimir interaction between a perfect conductor and graphene described by the Dirac model

We adopt the Dirac model for graphene and calculate the Casimir interaction energy between a plane suspended graphene sample and a parallel plane perfect conductor. This is done in two ways. First, we use the quantum-field-theory approach and evaluate the leading-order diagram in a theory with 2+1-dimensional fermions interacting with 3+1-dimensional photons. Next, we consider an effective theory for the electromagnetic field with matching conditions induced by quantum quasiparticles in graphene. The first approach turns out to be the leading order in the coupling constant of the second one. The Casimir interaction for this system appears to be rather weak. It exhibits a strong dependence on the mass of the quasiparticles in graphene.

Two complementary molecular energy decomposition schemes: the Mayer and Ziegler-Rauk methods in comparison

In the present paper we discuss and compare two different energy decomposition schemes: Mayer's Hartree-Fock energy decomposition into diatomic and monoatomic contributions [Chem. Phys. Lett. 382, 265 (2003)], and the Ziegler-Rauk dissociation energy decomposition [Inorg. Chem. 18, 1558 (1979)]. The Ziegler-Rauk scheme is based on a separation of a molecule into fragments, while Mayer's scheme can be used in the cases where a fragmentation of the system in clearly separable parts is not possible. In the Mayer scheme, the density of a free atom is deformed to give the one-atom Mulliken density that subsequently interacts to give rise to the diatomic interaction energy. We give a detailed analysis of the diatomic energy contributions in the Mayer scheme and a close look onto the one-atom Mulliken densities. The Mulliken density ρA has a single large maximum around the nuclear position of the atom A, but exhibits slightly negative values in the vicinity of neighboring atoms. The main connecting point between both analysis schemes is the electrostatic energy. Both decomposition schemes utilize the same electrostatic energy expression, but differ in how fragment densities are defined. In the Mayer scheme, the electrostatic component originates from the interaction of the Mulliken densities, while in the Ziegler-Rauk scheme, the undisturbed fragment densities interact. The values of the electrostatic energy resulting from the two schemes differ significantly but typically have the same order of magnitude. Both methods are useful and complementary since Mayer's decomposition focuses on the energy of the finally formed molecule, whereas the Ziegler-Rauk scheme describes the bond formation starting from undeformed fragment densities

Accurate calculations of intermolecular interaction energies using explicitly correlated coupled cluster wave functions and a dispersion-weighted MP2 method.

Explicitly correlated coupled-cluster calculations of intermolecular interaction energies for the S22 benchmark set of Jurecka, Sponer, Cerny, and Hobza (Chem. Phys. Phys. Chem. 2006, 8, 1985) are presented. Results obtained with the recently proposed CCSD(T)-F12a method and augmented double-zeta basis sets are found to be in very close agreement with basis set extrapolated conventional CCSD(T) results. Furthermore, we propose a dispersion-weighted MP2 (DW-MP2) approximation that combines the good accuracy of MP2 for complexes with predominately electrostatic bonding and SCS-MP2 for dispersion-dominated ones. The MP2-F12 and SCS-MP2-F12 correlation energies are weighted by a switching function that depends on the relative HF and correlation contributions to the interaction energy. For the S22 set, this yields a mean absolute deviation of 0.2 kcal/mol from the CCSD(T)-F12a results. The method, which allows obtaining accurate results at low cost, is also tested for a number of dimers that are not in the training set.

Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors.

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.

Two complementary molecular energy decomposition schemes: the Mayer and Ziegler-Rauk methods in comparison

In the present paper we discuss and compare two different energy decomposition schemes: Mayer's Hartree-Fock energy decomposition into diatomic and monoatomic contributions [Chem. Phys. Lett. 382, 265 (2003)], and the Ziegler-Rauk dissociation energy decomposition [Inorg. Chem. 18, 1558 (1979)]. The Ziegler-Rauk scheme is based on a separation of a molecule into fragments, while Mayer's scheme can be used in the cases where a fragmentation of the system in clearly separable parts is not possible. In the Mayer scheme, the density of a free atom is deformed to give the one-atom Mulliken density that subsequently interacts to give rise to the diatomic interaction energy. We give a detailed analysis of the diatomic energy contributions in the Mayer scheme and a close look onto the one-atom Mulliken densities. The Mulliken density ρA has a single large maximum around the nuclear position of the atom A, but exhibits slightly negative values in the vicinity of neighboring atoms. The main connecting point between both analysis schemes is the electrostatic energy. Both decomposition schemes utilize the same electrostatic energy expression, but differ in how fragment densities are defined. In the Mayer scheme, the electrostatic component originates from the interaction of the Mulliken densities, while in the Ziegler-Rauk scheme, the undisturbed fragment densities interact. The values of the electrostatic energy resulting from the two schemes differ significantly but typically have the same order of magnitude. Both methods are useful and complementary since Mayer's decomposition focuses on the energy of the finally formed molecule, whereas the Ziegler-Rauk scheme describes the bond formation starting from undeformed fragment densities

A PRELIMINARY ANALYSIS OF THE ENERGY TRANSFER BETWEEN THE DARK SECTORS OF THE UNIVERSE

We study the mutual interaction between the dark sectors (dark matter and dark energy) of the Universe by resorting to the extended thermodynamics of irreversible processes and constrain the former with supernova type Ia data. As a by-product, the present dark matter temperature results are not extremely small and can meet the independent estimate of the temperature of the gas of sterile neutrinos.

On energy transfer between vibrating systems under linear and nonlinear interactions

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Energy transfer between Pr3+ ions in a fluoroindate glass

The efficiency of energy transfer (ET) between Pr3+ ions in a fluoroindate glass is determined. ET rates, WET, were determined for dilute samples and the results show a dependence of WET on the Pr3+ concentration. ET processes which contribute to resonance fluorescence and frequency upconversion emission were studied. The origin of the interaction energy among the Pr3+ ions was determined to be dipole - dipole. © 1998 Elsevier Science B.V. All rights reserved.

Study of water and dimethylformamide interaction by computer simulation

Monte Carlo simulations of water-dimethylformamide (DMF) mixtures were performed in the isothermal and isobaric ensemble at 298.15 K and 1 atm. The intermolecular interaction energy was calculated using the classical 6-12 Lennard-Jones pairwise potential plus a Coulomb term. The TIP4P model was used for simulating water molecules, and a six-site model previously optimised by us was used to represent DMF. The potential energy for the water-DMF interaction was obtained via standard geometric combining rules using the original potential parameters for the pure liquids. The radial distribution functions calculated for water-DMF mixtures show well characterised hydrogen bonds between the oxygen site of DMF and hydrogen of water. A structureless correlation curve was observed for the interaction between the hydrogen site of the carbonyl group and the oxygen site of water. Hydration effects on the stabilisation of the DMF molecule in aqueous solution have been investigated using statistical perturbation theory. The results show that energetic changes involved in the hydration process are not strong enough to stabilise another configuration of DMF than the planar one.

Protein-protein interaction: a genetic selection for compensating mutations at the barnase-barstar interface.

Barnase and barstar are trivial names of the extracellular RNase and its intracellular inhibitor produced by Bacillus amyloliquefaciens. Inhibition involves the formation of a very tight one-to-one complex of the two proteins. With the crystallographic solution of the structure of the barnase-barstar complex and the development of methods for measuring the free energy of binding, the pair can be used to study protein-protein recognition in detail. In this report, we describe the isolation of suppressor mutations in barstar that compensate for the loss in interaction energy caused by a mutation in barnase. Our suppressor search is based on in vivo selection for barstar variants that are able to protect host cells against the RNAse activity of those barnase mutants not properly inhibited by wild-type barstar. This approach utilizes a plasmid system in which barnase expression is tightly controlled to keep the mutant barnase gene silent. When expression of barnase is turned on, failure to form a complex between the mutant barnase and barstar has a lethal effect on host cells unless overcome by substitution of the wild-type barstar by a functional suppressor derivative. A set of barstar suppressors has been identified for barnase mutants with substitutions in two amino acid positions (residues 102 and 59), which are critically involved in both RNase activity and barstar binding. The mutations selected as suppressors could not have been predicted on the basis of the known protein structures. The single barstar mutation with the highest information content for inhibition of barnase (H102K) has the substitution Y30W. The reduction in binding caused by the R59E mutation in barnase can be partly reversed by changing Glu-76 of barstar, which forms a salt bridge with the Arg-59 in the wild-type complex, to arginine, thus completing an interchange of the two charges.

Effect of surface-perturbed intermolecular interaction on adsorption of simple gases on a graphitized carbon surface

In this paper, we investigate the effect of the solid surface on the fluid-fluid intermolecular potential energy. This modified fluid-fluid interaction energy due to the inducement of a solid surface is used in the grand canonical Monte Carlo (GCMC) simulation of various noble gases, nitrogen, and methane on graphitized thermal carbon black. This effect is such that the effective interaction potential energy between two particles close to surface is less than the potential energy if the solid substrate is not present. With this modification the GCMC simulation results agree extremely well with the experimental data over a wide range of pressures while the simulation results with the unmodified potential energy give rise to a shoulder near the neighborhood of monolayer coverage and the significant overprediction of the second and higher layer coverages. In particular the unmodified GCMC results exhibit very sharp change in those higher layers while the experimental data have a much gradual change in the uptake. We will illustrate this theory with adsorption data of argon, xenon, neon, nitrogen, and methane on graphitized thermal carbon black.

Hydrogen bond interactions between acetone and supercritical water

Hydrogen bond interactions between acetone and supercritical water are investigated using a combined and sequential Monte Carlo/quantum mechanics (S-MC/QM) approach. Simulation results show a dominant presence of con. gurations with one hydrogen bond for different supercritical states, indicating that this specific interaction plays an important role on the solvation properties of acetone in supercritical water. Using QM MP2/aug-cc-pVDZ the calculated average interaction energy reveals that the hydrogen-bonded acetone-water complex is energetically more stable under supercritical conditions than ambient conditions and its stability is little affected by variations of temperature and/or pressure. All average results reported here are statistically converged.

Thermodynamic origin of cooperativity in actomyosin interactions: The coupling of short-range interactions with actin bending stiffness in an Ising-like model

We present Monte Carlo simulations for a molecular motor system found in virtually all eukaryotic cells, the acto-myosin motor system, composed of a group of organic macromolecules. Cell motors were mapped to an Ising-like model, where the interaction field is transmitted through a tropomyosin polymer chain. The presence of Ca(2+) induces tropomyosin to block or unblock binding sites of the myosin motor leading to its activation or deactivation. We used the Metropolis algorithm to find the transient and the equilibrium states of the acto-myosin system composed of solvent, actin, tropomyosin, troponin, Ca(2+), and myosin-S1 at a given temperature, including the spatial configuration of tropomyosin on the actin filament surface. Our model describes the short- and long-range cooperativity during actin-myosin binding which emerges from the bending stiffness of the tropomyosin complex. We found all transition rates between the states only using the interaction energy of the constituents. The agreement between our model and experimental data also supports the recent theory of flexible tropomyosin.

Nanogravimetric and voltammetric studies of a Pt-Rh alloy surface and its behavior for methanol oxidation

This paper describes the preparation of a Pt-Rh alloy surface electrodeposited on Pt electrodes and its electrocatalytic characterization for methanol oxidation. The X-ray photoelectronic spectroscopy ( XPS) results demonstrate that the surface composition is approximately 24 at-% Rh and 76 % Pt. The cyclic voltammetry (CV) and electrochemical quartz crystal (EQCN) results for the alloy were associated, for platinum, to the well known profile in acidic medium. For Rh, on the alloy, the generation of rhodium hydroxide species (Rh(OH)(3) and RhO(OH)(3)) was measured. During the successive oxidation-reduction cycles the mass returns to its original value, indicating the reversibility of the processes. It was not observed rhodium dissolution during the cycling. The 76/24 at % Pt-Rh alloy presented singular electrocatalytic activity for methanol electrooxidation, which started at more negative potentials compared to pure Pt (70 mV). During the sweep towards more negative potentials, there is only weak CO re-adsorption on both Rh and Pt-Rh alloy surfaces, which can be explained by considering the interaction energy between Rh and CO.