Molecular dynamics simulation and binding free energy studies of novel leads belonging to the benzofuran class inhibitors of Mycobacterium tuberculosis Polyketide Synthase 13

In this work, the binding mechanism of new Polyketide Synthase 13 (Pks13) inhibitors has been studied through molecular dynamics simulation and free energy calculations. The drug Tam1 and its analogs, belonging to the benzofuran class, were submitted to 100 ns simulations, and according to the results obtained for root mean square deviation, all the simulations converged from approximately 30 ns. For the analysis of backbone flotation, the root mean square fluctuations were plotted for the Cα atoms; analysis revealed that the greatest fluctuation occurred in the residues that are part of the protein lid domain. The binding free energy value (ΔG_bind) obtained for the Tam16 lead molecule was of ?51.43 kcal/mol. When comparing this result with the ΔG_bind values for the remaining analogs, the drug Tam16 was found to be the highest ranked: this result is in agreement with the experimental results obtained by Aggarwal and collaborators, where it was verified that the IC₅₀ for Tam16 is the smallest necessary to inhibit the Pks13 (IC₅₀ = 0.19 μM). The energy decomposition analysis suggested that the residues which most interact with inhibitors are: Ser1636, Tyr1637, Asn1640, Ala1667, Phe1670, and Tyr1674, from which the greatest energy contribution to Phe1670 was particularly notable. For the lead molecule Tam16, a hydrogen bond with the hydroxyl of the phenol not observed in the other analogs induced a more stable molecular structure. Aggarwal and colleagues reported this hydrogen bonding as being responsible for the stability of the molecule, optimizing its physic-chemical, toxicological, and pharmacokinetic properties.     

Ensemble-based virtual screening: identification of a potential allosteric inhibitor of Bcr-Abl

Ensemble-based virtual screening using different conformations of a target protein is gaining popularity, as it can leverage information from target flexibility for effective lead identification. In this paper, molecular dynamics simulation followed by RMSD-based clustering was employed to generate and choose distinct conformations of Bcr-Abl. Three representative structures from the most-populated clusters along with the crystal structure conformation (PDBID: 3K5V) were used to perform docking-based virtual screening of 14,400 compounds (in the Maybridge database) in order to identify potential allosteric site binders. Seven compounds found as hits in at least three of the four virtual screenings had higher Glide docking scores than the co-crystallized allosteric inhibitor GNF-2. Detailed computational analyses of the seven hits identified SEW02675 (ΔG _bind = ?164.92 kJ/mol with the wild-type (wt) Bcr-Abl and ?167.37 kJ/mol with the T334I Bcr-Abl mutant) as a better allosteric site binder with both the wt and the mutant Bcr-Abl protein than the reference allosteric inhibitor GNF-2 (ΔG _bind = ?103.12 with wt and ?142.96 kJ/mol with T334I). Moreover, the presence of SEW02675 in the allosteric site enhanced the binding of imatinib (ΔG _bind = ?367.58 with wt and ?294.56 kJ/mol with T334I) to the ATP sites of the wt and the mutant Bcr-Abl. However, when GNF-2 was present in the allosteric site, the binding of imatinib (ΔG _bind = ?351.76 with wt and ?273.94 kJ/mol with T334I) to the ATP site was weaker. The in silico findings suggest that SEW02675 could be used in combination with imatinib to treat chronic myeloid leukemia, and that it could help to overcome resistance due to T334I Bcr-Abl mutation.

Open image in new window

Graphical abstract Virtual screening strategy to identify allosteric inhbitors of Bcr-Abl for the treatment of Chronic myeloid leukemia.

     

Repurposing of FDA-approved drugs to target MurB and MurE enzymes in Mycobacterium tuberculosis

Abstract

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb) is one amongst the top 10 causes of death worldwide. The growing rise in antibiotic resistance compounded with slow and expensive drug discovery has further aggravated the situation. ‘Drug repurposing’ is a promising approach where known drugs are examined for a new indication. In the present study, we have attempted to identify drugs that could target MurB and MurE enzymes involved in the muramic acid synthesis pathway (Mur Pathway) in Mtb. FDA-approved drugs from two repositories i.e. Drug Bank (1932 drugs) and e-LEA3D (1852 drugs) were screened against these proteins. Several criteria were applied to study the protein-drug interactions and the consensus drugs were further studied by molecular dynamics (MD) simulation. Our study found Sulfadoxine (–7.3?kcal/mol) and Pyrimethamine (–7.8?kcal/mol) to show stable interaction with MurB while Lifitegrast (–10.5?kcal/mol) and Sildenafil (–9.1?kcal/mol) showed most reliable interaction with MurE. Furthermore, binding free energy (ΔG_bind), RMSD and RMSF data and the number of hydrogen bonds corroborated the stability of interactions and hence these drugs for repurposing should be explored further.

Communicated by Ramaswamy H. Sarma     

Binding affinity of pyrano[3, 2-f]quinoline and DNA: spectroscopic and docking approach

The interaction between pyrano[3, 2-f]quinoline (PQ) and calf thymus DNA (CTDNA) using spectroscopic and molecular modeling approach has been presented here. Apparent association constant (1.05×10⁵ L/mol) calculated from UV-vis specta, indicates a moderate complex formation between CTDNA and PQ. The quenching phenomena as obtained from emission spectra of ethidium bromide (EB)–CTDNA by PQ was found to be a dynamic one and the binding constants found to be 8.64, 9.25, 11.17, 12.03 × 10⁴ L/mol at 293, 300, 308, and 315 K. Thermodynamic parameter enthalpy change (ΔH) and entropy change (ΔS), indicates weak force like van der Walls force and hydrogen bonds having the key role in this binding process. The results of circular dichroism (CD) demonstrate that PQ has not induced characteristic changed in CTDNA. Results achieved from UV absorption and fluorescence spectroscopy indicating the binding mode of PQ with DNA seems to be a nonintercalative binding. The theoretical results as originating from molecular modeling showed that PQ possibly will bind into the hydrophobic region of DNA having docking binding energy = ?10.03 kcal/mol and the obtained results are in consonance with the inferences obtained from experimental data. This result is important for the better understanding of pharmaceutical aspects of binding affinity of PQ and CTDNA.     

Multi-drug resistance profile of PR20 HIV-1 protease is attributed to distorted conformational and drug binding landscape: molecular dynamics insights

The PR20 HIV-1 protease, a variant with 20 mutations, exhibits high levels of multi-drug resistance; however, to date, there has been no report detailing the impact of these 20 mutations on the conformational and drug binding landscape at a molecular level. In this report, we demonstrate the first account of a comprehensive study designed to elaborate on the impact of these mutations on the dynamic features as well as drug binding and resistance profile, using extensive molecular dynamics analyses. Comparative MD simulations for the wild-type and PR20 HIV proteases, starting from bound and unbound conformations in each case, were performed. Results showed that the apo conformation of the PR20 variant of the HIV protease displayed a tendency to remain in the open conformation for a longer period of time when compared to the wild type. This led to a phenomena in which the inhibitor seated at the active site of PR20 tends to diffuse away from the binding site leading to a significant change in inhibitor–protein association. Calculating the per-residue fluctuation (RMSF) and radius of gyration, further validated these findings. MM/GBSA showed that the occurrence of 20 mutations led to a drop in the calculated binding free energies (ΔG_bind) by ~25.17 kcal/mol and ~5 kcal/mol for p2-NC, a natural peptide substrate, and darunavir, respectively, when compared to wild type. Furthermore, the residue interaction network showed a diminished inter-residue hydrogen bond network and changes in inter-residue connections as a result of these mutations. The increased conformational flexibility in PR20 as a result of loss of intra- and inter-molecular hydrogen bond interactions and other prominent binding forces led to a loss of protease grip on ligand. It is interesting to note that the difference in conformational flexibility between PR20 and WT conformations was much higher in the case of substrate-bound conformation as compared to DRV. Thus, developing analogues of DRV by retaining its key pharmacophore features will be the way forward in the search for novel protease inhibitors against multi-drug resistant strains.     

Recognition of specific and nonspecific DNA by human lactoferrin

The general principles of recognition of nucleic acids by proteins are among the most exciting problems of molecular biology. Human lactoferrin (LF) is a remarkable protein possessing many independent biological functions, including interaction with DNA. In human milk, LF is a major DNase featuring two DNA‐binding sites with different affinities for DNA. The mechanism of DNA recognition by LF was studied here for the first time. Electrophoretic mobility shift assay and fluorescence measurements were used to probe for interactions of the high‐affinity DNA‐binding site of LF with a series of model‐specific and nonspecific DNA ligands, and the structural determinants of DNA recognition by LF were characterized quantitatively. The minimal ligands for this binding site were orthophosphate (K_i = 5 mM), deoxyribose 5'‐phosphate (K_i = 3 mM), and different dNMPs (K_i = 0.56–1.6 mM). LF interacted additionally with 9–12 nucleotides or nucleotide pairs of single‐ and double‐stranded ribo‐ and deoxyribooligonucleotides of different lengths and sequences, mainly through weak additive contacts with internucleoside phosphate groups. Such nonspecific interactions of LF with noncognate single‐ and double‐stranded d(pN)₁₀ provided ~6 to ~7.5 orders of magnitude of the enzyme affinity for any DNA. This corresponds to the Gibbs free energy of binding (ΔG⁰) of ?8.5 to ?10.0 kcal/mol. Formation of specific contacts between the LF and its cognate DNA results in an increase of the DNA affinity for the enzyme by approximately 1 order of magnitude (K_d = 10 nM; ΔG⁰ ≈ ?11.1 kcal/mol). A general function for the LF affinity for nonspecific d(pN)_n of different sequences and lengths was obtained, giving the K_d values comparable with the experimentally measured ones. A thermodynamic model was constructed to describe the interactions of LF with DNA. Copyright © 2013 John Wiley & Sons, Ltd.     

Recognition of a Flipped Base in a Hairpinloop DNA by a Small Peptide

Two tiny hairpin DNAs, CORE (dAGGCTTCGGCCT) and AP2 (dAGGCTXCGGCCT; X: abasic nucleotide), fold into almost the same tetraloop hairpin structure with one exception, that is, the sixth thymine (T6) of CORE is exposed to the solvent water (Kawakami, J. et al., Chem. Lett. 2001, 258–259). In the present study, we selected small peptides that bind to CORE or AP2 from a combinatorial pentapeptide library with 2.5 × 10⁶ variants. On the basis of the structural information, the selected peptide sequences should indicate the essential qualifications for recognition of the hairpin loop DNA with and without a flipped base. In the selected DNA binding peptides, aromatic amino acids such as histidine for CORE and glutamine/aspartic acid for AP2 were found to be abundant amino acids. This amino acid preference suggests that CORE-binding peptides use π–π stacking to recognize the target while hydrogen bonding is dominant for AP2-binding peptides. To investigate the binding properties of the selected peptide to the target, surface plasmon resonance was used. The binding constant of the interaction between CORE and a CORE-binding peptide (HWHHE) was about 1.1 × 10⁶ M^?1 at 25°C and the resulting binding free energy change at 25°C (ΔG°₂₅) was ?8.2 kcal mol^?1. The binding of the peptide to AP2 was also analyzed and the resulting binding constant and ΔG°₂₅ were about 4.2 × 10⁴ M^?1 and ?6.3 kcal mol^?1, respectively. The difference in the binding free energy changes (ΔΔG°₂₅) of 1.9 kcal mol^?1 was comparable to the values reported in other systems and was considered a consequence of the loss of π–π stacking. Moreover, the stabilization effect by stacking affected the dissociation step as well as the association step. Our results suggest that the existence of an aromatic ring (T6 base) produces new dominant interactions between peptides and nucleic acids, although hydrogen bonding is the preferable mode of interaction in the absence of the flipping base. These findings regarding CORE and AP2 recognition are expected to give useful information in the design of novel artificial DNA binding peptides.     

Binding profiles of cholesterol ester transfer protein with current inhibitors: a look at mechanism and drawback

Although the pharmacological inhibition of cholesterol ester transport protein (CETP) has been proposed as a method of preventing and treating cardiovascular disease (CVD), the adverse effects of current inhibitors have cast doubt on the interaction mechanisms of inhibitors and CETP. In response, a molecular dynamics simulation was used to investigate their interaction and shed light on the lipid exchange mechanism of CETP. Results showed that torcetrapib, anacetrapib, and evacetrapib can induce the incremental rigidity of CETP, yet decrease the stability of Helix X and the hydrophobic tunnel of CETP, with passable binding abilities (ΔG_bind, ?61.08, ?64.23, and ?61.57 kcal mol^?1). During their binding processes, Van der Waals components (ΔE_vdw + ΔG_SA) play a dominant role, and the inhibitory effects closely correlated with residues Cys13, Val198, Gln199, Ser230, His232, and Phe263, which could reduce the flexibility of N- and C- termini and Helix X, as well as the stability of hydrophobic tunnel, into which the three inhibitors could enter and promote the formation of intramolecular H-bonds such as Thr138–Asn192 and Arg37–Glu186. Additionally, the three inhibitors could restrain the formation of an opening at the CETP N-terminal, which given the other findings suggests the tunneling mechanism of CETP transfer. The paper closes with an explanation of conceivable causes of the insufficient efficacy of the inhibitors, and puts forward the rationality in targeting the CETP distal end for CVD therapies.     

Exploring additivity effects of double mutations on the binding affinity of protein‐protein complexes

Additivity in binding affinity of protein‐protein complexes refers to the change in free energy of binding (ΔΔG_bind) for double (or multiple) mutations which is approximately equal to the sum of their corresponding single mutation ΔΔG_bind values. In this study, we have explored the additivity effect of double mutants, which shows a linear relationship between the binding affinity of double and sum of single mutants with a correlation of 0.90. However, the comparison of ΔΔG_bind values showed a mean absolute deviation of 0.86 kcal/mol, and 25.6% of the double mutants show a deviation of more than 1 kcal/mol, which are identified as non‐additive. The additivity effects have been analyzed based on the influence of structural features such as accessible surface area, long range order, binding propensity change, surrounding hydrophobicity, flexibility, atomic contacts between the mutations and distance between the 2 mutations. We found that non‐additive mutations tend to be closer to each other and have more contacts. We have also used machine learning methods to discriminate additive and non‐additive mutations using structure‐based features, which showed the accuracies in the range of 0.77–0.92 for protein‐protein complexes belonging to different functions. Further, we have compared the additivity effects of protein stability along with binding affinity and explored the similarities and differences between them. The results obtained in this study provide insights into the effects of various structural features on binding affinity of double mutants, and will aid the development of accurate methods to predict the binding affinity of double mutants.     

Structural insights into conformational stability and binding of thiazolo-[2,3-b] quinazolinone derivatives with EGFR-TKD and in-vitro study

Heterocyclic molecules are well-known drugs against various diseases including cancer. Many tyrosine kinase inhibitors including erlotinib, osimertinib, and sunitinib were developed and approved but caused adverse effects among treated patients. Which prevents them from being used as cancer therapeutics. In this study, we strategically developed heterocyclic thiazolo-[2,3-b]quinazolinone derivatives by an organic synthesis approach. These synthesized molecules were assessed against the epidermal growth factor receptor tyrosine kinase domain (EGFR-TKD) by in silico methods. Molecular docking simulations unravel derivative 17 showed better binding energy scores and followed Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) properties. The binding affinity displayed by synthetic congener and reference molecule erlotinib was found to be ?8.26 ± 0.0033 kcal/mol and ?7.54 ± 0.1411 kcal/mol with the kinase domain. Further, molecular dynamic simulations were conducted thrice to validate the molecular docking study and achieved significant results. Both synthetic derivative and reference molecule attained stability in the active site of the TKD. The synthetic congener and erlotinib showed free energy binding (ΔG_bind) ?102.975 ± 3.714 kJ/mol and ?130.378 ± 0.355 kJ/mol computed by Molecular Mechanics Poison Boltzmann Surface Area (MM-PBSA) method. In addition, the motions of each sampled system including the Apo complex were determined by the principal component analysis and Gibbs energy landscape analysis. The in-vitro apoptosis study was performed using MCF-7 and H-1299 cancer cell lines. However, thiazolo-[2,3-]-quinazoline derivative 17 showed fair anti-proliferative activity against MCF-7 and H-1299. Further, the in-vivo study is necessary to determine the effectivity of the potent anti-proliferative, non-toxic molecule against TKD.     

Protoberberine Alkaloids Berberine,Palmatine, and Coralyne Binding to Poly(dT)⋅(Poly(dA)⋅Poly(dT)) Triplex: Comparative Structural Aspects and Energetics Profiles of the Interaction

The interaction of bioactive protoberberine alkaloids berberine, palmatine, and coralyne with the DNA triplex poly(dT)⋅(poly(dA)⋅poly(dT)) was studied using biophysical and calorimetric techniques. All three alkaloids bound the triplex cooperatively. Berberine and palmatine predominantly stabilized the triplex structure, while coralyne stabilized both triplex and duplex structures as inferred from optical thermal melting profiles. Fluorescence quenching, polarization, and viscometric studies hinted at an intercalative mode of binding for the alkaloids to the triplex, coralyne being more strongly intercalated compared to partial intercalation of berberine and palmatine. The overall affinity of coralyne was two order higher (2.29×10⁷ M ⁻¹) than that of berberine (3.43×10⁵ M ⁻¹) and palmatine (2.34×10⁵ M ⁻¹). Isothermal titration calorimetric studies revealed that the binding to the triplex was favored by negative enthalpy change (ΔH=−3.34 kcal/mol) with favorable entropy contribution (TΔS = 4.07 kcal/mol) for berberine, favored by almost equal negative enthalpy (ΔH =−3.88 kcal/mol) and entropy changes (TΔS = 3.37 kcal/mol) for palmatine, but driven by large enthalpy contributions (ΔH =−25.62 kcal/mol and TΔS =−15.21 kcal/mol) for coralyne. These results provide new insights on the binding of isoquinoline alkaloids to the DNA triplex structure.     

Evaluation of DNA,BSA binding,and antimicrobial activity of new synthesized neodymium complex containing 29-dimethyl 110-phenanthroline

In order to evaluate biological potential of a novel synthesized complex [Nd(dmp)₂Cl₃.OH₂] where dmp is 29-dimethyl 110-phenanthroline, the DNA-binding, cleavage, BSA binding, and antimicrobial activity properties of the complex are investigated by multispectroscopic techniques study in physiological buffer (pH 7.2).The intrinsic binding constant (K_b) for interaction of Nd(III) complex and FS–DNA is calculated by UV–Vis (K_b = 2.7 ± 0.07 × 10⁵) and fluorescence spectroscopy (K_b = 1.13 ± 0.03 × 10⁵). The Stern–Volmer constant (K_SV), thermodynamic parameters including free energy change (ΔG°), enthalpy change (?H°), and entropy change (?S°), are calculated by fluorescent data and Vant’ Hoff equation. The experimental results show that the complex can bind to FS–DNA and the major binding mode is groove binding. Meanwhile, the interaction of Nd(III) complex with protein, bovine serum albumin (BSA), has also been studied by using absorption and emission spectroscopic tools. The experimental results show that the complex exhibits good binding propensity to BSA. The positive ΔH° and ?S° values indicate that the hydrophobic interaction is main force in the binding of the Nd(III) complex to BSA, and the complex can quench the intrinsic fluorescence of BSA remarkably through a static quenching process. Also, DNA cleavage was investigated by agarose gel electrophoresis that according to the results cleavage of DNA increased with increasing of concentration of the complex. Antimicrobial screening test gives good results in the presence of Nd(III) complex system.     

The enthalpy of oxidation of flavin mononucleotide: Temperature dependence of in vitro bacterial luciferase bioluminescence

The enthalpy of the bioluminescent reaction

$\text{FMNH}{}_2\text{+ RCHO + O}{}_2\text{→}\text{luciferase}\text{FMN + RCOO + H}{}_3\text{O}{}^{\mathrm{+}}\text{+ hv}$

has been studied by direct calorimetric methods. Bacterial luciferase, isolated from Beneckea harveyi (formerly strain MAV) has been used to catalyze the oxidation of reduced flavin mononucleotide (FMNH₂) and a long chain aliphatic aldehyde (dodecanal, RCHO) by molecular oxygen to give the indicated products and blue-green light. The enthalpy measured for this process was found to be ΔH_L = ?338.9 k.J (mol FMN)^?1 (?81.0 kcal) at 25.00 °C and ?402.9 kJ (mol FMN)^?1 (?96.3 kcal) at 7.00 °C. Calculations based on redox electrode potentials indicate a corresponding value of the free energy change, ΔG_L = ?464.8 kJ (mol FMN)^?1 (?111.1 kcal), at 25 °C. Measurements were performed in 0.15 m phosphate buffer, pH 7.0 and the values were arrived at by correcting the observed heats for the heat associated with the autoxidation process: FMNH₂ + O₂ ? FMN + H₂O₂; ΔH_D = ?158.5 kJ (mol FMN)^?1 (?37.8). These data and a detailed thermodynamic analysis have demonstrated the need for two parameters, referred to as the intrinsic free energy, ΔG₁, and intrinsic enthalpy, ΔH₁, which are functionally defined by the relations ΔG_I = ΔG_L ? uhvΔH_I = ΔH_L ? uhv, where u is the quantum yield of the reaction expressed in einsteins mole^?1.These parameters reflect the thermochemistry of the bioluminescent reaction corrected for emitted photons. Thus, they are useful for comparing the thermochemistry of a chemiluminescent process. Their values for the bacterial luciferase system at 25 °C and pH 7.0 are ?391.6 and ?266.9 kJ (mol FMN)^?1 (?93.6 and ?63.8 kcal), respectively, assuming a value of 0.3 for the quantum yield. The calorimetric data also suggest the existence of a long-lived species which persists after photon emission.     

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.     

The standard Gibbs free energy change of hydrolysis of alpha-D-ribose 1-phosphate

The standard Gibbs free energy change of hydrolysis of α-d-ribose 1-phosphate has been measured at pH 7.0, ionic strength 0.1 m, and 25 °C by combining the corresponding values of the two following reactions: adenosine + H₂O ág adenine + ribose (ΔG^0′ = ?2.3 ± 0.1 kcal/mol), catalyzed by adenosine nucleosidase, and ribose 1-phosphate + adenine ág adenosine + P_i (ΔG^0′ = ?3.1 ± 0.1 kcal/mol), catalyzed by adenosine phosphorylase. The standard Gibbs free energy changes were calculated for both reactions from the equilibrium constant. A value of -5.4 ± 0.15 kcal/mol, comparable to that of other hemiacetal phosphoric esters, was obtained for the hydrolysis of ribose 1-phosphate.     

Isocyanurates with Planar Chirality: Design,Optical Resolution,and Isomerization

Designs and syntheses of isocyanurates ( 1–3 ) are described on the basis of a novel concept that two enantiotopic faces of C_s‐symmetric, prochiral planar molecules are differentiated with a location of groups at the top or bottom of the planar skeleton using a rigid linker. Such isocyanurates are atropisomeric. The planar‐chiral structures of 1 and 2 _anti (anti‐conformer of 2 ) were confirmed by single‐crystal X‐ray analyses, and the space groups were P1 (for 1 ) and P2₁/c (for 2 _anti), resulting that the crystals were racemates. Optical resolutions of 1–3 were successfully accomplished by using chiral high‐performance liquid chromatography technique in combination with circular dichroism, absorption, and nuclear magnetic resonance spectroscopies and mass spectrometry. Furthermore, the rotational barriers (ΔG^?s) related to isomerizations of 1–3 were estimated to be 27.2 (for 1 at 50 °C), 27.6 (for 2 _anti at 50 °C), and 40.6 (for 3 _syn at 150 °C) kcal/mol. The ΔG^?s of 2 and 3 were higher than that of 1 and, in particular, that of 3 was highest among them. This result indicates that an introduction of bulky substituents and an intramolecular bridging are effective for inhibitions of the isomerizations. Chirality 00:000–000, 2012. © 2012 Wiley Periodicals, Inc.     

Contribution of hydrophobic interactions to protein stability

Our goal was to gain a better understanding of the contribution of hydrophobic interactions to protein stability. We measured the change in conformational stability, Δ(ΔG), for hydrophobic mutants of four proteins: villin headpiece subdomain (VHP) with 36 residues, a surface protein from Borrelia burgdorferi (VlsE) with 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa and T1. We compared our results with those of previous studies and reached the following conclusions: (1) Hydrophobic interactions contribute less to the stability of a small protein, VHP (0.6 ± 0.3 kcal/mol per -CH₂- group), than to the stability of a large protein, VlsE (1.6 ± 0.3 kcal/mol per -CH₂- group). (2) Hydrophobic interactions make the major contribution to the stability of VHP (40 kcal/mol) and the major contributors are (in kilocalories per mole) Phe18 (3.9), Met13 (3.1), Phe7 (2.9), Phe11 (2.7), and Leu21 (2.7). (3) Based on the Δ(ΔG) values for 148 hydrophobic mutants in 13 proteins, burying a -CH₂- group on folding contributes, on average, 1.1 ± 0.5 kcal/mol to protein stability. (4) The experimental Δ(ΔG) values for aliphatic side chains (Ala, Val, Ile, and Leu) are in good agreement with their ΔG_tr values from water to cyclohexane. (5) For 22 proteins with 36 to 534 residues, hydrophobic interactions contribute 60 ± 4% and hydrogen bonds contribute 40 ± 4% to protein stability. (6) Conformational entropy contributes about 2.4 kcal/mol per residue to protein instability. The globular conformation of proteins is stabilized predominantly by hydrophobic interactions.     

Binding of 1-substituted carbazolyl-3,4-dihydro-β-carbolines with DNA: Molecular dynamics simulation and MM-GBSA analysis

Molecular Mechanics-Generalized Born-Solvent Accessibility free energy calculations were used to analyse DNA binding affinity of 1-substituted carbazolyl-3,4-dihydro-β-carboline molecules. In this study, DNA structure with sequence of d(CGATCG)2 was used for simulations. 15 ns molecular dynamics simulations of the studied complexes were performed. The calculated free energy was compared with experimental antitumor activity (IC₅₀). The predicted free energies decreased with the increase of IC₅₀ values. It was shown that molecules 1–6 bind to DNA via intercalation mode, while molecules 7–9 bind through groove binding mode. Also, it was found that the vdW energy term (ΔE_vdW) and the non-polar desolvation energy (ΔG_SA) are the favorable terms for binding energy, whereas net electrostatic energies (ΔE_ele + ΔG_GB) and conformational entropy energy (TΔS) are unfavorable ones.