Structural Basis and Mechanism of Chiral Benzedrine Molecules Interacting With Third Dopamine Receptor

In order to investigate the chiral benzedrine molecules corresponding to their different characteristics in biochemical systems, we studied their interaction with D₃R using the docking method, molecular dynamic simulation, and quantum chemistry. The obtained results indicate that the active residues for R‐benzedrine (RAT) bound with D₃R are Ala132, Asp133, and Tyr55, while Asn57, Asp133, Asp168, Cys172, Gly54, Trp24, and Vall136 act as the active residues for S‐benzedrine (SAT). The different active pockets are observed for ART or SAT because they possess different active residues. The binding energies between RAT and SAT with D₃R were determined to be ?44.0 kJ.mol^?1 and ?71.2 kJ.mol^?1, respectively. These results demonstrate that SAT within the studied pocket of D₃R has a stronger capability of binding with D₃R, while it is more feasible for RAT to leave from the interior positions of D₃R. In addition, the results suggest that the D₃R protein can recognize chiral benzedrine molecules and influence their different addictive and pharmacological effects in biochemical systems. Chirality 28:674–685, 2016. © 2016 Wiley Periodicals, Inc.     

Kinetics investigation of the hydrogen abstraction reaction between CH<Subscript>3</Subscript>SS and CN radicals

The reaction mechanisms and rates for the H abstraction reactions between CH₃SS and CN radicals in the gas phase were investigated with density functional theory (DFT) methods. The geometries, harmonic vibrational frequencies, and energies of all stationary points were obtained at B3PW91/6-311G(d,p) level of theory. Relationships between the reactants, intermediates, transition states and products were confirmed, with the frequency and the intrinsic reaction coordinate (IRC) analysis at the same theoretical level. High accurate energy information was provided by the G3(MP2) method combined with the standard statistical thermodynamics. Gibbs free energies at 298.15 K for all of the reaction steps were reported, and were used to describe the profile diagrams of the potential energy surface. The rate constants were evaluated with both the classical transition state theory and the canonical variational transition state theory, in which the small-curvature tunneling correction was included. A total number of 9 intermediates (IMs) and 17 transition states (TSs) were obtained. It is shown that IM1 is the most stable intermediate by the largest energy release, and the channel of CH₃SS?+?CN?→?IM3?→?TS10?→?P1(CH₂SS?+?HCN) is the dominant reaction with the lowest energy barrier of 144.7 kJ mol^?1. The fitted Arrhenius expressions of the calculated CVT/SCT rate constants for the rate-determining step of the favorable channel is k =7.73?×?10⁶? T ^1.40exp(?14,423.8/T) s^?1 in the temperature range of 200–2000 K. The apparent activation energy E _a(app.) for the main channel is ?102.5 kJ mol^?1, which is comparable with the G3(MP2) energy barrier of ?91.8 kJ mol^?1 of TS10 (relative to the reactants).     

Phase-transition properties of glycerol-monopalmitate lipid bilayers investigated by molecular dynamics simulation: influence of the system size and force-field parameters

Phase-transition properties of glycerol-1-monopalmitate (GMP) bilayers are investigated using explicit-solvent molecular dynamics (MD) simulations, initiated from structures appropriate for the gel (GL) or liquid crystal (LC) phases, and carried out at different hydration levels and temperatures. Building up on a previous study and based on 600 ns simulations, the influence of the system size and of the force field on the equilibrium thermodynamic and dynamic parameters of the bilayers in the GL and LC phases, as well as on the temperature T_m and properties of the GL ? LC phase transition, are analysed. Qualitatively speaking, the results agree with the available experimental data for the area per lipid in the two phases and for the phase-transition temperatures at the three hydration levels irrespective of the selected model parameters. They also suggest that the total number of hydrogen bonds formed between a lipid headgroup and its environment is essentially constant, amounting to about four in both the LC and the GL phases. Quantitatively speaking, the dependence of T_m on the hydration level is found to be non-systematic across the different combinations of model parameters. This results in part from a sensitivity of the results on the system size and force-field parameters but also from the limited accuracy of the bracketing approach employed here to estimate T_m. Finally, a simple kinetic model is proposed to account for the timescales of the transitions. This model involves enthalpy and entropy increases of about 26 kJ mol^? 1 and 83 J mol^? 1 K^? 1 per lipid, upon going from the GL to the LC phase. The transition state is associated with activation parameters corresponding to 13% and 11%, respectively, of these values along the GL → LC transition, resulting in an activation free energy of about 0.3 kJ mol^? 1 per lipid at T_m.     

The addition reactions of germylenoid H<Subscript>2</Subscript>GeAlCl<Subscript>3</Subscript> with ethylene: a theoretical investigation

Theoretical calculations using the M062X and QCISD methods were performed on the addition reactions of the aluminum germylenoid H₂GeAlCl₃ with ethylene. The most two stable structures of germylenoid H₂GeAlCl₃, i.e., the p-complex and three-membered ring structures, respectively, were employed as reactants. The calculated results indicate that, for the p-complex, H₂GeAlCl₃ there are two pathways, I and II, of which path I involves just one transition state, while path II involves two transition states between reactants and products. Comparing the reaction barrier heights of path I (44.6 kJ mol^?1) and II (37.6 kJ mol^?1), the two pathways are competitive, with similar barriers under the same conditions, while for the three-membered ring structure, another two pathways, III and IV, also exist. Path III has one transition state; however, in path IV, two transition states exist. By comparing their barrier heights, path III (barrier height 39.2 kJ mol^?1) could occur more easily than path IV (barrier height 92.8 kJ mol^?1). Considering solvent effects on these addition reactions, the PCM model and CH₂Cl₂ solvent were used in calculations, and the calculated results demonstrate that CH₂Cl₂ solvent is unfavorable for the reactions, except for path II. In CH₂Cl₂ solvent, paths II and III are more favorable than paths I and IV.     

Theoretical study of chemisorption of cyanuric fluoride and S-triazine on the surface of Al-doped graphene

We studied the adsorption of cyanuric fluoride (CF) and s-triazine (ST) molecules on the surface of pristine as well as Al-doped graphenes using density functional theory calculations. Our results reveal low adsorption on the surface of pristine graphene; but by modification of surface using aluminium, resulted Al-doped graphene becomes more reactive towards both CF and ST molecules. We aimed to focus on the adsorption energy, electronic structure, charge analysis, density of state and global indices of each system upon adsorption of CF and ST molecules on the above-mentioned surfaces. Our calculated adsorption energies for the most stable position configurations of CF and ST on Al-doped graphene were ?76.53 kJ mol^?1 (?57.45 kJ mol^?1 BSSE corrected energy) and ?115.55 kJ mol^?1 (?86.87 kJ mol^?1 BSSE corrected energy), respectively, which point to the chemisorption process. For each CF and ST molecule, upon adsorption on the surface of Al-doped graphene, the band gap of HOMO-LUMO was reduced considerably and it becomes a p-type semiconductor, whereas there is no hybridisation between the above-mentioned molecules and pristine graphene.     

Homologous yeast lipases/acyltransferases exhibit remarkable cold-active properties

Lipases/acyltransferases catalyse acyltransfer to various nucleophiles preferentially to hydrolysis even in aqueous media with high thermodynamic activity of water (a _w >0.9). Characterization of hydrolysis and acyltransfer activities in a large range of temperature (5 to 80 °C) of secreted recombinant homologous lipases of the Pseudozyma antarctica lipase A superfamily (CaLA) expressed in Pichia pastoris, enlighten the exceptional cold-activity of two remarkable lipases/acyltransferases: CpLIP2 from Candida parapsilosis and CtroL4 from Candida tropicalis. The activation energy of the reactions catalysed by CpLIP2 and CtroL4 was 18–23 kJ mol^?1 for hydrolysis and less than 15 kJ mol^?1 for transesterification between 5 and 35 °C, while it was respectively 43 and 47 kJ mol^?1 with the thermostable CaLA. A remarkable consequence is the high rate of the reactions catalysed by CpLIP2 and CtroL4 at very low temperatures, with CpLIP2 displaying at 5 °C 65 % of its alcoholysis activity and 45 % of its hydrolysis activity at 30 °C. These results suggest that, within the CaLA superfamily and its homologous subgroups, common structural determinants might allow both acyltransfer and cold-active properties. Such biocatalysts are of great interest for the efficient synthesis or functionalization of temperature-sensitive lipid derivatives, or more generally to lessen the environmental impact of biocatalytic processes.     

The role of CS<Subscript>2</Subscript> in CS<Subscript>2</Subscript>/NMP mixed solvent in weakening the hydrogen bond of OH–N in coal: a DFT investigation

The interaction processes of trace amounts of N-methyl-2-pyrrolidinone (NMP), CS₂/NMP (1:1 by volume) and pure NMP solvent with the hydrogen bond of OH?N in coal were constructed and simulated by density functional theory methods. The distances and bond orders between the main related atoms, and the hydrogen bond energy of OH?N were calculated. The calculated results show that pure NMP solvent does not weaken the hydrogen bond of OH?N in coal. However, trace amounts of NMP and CS₂/NMP (1:1 by volume) have a strong capacity to weaken the hydrogen bond of OH?N in coal. The H2–N3 distances are elongated from 1.87 Å to 3.80 Å and 3.44 Å, the bond orders of H2–N3 all disappear, and the corresponding hydrogen bond energies of OH?N in coal decrease from 45.72 kJ mol^?1 to 7.06 and 11.24 kJ mol^?1, respectively. These results show that CS₂ added to pure NMP solvent plays an important role in releasing the original capacity of NMP to weaken the hydrogen bond of OH?N in coal, in agreement with experimental observations.     

1,3‐Propanediol adsorption on a cation exchange resin: Adsorption isotherm,thermodynamics, and mechanistic studies

1,3‐Propanediol (1,3‐PD) is widely used in cosmetics, foods, antifreezes, and polyester. A low‐cost cation exchange resin, 001×7 H‐form resin, was selected for 1,3‐PD adsorption obtained from microbial fermentation of crude glycerol. The thermodynamics and kinetics of adsorption were studied. To identify the characteristics of the adsorption process, the influence of 1,3‐PD concentration, temperature, and resin particle diameters was studied. The temperature dependence of the adsorption equilibrium in the range of 288 to 318 K was observed to satisfy the Langmuir isotherm well. The thermodynamic parameters, adsorption enthalpy, entropy, and Gibbs free energy, were determined as 36.2 kJ·moL^?1, 160 mol^?1·K^?1, and ?11.2 kJ·moL^?1, respectively, which indicated that the adsorption was spontaneous and endothermic. The adsorption kinetics was accurately represented by the shell progressive model and indicated that the particle diffusion was the rate‐limiting step. Based on the kinetic simulation, the rate constant of exchange (k₀), order reaction (α), and the apparent activation energy reaction (E_a) were obtained as 3.11×10^?3, 0.644, and 11.5 kJ·moL^?1, respectively. This kinetic and thermodynamic analysis of 1,3‐PD recovery presented in this article is also applicable for the separation of other polyols by resin adsorption, which will promote value‐added utilization of glycerol.     

Concentration effect of cimetidine with POPC bilayer: a molecular dynamics simulation study

All-atom molecular dynamics simulations have been performed on cimetidine in the presence of a palmitoyloleoylphosphatidylcholine (POPC) bilayer. The free energy profile of a single cimetidine molecule passing across POPC bilayer displays a minimum at the interface of bilayer and water. Ten cimetidine molecules were inserted into POPC bilayer to obtain an 8 mol % drug model, and molecular dynamics results showed that cimetidine molecules reside at the polar region of POPC bilayer with sulphur atoms directing to the hydrophobic region. By comparing the one drug model with 8 mol % drug model, one can see that the central barrier to cross the membrane increases while the free energy in bulk water decreases, indicating that the ability of cimetidine passing across the POPC bilayer weakens at increased concentration. In addition, the free energy minimum shifts closer to the hydrophobic core. Our results indicate that with the increased drug concentration, it is more difficult for cimetidine to enter and pass across POPC bilayer.     

The dynamic motion of a M (M = Ca, Yb) atom inside the C74 (D 3h) cage: a relativistic DFT study

The interaction between M (M?=?Ca, Yb) atom and C₇₄ (D _3h) has been investigated by all electron relativistic density function theory. With the aid of the representative patch of C₇₄ (D _3h), we studied the interaction between C₇₄ (D _3h) and M (M?=?Ca, Yb) atom and obtained the interaction potential. Optimized structures show that there are three equivalent stable isomers and there is one transition state between every two stable isomers. According to the minimum energy pathway, the possible movement trajectory of M (M?=?Ca, Yb) atom in the C₇₄ (D _3h) cage is explored. The calculated energy barrier for Yb atoms moving from the stable isomer to the transition state is 10.4 kcal mol^?1 and the energy barrier for Ca atoms is 6.1 kcal mol^?1. The calculated NMR spectra of M@C₇₄ (M?=?Ca, Yb) are in good agreement with the experimental data. There are nine lines in the spectra: one 1/6 intensity signal, four half intensity signals and four full intensity signals.     

Thermodynamic Prediction of Glycine Polymerization as a Function of Temperature and pH Consistent with Experimentally Obtained Results

Prediction of the thermodynamic behaviors of biomolecules at high temperature and pressure is fundamental to understanding the role of hydrothermal systems in the origin and evolution of life on the primitive Earth. However, available thermodynamic dataset for amino acids, essential components for life, cannot represent experimentally observed polymerization behaviors of amino acids accurately under hydrothermal conditions. This report presents the thermodynamic data and the revised HKF parameters for the simplest amino acid “Gly” and its polymers (GlyGly, GlyGlyGly and DKP) based on experimental thermodynamic data from the literature. Values for the ionization states of Gly (Gly⁺ and Gly^?) and Gly peptides (GlyGly⁺, GlyGly^?, GlyGlyGly⁺, and GlyGlyGly^?) were also retrieved from reported experimental data by combining group additivity algorithms. The obtained dataset enables prediction of the polymerization behavior of Gly as a function of temperature and pH, consistent with experimentally obtained results in the literature. The revised thermodynamic data for zwitterionic Gly, GlyGly, and DKP were also used to estimate the energetics of amino acid polymerization into proteins. Results show that the Gibbs energy necessary to synthesize a mole of peptide bond is more than 10 kJ mol^?1 less than previously estimated over widely various temperatures (e.g., 28.3 kJ mol^?1 → 17.1 kJ mol^?1 at 25 °C and 1 bar). Protein synthesis under abiotic conditions might therefore be more feasible than earlier studies have shown.     

Conformational transitions of a dipeptide in water: Effects of imposed pathways using umbrella sampling techniques

The free energy difference between two states of a molecular system separated by an energy barrier can generally be computed using the technique of umbrella sampling along a chosen reaction coordinate or pathway. The effect of a particular choice of pathway upon the obtained free energy difference is investigated by molecular dynamics simulation of a model system consisting of a glycine dipeptide in aqueous solution. Two different reaction coordinates connecting the so-called C₅ and C₇ conformations, one involving intramolecular hydrogen bonds and the other involving the peptide ?, ψ angles, are considered. The Gibbs free energy differences ΔG(C₅ – C₇) are small in both cases, 1.5 ± 1 kJ mol^?1 and 2.2 ± 1 kJ mol ^?1, respectively. The two different reaction coordinates yield free energy differences that are identical to within their statistical error. It is found that the exchange of solute–solute, solute–water, and water–water hydrogen bonds involves free energy changes of less than k_BT, which points at the existence of a multitutde of low free energy pathways connecting the C₅ and C₇ dipeptide conformations. © 1994 John Wiley & Sons, Inc.     

Spectroscopic Studies on the Interaction of Fluorine Containing Triazole with Bovine Serum Albumin

The binding of one fluorine including triazole (C₁₀H₉FN₄S, FTZ) to bovine serum albumin (BSA) was studied by spectroscopic techniques including fluorescence spectroscopy, UV–Vis absorption, and circular dichroism (CD) spectroscopy under simulative physiological conditions. Fluorescence data revealed that the fluorescence quenching of BSA by FTZ was the result of forming a complex of BSA–FTZ, and the binding constants (K _a) at three different temperatures (298, 304, and 310 K) were 1.516?×?10⁴, 1.627?×?10⁴, and 1.711?×?10⁴?mol L^?1, respectively, according to the modified Stern–Volmer equation. The thermodynamic parameters ΔH and ΔS were estimated to be 7.752 kJ mol^?1 and 125.217 J?mol^?1?K^?1, respectively, indicating that hydrophobic interaction played a major role in stabilizing the BSA–FTZ complex. It was observed that site I was the main binding site for FTZ to BSA from the competitive experiments. The distance r between donor (BSA) and acceptor (FTZ) was calculated to be 7.42 nm based on the Förster theory of non-radioactive energy transfer. Furthermore, the analysis of fluorescence data and CD data revealed that the conformation of BSA changed upon the interaction with FTZ.     

Gene cloning and characterization of recombinant L-Asparaginase from <Emphasis Type="Italic">Bacillus subtilis</Emphasis> strain R5

L-asparaginase gene from Bacillus subtilis strain R5 (Asn-R5), comprising 990 nucleotides corresponding to a polypeptide of 329 amino acids, was cloned and expressed in Escherichia coli. Recombinant Asn-R5 was produced in soluble and active form exhibiting a specific activity of 223 μmol min^?1 mg^?1. The optimal temperature and pH for L-asparaginase activity of Asn-R5 were 35 °C and 9.0, respectively. Asn-R5 displayed a 50% activity with D-asparagine and 2% with L-glutamine compared to 100% with L-asparagine. No activity could be detected when D-glutamine was used as substrate. Half-life of the enzyme was 180 min at 35 °C and 40 min at 50 °C. There was no effect of metal ions and EDTA on the activity indicating that Asn-R5 enzyme activity is not metal ion dependent. The K_m and V_max values were 2.4 mM and 265 μmol min^?1 mg^?1, respectively. Activation energy for reaction catalyzed by Asn-R5 was 28 kJ mol^?1. High L-asparaginase activity and thermostability of recombinant Asn-R5 may be beneficial for industrial production and application.     

Kinetics and thermodynamics of catalysis and thermal inactivation of a novel α-amylase (Tp-AmyS) from Thermotoga petrophila

Abstract

This research is focussed on kinetic, thermodynamic and thermal inactivation of a novel thermostable recombinant α-amylase (Tp-AmyS) from Thermotoga petrophila. The amylase gene was cloned in pHIS-parallel1 expression vector and overexpressed in Escherichia coli. The steady-state kinetic parameters (V_max, K_m, k_cat and k_cat/K_m) for the hydrolysis of amylose (1.39?mg/min, 0.57?mg, 148.6?s^?1, 260.7), amylopectin (2.3?mg/min, 1.09?mg, 247.1?s^?1, 226.7), soluble starch (2.67?mg/min, 2.98?mg, 284.2?s^?1, 95.4) and raw starch (2.1?mg/min, 3.6?mg, 224.7?s^?1, 61.9) were determined. The activation energy (E_a), free energy (ΔG), enthalpy (ΔH) and entropy of activation (ΔS) at 98?°C were 42.9?kJ mol^?1, 74?kJ mol^?1, 39.9?kJ mol^?1 and ?92.3 J mol^?1 K^?1, respectively, for soluble starch hydrolysis. While ΔG of substrate binding (ΔG_E-S) and ΔG of transition state binding (ΔG_E-T) were 3.38 and ?14.1?kJ mol^?1, respectively. Whereas, EaD, Gibbs free energy (ΔG*), increase in the enthalpy (ΔH*) and activation entropy (ΔS*) for activation of the unfolding of transition state were 108, 107, 105?kJ mol^?1 and ?4.1 J mol^?1 K^?1. The thermodynamics of irreversible thermal inactivation of Tp-AmyS revealed that at high temperature the process involves the aggregation of the protein.     

Permeability of lipid bilayer to anthracycline derivatives. Role of the bilayer composition and of the temperature

The uptake of three anthracycline derivatives: doxorubicin, daunorubicin and pirarubicin, into large unilamellar vesicles (LUV) in response to a driving force provided by DNA encapsulated inside the LUV has been investigated as a function of the temperature and of the bilayers lipid composition. The kinetics of the decay of the anthracycline fluorescence in the presence of DNA-containing liposome was used to follow the diffusion of the drug through the membrane. For the three drugs, the permeability coefficient of the neutral form of the drug (P⁰) decreases as the amount of negatively charged phospholipid in the bilayers increases. This can be explained by the fact that the kinetics of passive diffusion of the drugs depends on the amount of neutral form embedded in the polar head group region, which decreases as the quantity of negatively charged phospholipids increases. P⁰ also decreases as the amount of cholesterol, that makes the bilayer more rigid, increases. The activation energies, E_a, for the passage of the neutral form of these anthracyclines through the bilayers lie within 100±15 kJ·mol⁻¹, except for pirarubicin and doxorubicin through anionic phospholipid-rich membranes (E_a=57 kJ·mol⁻¹) and cholesterol-rich membranes (E_a=167 kJ·mol⁻¹).     

Oxidation reactivity of zinc–cysteine clusters in metallothionein

Evaluating the reactivity of the metal–thiolate clusters in metallothionein (MT) is a key step in understanding the biological functions of this protein. The effects of the metal clustering and protein environment on the thiolate reactivity with hydrogen peroxide (H₂O₂) were investigated by performing quantum theory calculations with chemical accuracy at two levels of complexity. At the first level, the reactivity with H₂O₂ of a model system ([(Zn)₃(MeS)₉]^3?, MeS is methanethiolate) of the β domain cluster of MT was evaluated using density functional theory (DFT) with the mPW1PW91 functional. At the second level of complexity, the protein environment was included in the reactant system and the calculations were performed with the hybrid ONIOM method combining the DFT–mPW1PW91 and the semiempirical PM6 levels of theory. In these conditions, the energy barrier for the oxidation of the most reactive terminal thiolate was 21.5 kcal mol^?1. This is 3 kcal mol^?1 higher than that calculated for the terminal thiolate in the model system [(Zn)₃(MeS)₉]^3? and about 7 kcal mol^?1 higher than that obtained for the free thiolate. In spite of this rise of the energy barrier induced by the protein environment, the thiolate oxidation by H₂O₂ is confirmed as a possible way for metal release from MT. On the other hand, the results suggest that the antioxidant role of MT in the living cell cannot be as important as that of glutathione (which bears a free thiol).     

Thermal effects on motion and orientation of the cholestane spin label in planar multibilayers of dimyristoylphosphatidylcholine and dipalmitoylphosphatidylcholine

A detailed picture of the orientation and restricted motion of the cholestane spin label (3-spiro-doxyl-5α-cholestane) in planar multibilayers of dipalmitoylphosphatidylcholine and dimyristoylphosphatidylcholine has been recorded by simultaneous simulation of ESR spectra obtained with the magnetic field parallel and perpendicular to the bilayers (Shimoyama, Y., Eriksson, L.E.G. and Ehrenberg, A. (1978) Biochim. Biophys. Acta 508, 213–235). The analysis has been made over the temperature range ?30°C to 60°C on samples containing 20 to 22% water. At low temperatures the cholestane spin label is tilted with respect to the lipid bilayer normal by an angle of approx. 30° which disappears at the pretransition. In this low temperature range the restricted twisting motion has an activation energy of 5.5 kJ·mol^?1. Above the main transition the twisting motion is unrestricted and has the activation energy 20 kJ·mol^?1. From below the pretransition to above the main transition the velocity of the twisting motion increases by an order of magnitude. The amplitude of the wobbling motion increases abruptly from 0° to 35° at the main transition.