Anesthetic-protein interaction Random versus helix polylysine monolayers and interaction with 1-alkanols

Penetration of 1-alkanols into monolayers of hydrophobic polypeptides, poly(ε-benzyloxycarbonyl-l-lysine) and poly(ε-benzyloxycarbonyl-dl-lysine), was compared with their adsorption on the air/water interface in the absence of monolayers. The polypeptide prepared from l-lysine is generally considered to be in the α-helical form whereas dl-copolymer polypeptide contains random-coiled portions due to the structural incompatibility between the two isomers. The free energy of adsorption of 1-alkanols on the air/water interface at dilute concentrations was ?0.68 kcal·mol^?1 per methylene group and 0.15 kcal·mol^?1 for the hydroxyl group at 25°C. In the close-packed state, the surface area occupied by each molecule of 1-alkanols of varying carbon chain-lengths showed nearly a constant value of about 27.2 Ų, indicating perpendicular orientation of the alkanol molecules at the interface. About 75% of the water surface was covered by 1-butanol in this close-packed state. The mode of adsorption of 1-alkanols on the vacant air/water interface followed the Gibbs surface excess while the mode on the polypeptide membranes followed the Langmuir adsorption isotherm, indicating that the latter is characterized by the presence of a finite number of binding sites. The free energies of adsorption of 1-alkanols on the l-polymer monolayers were more negative than those on the vacant air/water interface and less negative than those on the dl-copolymer monolayers. Thus, the affinity of 1-alkanols to the interface was in the order of vacant air/water interface <l-polymer <dl-copolymer. The difference between the air/water interface and l-polymer was about 0.54 kcal·mol^?1 and that between l-polymer and dl-copolymer was 0.17 kcal·mol^?1 at 25°C: the adsorption of 1-alkanols to the dl-copolymer was favored compared to the l-polymer. The polar moieties of the backbone of the dl-copolymer may be exposed to the aqueous phase at the disordered portion. Dipole interaction between this portion and 1-alkanol molecules may account for the enhanced adsorption of the alkanols to the dl-copolymer.     

A new era of the prototropic tautomerism of the quercetin molecule: A QM/QTAIM computational advances

Abstract

In this study for the first time we have revealed and investigated in details 123 different prototropic tautomers of the most stable conformer of the quercetin molecule using quantum-mechanical calculations at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of QM theory. We have found that in the most energetically favorable prototropic tautomer mobile hydrogen atoms are localized at the О3, О3′, О4′, О5, and О7 exocyclic oxygen atoms. Molecular tautomers are in the range of the Gibbs free energies from 0.0 to 69.8?kcal·mol^?1, while zwitterionic ones – from 30.1 до 172.8?kcal·mol^?1 at normal conditions. It was also reliably established that the weakest point causing the decyclization of the molecule is its C ring – this reaction is launched by the transition of the proton from the C8H group to the endocyclic O1 oxygen atom. All prototropic tautomers, except two cases, are joined by the intramolecular cooperative specific interactions (from 1 to 5) – H-bonds and attractive van der Waals contacts, which have been revealed and characterized by QTAIM analysis.

Communicated by Ramaswamy H. Sarma     

A theoretical investigation on the conformation and the interaction of CHF2OCF2CHF2 (desflurane II) with one water molecule

The conformation and the interaction of CHF₂OCF₂CHF₂ (desflurane II) with one water molecule is investigated theoretically using the ab initio MP2/aug-cc-pvdz and DFT-based M062X/6-311++G(d,p) methods. The calculations include the optimized geometries, the harmonic frequencies of relevant vibrational modes along with a natural bond orbital (NBO) analysis including the NBO charges, the hybridization of the C atom and the intra- and intermolecular hyperconjugation energies. In the two most stable conformers, the CH bond of the F₂HCO- group occupies the gauche position. The hyperconjugation energies are about the same for both conformers and the conformational preference depends on the interaction between the non-bonded F and H atoms. The deprotonation enthalpies of the CH bonds are about the same for both conformers, the proton affinity of the less stable conformer being 3 kcal mol^?1 higher. Both conformers of desflurane II interact with water forming cyclic complexes characterized by CH…O and OH…F hydrogen bonds. The binding energies are moderate, ranging from ?2.4 to ?3.2 kcal mol^?1 at the MP2 level. The origin of the blue shifts of the ν(CH) vibrations is analyzed. In three of the complexes, the water molecule acts as an electron donor. Interestingly, in these cases a charge transfer is also directed to the non bonded OH group of the water molecule. This effect seems to be a property of polyfluorinated ethers.     

Kinetic and thermodynamic studies on nitrobenzylthioinosine binding to the nucleoside transporter of chinese hamster ovary cells

The binding of [G-³H]nitrobenzylthioinosine to intact Chinese hamster ovary cells has been studied kinetically and thermodynamically. The association of nitrobenzylthioinosine with cells is a second-order process which proceeds at 24°C with a rate constant of 2·10⁷ M^?1·s^?1. Dissociation of the complex was characterized as a simple first-order process with rate constant on the order of 7·10^?3 s^?1. The quotient of these is comparable to the dissociation constant as measured in equilibrium binding studies, 2.2·10^?10 M. The temperature dependence of the rate of association indicated an Arrhenius activation energy of 8.4 kcal·mol^?1, while that of the equilibrium constant for dissociation indicated a standard enthalpy change of 8.8 kcal·mol^?1. The large increase in affinity of nitrobenzylthioinosine as compared to natural nucleosides is attributable to an entropy-driven interaction with the binding site. Thymidine, dipyridamole and papaverine each decrease the apparent dissociation constant for the nitrobenzylthioinosine-cell complex; the latter, inhibitors of nucleoside transport, decrease the rate of dissociation of the complex.     

Base Pairing at the Stem-loop Junction in the SL1 Kissing Complex of HIV-1 RNA: A Thermodynamic Study Probed by Molecular Dynamics Simulation

Abstract

The thermodynamics of the opening/closure process of a GC base pair located at the stem-loop junction of the SL1 sequence from HIV-1_Lai genomic RNA was investigated in the context of a loop-loop homodimer (or kissing complex) using molecular dynamics simulation. The free energy, enthalpy and entropy changes for the closing reaction are 0 kcal·mol^?1, ?11 kcal·mol^?1and ?0.037 kcal·mol^?1-K^?1 at 300° K respectively. Furthermore it is found that the free energy change is the same for the formation of a 11 nucleotide loop closed with UG and for the formation of a 9 nucleotide loop closed with GC. Our study evidences the high flexibility of the nucleotides at the stem-loop junction explaining the weak stability of this structure.     

Conformational diversity of the quercetin molecule: a quantum-chemical view

Abstract

This paper focuses on the comprehensive conformational analysis of the quercetin molecule with a broad range of the therapeutic and biological actions. All possible conformers of these molecule, corresponding to the local minima on the potential energy hypersurface, have been obtained by the sequential rotation of its five hydroxyl groups and also by the rotation of its (A?+?C) and B rings relatively each other. Altogether, it was established 48 stable conformers, among which 24 conformers possess planar structure and 24 conformers – nonplanar structure. Their structural, symmetrical, energetical and polar characteristics have been investigated in details. Quantum-mechanical calculations indicate that conformers of the quercetin molecule are polar structures with a dipole moment, which varies within the range from 0.35 to 9.87 Debay for different conformers. Relative Gibbs free energies of these conformers are located within the range from 0.0 to 25.3?kcal·mol^?1 in vacuum under normal conditions. Impact of the continuum with ε?=?4 leads to the decreasing of the Gibbs free energies (–0.19–18.15?kcal·mol^?1) and increasing of the dipole moment (0.57–12.48?D). It was shown that conformers of the quercetin molecule differ from each other by the intramolecular specific contacts (two or three), stabilizing all possible conformers of the molecule – H-bonds (both classical ОН…О and so-called unusual С′Н…О and ОН…С′) and attractive van-der-Waals contacts О…О. Obtained conformational analysis for the quercetin molecule enables to provide deeper understanding of the ‘structure-function’ relationship and also to suggest its mechanisms of the therapeutic and biological actions.

Communicated by Ramaswamy H. Sarma     

Binding of NDGA and morin with phospholipase A2: experimental and computational evidences

The effects of morin and nordihydroguaiaretic acid (NDGA), two plant secondary metabolites, on porcine pancreatic phospholipase A₂ (PLA₂) were investigated by isothermal titration calorimetry (ITC) and in silico docking analyses. The binding energies obtained for NDGA and morin from the ITC studies are ? 6.36 and ? 5.91 kcal mol^? 1, respectively. Similarly, the glide scores obtained for NDGA and morin towards PLA₂ were ? 7.32 and ? 7.23 kcal mol^? 1, respectively. Further the docked complexes were subjected to MD simulation in the presence of explicit water molecules to check the binding stability of the ligands in the active site of PLA₂. The bound ligands make hydrogen bonds with the active site residues of the enzyme and coordinate bonds with catalytically important Ca²⁺ ion. The binding of ligands at the active site of PLA₂ may also contribute to the reported anti-inflammatory properties of NDGA and morin.     

The effect of temperature of the rate of photosynthetic electron transfer in chloroplasts of chilling-sensitive and chilling-resistant plants

1. Photochemical activities as a function of temperature have been compared in chloroplasts isolated from chilling-sensitive (below approximately 12 °C) and chilling-resistant plants.2. An Arrhenius plot of the photoreduction of NADP⁺ from water by chloroplasts isolated from tomato (Lycopersicon esculentum var. Gross Lisse), a chilling-sensitive plant, shows a change in slope at about 12 °C. Between 25 and 14 °C the activation energy for this reaction is 8.3 kcal·mole^?1. Between 11 and 3 °C the activation energy increases to 22 kcal·mole^?1. Photoreduction of NADP⁺ by chloroplasts from another chilling-sensitive plant, bean (Phaseolus vulgaris var. brown beauty), shows an increase in activation energy from 5.9 to 17.5 kcal·mole^?1 below about 12 °C.3. The photoreduction of NADP⁺ by chloroplasts isolated from two chilling-resistant plants, lettuce (Lactuca sativa var. winter lake) and pea (Pisum sativum var. greenfeast), shows constant activation energies of 5.4 and 8.0 kcal·mole^?1, respectively, over the temperature range 3–25 °C.4. The effect of temperature on photosynthetic electron transfer in the chloroplasts of chilling-sensitive plants is localized in Photosystem I region of photosynthesis. Both the photoreduction of NADP⁺ from reduced 2,6-dichlorophenol-indophenol and the ferredoxin-NADP⁺ reductase (EC 1.6.99.4) activity of choroplasts of chilling-sensitive plants show increases in activation energies at approximately 12 °C whereas Photosystem II activity of chloroplasts of chilling-sensitive plants shows a constant activation energy over the temperature range 3–25 °C. The photoreduction of Diquat (1,1′-ethylene-2,2′-dipyridylium dibromide) from water by bean chloroplasts, however, does not show a change in activation energy over the same temperature range. The activation energies of each of these reactions in chilling-resistant plants is constant between 3 and 25 °C.5. The effect of temperature on the activation energy of these reactions in chloroplasts from chilling-sensitive plants is reversible.6. In chilling-sensitive plants, the increased activation energies below approximately 12 °C, with consequent decreased rates of reaction for the photoreduction of NADP⁺, would result in impaired photosynthetic activity at chilling temperatures. This could explain the changes in chloroplast structure and function when chilling-sensitive plants are exposed to chilling temperatures.     

A theoretical investigation of the characteristics of hydrogen/halogen bonding interactions in dibromo-nitroaniline

In this work, computations of density functional theory (DFT) were carried out to investigate the nature of interactions in solid 2,6-dibromo-4-nitroaniline (DBNA). This system was selected to mimic the hydrogen/halogen bonding found within crystal structures as well as within biological molecules. DFT (M06-2X/6-311++G**) calculations indicated that the binding energies for different of interactions lie in the range between ?1.66 and ?9.77 kcal mol^?1. The quantum theory of atoms in molecules (QTAIM) was applied to provide more insight into the nature of these interactions. Symmetry-adapted perturbation theory (SAPT) analysis indicated that stability of the Br···Br halogen bonds is predicted to be attributable mainly to dispersion, while electrostatic forces, which have been widely believed to be responsible for these types of interactions, play a smaller role. Our results indicate that, for those nuclei participating in hydrogen/halogen bonding interactions, nuclear quadrupole resonance parameters exhibit considerable changes on going from the isolated molecule model to crystalline DBNA.

DPT tautomerization of the long A∙A* Watson-Crick base pair formed by the amino and imino tautomers of adenine: combined QM and QTAIM investigation

Combining quantum-mechanical (QM) calculations with quantum theory of atoms in molecules (QTAIM) and using the methodology of sweeps of the energetic, electron-topological, geometric and polar parameters, which describe the course of the tautomerization along the intrinsic reaction coordinate (IRC), we showed for the first time that the biologically important A?A* base pair (C_s symmetry) formed by the amino and imino tautomers of adenine (A) tautomerizes via asynchronous concerted double proton transfer (DPT) through a transition state (TS), which is the A⁺?A^? zwitterion with the separated charge, with C_s symmetry. The nine key points, which can be considered as electron-topological “fingerprints” of the asynchronous concerted A?A*?A*?A tautomerization process via the DPT, were detected and completely investigated along the IRC of the A?A*?A*?A tautomerization. Based on the sweeps of the H-bond energies, it was found that intermolecular antiparallel N6Н?N6 (7.01 kcal mol^?1) and N1H?N1 (6.88 kcal mol^?1) H-bonds are significantly cooperative and mutually reinforce each other. It was shown for the first time that the A?A*?A*?A tautomerization is assisted by the third C2H?HC2 dihydrogen bond (DHB), which, in contrast to the two others N6H?N6 and N1H?N1 H-bonds, exists within the IRC range from ?2.92 to 2.92 Å. The DHB cooperatively strengthens, reaching its maximum energy 0.42 kcal mol^?1 at IRC?=??0.52 Å and minimum energy 0.25 kcal mol^?1 at IRC?=??2.92 Å, and is accompanied by strengthening of the two other aforementioned classical H-bonds. We established that the C2H?HC2 DHB completely satisfies the electron-topological criteria for H-bonding, in particular Bader’s and all eight “two-molecule” Koch and Popelier’s criteria. The positive value of the Grunenberg’s compliance constant (5.203 Å/mdyn) at the TS_A?A*?A*?A proves that the C2H?HC2 DHB is a stabilizing interaction. NBO analysis predicts transfer of charge from σ(C2–H) bonding orbital to σ*(H–C2) anti-bonding orbital; at this point, the stabilization energy E⁽²⁾ is equal to 0.19 kcal mol^?1 at the TS_A?A*?A*?A.     

Conformational study of trinucleoside tetraphosphate d(pCpGpCp): Transition of right-handed form to left-handed form

Using the semiempirical potential functions, conformational energies of the model compounds DMP^?, d(pCp), d(pGp), and d(pCpGpCp) are calculated, and the B → Z transition is discussed along the pseudorotational path of the sugar ring. For dimethylmonophosphate anion, DMP^?, the energy contour map is presented and the stabilities of the phosphodiester conformations discussed. For the sugar ring without the base attached, the minimum energies for each sugar-puckering form are calculated along the pseudorotational path. The energy barrier of the interconversion between the C(3′)-endo form and the C(2′)-endo form is calculated to be about 2.0 kcal/mol. From the conformational energy calculations of the interconversions of mononucleoside diphosphates, d(pCp) and d(pGp), between the C(2′)-endo conformer and the C(3′)-endo conformer, the purine sugar segment is known to be more convertible than the pyrimidine sugar segment. The results also support the finding that the pseudorotational transition occurred with the O(1′)-endo form more easily than with the O(1′)-exo form. Based on the results of conformational studies of DMP^?, d(pCp), and d(pGp), a topological transition of the handedness of the model compound, d(pCpGpCp), is studied. The left-handed Z-form is found to be less stable by about 8.5 kcal/mol than is the right-handed B-form. The energy barrier of the Z → B transition is calculated to be about 17.4 kcal/mol. The contributions of the electrostatic and nonbonded energies to the energy barrier are discussed in connection with the conformation changes of the model compound, d(pCpGpCp).     

DFT comparison of the OH-initiated degradation mechanisms for five chlorophenoxy herbicides

To compare the OH-initiated reaction mechanisms of five chlorophenoxy herbicides, density functional theory (DFT) calculations of reactions in which ·OH attacks one of three active positions on each herbicide were carried out at the MPWB1K/6-311 + G(3df,2p)//MPWB1K/6-31 + G(d,p) level. For each herbicide, the calculation results show that ·OH addition to the C1 atom, which is the nexus between the benzene ring and the side group, possesses the lowest energy barrier among the three kinds of reactions, indicating that ·OH addition–substitution of the side chain is the most energetically and kinetically favorable reaction mechanism. Comparisons among the herbicides show that the mechanisms are affected by the steric hindrance and the electronegativities of the –CH₃ and –Cl groups. When comparing the addition of ·OH to the C1 site among the five herbicides, the activation energy for the reaction of ·OH with DCPP reaction is the lowest (3.61 kcal mol^?1), while that for the ·OH and 4-CPA reaction was the highest (5.91 kcal mol^?1). ·OH addition to the C4 site presents the highest energy barriers among the three kinds of reactions, indicating that the para Cl is difficult to break down. When comparing the H-atom abstraction reactions of the five herbicides, the H atoms in the –CH₂– group of 2,4-D are the easiest for ·OH to abstract, whereas those of DCPP and MCPP are more difficult to abstract, due to the steric hindrance of the –CH₃ group. Additionally, the results obtained from the PCM calculations reveal that most of the reactions occur more easily in water than in gas, though the mechanisms involved are the same as those discussed above.     

DFT investigation on the decarboxylation mechanism of <Emphasis Type="Italic">ortho</Emphasis> hydroxy benzoic acids with acid catalysis

A density functional theory (DFT) study was performed to explore the mechanisms of the acid-catalyzed decarboxylation reaction of salicylic acids using the B3LYP method with 6-31++G(d,p) basis set in both gas phase and aqueous environment. The α-protonated cation of carboxylate acid was formed during the decarboxylation process in acidic conditions, and the presence of hydrogen ions promotes decarboxylation greatly by significantly decreasing the overall reaction energy barriers to 20.98 kcal mol^?1 in gas phase and 20.93 kcal mol^?1 in water, respectively. The hydrogen in the α-carbon came directly from the acid rather than from the carboxyl group in neutral state. Compared with the reaction in gas phase, water in aqueous state causes the reaction to occur more easily. Substituents of methyl group, chlorine and fluorine at the ortho-position to the carboxyl of salicylic acid could further lower the decarboxylation energy barriers and facilitate the reaction.     

Hybrid MC/QC simulations of water-assisted proton transfer in nucleosides. Guanosine and its analog acyclovir

To provide an in-depth insight into the molecular basis of spontaneous tautomerism in DNA and RNA base pairs, a hybrid Monte Carlo (MC)–quantum chemical (QC) methodology is implemented to map two-dimensional potential energy surfaces along the reaction coordinates of solvent-assisted proton transfer processes in guanosine and its analog acyclovir in aqueous solution. The solvent effects were simulated by explicit inclusion of water molecules that model the relevant part of the first hydration shell around the solute. The position of these water molecules was estimated by carrying out a classical Metropolis Monte Carlo simulation of dilute water solutions of the guanosine (Gs) and acyclovir (ACV) and subsequently analyzing solute–solvent intermolecular interactions in the statistically-independent MC-generated configurations. The solvent-assisted proton transfer processes were further investigated using two different ab initio MP2 quantum chemical approaches. In the first one, potential energy surfaces of the ‘bare’ finite solute–solvent clusters containing Gs/ACV and four water molecules (MP2/6-31+G(d,p) level) were explored, while within the second approach, these clusters were embedded in ‘bulk’ solvent treated as polarizable continuum (C-PCM/MP2/6-31+G(d,p) level of theory). It was found that in the gas phase and in water solution, the most stable tautomer for guanosine and acyclovir is the 1H-2-amino-6-oxo form followed by the 2-amino-6-(sZ)-hydroxy form. The energy barriers of the water-assisted proton transfer reaction in guanosine and in acyclovir are found to be very similar – 11.74 kcal mol^?1 for guanosine and 11.16 kcal mol^?1 for acyclovir, and the respective rate constants (k = 1.5?×?10¹ s^?1, guanosine and k = 4.09?×?10¹ s^?1, acyclovir), are sufficiently large to generate the 2-amino-6-(sZ)-hydroxy tautomer. The analysis of the reaction profiles in both compounds shows that the proton transfer processes occur through the asynchronous concerted mechanism.     

Effect of stepwise microhydration on the methylammonium···phenol and ammonium···phenol interaction

A computational study has been performed for studying the characteristics of the interaction of phenol with ammonium and methylammonium cations. The effect of the presence of water molecules has also been considered by microhydrating the clusters with up to three water molecules. Clusters of phenol with ammonium and methylammonium cations present similar characteristics, though ammonium complexes have been found to be more stable than the methylammonium ones. The first water molecule included in the complexes interacts with a N-H group of ammoniun cations and simultaneously with the hydroxyl oxygen atom of phenol (or the aromatic ring). This first water molecule is more tightly bound in the complex, so the stability gain as more water molecules are included drops significantly by 2-3 kcal?mol^?1 with respect to the first one. As more water molecules are included, the differences between favorable coordination sites (the cation, the hydroxyl group or a previous water molecule) decrease. As a consequence, several of the most stable complexes located including three water molecules already exhibit hydrogen bonds between the hydroxyl group and one water molecule. The results indicate that a cyclic pattern formed by a series of hydrogen bonds: π···H-N-H···O-H···O-?, is characteristic of the most stable minima, being kept as more water molecules are included in the system. Therefore, this pattern can be expected to be crucial in ammonium cations···phenol interaction if exposed to the solvent to any degree.     

Molecular dynamics simulations of void defects in the energetic material HMX

A molecular dynamics (MD) simulation was carried out to characterize the dynamic evolution of void defects in crystalline octahydro-1, 3, 5, 7-tetranitro-1, 3, 5, 7-tetrazocine (HMX). Different models were constructed with the same concentration of vacancies (10 %) to discuss the size effects of void. Energetic ground state properties were determined by annealing simulations. The void formation energy per molecule removed was found to be 55–63 kcal/mol^?1, and the average binding energy per molecule was between 32 and 34 kcal/mol^?1 according to the change in void size. Voids with larger size had lower formation energy. Local binding energies for molecules directly on the void surface decreased greatly compared to those in defect-free lattice, and then gradually increased until the distance away from the void surface was around 10 Å. Analysis of 1 ns MD simulations revealed that the larger the void size, the easier is void collapse. Mean square displacements (MSDs) showed that HMX molecules that had collapsed into void present liquid structure characteristics. Four unique low-energy conformers were found for HMX molecules in void: two whose conformational geometries corresponded closely to those found in HMX polymorphs and two, additional, lower energy conformers that were not seen in the crystalline phases. The ratio of different conformers changed with the simulated temperature, in that the ratio of α conformer increased with the increase in temperature.     

Theoretical studies of the interaction between enflurane and water

Increase of the atmospheric concentration of halogenated organic compounds is partially responsible for a change of the global climate. In this work we have investigated the interaction between halogenated ether and water, which is one of the most important constituent of the atmosphere. The structures of the complexes formed by the two most stable conformers of enflurane (a volatile anaesthetic) with one and two water molecules were calculated by means of the counterpoise CP-corrected gradient optimization at the MP2/6–311++G(d,p) level. In these complexes the CH…O_w hydrogen bonds are formed, with the H…O_w distances varying between 2.23 and 2.32 Å. A small contraction of the CH bonds and the blue shifts of the ν(CH) stretching vibrations are predicted. There is also a weak interaction between one of the F atoms and the H atom of water, with the H_w…F distances between 2.41 and 2.87 Å. The CCSD(T)/CBS calculated stabilization energies in these complexes are between ?5.89 and ?4.66 kcal?mol^?1, while the enthalpies of formation are between ?4.35 and ?3.22 kcal?mol^?1. The Cl halogen bonding between enflurane and water has been found in two complexes. The intermolecular (Cl···O) distance is smaller than the sum of the corresponding van der Waals radii. The CCSD(T)/CBS stabilization energies for these complexes are about ?2 kcal?mol^?1.

Insight into the process of product expulsion in cellobiohydrolase Cel6A from Trichoderma reesei by computational modeling

Glycoside hydrolase cellulase family 6 from Trichoderma reesei (TrCel6A) is an important cellobiohydrolase to hydrolyze cellooligosaccharide into cellobiose. The knowledge of enzymatic mechanisms is critical for improving the conversion efficiency of cellulose into ethanol or other chemicals. However, the process of product expulsion, a key component of enzymatic depolymerization, from TrCel6A has not yet been described in detail. Here, conventional molecular dynamics and steered molecular dynamics (SMD) were applied to study product expulsion from TrCel6A. Tyr103 may be a crucial residue in product expulsion given that it exhibits two different posthydrolytic conformations. In one conformation, Tyr103 rotates to open the ?3 subsite. However, Tyr103 does not rotate in the other conformation. Three different routes for product expulsion were proposed on the basis of the two different conformations. The total energy barriers of the three routes were calculated through SMD simulations. The total energy barrier of product expulsion through Route 1, in which Tyr103 does not rotate, was 22.2 kcal·mol^?1. The total energy barriers of product expulsion through Routes 2 and 3, in which Tyr103 rotates to open the ?3 subsite, were 10.3 and 14.4 kcal·mol^?1, respectively. Therefore, Routes 2 and 3 have lower energy barriers than Route 1, and Route 2 is the thermodynamically optimal route for product expulsion. Consequently, the rotation of Tyr103 may be crucial for product release from TrCel6A. Results of this work have potential applications in cellulase engineering.     

A DFT study of adsorption and decomposition of hexahydro-1,3,5-trinitro-1,3,5-triazine on Mg(0001) surface

The adsorption and decomposition of hexogen (RDX) molecule on the Mg(0001) surface were investigated by the generalized gradient approximation (GGA) of density functional theory (DFT). The calculations employed a supercell (4?×?4?×?4) slab model and three-dimensional periodic boundary conditions. The strong attractive forces between RDX molecule and magnesium atoms induce the RDX’s N???O bond breaking. Subsequently, the dissociated oxygen atoms and radical fragment of RDX oxidize the Mg surface. The largest adsorption energy is ?2104.0 kJ mol^-1. We also investigated the decomposition mechanism of RDX molecule on the Mg(0001) surface. The activation energy for the dissociation step of configuration V4 is as small as 2.5 kJ mol^-1, while activation energies of other configurations are much larger, in the range of 964.9–1375.1 kJ mol^-1. Mg powder is more active than Al powder, and Mg powder performs better in increasing the combustion exothermicity of RDX as well.     

Interfacial adsorption of an inhalation anesthetic onto ionic surfactant micelles and its desorption by high pressure

The effects of pressure and temperature on the critical micelle concentration (CMC) of sodium dodecylsulfate (SDS) were measured in the presence of various concentrations of an inhalation anesthetic, methoxyflurane. The change in the partial molal volume of SDS on micellization, $\text{Δ}\text{V}{}_{\mathrm{<mtext>m</mtext>}}$, increased with the increase in the concentration of methoxyflurane. The CMC-decreasing power, which is defined as the slope of the linear plot between ln(CMC) vs. mole fraction of anesthetic, was determined as a function of pressure and temperature. Since the CMC-decreasing power is correlated to the micelle/water partition coefficient of anesthetic, the volume change of the transfer (ΔV_p^o) of methoxyflurane from water to the micelle can be determined from the pressure dependence of the CMC-decreasing power. The value of ΔV_p^o amounts 6.5±1.8 cm³·mol^?1, which is in reasonable agreement with the volume change determined directly from the density data, 5.5±0.6 cm³ · mol^?1. Under the convention of thermodynamics, this indicates that the application of pressure squeezes out anesthetic molecules from the micelle. The transfer enthalpy of anesthetic from water to the micelle is slightly endothermic. The partial molal volume of methoxyflurane in the micelle (112.0 cm³·mol^?1) is smaller than that in decane (120.5 cm³·mol^?1) and is larger than that in water (108.0 cm³·mol^?1). This indicates that the anesthetic molecules are incorporated into the micellar surface region, i.e., the palisade layer of the micelle in contact with water molecules, rather than into the micelle core.