Calculations of the free energy of interaction of the c-Fos−c-Jun coiled coil: Effects of the solvation model and the inclusion of polarization effects

The leucine zipper region of activator protein-1 (AP-1) comprises the c-Jun and c-Fos proteins and constitutes a well-known coiled coil protein−protein interaction motif. We have used molecular dynamics (MD) simulations in conjunction with the molecular mechanics/Poisson−Boltzmann generalized-Born surface area [MM/PB(GB)SA] methods to predict the free energy of interaction of these proteins. In particular, the influence of the choice of solvation model, protein force field, and water potential on the stability and dynamic properties of the c-Fos−c-Jun complex were investigated. Use of the AMBER polarizable force field ff02 in combination with the polarizable POL3 water potential was found to result in increased stability of the c-Fos−c-Jun complex. MM/PB(GB)SA calculations revealed that MD simulations using the POL3 water potential give the lowest predicted free energies of interaction compared to other nonpolarizable water potentials. In addition, the calculated absolute free energy of binding was predicted to be closest to the experimental value using the MM/GBSA method with independent MD simulation trajectories using the POL3 water potential and the polarizable ff02 force field, while all other binding affinities were overestimated.

Theoretical calculations of base-base interactions in nucleic acids: II Stacking interactions in polynucleotides

Base-base interactions were computed for single- and double- stranded polynucleotides, for all possible base sequences. In each case, both right and left stacking arrangements are energetically possible. The preference of one over the other depends upon the base-sequence and the orientation of the bases with respect to helix-axis. Inverted stacking arrangement is also energetically possible for both single- and double-stranded polynucleotides. Finally, interaction energies of a regular duplex and the alternative structures3 were compared. It was found that the type II model3 is energetically more favourable than the rest.

Studies on the conformation of amino acids XII. Energy calculations on prolyl residue

The favoured conformations of the prolyl residue have been obtained by calculating their potential energies arising from bond-angle strain, torsion-angle strain, non-bonded and electrostatic interatomic energies. In addition to the five membered ring, the peptide unit at the amino end (with ω = 180°) and the C′ atom at the carboxyl end have been taken into account. It is found that there are two local minima in the configurational space of the parameters defining the conformation, as is actually observed-one (denoted by B) with Cγ displaced on the same side as C′, which is lower in energy than the other (denoted by A) with Cγ displaced on the opposite side of C′. The other four atoms Cδ, N, Cα, Cβ are nearly in a plane. The conformations of minimum energy (for both A and B) have bond angles very close to the mean observed values while the torsion angles are well within the range observed in various structures for each type. Taking into account the fact that the influence of neighbouring molecules in a crystal structure may make the conformation of a molecule different from the minimal one, the ranges of the conformational parameters for which the energy is within 0.6 kcal/mole above the minimum value (called the "most probable range") and within 1.2 kcal/mole (called the "probable range") have been determined. The ranges thus obtained, agree well with observation, and most of the observed data lie within the most probable ranges, although differing appreciably from the conformation of minimum energy. The study has been extended, in a limited way, to the conformation of the ring in the amino acid proline. Since the nitrogen is tetrahedral in this (as contrasted with being planar in the prolyl residue), it is found that any one of the five atoms can be out of plane (either way), with the other four lying nearly in a plane. These correspond to low energy conformations (up to 1.2 kcal/mole above the minimum). One such example, in which the Cα atom is out of plane is known for dl-proline · HCl. It is also shown that in these calculations energies due to bond length distortions can be neglected to a good degree of approximation, provided the 'best' values of the bond lengths for the particular compound are used in the theoretical calculations.

A Thermodynamic Proposal for the Exposition of Critically and other Subtle Features in Self-Deflagrating Ammonium Perchlorate Systems

A first comprehensive investigation on the deflagration of ammonium perchlorate (AP) in the subcritical regime, below the low pressure deflagration limit (LPL, 2.03 MPa) christened as regime I$^{\prime}$, is discussed by using an elegant thermodynamic approach. In this regime, deflagration was effected by augmenting the initial temperature (T$_{0}$) of the AP strand and by adding fuels like aliphatic dicarboxylic acids or polymers like carboxy terminated polybutadiene (CTPB). From this thermodynamic model, considering the dependence of burning rate ($\dot{r}$) on pressure (P) and T$_{0}$, the true condensed (E$_{\text{s,c}}$) and gas phase (E$_{\text{s,g}}$) activation energies, just below and above the surface respectively, have been obtained and the data clearly distinguishes the deflagration mechanisms in regime I$^{\prime}$ and I (2.03-6.08 MPa). Substantial reduction in the E$_{\text{s,c}}$ of regime I$^{\prime}$, compared to that of regime I, is attributed to HClO$_{4}$ catalysed decomposition of AP. HClO$_{4}$ formation, which occurs only in regime I$^{\prime}$, promotes dent formation on the surface as revealed by the reflectance photomicrographs, in contrast to the smooth surface in regime I. The HClO$_{4}$ vapours, in regime I$^{\prime}$, also catalyse the gas phase reactions and thus bring down the E$_{\text{s,g}}$ too. The excess heat transferred on to the surface from the gas phase is used to melt AP and hence E$_{\text{s,c}}$, in regime I, corresponds to the melt AP decomposition. It is consistent with the similar variation observed for both the melt layer thickness and $\dot{r}$ as a function of P. Thermochemical calculations of the surface heat release support the thermodynamic model and reveal that the AP sublimation reduces the required critical exothermicity of 1108.8 kJ kg$^{-1}$ at the surface. It accounts for the AP not sustaining combustion in the subcritical regime I$^{\prime}$. Further support for the model comes from the temperature-time profiles of the combustion train of AP. The gas and condensed phase enthalpies, derived from the profile, give excellent agreement with those computed thermochemically. The $\sigma _{\text{p}}$ expressions derived from this model establish the mechanistic distinction of regime I$^{\prime}$ and I and thus lend support to the thermodynamic model. On comparing the deflagration of strand against powder AP, the proposed thermodynamic model correctly predicts that the total enthalpy of the condensed and gas phases remains unaltered. However, 16% of AP particles undergo buoyant lifting into the gas phase in the `free board region' (FBR) and this renders the demarcation of the true surface difficult. It is found that T$_{\text{s}}$ lies in the FBR and due to this, in regime I$^{\prime}$, the E$_{\text{s,c}}$ of powder AP matches with the E$_{\text{s,g}}$ of the pellet. The model was extended to AP/dicarboxylic acids and AP/CTPB mixture. The condensed ($\Delta $H$_{1}$) and gas phase ($\Delta $H$_{2}$) enthalpies were obtained from the temperature profile analyses which fit well with those computed thermochemically. The $\Delta $H$_{1}$ of the AP/succinic acid mixture was found just at the threshold of sustaining combustion. Indeed the lower homologue malonic acid, as predicted, does not sustain combustion. In vaporizable fuels like sebacic acid the E$_{\text{s,c}}$ in regime I$^{\prime}$, understandably, conforms to the AP decomposition. However, the E$_{\text{s,c}}$ in AP/CTPB system corresponds to the softening of the polymer which covers AP particles to promote extensive condensed phase reactions. The proposed thermodynamic model also satisfactorily explains certain unique features like intermittent, plateau and flameless combustion in AP/ polymeric fuel systems.

Size dependence of the bulk modulus of semiconductor nanocrystals from first-principles calculations

The variation in the bulk modulus of semiconductor nanoparticles has been studied within first-principles electronic-structure calculations using the local density approximation (LDA) for the exchange correlation. Quantum Monte Carlo calculations carried out for a silicon nanocrystal Si87H76 provided reasonable agreement with the LDA results. An enhancement was observed in the bulk modulus as the size of the nanoparticle was decreased, with modest enhancements being predicted for the largest nanoparticles studied here, a size just accessible in experiments. To access larger sizes, we fit our calculated bulk moduli to the same empirical law for all materials, the asymptote of which is the bulk value of the modulus. This was found to be within 2-10% of the independently calculated value. The origin of the enhancement has been discussed in terms of Cohen's empirical law M.L. Cohen, Phys. Rev. B 32, 7988 (1985)] as well as other possible scenarios.

Ab initio molecular orbital calculations on ion pair-water complexes of metal halides and oxides

Ab initio MO calculations are performed on a series of ion-molecular and ion pair-molecular complexes of H2O + MX (MX = LiF, LiCl, NaCl, BeO and MgO) systems. BSSE-corrected stabilization energies, optimized geometrical parameters, internal force constants and harmonic vibrational frequencies have been evaluated for all the structures of interest. The trends observed in the geometrical parameters and other properties calculated for the mono-hydrated contact ion pair complexes parallel those computed for the complexes of the individual ions. The bifurcated structures are found to be saddle points with an imaginary frequency corresponding to the rocking mode of water molecules. The solvent-shared ion pair complexes have high interaction energies. Trends in the internal force constant and harmonic frequency values are discussed in terms of ion-molecular and ion-pair molecular interactions.

Ab initio, DFT and transition state theory calculations on 1,2-HF, HCI and CIF elimination reactions from CH2F-CH2Cl

This paper reports ab intio, DFT and transition state theory (TST) calculations on HF, HCI and CIF elimination reactions from CH2Cl-CH2F molecule. Both the ground state and the transition state for HX elimination reactions have been optimized at HF, MP2 and DFT calculations with 6-31G*, 6-31G** and 6-311++G** basis sets. In addition, CCSD(T) single point calculations were carried out with MP2/6-311++G** optimized geometry for more accurate determination of the energies of the minima and transition state, compared to the other methods employed here. Classical barriers are converted to Arrhenius activation energy by TST calculations for comparisons with experimental results. The pre-exponential factors, A, calculated at all levels of theory are significantly larger than the experimental values. For activation energy, E-a DFT gives good results for HF elimination, within 4-8 W mol(-1) from experimental values. None of the methods employed, including CCSD(T), give comparable results for HCI elimination reactions. However, rate constants calculated by CCSD(T) method are in very good agreement with experiment for HCI elimination and they are in reasonable agreement for HF elimination reactions. Due to the strong correlation between A and E., the rate constants could be fit to a lower A and E-a (as given by experimental fitting, corresponding to a tight TS) or to larger A and E-a (as given by high level ab initio calculations, corresponding to a loose TS). The barrier for CIF elimination is determined to be 607 U mol(-1) at HF level and it is unlikely to be important for CH2FCH2Cl. Results for other CH2X-CH2Y (X,Y = F/Cl) are included for comparison.

Computational Study of the Migration of Rhenium from One Enantioface of an Olefin to the Other Facilitated by (C-H)center dot center dot center dot Re Interactions

The migration of a metal atom in a metal olefin complex from one pi face of the olefin to the opposite pi face has been rarely documented. Gladysz and co-workers showed that such a movement is indeed possible in monosubstituted chiral Re olefin complexes, resulting in diastereomerization. Interestingly, this isomerization occurred without dissociation, and on the basis of kinetic isotope effects, the involvement of a trans C-H bond was indicated. Either oxidative addition or an agostic interaction of the vinylic C-H(D) bond with the metal could account for the experimentally observed kinetic isotope effect. In this study we compute the free energy of activation for the migration of Re from one enantioface of the olefin to the other through various pathways. On the basis of DFT calculations at the B3LYP level we show that a trans (C-H)center dot center dot center dot Re interaction and trans C-H oxidative addition provide a nondissociative path for the diastereomerization. The trans (C-H)center dot center dot center dot Re interaction path is computed to be more favorable by 2.3 kcal mol(-1) than the oxidative addition path. While direct experimental evidence was not able to discount the migration of the metal through the formation of a eta(2)-arene complex (conducted tour mechanism), computational results at the B3LYP level show that it is energetically more expensive. Surprisingly, a similar analysis carried out at the M06 level computes a lower energy path for the conducted tour mechanism and is not consistent with the experimental isotope effects observed. Metal-(C-H) interactions and oxidative additions of the metal into C-H bonds are closely separated in energy and might contribute to unusual fluxional processes such as this diastereomerization.

1. Ultrasonic studies of binary liquid structure in the critical region. Theory and experiment for the 2,6-lutidine/water system. 2. Hartree-Fock calculations of electric polarizabilities of some simple atoms and molecules, and their practicality. 3. Calculation of vibrational transition probabilities in collinear atom-diatom and diatom-diatom collisions with Lennard-Jones interaction

Part 1. Many interesting visual and mechanical phenomena occur in the critical region of fluids, both for the gas-liquid and liquid-liquid transitions. The precise thermodynamic and transport behavior here has some broad consequences for the molecular theory of liquids. Previous studies in this laboratory on a liquid-liquid critical mixture via ultrasonics supported a basically classical analysis of fluid behavior by M. Fixman (e. g., the free energy is assumed analytic in intensive variables in the thermodynamics)--at least when the fluid is not too close to critical. A breakdown in classical concepts is evidenced close to critical, in some well-defined ways. We have studied herein a liquid-liquid critical system of complementary nature (possessing a lower critical mixing or consolute temperature) to all previous mixtures, to look for new qualitative critical behavior. We did not find such new behavior in the ultrasonic absorption ascribable to the critical fluctuations, but we did find extra absorption due to chemical processes (yet these are related to the mixing behavior generating the lower consolute point). We rederived, corrected, and extended Fixman's analysis to interpret our experimental results in these more complex circumstances. The entire account of theory and experiment is prefaced by an extensive introduction recounting the general status of liquid state theory. The introduction provides a context for our present work, and also points out problems deserving attention. Interest in these problems was stimulated by this work but also by work in Part 3.

Part 2. Among variational theories of electronic structure, the Hartree-Fock theory has proved particularly valuable for a practical understanding of such properties as chemical binding, electric multipole moments, and X-ray scattering intensity. It also provides the most tractable method of calculating first-order properties under external or internal one-electron perturbations, either developed explicitly in orders of perturbation theory or in the fully self-consistent method. The accuracy and consistency of first-order properties are poorer than those of zero-order properties, but this is most often due to the use of explicit approximations in solving the perturbed equations, or to inadequacy of the variational basis in size or composition. We have calculated the electric polarizabilities of H₂, He, Li, Be, LiH, and N₂ by Hartree-Fock theory, using exact perturbation theory or the fully self-consistent method, as dictated by convenience. By careful studies on total basis set composition, we obtained good approximations to limiting Hartree-Fock values of polarizabilities with bases of reasonable size. The values for all species, and for each direction in the molecular cases, are within 8% of experiment, or of best theoretical values in the absence of the former. Our results support the use of unadorned Hartree-Pock theory for static polarizabilities needed in interpreting electron-molecule scattering data, collision-induced light scattering experiments, and other phenomena involving experimentally inaccessible polarizabilities.

Part 3. Numerical integration of the close-coupled scattering equations has been carried out to obtain vibrational transition probabilities for some models of the electronically adiabatic H₂-H₂ collision. All the models use a Lennard-Jones interaction potential between nearest atoms in the collision partners. We have analyzed the results for some insight into the vibrational excitation process in its dependence on the energy of collision, the nature of the vibrational binding potential, and other factors. We conclude also that replacement of earlier, simpler models of the interaction potential by the Lennard-Jones form adds very little realism for all the complication it introduces. A brief introduction precedes the presentation of our work and places it in the context of attempts to understand the collisional activation process in chemical reactions as well as some other chemical dynamics.

Studies of magnetic interactions in Ni-doped ZnO from first-principles calculations

The magnetic interactions in Ni-doped ZnO are calculated using GGA and GGA + U method of density functional theory. The following three cases: (i) Ni-doped ZnO, (ii) (Ni, Al)-codoped ZnO, and (iii) (Ni, Li)-codoped ZnO are studied. The ferromagnetic ordering is always favorable for the three cases within GGA method. However, the ferromagnetic state is sometimes favorable after treating within the method of GGA + U. The GGA underestimates the correlated interactions especially when the Ni ions align directly to each other. (C) 2007 Elsevier B.V. All rights reserved.

ENERGY-LEVEL CALCULATIONS OF THE RECONSTRUCTED 90-DEGREES PARTIAL DISLOCATION IN SILICON

An LCAO scheme taking into account 10 atomic orbitals (s-, p-, and d-type) applied to a supercell containing 256 atoms is used to calculate the bound states of the reconstructed 90-degrees partial dislocation in Si. The results differ significantly from our earlier calculations on the unreconstructed 90-degrees partial using the same method. We find two bands separate from each other in the entire Brillouin zone and the upper band penetrates deep into the indirect band gap which is in contradiction with the general opinion that core reconstruction clears the band gap of dislocation states.

1ST-PRINCIPLES CALCULATIONS FOR QUASI-PARTICLE ENERGIES OF GAP AND GAAS

We have applied the Green-function method in the GW approximation to calculate quasiparticle energies for the semiconductors GaP and GaAs. Good agreement between the calculated excitation energies and the experimental results was achieved. We obtained calculated direct band gaps of GaP and GaAs of 2.93 and 1.42 eV, respectively, in comparison with the experimental values of 2.90 and 1.52 eV, respectively. An ab initio pseudopotential method has been used to generate basis wave functions and charge densities for calculating the dielectric matrix elements and self-enegies. To evaluate the dynamical effects of the screened interaction, the generalized-plasma-pole model has been utilized to extend the dielectric matrix elements from static results to finite frequencies. We presen the calculated quasiparticle energies at various high-symmetry points of the Brillouin zone and compare them with the experimental results and other calculations.

Ab initio calculations of differential cross sections for single charge transfer in He-3(2+)+He-4 collisions

The single charge transfer process in He-3(2+)+He-4 collisions is investigated using the quantum-mechanical molecular-orbital close-coupling method, in which the adiabatic potentials and radial couplings are calculated by using the ab initio multireference single- and double-excitation configuration interaction methods. The differential cross sections for the single charge transfer are presented at the laboratorial energies E = 6 keV and 10 keV for the projectile He-3(2+). Comparison with the existing data shows that the present results are better in agreement with the experimental measurements than other calculations in the dominant small angle scattering, which is attributed to the accurate calculations of the adiabatic potentials and the radial couplings.

Ab initio calculations of group 4 metallocene reaction mechanisms: atomic layer deposition and bond activation catalysis

Thin film dielectrics based on titanium, zirconium or hafnium oxides are being introduced to increase the permittivity of insulating layers in transistors for micro/nanoelectronics and memory devices. Atomic layer deposition (ALD) is the process of choice for fabricating these films, as it allows for high control of composition and thickness in thin, conformal films which can be deposited on substrates with high aspect-ratio features. The success of this method depends crucially on the chemical properties of the precursor molecules. A successful ALD precursor should be volatile, stable in the gas-phase, but reactive on the substrate and growing surface, leading to inert by-products. In recent years, many different ALD precursors for metal oxides have been developed, but many of them suffer from low thermal stability. Much promise is shown by group 4 metal precursors that contain cyclopentadienyl (Cp = C5H5-xRx) ligands. One of the main advantages of Cp precursors is their thermal stability. In this work ab initio calculations were carried out at the level of density functional theory (DFT) on a range of heteroleptic metallocenes [M(Cp)4-n(L)n], M = Hf/Zr/Ti, L = Me and OMe, in order to find mechanistic reasons for their observed behaviour during ALD. Based on optimized monomer structures, reactivity is analyzed with respect to ligand elimination. The order in which different ligands are eliminated during ALD follows their energetics which was in agreement with experimental measurements. Titanocene-derived precursors, TiCp*(OMe)3, do not yield TiO2 films in atomic layer deposition (ALD) with water, while Ti(OMe)4 does. DFT was used to model the ALD reaction sequence and find the reason for the difference in growth behaviour. Both precursors adsorb initially via hydrogen-bonding. The simulations reveal that the Cp* ligand of TiCp*(OMe)3 lowers the Lewis acidity of the Ti centre and prevents its coordination to surface O (densification) during both of the ALD pulses. Blocking this step hindered further ALD reactions and for that reason no ALD growth is observed from TiCp*(OMe)3 and water. The thermal stability in the gas phase of Ti, Zr and Hf precursors that contain cyclopentadienyl ligands was also considered. The reaction that was found using DFT is an intramolecular α-H transfer that produces an alkylidene complex. The analysis shows that thermal stabilities of complexes of the type MCp2(CH3)2 increase down group 4 (M = Ti, Zr and Hf) due to an increase in the HOMO-LUMO band gap of the reactants, which itself increases with the electrophilicity of the metal. The reverse reaction of α-hydrogen abstraction in ZrCp2Me2 is 1,2-addition reaction of a C-H bond to a Zr=C bond. The same mechanism is investigated to determine if it operates for 1,2 addition of the tBu C-H across Hf=N in a corresponding Hf dimer complex. The aim of this work is to understand orbital interactions, how bonds break and how new bonds form, and in what state hydrogen is transferred during the reaction. Calculations reveal two synchronous and concerted electron transfers within a four-membered cyclic transition state in the plane between the cyclopentadienyl rings, one π(M=X)-to-σ(M-C) involving metal d orbitals and the other σ(C-H)-to-σ(X-H) mediating the transfer of neutral H, where X = C or N. The reaction of the hafnium dimer complex with CO that was studied for the purpose of understanding C-H bond activation has another interesting application, namely the cleavage of an N-N bond and resulting N-C bond formation. Analysis of the orbital plots reveals repulsion between the occupied orbitals on CO and the N-N unit where CO approaches along the N-N axis. The repulsions along the N-N axis are minimized by instead forming an asymmetrical intermediate in which CO first coordinates to one Hf and then to N. This breaks the symmetry of the N-N unit and the resultant mixing of MOs allows σ(NN) to be polarized, localizing electrons on the more distant N. This allowed σ(CO) and π(CO) donation to N and back-donation of π*(Hf2N2) to CO. Improved understanding of the chemistry of metal complexes can be gained from atomic-scale modelling and this provides valuable information for the design of new ALD precursors. The information gained from the model decomposition pathway can be additionally used to understand the chemistry of molecules in the ALD process as well as in catalytic systems.

Drug calculations part 1: a critique of the formula used by nurses

The role of mathematics is integral to nursing practice, and careful and accurate calculations are important to help prevent medication errors. This two-part article examines different methods for solving drug calculation problems. The first part critiques the commonly taught nursing drug calculation formula. Part 2, to be published next week, explores a range of other methods that are used in practice to calculate drug dosages.