Optimal specificity and function for flexible biomolecular recognition

Biomolecular associations often accompanied by large conformational changes, sometimes folding and unfolding. By exploring an exactly solvable model, we constructed the free energy landscape and established a general framework for studying the biomolecular flexible binding process. We derived an optimal criterion for the specificity and function for flexible biomolecular binding where the binding and conformational folding are coupled.

Large-area nanopatterning of self-assembled monolayers of alkanethiolates by interferometric lithography.

We demonstrate that interferometric lithography provides a fast, simple approach to the production of patterns in self-assembled monolayers (SAMs) with high resolution over square centimeter areas. As a proof of principle, two-beam interference patterns, formed using light from a frequency-doubled argon ion laser (244 nm), were used to pattern methyl-terminated SAMs on gold, facilitating the introduction of hydroxyl-terminated adsorbates and yielding patterns of surface free energy with a pitch of ca. 200 nm. The photopatterning of SAMs on Pd has been demonstrated for the first time, with interferometric exposure yielding patterns of surface free energy with similar features sizes to those obtained on gold. Gold nanostructures were formed by exposing SAMs to UV interference patterns and then immersing the samples in an ethanolic solution of mercaptoethylamine, which etched the metal substrate in exposed areas while unoxidized thiols acted as a resist and protected the metal from dissolution. Macroscopically extended gold nanowires were fabricated using single exposures and arrays of 66 nm gold dots at 180 nm centers were formed using orthogonal exposures in a fast, simple process. Exposure of oligo(ethylene glycol)-terminated SAMs to UV light caused photodegradation of the protein-resistant tail groups in a substrate-independent process. In contrast to many protein patterning methods, which utilize multiple steps to control surface binding, this single step process introduced aldehyde functional groups to the SAM surface at exposures as low as 0.3 J cm(-2), significantly less than the exposure required for oxidation of the thiol headgroup. Although interferometric methods rely upon a continuous gradient of exposure, it was possible to fabricate well-defined protein nanostructures by the introduction of aldehyde groups and removal of protein resistance in nanoscopic regions. Macroscopically extended, nanostructured assemblies of streptavidin were formed. Retention of functionality in the patterned materials was demonstrated by binding of biotinylated proteins.

Electrostatic Energetics of Bacillus subtilis Ribonuclease P Protein Determined by Nuclear Magnetic Resonance-Based Histidine pKa Measurements.

The pKa values of ionizable groups in proteins report the free energy of site-specific proton binding and provide a direct means of studying pH-dependent stability. We measured histidine pKa values (H3, H22, and H105) in the unfolded (U), intermediate (I), and sulfate-bound folded (F) states of RNase P protein, using an efficient and accurate nuclear magnetic resonance-monitored titration approach that utilizes internal reference compounds and a parametric fitting method. The three histidines in the sulfate-bound folded protein have pKa values depressed by 0.21 ± 0.01, 0.49 ± 0.01, and 1.00 ± 0.01 units, respectively, relative to that of the model compound N-acetyl-l-histidine methylamide. In the unliganded and unfolded protein, the pKa values are depressed relative to that of the model compound by 0.73 ± 0.02, 0.45 ± 0.02, and 0.68 ± 0.02 units, respectively. Above pH 5.5, H22 displays a separate resonance, which we have assigned to I, whose apparent pKa value is depressed by 1.03 ± 0.25 units, which is ∼0.5 units more than in either U or F. The depressed pKa values we observe are consistent with repulsive interactions between protonated histidine side chains and the net positive charge of the protein. However, the pKa differences between F and U are small for all three histidines, and they have little ionic strength dependence in F. Taken together, these observations suggest that unfavorable electrostatics alone do not account for the fact that RNase P protein is intrinsically unfolded in the absence of ligand. Multiple factors encoded in the P protein sequence account for its IUP property, which may play an important role in its function.

Differential scanning calorimetry in life science: thermodynamics, stability, molecular recognition and application in drug design

All biological phenomena depend on molecular recognition, which is either intermolecular like in ligand binding to a macromolecule or intramolecular like in protein folding. As a result, understanding the relationship between the structure of proteins and the energetics of their stability and binding with others (bio)molecules is a very interesting point in biochemistry and biotechnology. It is essential to the engineering of stable proteins and to the structure-based design of pharmaceutical ligands. The parameter generally used to characterize the stability of a system (the folded and unfolded state of the protein for example) is the equilibrium constant (K) or the free energy (deltaG(o)), which is the sum of enthalpic (deltaH(o)) and entropic (deltaS(o)) terms. These parameters are temperature dependent through the heat capacity change (deltaCp). The thermodynamic parameters deltaH(o) and deltaCp can be derived from spectroscopic experiments, using the van't Hoff method, or measured directly using calorimetry. Along with isothermal titration calorimetry (ITC), differential scanning calorimetry (DSC) is a powerful method, less described than ITC, for measuring directly the thermodynamic parameters which characterize biomolecules. In this article, we summarize the principal thermodynamics parameters, describe the DSC approach and review some systems to which it has been applied. DSC is much used for the study of the stability and the folding of biomolecules, but it can also be applied in order to understand biomolecular interactions and can thus be an interesting technique in the process of drug design.

Protection of DNA against low-energy electrons by amino acids: A first-principles molecular dynamics study

Using first-principles molecular dynamics simulations, we have investigated the notion that amino acids can play a protective role when DNA is exposed to excess electrons produced by ionizing radiation. In this study we focus on the interaction of glycine with the DNA nucleobase thymine. We studied thymine-glycine dimers and a condensed phase model consisting of one thymine molecule solvated in amorphous glycine. Our results show that the amino acid acts as a protective agent for the nucleobase in two ways. If the excess electron is initially captured by the thymine, then a proton is transferred in a barrier-less way from a neighboring hydrogen-bonded glycine. This stabilizes the excess electron by reducing the net partial charge on the thymine. In the second mechanism the excess electron is captured by a glycine, which acts as a electron scavenger that prevents electron localization in DNA. Both these mechanisms introduce obstacles to further reactions of the excess electron within a DNA strand, e.g. by raising the free energy barrier associated with strand breaks.

Understanding the Interaction between Low-Energy Electrons and DNA Nucleotides in Aqueous Solution

Reactions that can damage DNA have been simulated using a combination of molecular dynamics and density functional theory. In particular, the damage caused by the attachment of a low energy electron to the nucleobase. Simulations of anionic single nucleotides of DNA in an aqueous environment that was modeled explicitly have been performed. This has allowed us to examine the role played by the water molecules that surround the DNA in radiation damage mechanisms. Our simulations show that hydrogen bonding and protonation of the nucleotide by the water can have a significant effect on the barriers to strand breaking reactions. Furthermore, these effects are not the same for all four of the bases.

Mechanism of tocopherol transfer by human α-tocopherol transfer protein (α-hTTP)

Vitamin E is a well known fat soluble chain breaking antioxidant. It is a general tenn used to describe a family of eight stereoisomers of tocopherols. Selective retention of a-tocopherol in the human circulation system is regulated by the a -Tocopherol Transfer Protein (a-TIP). Using a fluorescently labelled a-tocopherol (NBD-a-Toc) synthesized in our laboratory, a fluorescence resonance energy transfer (FRET) assay was developed to monitor the kinetics of ligand transfer by a-hTTP in lipid vesicles. Preliminary results implied that NBD-a-Toe simply diffused from 6-His-a-hTTP to acceptor membranes since the kinetics of transfer were not responsive to a variety of conditions tested. After a series of trouble shooting experiments, we identified a minor contaminant, E coli. outer membrane porin F (OmpF) that co-purified with 6-His-a-hTTP from the metal affinity column as the source of the problem. In order to completely avoid OmpF contamination, a GST -a-hTTP fusion protein was purified from a glutathione agarose column followed by an on-column thrombin digestion to remove the GST tag. We then demonstrated that a-hTTP utilizes a collisional mechanism to deliver its ligand. Furthennore, a higher rate of a-tocopherol transfer to small unilamellar vesicles (SUV s) versus large unilamellar vesicles (LUV s) indicated that transfer is sensitive to membrane curvature. These findings suggest that ahTTP mediated a-Toc transfer is dominated by the hydrophobic nature of a-hTTP and the packing density of phospholipid head groups within acceptor membranes. Based on the calculated free energy change (dG) when a protein is transferred from water to the lipid bilayer, a model was generated to predict the orientation of a-hTTP when it interacts with lipid membranes. Guided by this model, several hydrophobic residues expected to penetrate deeply into the bilayer hydrophobic core, were mutated to either aspartate or alanine. Utilizing dual polarization interferometry and size exclusion vesicle binding assays, we identified the key residues for membrane binding to be F 165, F 169 and 1202. In addition, the rates of ligand transfer of the u-TTP mutants were directly correlated to their membrane binding capabilities, indicating that membrane binding was likely the rate limiting step in u-TTP mediated transfer of u-Toc. The propensity of u-TTP for highly curved membrane provides a connection to its colocalization with u-Toc in late endosomes.

Design and study of self-assembled functional organic and hybrid systems for biological applications

The focus of self-assembly as a strategy for the synthesis has been confined largely to molecules, because of the importance of manipulating the structure of matter at the molecular scale. We have investigated the influence of temperature and pH, in addition to the concentration of the capping agent used for the formation of the nano-bio conjugates. For example, the formation of the narrower size distribution of the nanoparticles was observed with the increase in the concentration of the protein, which supports the fact that γ-globulin acts both as a controller of nucleation as well as stabiliser. As analyzed through various photophysical, biophysical and microscopic techniques such as TEM, AFM, C-AFM, SEM, DLS, OPM, CD and FTIR, we observed that the initial photoactivation of γ-globulin at pH 12 for 3 h resulted in small protein fibres of ca. Further irradiation for 24 h, led to the formation of selfassembled long fibres of the protein of ca. 5-6 nm and observation of surface plasmon resonance band at around 520 nm with the concomitant quenching of luminescence intensity at 680 nm. The observation of light triggered self-assembly of the protein and its effect on controlling the fate of the anchored nanoparticles can be compared with the naturally occurring process such as photomorphogenesis.Furthermore,our approach offers a way to understand the role played by the self-assembly of the protein in ordering and knock out of the metal nanoparticles and also in the design of nano-biohybrid materials for medicinal and optoelectronic applications. Investigation of the potential applications of NIR absorbing and water soluble squaraine dyes 1-3 for protein labeling and anti-amyloid agents forms the subject matter of the third chapter of the thesis. The study of their interactions with various proteins revealed that 1-3 showed unique interactions towards serum albumins as well as lysozyme. 69%, 71% and 49% in the absorption spectra as well as significant quenching in the fluorescence intensity of the dyes 1-3, respectively. Half-reciprocal analysis of the absorption data and isothermal titration calorimetric (ITC) analysis of the titration experiments gave a 1:1 stoichiometry for the complexes formed between the lysozyme and squaraine dyes with association constants (Kass) in the range 104-105 M-1. We have determined the changes in the free energy (ΔG) for the complex formation and the values are found to be -30.78, -32.31 and -28.58 kJmol-1, respectively for the dyes 1, 2 and 3. Furthermore, we have observed a strong induced CD (ICD) signal corresponding to the squaraine chromophore in the case of the halogenated squaraine dyes 2 and 3 at 636 and 637 nm confirming the complex formation in these cases. To understand the nature of interaction of the squaraine dyes 1-3 with lysozyme, we have investigated the interaction of dyes 1-3 with different amino acids. These results indicated that the dyes 1-3 showed significant interactions with cysteine and glutamic acid which are present in the side chains of lysozyme. In addition the temperature dependent studies have revealed that the interaction of the dye and the lysozyme are irreversible. Furthermore, we have investigated the interactions of these NIR dyes 1-3 with β- amyloid fibres derived from lysozyme to evaluate their potential as inhibitors of this biologically important protein aggregation. These β-amyloid fibrils were insoluble protein aggregates that have been associated with a range of neurodegenerative diseases, including Huntington, Alzheimer’s, Parkinson’s, and Creutzfeldt-Jakob diseases. We have synthesized amyloid fibres from lysozyme through its incubation in acidic solution below pH 4 and by allowing to form amyloid fibres at elevated temperature. To quantify the binding affinities of the squaraine dyes 1-3 with β-amyloids, we have carried out the isothermal titration calorimetric (ITC) measurements. The association constants were determined and are found to be 1.2 × 105, 3.6× 105 and 3.2 × 105 M-1 for the dyes, 1-3, respectively. To gain more insights into the amyloid inhibiting nature of the squaraine dyes under investigations, we have carried out thioflavin assay, CD, isothermal titration calorimetry and microscopic analysis. The addition of the dyes 1-3 (5μM) led to the complete quenching in the apparent thioflavin fluorescence, thereby indicating the destabilization of β-amyloid fibres in the presence of the squaraine dyes. Further, the inhibition of the amyloid fibres by the squaraine dyes 1-3, has been evidenced though the DLS, TEM AFM and SAED, wherein we observed the complete destabilization of the amyloid fibre and transformation of the fibre into spherical particles of ca. These results demonstrate the fact that the squaraine dyes 1-3 can act as protein labeling agents as well as the inhibitors of the protein amyloidogenesis. The last chapter of the thesis describes the synthesis and investigation of selfassembly as well as bio-imaging aspects of a few novel tetraphenylethene conjugates 4-6.Expectedly, these conjugates showed significant solvatochromism and exhibited a hypsochromic shift (negative solvatochromism) as the solvent polarity increased, and these observations were justified though theoretical studies employing the B3LYP/6-31g method. We have investigated the self-assembly properties of these D-A conjugates though variation in the percentage of water in acetonitrile solution due to the formation of nanoaggregates. Further the contour map of the observed fluorescence intensity as a function of the fluorescence excitation and emission wavelength confirmed the formation of J-type aggregates in these cases. To have a better understanding of the type of self-assemblies formed from the TPE conjugates 4-6, we have carried out the morphological analysis through various microscopic techniques such as DLS, SEM and TEM. 70%, we observed rod shape architectures having ~ 780 nm in diameter and ~ 12 μM in length as evidenced through TEM and SEM analysis. We have made similar observations with the dodecyl conjugate 5 at ca. 70% and 50% water/acetonitrile mixtures, the aggregates formed from 4 and 5 were found to be highly crystalline and such structures were transformed to amorphous nature as the water fraction was increased to 99%. To evaluate the potential of the conjugate as bio-imaging agents, we have carried out their in vitro cytotoxicity and cellular uptake studies though MTT assay, flow cytometric and confocal laser scanning microscopic techniques. Thus nanoparticle of these conjugates which exhibited efficient emission, large stoke shift, good stability, biocompatibility and excellent cellular imaging properties can have potential applications for tracking cells as well as in cell-based therapies. In summary we have synthesized novel functional organic chromophores and have studied systematic investigation of self-assembly of these synthetic and biological building blocks under a variety of conditions. The investigation of interaction of water soluble NIR squaraine dyes with lysozyme indicates that these dyes can act as the protein labeling agents and the efficiency of inhibition of β-amyloid indicate, thereby their potential as anti-amyloid agents.

Resolving the energy paradox of chaperone/usher-mediated fibre assembly

Periplasmic chaperone/usher machineries are used for assembly of filamentous adhesion organelles of Gram-negative pathogens in a process that has been suggested to be driven by folding energy. Structures of mutant chaperone-subunit complexes revealed a final folding transition (condensation of the subunit hydrophobic core) on the release of organelle subunit from the chaperone-subunit pre-assembly complex and incorporation into the final fibre structure. However, in view of the large interface between chaperone and subunit in the pre-assembly complex and the reported stability of this complex, it is difficult to understand how final folding could release sufficient energy to drive assembly. In the present paper, we show the X-ray structure for a native chaperone-fibre complex that, together with thermodynamic data, shows that the final folding step is indeed an essential component of the assembly process. We show that completion of the hydrophobic core and incorporation into the fibre results in an exceptionally stable module, whereas the chaperone-subunit preassembly complex is greatly destabilized by the high-energy conformation of the bound subunit. This difference in stabilities creates a free energy potential that drives fibre formation.

Binding studies of a protonated dioxatetraazamacrocycle with carboxylate substrates

The tetraprotonated form of the dioxatetraazamacrocycle, 6,19-dioxa-3,9,16,22-tetraaza[22.2.2.2(11,14)]-triaconta-1(26),11,13,24, 27,29-hexaene, (H4L1)(4+), was used as the receptor for binding studies with carboxylate anionic substrates of different shapes, sizes, and charges [succinate (suc(2-)), cyclo- hexanetricarboxylate (cta(3-)), phthalate (ph(2-)), isophthalate (iph(2-)), terephthalate (tph(2-)), and benezenetricarboxylate (btc(3-))]. Association constants were determined by potentiometry in aqueous solution at 298.2 K and 0.10 M KCl and by H-1 NMR titration in D2O. The strongest association was found for the btc3- anion at 5-7 pH region. From both techniques it was possible to establish the binding preference trend of the receptor for the different substrates, and the H-1 NMR spectroscopy gave important suggestions about the type of interactions between partners and the location of the substrates in the supramolecular entities formed. The effective binding constants at pH 6 follow the order: btc(3-)>iph(2-)>cta(3-) =ph(2-)>tph(2-)>suc(2-). All the studies suggest that the anionic substrates bind to the receptor via N-H center dot center dot center dot O = C hydrogen bonds and electrostatic interactions, and the aromatic substrates can also establish pi-pi stacking interactions. The crystal structures of (H4L1)(4+) and its supramolecular assemblies with ph(2-) and tph(2-) were determined by X-ray diffraction. The last two structures showed that the association process in solid state occurs via multiple N-H center dot center dot center dot O = C hydrogen bonds with the anionic substrate located outside the macrocyclic cavity of the receptor. Molecular dynamics simulations carried out for the association of (H4L1)(4+) with tph(2-) and btC(3-) in water solution established at atomic level the existence of all interactions suggested by the experimental studies, which act cooperatively in the binding process. Furthermore, the binding free energies were estimated and the values are in agreement with the experimental ones, indicating that the binding of these two anionic substrates occurs into the receptor cavity. However, the tph(2-) has also propensity to leave the macrocyclic cavity and its molecular recognition can also happen at the top of the receptor. (C) 2008 Elsevier Ltd. All rights reserved.

Selectivity of bis(calix[4]diquinone) ionophores towards metal ions in solvent dimethylsulfoxide: a molecular mechanics and molecular dynamics study

Molecular modelling studies have been carried out on two bis(calix[4]diqu(inone) ionophores, each created from two (calix[4]diquinone)arenes bridged at their bottom rims via alkyl chains (CH2)(n), 1: n = 3, 2; n = 4, in order to understand the reported selectivity of these ligands towards different sized metal ions such as Na+, K+, Rb+, and Cs+ in dmso solution. Conformational. analyses have been carried out which show that in the lowest energy conformations of the two macrocycles, the individual calix[4]diquinones exhibit a combination of partial cone, 1,3-alternate and cone conformations. The interactions of these alkali metals with the macrocycles have been studied in the gas phase and in a periodic box of solvent dmso by molecular mechanics and molecular dynamics calculations. Molecular mechanics calculations have been carried out on the mode of entry of the ions into the macrocycles and suggest that this is likely to occur from the side of the central cavity, rather than through the main axis of the calix[4]diquinones. There are energy barriers of ca. 19 kcal mol(-1) for this entry path in the gas phase, but in solution no energy barrier is found. Molecular dynamics simulations show that in both 1 and 2, though particularly in the latter macrocycle, one or two solvent molecules are bonded to the metal throughout the course of the simulation, often to the exclusion, of one or more of the ether oxygen atoms. By contrast the carbonyl oxygen atoms remain bonded to the metal atoms throughout with bond lengths that remain significantly less than those to the ether oxygen atoms. Free energy perturbation studies have been carried out in dmso and indicate that for 1, the selectivity follows the order Rb+ approximate to K+ > Cs+ >> Na+, which is partially in agreement with the experimental results. The energy differences are small and indeed the ratio between stability constants found for Cs+ and K+ complexes is only 0.60, showing that 1 has only a slight preference for K+. For the larger receptor 2, which is better suited to metal complexation, the binding affinity follows the pattern Cs+ >> Rb+ >> K+ >> Na+, with energy differences of 5.75, 2.61, 2.78 kcal mol(-1) which is perfectly consistent with experimental results.

Homoleptic copper(II) and copper(I) complexes of 5,6-dihydro-5,6-epoxy-1,10-phenanthroline. Six-coordinate copper(I) in solution

Reaction of 5,6-dihydro-5,6-epoxy-1,10-phenanthroline (L) with Cu(ClO(4))(2)center dot 6H(2)O in methanol in 3:1 M ratio at room temperature yields light green [CuL(3)](ClO(4))(2)center dot H(2)O (1). The X-ray crystal structure of the hemi acetonitrile solvate [CuL(3)](ClO(4))(2)center dot 0.5CH(3)CN has been determined which shows Jahn-Teller distortion in the CuN(6) core present in the cation [CuL(3)](2+). Complex 1 gives an axial EPR spectrum in acetonitrile-toluene glass with g(parallel to) = 2.262 (A(parallel to) = 169 x 10 (4) cm (1)) and g(perpendicular to) = 2.069. The Cu(II/I) potential in 1 in CH(2)Cl(2) at a glassy carbon electrode is 0.32 V versus NHE. This potential does not change with the addition of extra L in the medium implicating generation of a six-coordinate copper(I) species [CuL(3)](+) in solution. B3LYP/LanL2DZ calculations show that the six Cu-N bond distances in [CuL(3)](+) are 2.33, 2.25, 2.32, 2.25, 2.28 and 2.25 angstrom while the ideal Cu(I)-N bond length in a symmetric Cu(I)N(6) moiety is estimated as 2.25 angstrom. Reaction of L with Cu(CH(3)CN)(4)ClO(4) in dehydrated methanol at room temperature even in 4:1 M proportion yields [CuL(2)]ClO(4) (2). Its (1)H NMR spectrum indicates that the metal in [CuL(2)](+) is tetrahedral. The Cu(II/I) potential in 2 is found to be 0.68 V versus NHE in CH(2)Cl(2) at a glassy carbon electrode. In presence of excess L, 2 yields the cyclic voltammogram of 1. From (1)H NMR titration, the free energy of binding of L to [CuL(2)](+) to produce [CuL(3)](+) in CD(2)Cl(2) at 298 K is estimated as -11.7 (+/-0.2) kJ mol (1).

Electronic Structure and Spectro-Structural Correlations of Fe(III)Zn(II) Biomimetics for Purple Acid Phosphatases: Relevance to DNA Cleavage and Cytotoxic Activity

Purple acid phosphatases (PAPs) are a group of metallohydrolases that contain a dinuclear Fe(II)M(II) center (M(II) = Fe, Mn, Zn) in the active site and are able to catalyze the hydrolysis of a variety of phosphoric acid esters. The dinuclear complex [(H(2)O)Fe(III)(mu-OH)Zn(II)(L-H)](CIO(4))(2) (2) with the ligand 2-[N-bis(2-pyridylmethyl)aminomethyl]-4-methyl-6-[N-(2-pyridylmethyl)(2-hydroxybenzyl) aminomethyl]phenol (H(2)L-H) has recently been prepared and is found to closely mimic the coordination environment of the Fe(III)Zn(II) active site found in red kidney bean PAP (Neves et al. J. Am. Chem. Soc. 2007, 129, 7486). The biomimetic shows significant catalytic activity in hydrolytic reactions. By using a variety of structural, spectroscopic, and computational techniques the electronic structure of the Fe(III) center of this biomimetic complex was determined. In the solid state the electronic ground state reflects the rhombically distorted Fe(III)N(2)O(4) octahedron with a dominant tetragonal compression align ad along the mu-OH-Fe-O(phenolate) direction. To probe the role of the Fe-O(phenolate) bond, the phenolate moiety was modified to contain electron-donating or -withdrawing groups (-CH(3), -H, -Br, -NO(2)) in the 5-position. Tie effects of the substituents on the electronic properties of the biomimetic complexes were studied with a range of experimental and computational techniques. This study establishes benchmarks against accurate crystallographic struck ral information using spectroscopic techniques that are not restricted to single crystals. Kinetic studies on the hydrolysis reaction revealed that the phosphodiesterase activity increases in the order -NO(2)<- Br <- H <- CH(3) when 2,4-bis(dinitrophenyl)phosphate (2,4-bdnpp) was used as substrate, and a linear free energy relationship is found when log(k(cat)/k(0)) is plotted against the Hammett parameter a. However, nuclease activity measurements in the cleavage of double stranded DNA showed that the complexes containing the electron-withdrawing -NO(2) and electron-donating CH3 groups are the most active while the cytotoxic activity of the biomimetics on leukemia and lung tumoral cells is highest for complexes with electron-donating groups.

Gaining ligand selectivity in thyroid hormone receptors via entropy

Nuclear receptors are important targets for pharmaceuticals, but similarities between family members cause difficulties in obtaining highly selective compounds. Synthetic ligands that are selective for thyroid hormone (TH) receptor beta (TR beta) vs. TR alpha reduce cholesterol and fat without effects on heart rate; thus, it is important to understand TR beta-selective binding. Binding of 3 selective ligands (GC-1, KB141, and GC-24) is characterized at the atomic level; preferential binding depends on a nonconserved residue (Asn-331 beta) in the TR beta ligand-binding cavity (LBC), and GC-24 gains extra selectivity from insertion of a bulky side group into an extension of the LBC that only opens up with this ligand. Here we report that the natural TH 3,5,3`-triodothyroacetic acid (Triac) exhibits a previously unrecognized mechanism of TR beta selectivity. TR x-ray structures reveal better fit of ligand with the TR alpha LBC. The TR beta LBC, however, expands relative to TR alpha in the presence of Triac (549 angstrom(3) vs. 461 angstrom(3)), and molecular dynamics simulations reveal that water occupies the extra space. Increased solvation compensates for weaker interactions of ligand with TR beta and permits greater flexibility of the Triac carboxylate group in TR beta than in TR alpha. We propose that this effect results in lower entropic restraint and decreases free energy of interactions between Triac and TR beta, explaining subtype-selective binding. Similar effects could potentially be exploited in nuclear receptor drug design.

Ab Initio Analysis of Monomers and Dimers of Trialkylphosphine Oxides: Structural and Thermodynamic Stability

Structural and thermodynamic stabilities of monomers and dimers of trialkylphosphine oxides (TRPO) were Studied using quantum chemistry calculations. Density functional theory calculations were carried Out and the structures Of four TRPO have been determined: TMPO (methyl; R = CH(3)), TEPO (ethyl; R = CH(3)CH(2)), TBPO (n-butyl; R = CH(3)(CH(2))(3)), and TOPO (n-octyl; R = CH(3)(CH(2))(7)). TRPO homodimers were investigated considering two isomeric possibilities for each dimer. Relative binding energies and the enthalpic and entropic contributions to the Gibbs free energy were Calculated for all dimers. The formation of dimers from the individual monomeric TRPO species as a function of temperature was also analyzed. (C) 2008 Wiley Periodicals, Inc. Int J Quantum Chem 109: 250-258, 2009