Potential energy surface and moleculardynamics of MbNO: existence of an unsuspected FeON minimum

Ligands such as CO, O2, or NO are involved in the biological function of myoglobin. Here we investigate the energetics and dynamics of NO interacting with the Fe(II) heme group in native myoglobin using ab initio and molecular dynamics simulations. At the global minimum of the ab initio potential energy surface (PES), the binding energy of 23.4 kcal/mol and the Fe-NO structure compare well with the experimental results. Interestingly, the PES is found to exhibit two minima: There exists a metastable, linear Fe-O-N minimum in addition to the known, bent Fe-N-O global minimum conformation. Moreover, the T-shaped configuration is found to be a saddle point, in contrast to the corresponding minimum for NO interacting with Fe(III). To use the ab initio results for finite temperature molecular dynamics simulations, an analytical function was fitted to represent the Fe-NO interaction. The simulations show that the secondary minimum is dynamically stable up to 250 K and has a lifetime of several hundred picoseconds at 300 K. The difference in the topology of the heme-NO PES from that assumed previously (one deep, single Fe-NO minimum) suggests that it is important to use the full PES for a quantitative understanding of this system. Why the metastable state has not been observed in the many spectroscopic studies of myoglobin interacting with NO is discussed, and possible approaches to finding it are outlined.

Dopant-vacancy binding effects in Li-doped magnesium hydride

We use a combination of ab initio calculations and statistical mechanics to investigate the substitution of Li+ for Mg2+ in magnesium hydride (MgH2) accompanied by the formation of hydrogen vacancies with positive charge (with respect to the original ion at the site). We show that the binding energy between dopants and vacancy defects leads to a significant fraction of trapped vacancies and therefore a dramatic reduction in the number of free vacancies available for diffusion. The concentration of free vacancies initially increases with dopant concentration but reaches a maximum at around 1 mol % Li doping and slowly decreases with further doping. At the optimal level of doping, the corresponding concentration of free vacancies is much higher than the equilibrium concentrations of charged and neutral vacancies in pure MgH2 at typical hydrogen storage conditions. We also show that Li-doped MgH2 is thermodynamically metastable with respect to phase separation into pure magnesium and lithium hydrides at any significant Li concentration, even after considering the stabilization provided by dopant-vacancy interactions and configurational entropic effects. Our results suggest that lithium doping may enhance hydrogen diffusion hydride but only to a limited extent determined by an optimal dopant concentration and conditioned to the stability of the doped phase.

Many-body energy decomposition of hydrogen-bonded glycine clusters in gas-phase

A detailed analysis of the many-body contribution to the interaction energies of the gas-phase hydrogen-bonded glycine clusters, (Gly)(N), N = 1-4 is presented. The energetics of the hydrogen-bonded dimer, trimer and tetramer complexes have been analyzed using density-functional theory. The magnitude of the two-through four-body energy terms have been calculated and compared. The relaxation energy and the two-body energy terms are the principal contributors to the total binding energy. Four-body contribution is negligible. However, the three-body contribution is found to be sizable and the formation of the cyclic glycine trimer presents geometric strains that make it less favorable. (C) 2010 Elsevier B.V. All rights reserved.

Dimensional Compactification and Two-Particle Binding

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

Binding and structure of tetramers in the scaling limit

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

Low-energy correlations in the positronium-hydrogen-atom system

Employing a nonlocal model potential for electron exchange we study positronium-hydrogen-atom (Ps-H) scattering using a five-state coupled-channel model allowing for Ps(2s,2p)H(1s) and Ps(1s)H(2s,2p) excitations. We find remarkable correlations among S-wave Ps-H binding energy, scattering length, effective range, and resonance energy in the electronic singlet state. Using these correlations we predict fairly accurate values of singlet Ps-H scattering length (3.50a0) and effective range (1.65a0).

Interaction of Bacillus thuringiensis Cry1 and Vip3A proteins with Spodoptera frugiperda midgut binding sites

Vip3Aa, Vip3Af, Cry1Ab, and Cry1Fa were tested for their toxicities and binding interactions. Vip3A proteins were more toxic than Cry1 proteins. Binding assays showed independent specific binding sites for Cry1 and Vip3A proteins. Cry1Ab and Cry1Fa competed for the same binding sites, whereas Vip3Aa competed for those of Vip3Af. Copyright © 2009, American Society for Microbiology. All Rights Reserved.

Search for new physics in final states with a lepton and missing transverse energy in pp collisions at the LHC

This Letter describes the search for an enhanced production rate of events with a charged lepton and a neutrino in high-energy pp collisions at the LHC. The analysis uses data collected with the CMS detector, with an integrated luminosity of 5.0 fb-1 at √s=7 TeV, and a further 3.7 fb -1 at √s=8 TeV. No evidence is found for an excess. The results are interpreted in terms of limits on a heavy charged gauge boson (W ′) in the sequential standard model, a split universal extra dimension model, and contact interactions in the helicity-nonconserving model. For the last, values of the binding energy below 10.5 (8.8) TeV in the electron (muon) channel are excluded at a 95% confidence level. Interpreting the ℓν final state in terms of a heavy W′ with standard model couplings, masses below 2.90 TeV are excluded. © 2013 CERN.

Molecular design of new P3HT derivatives: Adjusting electronic energy levels for blends with PCBM

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

Theoretical prevision for the low-energy 3S1-3D1 mixing parameters

Using a recent shape-independent approximation for the 3S1-3D1 mixing parameter, theoretical prevision for the low-energy mixing parameters is made. The present prevision is consistent with the deuteron binding energy, its asymptotic D-state to S-state ratio, ηd, the triplet-scattering length, and the meson exchange tail of the tensor nucleon-nucleon potential. The theoretical prevision up to an incident laboratory energy of 25 MeV is consistent with the recent multi-energy determination of mixing parameters, but is much higher than many single-energy determinations of the same. The low single-energy values of the mixing parameter could be reproduced by meson-theoretical potentials only with a substantially reduced ηd. © 1994 The American Physical Society.

Metal binding selectivity of oxa-aza macrocyclic ligand: a DFT study of first- and second-row transition metal for four coordination systems

A detailed theoretical study of the 1,7,1l,17-tetraoxa-2,6,12,16-tetraaza-cycloeicosane ligand ([20]AneN(4)O(4)) coordinated to Fe2+, Co2+, Ni2+, Ru2+, Rh2+, and Pd2+ transition metal ions was carried out with the B3LYP method. Two different cases were performed: when nitrogen is the donor atom (1a (q) ) and also with the oxygen as the donor atom (1b (q) ). For all the cases performed in this study 1a (q) structures were always more stable than the 1b (q) ones. Considering each row is possible to see that the energy increases with the increase of the atomic number. The M2+ cation binding energies for the 1a (q) complexes increase with the following order: Fe2+ < Ru2+ < Co2+ < Ni2+ < Rh2+ < Pd2+.

Benchmark Structures and Binding Energies of Small Water Clusters with Anharmonicity Corrections

For (H2O)n where n = 1–10, we used a scheme combining molecular dynamics sampling with high level ab initio calculations to locate the global and many low lying local minima for each cluster. For each isomer, we extrapolated the RI-MP2 energies to their complete basis set limit, included a CCSD(T) correction using a smaller basis set and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled and unscaled harmonic vibrational frequencies. The vibrational scaling factors were determined specifically for water clusters by comparing harmonic frequencies with VPT2 fundamental frequencies. We find the CCSD(T) correction to the RI-MP2 binding energy to be small (<1%) but still important in determining accurate conformational energies. Anharmonic corrections are found to be non-negligble; they do not alter the energetic ordering of isomers, but they do lower the free energies of formation of the water clusters by as much as 4 kcal/mol at 298.15 K.

Benchmark Structures and Binding Energies of Small Water Clusters with Anharmonicity Corrections

For (H2O)n where n = 1–10, we used a scheme combining molecular dynamics sampling with high level ab initio calculations to locate the global and many low lying local minima for each cluster. For each isomer, we extrapolated the RI-MP2 energies to their complete basis set limit, included a CCSD(T) correction using a smaller basis set and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled and unscaled harmonic vibrational frequencies. The vibrational scaling factors were determined specifically for water clusters by comparing harmonic frequencies with VPT2 fundamental frequencies. We find the CCSD(T) correction to the RI-MP2 binding energy to be small (

Unraveling principles of lead discovery: from unfrustrated energy landscapes to novel molecular anchors.

The search for novel leads is a critical step in the drug discovery process. Computational approaches to identify new lead molecules have focused on discovering complete ligands by evaluating the binding affinity of a large number of candidates, a task of considerable complexity. A new computational method is introduced in this work based on the premise that the primary molecular recognition event in the protein binding site may be accomplished by small core fragments that serve as molecular anchors, providing a structurally stable platform that can be subsequently tailored into complete ligands. To fulfill its role, we show that an effective molecular anchor must meet both the thermodynamic requirement of relative energetic stability of a single binding mode and its consistent kinetic accessibility, which may be measured by the structural consensus of multiple docking simulations. From a large number of candidates, this technique is able to identify known core fragments responsible for primary recognition by the FK506 binding protein (FKBP-12), along with a diverse repertoire of novel molecular cores. By contrast, absolute energetic criteria for selecting molecular anchors are found to be promiscuous. A relationship between a minimum frustration principle of binding energy landscapes and receptor-specific molecular anchors in their role as "recognition nuclei" is established, thereby unraveling a mechanism of lead discovery and providing a practical route to receptor-biased computational combinatorial chemistry.

A mean field model of ligand-protein interactions: implications for the structural assessment of human immunodeficiency virus type 1 protease complexes and receptor-specific binding.

We propose a general mean field model of ligand-protein interactions to determine the thermodynamic equilibrium of a system at finite temperature. The method is employed in structural assessments of two human immuno-deficiency virus type 1 protease complexes where the gross effects of protein flexibility are incorporated by utilizing a data base of crystal structures. Analysis of the energy spectra for these complexes has revealed that structural and thermo-dynamic aspects of molecular recognition can be rationalized on the basis of the extent of frustration in the binding energy landscape. In particular, the relationship between receptor-specific binding of these ligands to human immunodeficiency virus type 1 protease and a minimal frustration principle is analyzed.