Theoretical investigation of the structures and properties of CL-20/DNB cocrystal and associated PBXs by molecular dynamics simulation

In this work, a CL-20/DNB cocrystal explosive model was established and six different kinds of fluoropolymers, i.e., PVDF, PCTFE, F₂₃₁₁, F₂₃₁₂, F₂₃₁₃ and F₂₃₁₄ were added into the (1 0 0), (0 1 0), (0 0 1) crystal orientations to obtain the corresponding polymer bonded explosives (PBXs). The influence of fluoropolymers on PBX properties (energetic property, stability and mechanical properties) was investigated and evaluated using molecular dynamics (MD) methods. The results reveal a decrease in engineering moduli, an increase in Cauchy pressure (i.e., rigidity and stiffness is lessened), and an increase in plastic properties and ductility, thus indicating that the fluoropolymers have a beneficial influence on the mechanical properties of PBXs. Of all the PBXs models tested, the mechanical properties of CL-20/DNB/F₂₃₁₁ were the best. Binding energies show that CL-20/DNB/F₂₃₁₁ has the highest intermolecular interaction energy and best compatibility and stability. Therefore, F₂₃₁₁ is the most suitable fluoropolymer for PBXs. The mechanical properties and binding energies of the three crystal orientations vary in the order (0 1 0)?>?(0 0 1)?>?(1 0 0), i.e., the mechanical properties of the (0 1 0) crystal orientation are best, and this is the most stable crystal orientation. Detonation performance results show that the density and detonation parameters of PBXs are lower than those of the pure CL-20 and CL-20/DNB cocrystal explosive. The power and energetic performance of PBXs are thus weakened; however, these PBXs still have excellent detonation performance and are very promising. The results and conclusions provide some helpful guidance and novel instructions for the design and manufacture of PBXs.     

Differences in the binding of aromatic substrates to horseradish peroxidase revealed by fluorescence line narrowing

The heme in horseradish peroxidase (HRP) isoenzyme C was replaced by mesoporphyrin (MP), and the binding effect of the aromatic substrates benzo-and naphthohydroxamic acid (BHA, NHA), resorcinol (RE), isomeric resorcylic acids (alpha-, beta-, gamma-RE), and hydroquinone (HQ) was studied at pH 5 by conventional and laser-excited fluorescence spectroscopy on the basis of the signal of the porphyrin. Under laser excitation at cryogenic temperatures site selection was demonstrated, and the fluorescence line narrowing data were used to characterize the HRP/substrate complexes by the inhomogeneous distribution function for the S0----S1 (0----0) transition energy and the vibrational energies in the S1 electronic state. A comparison with ground-state vibrational energies for MP in chloroform/ether showed a downward shift in vibrational energies for S1 by approximately 20 cm-1. The association characteristics of the substrates were in accordance with previous literature data indicating NHA to be of the strongest binding affinity. For BHA, spectral evidence was obtained for a second type of binding site where hydrophobic interactions with the porphyrin ring may be possible. The effect of the RE's was similar to each other, but only beta-RE showed saturation. Complexation in every case caused the strong reduction of the splitting in the 0----0 transition energy for the tautomeric forms of MP and an increase in the 0----0 energy by 100-200 cm-1 depending on the substrate. The substrate binding also affected the phonon coupling of vibronic transitions exciting into the delta v = 927- and 976-cm-1 modes; in the latter case, the vibrational energy was also increased to 983 cm-1 for beta-RE. In the same energy range, however, the transition into the delta nu = 958-960-cm-1 mode was not affected by binding. Both the magnitude of the energy shifting and the change in the strength of phonon coupling gave the same relation, BHA less than NHA less than HQ less than RE's, indicating a common conformational origin. A reduction of the fluctuational freedom of the protein chain at room temperature within the heme pocket was suggested on the basis of the reduction of the width of the inhomogeneous distribution of 0----0 energies (from 60-70 to approximately 30 cm-1 in case of HRP/HQ) upon substrate binding. Ways to relate the transition energy splitting and shifting effects to conformational changes are discussed by invoking the Jahn-Teller effect.     

Quantum and molecular dynamics study for binding of macrocyclic inhibitors to human alpha-thrombin

Molecular dynamics simulations followed by quantum mechanical calculation and Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) analysis have been carried out to study binding of proline- and pyrazinone-based macrocyclic inhibitors (L86 and T76) to human alpha-thrombin. Detailed binding interaction energies between these inhibitors and individual protein fragments are calculated using DFT method based on a new quantum mechanical approach for computing protein-ligand interaction energy. The analysis of detailed interaction energies provides insight on the protein-ligand binding mechanism. Study shows that T76 and L86 bind to thrombin in a very similar "inhibition mode" except that T76 has relatively weaker binding interaction with Glu(217). The analysis from quantum calculation of binding interaction is consistent with the MM-PBSA calculation of binding free energy, and the calculated free energies for L86/T76-thrombin binding agree well with the experimental data.     

Model assembly study of the ligand binding by p-hydroxybenzoate hydroxylase: Correlation between the calculated binding energies and the experimental dissociation constants

The energies of binding of seven ligands by p-hydroxybenzoate hydroxylase (PHBH) were calculated theoretically. Direct enzyme–ligand interaction energies were calculated using the ab initio quantum mechanical model assembly of the active site at the 3-21G level. Solvation energies of the ligands needed in the evaluation of the binding energies were calculated with the semiempirical AM1–SM2 method and the long-range electrostatic interaction energies between the ligands and the protein matrix classically using the static charge distributions of the ligands and the protein. Energies for proton-transfer between the ligands OH or SH substituent at position 4 and the active-site tyrosine within the ab initio model assemblies were calculated and compared to the corresponding pK_as in aqueous solution. Excluding 3,4-dihydroxybenzoate, the natural product of PHBH, a linear relationship between the calculated binding energies and the experimental binding free energies was found with a correlation coefficient of 0.90. Contributions of the direct enzyme–ligand interaction energies, solvation energies and the long-range electrostatic interaction energies to the calculated binding energies were analyzed. The proton-transfer energies of the ligands with substituents ortho to the ionized OH were found to be perturbed less in the model calculations than the energies of their meta isomers as deduced from the corresponding pK_as. © 1995 Wiley-Liss, Inc.     

The Influence of Secondary Processing on the Structural Relaxation Dynamics of Fluticasone Propionate

This study investigated the structural relaxation of micronized fluticasone propionate (FP) under different lagering conditions and its influence on aerodynamic particle size distribution (APSD) of binary and tertiary carrier-based dry powder inhaler (DPI) formulations. Micronized FP was lagered under low humidity (LH 25 C, 33% RH [relative humidity]), high humidity (HH 25°C, 75% RH) for 30, 60, and 90 days, respectively, and high temperature (HT 60°C, 44% RH) for 14 days. Physicochemical, surface interfacial properties via cohesive-adhesive balance (CAB) measurements and amorphous disorder levels of the FP samples were characterized. Particle size, surface area, and rugosity suggested minimal morphological changes of the lagered FP samples, with the exception of the 90-day HH (HH90) sample. HH90 FP samples appeared to undergo surface reconstruction with a reduction in surface rugosity. LH and HH lagering reduced the levels of amorphous content over 90-day exposure, which influenced the CAB measurements with lactose monohydrate and salmeterol xinafoate (SX). CAB analysis suggested that LH and HH lagering led to different interfacial interactions with lactose monohydrate but an increasing adhesive affinity with SX. HT lagering led to no detectable levels of the amorphous disorder, resulting in an increase in the adhesive interaction with lactose monohydrate. APSD analysis suggested that the fine particle mass of FP and SX was affected by the lagering of the FP. In conclusion, environmental conditions during the lagering of FP may have a profound effect on physicochemical and interfacial properties as well as product performance of binary and tertiary carrier-based DPI formulations.

Electronic supplementary material

The online version of this article (doi:10.1208/s12249-014-0222-8) contains supplementary material, which is available to authorized users.KEY WORDS: cohesive-adhesive balance, laagering, mechanical activation, particle adhesion, process-induced structural disorder     

Comparison of the stability of fish light meromyosins by guanidine hydrochloride denaturation

The guanidine hydrochloride denaturation of light meromyosins (LMMs) of fish (carp, sardine and greenling) and rabbit was investigated to determine their structural stability quantitatively. The circular dichroism (CD) and fluorescence spectroscopies were applied to monitor denaturation. The CD results indicate that the LMM α-helix undergoes a two-step unfolding. The free energy of denaturation was calculated based on the linear extrapolation method and the denaturant binding model. Total free energies of the two-step unfolding of the α-helix are related to the water temperatures in which the fish live and the body temperature of rabbit. The stability of α-helical structure of LMM was in the following descending order: rabbit>carp>sardine>greenling. The free energies of denaturation obtained by tryptophan fluorescence differ from the free energies of the unfolding α-helix. The data from the two spectroscopic measurements are discussed along with the conformational changes of LMMs.     

Theoretical study on the stability,electronic, mechanical,vibrational and thermodynamic properties of rare-earth intermetallic compound Rh3Ce

The structural, stability, electronic, mechanical, vibrational and thermodynamic properties of rare-earth intermetallic compound Rh₃Ce have been explored systematically by using first-principle calculations. The evaluation of the equilibrium lattice parameters were obtained firstly. Remarkably, the result of calculated unit cell volume, derived by the total energies as a function of volume, is consistent with other results. Next, the values of cohesive energy (E_c), formation enthalpy (ΔH) have verified that Rh₃Ce is a stable compound. In addition, the band structure and the total density of states indicate a metallic behaviour. Furthermore, the Mulliken charges were calculated to understand the bonding in Rh₃Ce compound. Otherwise, the elastic constants(C_ij) as well as other modulus were also calculated to evaluated the mechanical properties of Rh₃Ce. Phonon dispersion curves for Rh₃Ce were depicted to access the vibrational properties. Finally, the thermodynamic properties of Rh₃Ce were summarised range from 0 to 60?GPa, 0 to 1800?K, respectively. We also pointed out that the thermal expansion(α), heat capacity(C_v), entropy(S), Debye temperature(Θ) and Güneisen parameter (γ) change under pressure and temperature.     

Thermodynamics of the Ca2+ binding to bovine alpha-lactalbumin

Bovine alpha-lactalbumin contains one strong Ca2+-binding site. The free energy (delta G0), enthalpy (delta H0), and entropy (delta S0) of binding of Ca2+ to this site have been calculated from microcalorimetric experiments. The enthalpy of binding was dependent on the metal-free bovine alpha-lactalbumin concentration. At 0.8 mg ml-1, metal-free bovine alpha-lactalbumin delta H0 was -110 +/- 6 kJ mol-1. At this concentration the binding constant was estimated from a mathematical analysis of the titration curves to be greater than 10(7) M-1. This means that delta G0 is smaller than -40 kJ mol-1 and delta S0 is less negative than -235 J.K-1 mol-1. The binding of Ca2+ is therefore enthalpy-driven. From binding experiments as a function of temperature, a delta Cp value of -4.1 kJ.K-1 mol-1 was calculated. This value is dependent on the protein concentration. A tentative explanation for this large value is given.     

Empirical calculation of the relative free energies of peptide binding to the molecular chaperone DnaK

We describe a methodology to calculate the relative free energies of protein-peptide complex formation. The interaction energy was decomposed into nonpolar, electrostatic and entropic contributions. A free energy-surface area relationship served to calculate the nonpolar free energy term. The electrostatic free energy was calculated with the finite difference Poisson-Boltzmann method and the entropic contribution was estimated from the loss in the conformational entropy of the peptide side chains. We applied this methodology to a series of DnaK*peptide complexes. On the basis of the single known crystal structure of the peptide-binding domain of DnaK with a bound heptapeptide, we modeled ten other DnaK*heptapeptide complexes with experimentally measured K(d) values from 0.06 microM to 11 microM, using molecular dynamics to refine the structures of the complexes. Molecular dynamic trajectories, after equilibration, were used for calculating the energies with greater accuracy. The calculated relative binding free energies were compared with the experimentally determined free energies. Linear scaling of the calculated terms was applied to fit them to the experimental values. The calculated binding free energies were between -7.1 kcal/mol and - 9.4 kcal/mol with a correlation coefficient of 0.86. The calculated nonpolar contributions are mainly due to the central hydrophobic binding pocket of DnaK for three amino acid residues. Negative electrostatic fields generated by the protein increase the binding affinity for basic residues flanking the hydrophobic core of the peptide ligand. Analysis of the individual energy contributions indicated that the nonpolar contributions are predominant compared to the other energy terms even for peptides with low affinity and that inclusion of the change in conformational entropy of the peptide side chains does not improve the discriminative power of the calculation. The method seems to be useful for predicting relative binding energies of peptide ligands of DnaK and might be applicable to other protein-peptide systems, particularly if only the structure of one protein-ligand complex is available.     

Energetic analysis of binding of progesterone and 5 beta-androstane-3,17-dione to anti-progesterone antibody DB3 using molecular dynamics and free energy calculations.

Molecular dynamics simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) free energy calculations were used to study the energetics of the binding of progesterone (PRG) and 5 beta-androstane-3,17-dione (5AD) to anti-PRG antibody DB3. Although the two steroids bind to DB3 in different orientations, their binding affinities are of the same magnitude, 1 nM for PRG and 8 nM for 5AD. The calculated relative binding free energy of the steroids, 8.8 kJ/mol, is in fair agreement with the experimental energy, 5.4 kJ/mol. In addition, computational alanine scanning was applied to study the role of selected amino acid residues of the ligand-binding site on the steroid cross-reactivity. The electrostatic and van der Waals components of the total binding free energies were found to favour more the binding of PRG, whereas solvation energies were more favourable for the binding of 5AD. The differences in the free energy components are due to the binding of the A rings of the steroids to different binding pockets: PRG is bound to a pocket in which electrostatic antibody-steroid interactions are dominating, whereas 5AD is bound to a pocket in which van der Waals and hydrophobic interactions dominate.     

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.     

Ligand-induced conformational changes: improved predictions of ligand binding conformations and affinities

A computational docking strategy using multiple conformations of the target protein is discussed and evaluated. A series of low molecular weight, competitive, nonpeptide protein tyrosine phosphatase inhibitors are considered for which the x-ray crystallographic structures in complex with protein tyrosine phosphatase 1B (PTP1B) are known. To obtain a quantitative measure of the impact of conformational changes induced by the inhibitors, these were docked to the active site region of various structures of PTP1B using the docking program FlexX. Firstly, the inhibitors were docked to a PTP1B crystal structure cocrystallized with a hexapeptide. The estimated binding energies for various docking modes as well as the RMS differences between the docked compounds and the crystallographic structure were calculated. In this scenario the estimated binding energies were not predictive inasmuch as docking modes with low estimated binding energies corresponded to relatively large RMS differences when aligned with the corresponding crystal structure. Secondly, the inhibitors were docked to their parent protein structures in which they were cocrystallized. In this case, there was a good correlation between low predicted binding energy and a correct docking mode. Thirdly, to improve the predictability of the docking procedure in the general case, where only a single target protein structure is known, we evaluate an approach which takes possible protein side-chain conformational changes into account. Here, side chains exposed to the active site were considered in their allowed rotamer conformations and protein models containing all possible combinations of side-chain rotamers were generated. To evaluate which of these modeled active sites is the most likely binding site conformation for a certain inhibitor, the inhibitors were docked against all active site models. The receptor rotamer model corresponding to the lowest estimated binding energy is taken as the top candidate. Using this protocol, correct inhibitor binding modes could successfully be discriminated from proposed incorrect binding modes. Moreover, the ranking of the estimated ligand binding energies was in good agreement with experimentally observed binding affinities.     

Fluorescent derivatives of the GFP chromophore give a new insight into the GFP fluorescence process

The photophysical properties of synthetic compounds derived from the imidazolidinone chromophore of the green fluorescent protein were determined. Various electron-withdrawing or electron-donating substituents were introduced to mimic the effect of the chromophore surroundings in the protein. The absorption and emission spectra as well as the fluorescence quantum yields in dioxane and glycerol were shown to be highly dependent on the electronic properties of the substituents. We propose a kinetic scheme that takes into account the temperature-dependent twisting of the excited molecule. If the activation energy is low, the molecule most often undergoes an excited-state intramolecular twisting that leads it to the ground state through an avoided crossing between the S(1) and S(0) energy surfaces. For a high activation energy, the torsional motion within the compounds is limited and the ground-state recovery will occur preferentially by fluorescence emission. The excellent correlation between the fluorescence quantum yields and the calculated activation energies to torsion points to the above-mentioned avoided crossing as the main nonradiative deactivation channel in these compounds. Finally, our results are discussed with regard to the chromophore in green fluorescent protein and some of its mutants.     

Free energy simulations and MM-PBSA analyses on the affinity and specificity of steroid binding to antiestradiol antibody

Antiestradiol antibody 57-2 binds 17beta-estradiol (E2) with moderately high affinity (K(a) = 5 x 10(8) M(-1)). The structurally related natural estrogens estrone and estriol as well synthetic 17-deoxy-estradiol and 17alpha-estradiol are bound to the antibody with 3.7-4.9 kcal mol(-1) lower binding free energies than E2. Free energy perturbation (FEP) simulations and the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method were applied to investigate the factors responsible for the relatively low cross-reactivity of the antibody with these four steroids, differing from E2 by the substituents of the steroid D-ring. In addition, computational alanine scanning of the binding site residues was carried out with the MM-PBSA method. Both the FEP and MM-PBSA methods reproduced the experimental relative affinities of the five steroids in good agreement with experiment. On the basis of FEP simulations, the number of hydrogen bonds formed between the antibody and steroids, which varied from 0 to 3 in the steroids studied, determined directly the magnitude of the steroid-antibody interaction free energies. One hydrogen bond was calculated to contribute about 3 kcal mol(-1) to the interaction energy. Because the relative binding free energies of estrone (two antibody-steroid hydrogen bonds), estriol (three hydrogen bonds), 17-deoxy-estradiol (no hydrogen bonds), and 17alpha-estradiol (two hydrogen bonds) are close to each other and clearly lower than that of E2 (three hydrogen bonds), the water-steroid interactions lost upon binding to the antibody make an important contribution to the binding free energies. The MM-PBSA calculations showed that the binding of steroids to the antiestradiol antibody is driven by van der Waals interactions, whereas specificity is solely due to electrostatic interactions. In addition, binding of steroids to the antiestradiol antibody 57-2 was compared to the binding to the antiprogesterone antibody DB3 and antitestosterone antibody 3-C4F5, studied earlier with the MM-PBSA method.     

Computational Analysis of Phosphopeptide Binding to the Polo-Box Domain of the Mitotic Kinase PLK1 Using Molecular Dynamics Simulation

The Polo-Like Kinase 1 (PLK1) acts as a central regulator of mitosis and is over-expressed in a wide range of human tumours where high levels of expression correlate with a poor prognosis. PLK1 comprises two structural elements, a kinase domain and a polo-box domain (PBD). The PBD binds phosphorylated substrates to control substrate phosphorylation by the kinase domain. Although the PBD preferentially binds to phosphopeptides, it has a relatively broad sequence specificity in comparison with other phosphopeptide binding domains. We analysed the molecular determinants of recognition by performing molecular dynamics simulations of the PBD with one of its natural substrates, CDC25c. Predicted binding free energies were calculated using a molecular mechanics, Poisson-Boltzmann surface area approach. We calculated the per-residue contributions to the binding free energy change, showing that the phosphothreonine residue and the mainchain account for the vast majority of the interaction energy. This explains the very broad sequence specificity with respect to other sidechain residues. Finally, we considered the key role of bridging water molecules at the binding interface. We employed inhomogeneous fluid solvation theory to consider the free energy of water molecules on the protein surface with respect to bulk water molecules. Such an analysis highlights binding hotspots created by elimination of water molecules from hydrophobic surfaces. It also predicts that a number of water molecules are stabilized by the presence of the charged phosphate group, and that this will have a significant effect on the binding affinity. Our findings suggest a molecular rationale for the promiscuous binding of the PBD and highlight a role for bridging water molecules at the interface. We expect that this method of analysis will be very useful for probing other protein surfaces to identify binding hotspots for natural binding partners and small molecule inhibitors.     

Energetic basis of molecular recognition in a DNA aptamer

The thermal stability and ligand binding properties of the L-argininamide-binding DNA aptamer (5'-GATCGAAACGTAGCGCCTTCGATC-3') were studied by spectroscopic and calorimetric methods. Differential calorimetric studies showed that the uncomplexed aptamer melted in a two-state reaction with a melting temperature T(m)=50.2+/-0.2 degrees C and a folding enthalpy DeltaH(0)(fold)=-49.0+/-2.1 kcal mol(-1). These values agree with values of T(m)=49.6 degrees C and DeltaH(0)(fold)=-51.2 kcal mol(-1) predicted for a simple hairpin structure. Melting of the uncomplexed aptamer was dependent upon salt concentration, but independent of strand concentration. The T(m) of aptamer melting was found to increase as L-argininamide concentrations increased. Analysis of circular dichroism titration data using a single-site binding model resulted in the determination of a binding free energy DeltaG(0)(bind)=-5.1 kcal mol(-1). Isothermal titration calorimetry studies revealed an exothermic binding reaction with DeltaH(0)(bind)=-8.7 kcal mol(-1). Combination of enthalpy and free energy produce an unfavorable entropy of -TDeltaS(0)=+3.6 kcal mol(-1). A molar heat capacity change of -116 cal mol(-1) K(-1) was determined from calorimetric measurements at four temperatures over the range of 15-40 degrees C. Molecular dynamics simulations were used to explore the structures of the unligated and ligated aptamer structures. From the calculated changes in solvent accessible surface areas of these structures a molar heat capacity change of -125 cal mol(-1) K(-1) was calculated, a value in excellent agreement with the experimental value. The thermodynamic signature, along with the coupled CD spectral changes, suggest that the binding of L-argininamide to its DNA aptamer is an induced-fit process in which the binding of the ligand is thermodynamically coupled to a conformational ordering of the nucleic acid.     

Adhesion energy of receptor-mediated interaction measured by elastic deformation

We investigated the role of receptor binding affinity in surface adhesion. A sensitive technique was developed to measure the surface energy of receptor-mediated adhesion. The experimental system involved a functionalized elastic agarose bead resting on a functionalized glass coverslip. Attractive intersurface forces pulled the two surfaces together, deforming the bead to produce an enlarged contact area. The Johnson-Kendall-Roberts (JKR) model was used to relate the surface energy of the interaction to the elasticity of the bead and the area of contact. The surface energies for different combinations of modified surfaces in solution were obtained from reflection interference contrast microscopy (RICM) measurements of the contact area formed by the bead and the coverslip. Studies with surfaces functionalized with ligand-receptor pairs showed that the relationship between surface energy and the association constant of the ligand binding has two regimes. At low binding affinity, surface energy increased linearly with the association constant, while surface energy increased logarithmically with the association constant in the high affinity regime.     

Energy of binding of Aspergillus oryzae beta-glucosidase with the substrate, and the mechanism of its enzymic action

Several beta-D-glucopyranosides (p-nitrophenyl, phenyl, and ethyl), 1-thio-beta-D-glucopyranosides, and phenyl 2-deoxy, 3-deoxy, 4-deoxy, and 6-deoxy beta-D-glucopyranosides were synthesized and used to study the mechanism of the enzymatic action of Taka-beta-glucosidase [EC 3.2.1.21 Aspergillus oryzae]. Kinetic constants of the enzyme for these glycosides were determined from S/V-S or 1/V-1/S plots, and the hydrolysis rates of these compounds with the enzyme, acid (3 N HCl) and alkali (3 N NaOH) were compared. Inhibition of the enzyme by 1,5-anhydroglucitol, glucal, dihydroglucal, and 1,6-anhydroglucopyranose was also examined. Glucal and 1,5-anhydroglucitol showed strong competitive inhibition. Free energy of binding of each hydroxyl group of glucosidic glucose with the enzyme was estimated from Kms of phenyl beta-glucoside and its deoxy analogues, and also Ki values of some inhibitors. The free energies of binding of 2-OH, 3-OH, 4-OH, and 6-OH were calculated to be 1.1, 2.4, 0.7, and 1.8 kcal/mol, respectively. The free energy of binding of phenoxide at C-1 (0.3 kcal/mol) was calculated from the Km of Ph-beta-Glc and Ki of 1,5-anhydroglucitol. The energy of binding of 5-CH2OH (2.3 kcal/mol) was obtained from the Km of Ph-beta-Glc and that of Ph-beta-Xyl. The sum (6.8 kcal/mol) of each partial binding free energy was close to the value of binding free energy of Ph-beta-Glc (7.0 kcal/mol) calculated by the equation; -delta Gbind = -RT ln Km-T delta Smix, showing that the methods of estimation of each binding energy used in the present study seemed reasonable. Glucal, having a pyranose form distorted slightly, showed strong competitive inhibition and the Ki of this inhibitor was smaller than the Km of Ph-beta-Glc, suggesting that the sugar ring bound to the active site was distored to a half chair form which is labile to acid hydrolysis.     

Zinc finger protein binding to DNA: an energy perspective using molecular dynamics simulation and free energy calculations on mutants of both zinc finger domains and their specific DNA bases

Energy calculations based on MM-GBSA were employed to study various zinc finger protein (ZF) motifs binding to DNA. Mutants of both the DNA bound to their specific amino acids were studied. Calculated energies gave evidence for a relationship between binding energy and affinity of ZF motifs to their sites on DNA. ΔG values were ?15.82(12), ?3.66(12), and ?12.14(11.6) kcal/mol for finger one, finger two, and finger three, respectively. The mutations in the DNA bases reduced the value of the negative energies of binding (maximum value for ΔΔG = 42Kcal/mol for F1 when GCG mutated to GGG, and ΔΔG = 22 kcal/mol for F2, the loss in total energy of binding originated in the loss in electrostatic energies upon mutation (r = .98). The mutations in key amino acids in the ZF motif in positions-1, 2, 3, and 6 showed reduced binding energies to DNA with correlation coefficients between total free energy and electrostatic was .99 and with Van der Waal was .93. Results agree with experimentally found selectivity which showed that Arginine in position-1 is specific to G, while Aspartic acid (D) in position 2 plays a complicated role in binding. There is a correlation between the MD calculated free energies of binding and those obtained experimentally for prepared ZF motifs bound to triplet bases in other reports (), our results may help in the design of ZF motifs based on the established recognition codes based on energies and contributing energies to the total energy.