Flanking DNA-sequences contribute to the specific binding of cI-repressor and OR1.

The binding of cI-repressor to a series of mutant operators containing OR1 of the right operator of bacteriophage lambda was investigated. Sites OR2 and/or OR3 were inactivated by either point or deletion mutations. The free energy of binding repressor to OR1 in the wildtype operator, delta G1, is -13.7 +/- 0.3 kcal/mol. delta G1 determined for an OR2- operator created by a single point mutation in OR2 is -13.6 +/- 0.2 kcal/mol. In contrast, delta G1 for the binding of repressor to a cloned synthetic OR1 operator containing only 24 bp of lambda sequence is -12.2 +/- 0.1 kcal/mol. When sequence 5' to OR1 is present, the binding affinity increases to -13.0 +/- 0.1 kcal/mol. In addition, the proximity of OR1 to a fragment-end decreases delta G1 from -13.7 to -12.3 +/- 0.1 kcal/mol. These results suggest that the DNA sequence outside the 17 bp OR1 binding-site contributes to the specific binding of cI-repressor.     

Cooperative binding of R17 coat protein to RNA.

The binding of the R17 coat protein to synthetic RNAs containing one or two coat protein binding sites was characterized by using nitrocellulose filter and gel-retention assays. RNAs with two available sites bound coat protein in a cooperative manner, resulting in a higher affinity and reduced sensitivity to pH, ionic strength, and temperature when compared with RNAs containing only a single site. The cooperativity can contribute up to -5 kcal/mol to the overall binding affinity with the greatest cooperativity found at low pH, high ionic strength, and high temperatures. Similar solution properties for the encapsidation of the related fr and f2 phage suggest that the cooperativity is due to favorable interactions between the two coat proteins bound to the RNA. This system therefore resembles an intermediate state of phage assembly. No cooperative binding was observed for RNAs containing a single site and a 5' or 3' extension of nonspecific sequence, indicating that R17 coat protein has a very low nonspecific binding affinity. Unexpectedly weak binding was observed for several RNAs due to the presence of alternative conformational states of the RNA.     

Investigation of phosphotyrosine recognition by the SH2 domain of the Src kinase.

The binding of tyrosine phosphorylated targets by SH2 domains is required for propagation of many cellular signals in higher eukaryotes; however, the determinants of phosphotyrosine (pTyr) recognition by SH2 domains are not well understood. In order to identify the attributes of pTyr required for high affinity interaction with SH2 domains, the binding of the SH2 domain of the Src kinase (Src SH2 domain) to a dephosphorylated peptide, a phosphoserine-containing peptide, and the amino acid pTyr was studied using titration calorimetry and compared with the binding of a high affinity tyrosyl phosphopeptide. The dephosphorylated peptide and the phosphoserine containing peptide both bind extremely weakly to the Src SH2 domain (DeltaGo (dephosphorylated)=-3.6 kcal/mol, DeltaGo (phosphoserine) >-3.7 kcal/mol); however, the DeltaGo value of pTyr binding is more favorable (-4.7 kcal/mol, or 50 % of the entire binding free energy of a high affinity tyrosyl phosphopeptide). These results indicate that both the phosphate and the tyrosine ring of the pTyr are critical determinants of high affinity binding. Alanine mutagenesis was also used to evaluate the energetic contribution to binding of ten residues located in the pTyr-binding site. Mutation of the strictly conserved Arg betaB5 resulted in a large increase in DeltaGo (DeltaDeltaGo=3.2 kcal/mol) while elimination of the other examined residues each resulted in a significantly smaller (DeltaDeltaGo<1.4 kcal/mol) reduction in affinity, indicating that Arg betaB5 is the single most important determinant of pTyr recognition. However, mutation of Cys betaC3, a residue unique to the Src SH2 domain, surprisingly increased affinity by eightfold (DeltaDeltaGo=-1.1 kcal/mol). Using a double mutant cycle analysis, it was revealed that residues of the pTyr-binding pocket are not coupled to the peptide residues C-terminal to the pTyr. In addition, comparison of each residue's DeltaDeltaGo value upon mutation with that residue's sequence conservation among SH2 domains revealed only a modest correlation between a residue's energetic contribution to pTyr recognition and its conservation throughout evolution. The results of this investigation highlight the importance of a single critical interaction, the buried ionic bond between the phosphate of the pTyr and Arg betaB5 of the SH2 domain, driving the binding of SH2 domains to tyrosine phosphorylated targets.     

Theoretical study of the ligand-CYP2B4 complexes: effect of structure on binding free energies and heme spin state

The molecular origins of temperature-dependent ligand-binding affinities and ligand-induced heme spin state conversion have been investigated using free energy analysis and DFT calculations for substrates and inhibitors of cytochrome P450 2B4 (CYP2B4), employing models of CYP2B4 based on CYP2C5(3LVdH)/CYP2C9 crystal structures, and the results compared with experiment. DFT calculations indicate that large heme-ligand interactions (ca. -15 kcal/mol) are required for inducing a high to low spin heme transition, which is correlated with large molecular electrostatic potentials (approximately -45 kcal/mol) at the ligand heteroatom. While type II ligands often contain oxygen and nitrogen heteroatoms that ligate heme iron, DFT results indicate that BP and MF heme complexes, with weak substrate-heme interactions (ca. -2 kcal/mol), and modest MEPS minima (>-35 kcal/mol) are high spin. In contrast, heme complexes of the CYP2B4 inhibitor, 4PI, the product of benzphetamine metabolism, DMBP, and water are low spin, have substantial heme-ligand interaction energies (<-15 kcal/mol) and deep MEPS minima (<-45 kcal/mol) near their heteroatoms. MMPBSA analysis of MD trajectories were made to estimate binding free energies of these ligands at the heme binding site of CYP2B4. In order to initially assess the realism of this approach, the binding free energy of 4PI inhibitor was computed and found to be a reasonable agreement with experiment: -7.7 kcal/mol [-7.2 kcal/mol (experiment)]. BP was determined to be a good substrate [-6.3 kcal/mol (with heme-ligand water), -7.3 kcal/mol (without ligand water)/-5.8 kcal/mol (experiment)], whereas the binding of MF was negligible, with only marginal binding binding free energy of -1.7 kcal/mol with 2-MF bound [-3.8 kcal/mol (experiment)], both with and without retained heme-ligand water. Analysis of the free energy components reveal that hydrophobic/nonpolar contributions account for approximately 90% of the total binding free energy of these substrates and are the source of their differential and temperature-dependent CYP2B4 binding. The results indicate the underlying origins of the experimentally observed differential binding affinities of BP and MF, and indicate the plausibility of the use of models derived from moderate sequence identity templates in conjunction with approximate free energy methods in the estimation of ligand-P450 binding affinities.     

Site-bound water and the shortcomings of a less than perfect transition state analogue

Kinetic measurements have shown that substantial enthalpy changes accompany substrate binding by cytidine deaminase, increasing markedly as the reaction proceeds from the ground state (1/K(m), DeltaH = -13 kcal/mol) to the transition state (1/K(tx), DeltaH = -20 kcal/mol) [Snider, M. J., et al. (2000) Biochemistry 39, 9746-9753]. In the present work, we determined the thermodynamic changes associated with the equilibrium binding of inhibitors by cytidine deaminase by isothermal titration calorimetry and van't Hoff analysis of the temperature dependence of their inhibition constants. The results indicate that the binding of the transition state analogue 3,4-dihydrouridine DeltaH = -21 kcal/mol), like that of the transition state itself (DeltaH = -20 kcal/mol), is associated with a large favorable change in enthalpy. The significantly smaller enthalpy change that accompanies the binding of 3,4-dihydrozebularine (DeltaH = -10 kcal/mol), an analogue of 3,4-dihydrouridine in which a hydrogen atom replaces this inhibitor's 4-OH group, is consistent with the view that polar interactions with the substrate at the site of its chemical transformation play a critical role in reducing the enthalpy of activation for substrate hydrolysis. The entropic shortcomings of 3,4-dihydrouridine, in capturing all of the free energy involved in binding the actual transition state, may arise from its inability to displace a water molecule that occupies the binding site normally occupied by product ammonia.     

Interacting binding sites of isoleucyl-tRNA synthetase from Escherichia coli studied by equilibrium partition.

The binding of tRNAIIe to isoleucyl-tRNA synthetase in the presence of isoleucine or ATP was investigated using the equilibrium partition method. Isoleucine decreased the affinity of tRNAIIe for the enzyme by a factor of about 5. For the free standard energy of interaction a value of about 1 kcal/mol (4.2 kJ/mol) was calculated. ATP exhibits qualitatively the same effect as isoleucine. A binding of two molecules isoleucine per molecule of enzyme could not be demonstrated even in the presence of ATP and pyrophosphatase.     

A transition state in pieces: major contributions of entropic effects to ligand binding by adenosine deaminase.

Nebularine undergoes hydration at the active site of adenosine deaminase, in a reaction analogous to a partial reaction in the displacement of ammonia from adenosine by water, to generate an inhibitory complex that captures much of the binding affinity expected of an ideal transition-state analogue. Enzyme affinities of several compounds related to nebularine 1,6-hydrate, and to its stable analog 2'-deoxycoformycin, were compared in an effort to identify the structural origins of strong binding. Binding of the stable transition-state analog inhibitor 2'-deoxycoformycin was rendered 9.8 kcal/mol less favorable by removal of substituent ribose, 9.7 kcal/mol less favorable by inversion of the 8-hydroxyl substituent of the diazepine ring, and 10.0 kcal/mol less favorable by removal of atoms 4-6 of the diazepine ring. Binding of the unstable transition-state analog nebularine hydrate was rendered at least 9.9 kcal/mol less favorable by removal of the 6-hydroxyl group and 10.2 kcal/mol less favorable by removal of atoms 1-3 of the pyrimidine ring. In each case, the enzyme exhibited only modest affinity (Kd greater than or equal to 10(-2) M) for the "missing piece", indicating that incorporation of 2 binding determinants within a single molecule permits an additional 7-12 kcal/mol of intrinsic binding energy to be manifested as observed binding energy. These results are consistent with earlier indications that adenosine deaminase may use 10.5 kcal/mol of the intrinsic free energy of binding of the two substrates to place them in positions appropriate for reaction at the active site, overcoming the unfavorable entropy change of -35 eu for the equilibrium of 1,6-hydration of purine ribonucleoside and reducing the equilibrium constant for attainment of the transition state in deamination of adenosine. Thus, adenosine deaminase may achieve up to 8 orders of magnitude of its catalytic power by converting the nonenzymatic, bimolecular, hydration reaction to a monomolecular reaction at its active site. Several new 6-substituted 1,6-dihydropurine ribonucleosides, prepared by photoaddition of formate and by low-temperature addition of organolithium reagents to a derivative of purine ribonucleoside, exhibited Ki values of 9-1400 microM against adenosine deaminase, in accord with the active site's considerable tolerance of bulky leaving groups in substrates. Inhibition by one diastereomer of 6-carboxy-1,6-dihydropurine ribonucleoside was found to be time-dependent, progressing from a weakly bound to a more strongly bound complex.     

Energetic contribution of non-essential 5' sequence to catalysis in a hepatitis delta virus ribozyme

Hepatitis delta virus (HDV) ribozymes employ multiple catalytic strategies to achieve overall rate enhancement of RNA cleavage. These strategies include general acid-base catalysis by a cytosine side chain and involvement of divalent metal ions. Here we used a trans-acting form of the antigenomic ribozyme to examine the contribution of the 5' sequence in the substrate to HDV ribozyme catalysis. The cleavage rate constants increased for substrates with 5' sequence alterations that reduced ground-state binding to the ribozyme. Quantitatively, a plot of activation free energy of chemical conversion versus Gibb's free energy of substrate binding revealed a linear relationship with a slope of -1. This relationship is consistent with a model in which components of the substrate immediately 5' to the cleavage site in the HDV ribozyme-substrate complex destabilize ground-state binding. The intrinsic binding energy derived from the ground-state destabilization could contribute up to 2 kcal/mol toward the total 8.5 kcal/mol reduction in activation free energy for RNA cleavage catalyzed by the HDV ribozyme.     

Catalysis by entropic effects: the action of cytidine deaminase on 5,6-dihydrocytidine

In neutral solution, 5,6-dihydrocytidine undergoes spontaneous deamination (k25 approximately 3.2 x 10(-5) s(-1)) much more rapidly than does cytidine (k25 approximately 3.0 x 10(-10) s(-1)), with a more favorable enthalpy of activation (DeltaDeltaH# = -8.7 kcal/mol) compensated by a less favorable entropy of activation (TDeltaDeltaS# = -1.8 kcal/mol at 25 degrees C). E. coli cytidine deaminase enhances the rate of deamination of 5,6-dihydrocytidine (kcat/k(non) = 4.4 x 10(5)) by enhancing the entropy of activation (DeltaDeltaH# = 0 kcal/mol; TDeltaDeltaS# = +7.6 kcal/mol, at 25 degrees C). Binding of the competitive inhibitor 3,4,5,6-tetrahydrouridine (THU), a stable analogue of 5,6-dihydrocytidine in the transition state for its deamination, is accompanied by a release of enthalpy (DeltaH = -7.1 kcal/mol, TDeltaDeltaS = +2.2 kcal/mol) that approaches the estimated enthalpy of binding of the actual substrate in the transition state for deamination of 5,6-dihydrocytidine (DeltaH = -8.1 kcal/mol, TDeltaDeltaS = +6.0 kcal/mol). Thus, the shortcomings of THU in capturing all of the binding affinity expected of an ideal transition-state analogue reflect a less favorable entropy of association. That difference may arise from the analogue's inability to displace a water molecule from the "leaving group site" at which ammonia is generated in the normal reaction. The effect on binding of removing the 4-OH group from the transition-state analogue THU, to form 3,4,5,6-tetrahydrozebularine (THZ) (DeltaDeltaH = -2.1 kcal/mol, TDeltaDeltaS = -4.4 kcal/mol), is mainly entropic, consistent with the inability of THZ to displace water from the "attacking group site". These results are consistent with earlier indications [Snider, M. J., and Wolfenden, R. (2001) Biochemistry 40, 11364] that site-bound water plays a prominent role in substrate activation and inhibitor binding by cytidine deaminase.     

Interaction of R17 coat protein with its RNA binding site for translational repression

The interaction between bacteriophage R17 coat protein and its RNA binding site for translational repression was studied as an example of a sequence-specific RNA-protein interaction. A nitrocellulose filter retention assay is used to demonstrate equimolar binding between the coat protein and a synthetic 21 nucleotide RNA fragment. The Kd at 2 degrees C in a buffer containing 0.19 M salt is about 1 nM. The relatively weak ionic strength dependence of Ka and a delta H = -19 kcal/mole indicates that most of the binding free energy is due to non-electrostatic interactions. Since a variety of RNAs failed to compete with the 21 nucleotide fragment for coat protein binding, the interaction appears highly sequence specific. We have synthesized more than 30 different variants of the binding site sequence in order to identify the portions of the RNA molecule which are important for protein binding. Out of the five single stranded residues examined, four were essential for protein binding whereas the fifth could be replaced by any nucleotide. One variant was found to bind better than the wild type sequence. Substitution of nucleotides which disrupted the secondary structure of the binding fragment resulted in very poor binding to the protein. These data indicated that there are several points of contact between the RNA and the protein and the correct hairpin secondary structure of the RNA is essential for protein binding.     

A compact RNA tertiary structure contains a buried backbone-K+ complex

The structure of a 58 nucleotide ribosomal RNA fragment buries several phosphate groups of a hairpin loop within a large tertiary core. During refinement of an X-ray crystal structure containing this RNA, a potassium ion was found to be contacted by six oxygen atoms from the buried phosphate groups; the ion is contained completely within the solvent-accessible surface of the RNA. The electrostatic potential at the ion chelation site is unusually large, and more than compensates for the substantial energetic penalties associated with partial dehydration of the ion and displacement of delocalized ions. The very large predicted binding free energy, approximately -30 kcal/mol, implies that the site must be occupied for the RNA to fold. These findings agree with previous studies of the ion-dependent folding of tertiary structure in this RNA, which concluded that a monovalent ion was bound in a partially dehydrated environment where Mg2+ could not easily compete for binding. By compensating the unfavorable free energy of buried phosphate groups with a chelated ion, the RNA is able to create a larger and more complex tertiary fold than would be possible otherwise.     

A 62,000 molecular weight spliceosome protein crosslinks to the intron polypyrimidine tract

Incubation in HeLa nuclear extract of a 32P-labeled 61 nucleotide-long RNA corresponding to the lariat branch site/polypyrimidine tract/3' splice site of the first intron of human beta-globin pre-mRNA led to the crosslinking of a single protein of approximately 62,000 mol. wt. (p62). p62 corresponds to a polypyrimidine tract-binding protein recently described by Garcia-Blanco et al. (Genes & Dev. 3: 1874-1886, 1989). Crosslinking of p62 to the 61 nt RNA was highly sequence specific. No p62 crosslinking was observed with a 60 nt pGEM vector RNA, a 63 nt RNA antisense to the 61-mer or a 72 nt U2 RNA sequence. p62 crosslinking to the 61 nt RNA was competed by unlabeled 61 nt RNA, by beta-globin pre-mRNA containing intron 1, and by poly(U) and poly(C), but was competed to a lesser extent or not at all by pGEM RNA, a beta-globin RNA lacking intron 1, or poly(A). Experiments with mutated RNAs revealed that neither the lariat branch site adenosine nor the 3' splice site were required for p62 crosslinking to polypyrimidine tract-containing RNA. Elimination of the polypyrimidine tract reduced p62 crosslinking, as did mutation of a polypyrimidine tract UU dinucleotide to GA. However, replacement of a pyrimidine-rich tract immediately adjacent (3') to the lariat branch site with a 57% A + G pGEM vector RNA sequence also significantly reduced p62 crosslinking, indicating the involvement of both this pyrimidine-rich region and the classical polypyrimidine tract adjacent to the 3' splice site. The sites of protein interaction were further defined by RNase H protection experiments, the results of which confirmed the patterns of p62 crosslinking to mutant RNAs. p62 crosslinking was efficiently competed by a DNA oligonucleotide having the same sequence as the 61 nt RNA, showing that p62 requires neither ribose 2' OH groups nor uracil bases for its interaction with the polypyrimidine tract. p62 was not crosslinked to double-stranded 61 nt RNA. Q-Sepharose chromatography of HeLa nuclear extract yielded an unbound fraction (QU) in which p62 was the only polypyrimidine tract-crosslinkable protein and a bound fraction (QB) in which, surprisingly, several non-p62 proteins were crosslinkable to the polypyrimidine tract RNA. Yet, when the two Q-Sepharose fractions were combined, p62 strongly out-competed the otherwise-crosslinkable QB proteins for polypyrimidine tract RNA crosslinking. This indicates that p62 may have the highest affinity and/or crosslinking efficiency for the intron polypyrimidine tract of any HeLa nuclear protein.     

In silico investigation of potential inhibitors to main protease and spike protein of SARS-CoV-2 in propolis

Docking analysis of propolis's natural compound was successfully performed against SARS-CoV-2 main protease (Mpro) and spike protein subunit 2 (S2). Initially, the propolis's protein was screened using chromatography analysis and successfully identified 22 compounds in the propolis. Four compounds were further investigated, i.e., neoblavaisoflavone, methylophiopogonone A, 3′-Methoxydaidzin, and genistin. The binding affinity of 3′-Methoxydaidzin was ?7.7 kcal/mol, which is similar to nelfinavir (control), while the others were ?7.6 kcal/mol. However, we found the key residue of Glu A:166 in the methylophiopogonone A and genistin, even though the predicted binding energy slightly higher than nelfinavir. In contrast, the predicted binding affinity of neoblavaisoflavone, methylophiopogonone A, 3′-Methoxydaidzin, and genistin against S2 were ?8.1, ?8.2, ?8.3, and ?8.3 kcal/mol, respectively, which is far below of the control (pravastatin, ?7.3 kcal/mol). Instead of conventional hydrogen bonding, the π bonding influenced the binding affinity against S2. The results reveal that this is the first report about methylophiopogonone A, 3′-Methoxydaidzin, and genistin as candidates for anti-viral agents. Those compounds can then be further explored and used as a parent backbone molecule to develop a new supplementation for preventing SARS-CoV-2 infections during COVID-19 outbreaks.     

Thermodynamic characterization of single mismatches found in naturally occurring RNA

Many naturally occurring RNA structures contain single mismatches. However, the algorithms currently used to predict RNA structure from sequence rely on a minimal set of data for single mismatches, most of which occur rather infrequently in nature. As a result, several approximations and assumptions are used to predict the stability of RNA duplexes containing the most common single mismatches. Therefore, the relative frequency of single mismatches was determined by compiling and searching a database of 955 RNA secondary structures. Thermodynamic parameters for duplex formation, derived from optical melting experiments, are reported for 28 oligoribonucleotides containing frequently occurring single mismatches. These data were then combined with previous data to construct a dataset of 64 single mismatches, including the 30 most common in the database. Because of this increase in experimental thermodynamic parameters for single mismatches that occur frequently in nature, more accurate free energy calculations have resulted. To improve the prediction of the thermodynamic parameters for duplexes containing single mismatches that have not been experimentally measured, single mismatch-specific nearest neighbor parameters were derived. The free energy of an RNA duplex containing a single mismatch that has not been thermodynamically characterized can be calculated by: DeltaG degrees 37,single mismatch = DeltaG degrees 37,mismatch nt + DeltaG degrees 37,mismatch-NN interaction + DeltaG degrees 37,AU/GU. Here, DeltaG degrees 37,mismatch is -0.4, -2.1, and -0.3 kcal/mol for A.G, G.G, and U.U mismatches, respectively; DeltaG degrees 37,mismatch-NN interaction is 0.7, -0.5, 0.4, -0.4, and -1.0 kcal/mol for 5'YRR3'/3'RRY5', 5'RYY3'/3'YYR5', 5'YYR3'/3'RYY5', 5'YRY3'/3'RYR5', and 5'RRY3'/3'YYR5' mismatch-nearest neighbor combinations, respectively, when A and G are categorized as purines (R) and C and U are categorized as pyrimidines (Y); and DeltaG degrees 37,AU/GU is a penalty of 1.2 kcal/mol for replacing a G-C base pair with either an A-U or G-U base pair. Similar predictive models were also derived for DeltaH degrees single mismatch and DeltaS degrees single mismatch. These new predictive models, in conjunction with the reported thermodynamics for frequently occurring single mismatches, should allow for more accurate calculations of the free energy of RNA duplexes containing single mismatches and, furthermore, allow for improved prediction of secondary structure from sequence.     

Cooperativity of papain-substrate interaction energies in the S2 to S2' subsites

Enzyme-substrate contacts in the hydrolysis of ester substrates by the cysteine protease papain were investigated by systematically altering backbone hydrogen-bonding and side-chain hydrophobic contacts in the substrate and determining each substrate's kinetic constants. The observed specificity energies [defined as delta delta G obs = -RT ln [(kcat/KM)first/(kcat/KM)second)]] of the substrate backbone hydrogen bonds were -2.7 kcal/mol for the P2 NH and -2.6 kcal/mol for the P1 NH when compared against substrates containing esters at those sites. The observed binding energies were -4.0 kcal/mol for the P2 Phe side chain, -1.0 kcal/mol for the P1' C=O, and -2.3 kcal/mol for the P2' NH. The latter three values probably all significantly underestimate the incremental binding energies. The P2 NH, P2 Phe side-chain, and P1 NH contacts display a strong interdependence, or cooperativity, of interaction energies that is characteristic of enzyme-substrate interactions. This interdependence arises largely from the entropic cost of forming the enzyme-substrate transition state. As favorable contacts are added successively to a substrate, the entropic penalty associated with each decreases and the free energy expressed approaches the incremental interaction energy. This is the first report of a graded cooperative effect. Elucidation of favorable enzyme-substrate contacts remote from the catalytic site will assist in the design of highly specific cysteine protease inhibitors.     

Binding mechanism of spinosine and venenatine molecules with p300 HAT enzyme: Molecular screening,molecular dynamics and free-energy analysis

The chromatin modification is regulated by the histone acetyltransferase (HAT) and histone deacetyltransferase (HDAC) enzymes; abnormal function of these enzymes leads to several malignant diseases. The inhibition of these enzymes using natural ligand molecules is an emerging technique to cure these diseases. The in vitro analysis of natural molecules, venenatine, spinosine, palmatine and taxodione are giving the best inhibition rate against p300 HAT enzyme. However, the detailed understanding of binding and the stability of these molecules with p300 HAT is not yet known. The aim of the present study is focused to determine the binding strength of the molecules from molecular dynamics simulation analysis. The docking analysis confirms that, the venenatine (−6.97 kcal/mol - conformer 8), spinosine (−6.52 kcal/mol conformer −10), palmatine (−5.72 kcal/mol conformer-3) and taxodione (−4.99 kcal/mol conformer-4) molecules form strong hydrogen bonding interactions with the key amino acid residues (Arg1410, Thr1411 and Trp1466) present in the active site of p300. In the molecular dynamics (MD) simulation, the spinosine retain these key interactions with the active site amino acid residues (Arg1410, Thr1411, and Trp1466) than venenatine and are stable throughout the simulation. The RMSD value of spinosine (0.5 to 1.3 Å) and venenatine (0.3 to 1.3 Å) are almost equal during the MD simulation. However, during the MD simulation, the intermolecular interaction between venenatine and the active site amino acid residues (Arg1410, Thr1411, and Trp1466) decreased on comparing with the spinosine-p300 interaction. The binding free energy of the spinosine (−15.30 kcal/mol) is relatively higher than the venenatine (−11.8 kcal/mol); this increment is attributed to the strong hydrogen bonding interactions of spinosine molecule with the active site amino acid residues of p300.     

Solvation and the hidden thermodynamics of a zinc finger probed by nonstandard repair of a protein crevice

The classical Zn finger contains a phenylalanine at the crux of its three architectural elements: a beta-hairpin, an alpha-helix, and a Zn(2+)-binding site. Surprisingly, phenylalanine is not required for high-affinity Zn2+ binding, but instead contributes to the specification of a precise DNA-binding surface. Substitution of phenylalanine by leucine leads to a floppy but native-like structure whose Zn affinity is maintained by marked entropy-enthalpy compensation (DeltaDeltaH -8.3 kcal/mol and -TDeltaDeltaS 7.7 kcal/mol). Phenylalanine and leucine differ in shape, size, and aromaticity. To distinguish which features correlate with dynamic stability, we have investigated a nonstandard finger containing cyclohexanylalanine at this site. The structure of the nonstandard finger is similar to that of the native domain. The cyclohexanyl ring assumes a chair conformation, and conformational fluctuations characteristic of the leucine variant are damped. Although the nonstandard finger exhibits a lower affinity for Zn2+ than does the native domain (DeltaDeltaG -1.2 kcal/mol), leucine-associated perturbations in enthalpy and entropy are almost completely attenuated (DeltaDeltaH -0.7 kcal/mol and -TDeltaDeltaS -0.5 kcal/mol). Strikingly, global changes in entropy (as inferred from calorimetry) are in each case opposite in sign from changes in configurational entropy (as inferred from NMR). This seeming paradox suggests that enthalpy-entropy compensation is dominated by solvent reorganization rather than nominal molecular properties. Together, these results demonstrate that dynamic and thermodynamic perturbations correlate with formation or repair of a solvated packing defect rather than type of physical interaction (aromatic or aliphatic) within the core.