Antithyroid drugs and their analogues: synthesis, structure, and mechanism of action

Thyroid hormones are essential for the development and differentiation of all cells of the human body. They regulate protein, fat, and carbohydrate metabolism. In this Account, we discuss the synthesis, structure, and mechanism of action of thyroid hormones and their analogues. The prohormone thyroxine (14) is synthesized on thyroglobulin by thyroid peroxidase (TPO), a heme enzyme that uses iodide and hydrogen peroxide to perform iodination and phenolic coupling reactions. The monodeiodination of T4 to 3,3',5-triiodothyronine (13) by selenium-containing deiodinases (ID-1, ID-2) is a key step in the activation of thyroid hormones. The type 3 deiodinase (ID-3) catalyzes the deactivation of thyroid hormone in a process that removes iodine selectively from the tyrosyl ring of T4 to produce 3,3',5'-triiodothyronine (rT3). Several physiological and pathological stimuli influence thyroid hormone synthesis. The overproduction of thyroid hormones leads to hyperthyroidism, which is treated by antithyroid drugs that either inhibit the thyroid hormone biosynthesis and/or decrease the conversion of T4 to T3. Antithyroid drugs are thiourea-based compounds, which indude propylthiouracil (PTU), methimazole (MM I), and carbimazole (CBZ). The thyroid gland actively concentrates these heterocyclic compounds against a concentration gradient Recently, the selenium analogues of PTU, MMI, and CBZ attracted significant attention because the selenium moiety in these compounds has a higher nucleophilicity than that of the sulfur moiety. Researchers have developed new methods for the synthesis of the selenium compounds. Several experimental and theoretical investigations revealed that the selone (C=Se) in the selenium analogues is more polarized than the thione (C=S) in the sulfur compounds, and the selones exist predominantly in their zwitterionic forms. Although the thionamide-based antithyroid drugs have been used for almost 70 years, the mechanism of their action is not completely understood. Most investigations have revealed that MMI and PTU irreversibly inhibit TPO. PTU, MTU, and their selenium analogues also inhibit ID-1, most likely by reacting with the selenenyl iodide intermediate. The good ID-1 inhibitory activity of Pill and its analogues can be ascribed to the presence of the -N(H)-C(=O)- functionality that can form hydrogen bonds with nearby amino add residues in the selenenyl sulfide state. In addition to the TPO and ID-1 inhibition, the selenium analogues are very good antioxidants. In the presence of cellular reducing agents such as GSH, these compounds catalytically reduce hydrogen peroxide. They can also efficiently scavenge peroxynitrite, a potent biological oxidant and nitrating agent.

Molecular Dynamics Perspective on the Protein Thermal Stability: A Case Study Using SAICAR Synthetase

The enzyme SAICAR synthetase ligates aspartate with CAIR (5'-phosphoribosyl-4-carboxy-5-aminoimidazole) forming SAICAR (5-amino-4-imidazole-N-succinocarboxamide ribonucleotide) in the presence of ATP. In continuation with our previous study on the thermostability of this enzyme in hyper-/thermophiles based on the structural aspects, here, we present the dynamic aspects that differentiate the mesophilic (E. coli, E. chaffeensis), thermophilic (G. kaustophilus), and hyperthermophilic (M. jannaschii, P. horikoshii) SAICAR synthetases by carrying out a total of 11 simulations. The five functional dimers from the above organisms were simulated using molecular dynamics for a period of 50 ns each at 300 K, 363 K, and an additional simulation at 333 K for the thermophilic protein. The basic features like root-mean-square deviations, root-mean-square fluctuations, surface accessibility, and radius of gyration revealed the instability of mesophiles at 363 K. Mean square displacements establish the reduced flexibility of hyper-/thermophiles at all temperatures. At the simulations time scale considered here, the long-distance networks are considerably affected in mesophilic structures at 363 K. In mesophiles, a comparatively higher number of short-lived (having less percent existence time) C alpha, hydrogen bonds, hydrophobic interactions are formed, and long-lived (with higher percentage existence time) contacts are lost. The number of time-averaged salt-bridges is at least 2-fold higher in hyperthermophiles at 363 K. The change in surface accessibility of salt-bridges at 363 K from 300 K is nearly doubled in mesophilic protein compared to proteins from other temperature classes.

Sensitivity of Water Dynamics to Biologically Significant Surfaces of Monomeric Insulin: Role of Topology and Electrostatic Interactions

In addition to the biologically active monomer of the protein insulin circulating in human blood, the molecule also exists in dimeric and hexameric forms that are used as storage. The insulin monomer contains two distinct surfaces, namely, the dimer forming surface (DFS) and the hexamer forming surface (HFS), that are specifically designed to facilitate the formation of the dimer and the hexamer, respectively. In order to characterize the structural and dynamical behavior of interfacial water molecules near these two surfaces (DFS and HFS), we performed atomistic molecular dynamics simulations of insulin with explicit water. Dynamical characterization reveals that the structural relaxation of the hydrogen bonds formed between the residues of DFS and the interfacial water molecules is faster than those formed between water and that of the HFS. Furthermore, the residence times of water molecules in the protein hydration layer for both the DFS and HFS are found to be significantly higher than those for some of the other proteins studied so far, such as HP-36 and lysozyme. In particular, we find that more structured water molecules, with higher residence times (similar to 300-500 ps), are present near HFS than those near DFS. A significant slowing down is observed in the decay of associated rotational auto time correlation functions of O-H bond vector of water in the vicinity of HFS. The surface topography and the arrangement of amino acid residues work together to organize the water molecules in the hydration layer in order to provide them with a preferred orientation. HFS having a large polar solvent accessible surface area and a convex extensive nonpolar region, drives the surrounding water molecules to acquire predominantly an outward H-atoms directed, clathrate-like structure. In contrast, near the DFS, the surrounding water molecules acquire an inward H-atoms directed orientation owing to the flat curvature of hydrophobic surface and the interrupted hydrophilic residual alignment. We have followed escape trajectory of several such quasi-bound water molecules from both the surfaces that reveal the significant differences between the two hydration layers.

Intermolecular interactions, charge-density distribution and the electrostatic properties of pyrazinamide anti-TB drug molecule: an experimental and theoretical charge-density study

An experimental charge-density analysis of pyrazinamide (a first line antitubercular drug) was performed using high-resolution X-ray diffraction data (sin theta/lambda)(max) = 1.1 angstrom(-1)] measured at 100 (2) K. The structure was solved by direct methods using SHELXS97 and refined by SHELXL97. The total electron density of the pyrazinamide molecule was modeled using the Hansen-Coppens multipole formalism implemented in the XD software. The topological properties of electron density determined from the experiment were compared with the theoretical results obtained from CRYSTAL09 at the B3LYP/6-31G** level of theory. The crystal structure was stabilized by N-H center dot center dot center dot N and N-H center dot center dot center dot O hydrogen bonds, in which the N3-H3B center dot center dot center dot N1 and N3-H3A center dot center dot center dot O1 interactions form two types of dimers in the crystal. Hirshfeld surface analysis was carried out to analyze the intermolecular interactions. The fingerprint plot reveals that the N center dot center dot center dot H and O center dot center dot center dot H hydrogen-bonding interactions contribute 26.1 and 18.4%, respectively, of the total Hirshfeld surface. The lattice energy of the molecule was calculated using density functional theory (B3LYP) methods with the 6-31G** basis set. The molecular electrostatic potential of the pyrazinamide molecule exhibits extended electronegative regions around O1, N1 and N2. The existence of a negative electrostatic potential (ESP) region just above the upper and lower surfaces of the pyrazine ring confirm the pi-electron cloud.

PROPENSITIES OF AMINO ACID RESIDUES IN PROTEINS FOR DIFFERENT REGIONS OF THE RAMACHANDRAN MAP

Using a dataset of 1164 crystal structures of largely non-homologous proteins defined at a resolution of 1.5 angstrom or better, we have investigated the (phi,psi) preferences of 20 residue types by considering the residues which occur in loops. Propensities of residue types to occur in the loops with (phi,psi) values in the aa region of the Ramachandran map has a poor correlation coefficient of 0.48 to the Chou-Fasman propensities of the residue types to occur in the a-helical segments. However the correlation coefficient between propensities of residues in loops to adopt beta conformations and those in beta-sheet is much higher (0.95). These observations suggest that a-helix formation is well influenced by the local amino acid sequence while intrinsic preference of residue types for beta-sheet plays a major role in the formation of beta-sheet. The main chain polar groups of residues in loops, that can affect the (phi,psi) values, can be involved in intra-molecular hydrogen bonding. Therefore we investigated further by considering subset of residues in loops with low (0 to 2) number of intra-molecular hydrogen bonds per residue involving main chain polar atoms. For this subset, the correlation coefficients between propensities for alpha-helix and alpha(R) region and between beta-sheet and beta-region are 0.26 and 0.64 respectively. This reiterates higher intrinsic tendency of beta-region favouring residues to adopt beta-sheet than alpha(R) region favouring residues to adopt alpha-helical structure.

INDUCTION OF FOLDED STRUCTURES IN DESIGNED PEPTIDES USING CONFORMATIONALLY CONSTRAINED RESIDUES

Folding into compact globular structures, with well-defined modules of secondary structure, appears to be a characteristic of long polypeptide chains, with a specific patterning of coded amino acid residues along the length of sequence. Cooperative hydrogen bond driven secondary structure formation and solvent forces, which contribute favorably to the entropy of folding, by promoting compaction of the polymeric chain, have long been discussed as major determinants of the folding process. First principles design approaches, which use non-coded amino acids, employ an alternative structure directing strategy, by using amino acid residues which exhibit a strong conformational bias for specific regions of the Ramachandran map. This overview of ongoing studies in the authors' laboratory, attempts to explore the use of conformationally restricted amino acid residues in the design of peptides with well-defined secondary structures. Short peptides composed of 20 genetically coded amino acids usually exist in solution as an ensemble of equilibrating conformations. Apolar peptide sequences, which are readily soluble in organic solvents like chloroform and methanol, facilitate formation of structures which are predominately driven by intramolecular hydrogen bond formation. The choice of sequences containing residues with a limited range of conformational choices strongly favors formation of local turn structures, stabilized by short range intramolecular hydrogen bonds. Two residue beta-turns can nucleate either helical or hairpin folding, depending on the precise conformation of the turn segment Restriction of the conformational space available to amino acid residues is easily achieved by introduction of an additional alkyl group at the C alpha carbon atom or by side chain backbone cyclization, as in proline. Studies of synthetic sequences incorporating two prototype residues alpha-aminoisobutyric acid (Aib) and D-proline (DPro) illustrate the utility of the strategy in construction of helices and hairpins. Extensions to the design of conformationally switchable sequences and structurally defined hybrid peptides containing backbone homologated residues are also surveyed.

Synthesis of 5,6-dihydro-4H-1,3-thiazine containing peptide mimics from N-(3-hydroxypropyl)thioamides and epimerization studies

The first synthesis of 1,3-thiazine fused peptide mimics is described from N-(3-hydroxypropyl)thioamides under MsCl/NEt3 conditions. The method is amenable to oligopeptidomimics with polar and apolar side chains. Substantial epimerization occurs at chiral C(2) exo methines in non-Pro fused mimics even under neutral conditions. H-1 NMR and crystal structure analyses indicate that the Thi analogues primarily associate with each other through intermolecular hydrogen bonds, involving the nitrogen of 1,3-thiazine and the N-H of the fused residue, which may be the basis for enamination-racemization process in these peptide mimics. (C) 2014 Elsevier Ltd. All rights reserved.

Component Selection in the Self-Assembly of Palladium(II) Nanocages and Cage-to-Cage Transformations

Dynamic supramolecular systems involving a tetratopic palladium(II) acceptor and three different pyridine-and imidazole-based donors have been used for self-selection by a synergistic effect of morphological information and coordination ability of ligands through specific coordination interactions. Three different cages were first synthesized by two-component self-assembly of individual donor and acceptor. When all four components were allowed to interact in a reaction mixture, only one out of three cages was isolated. The preferential binding affinity towards a particular partner was also established by transforming a non-preferred cage into a preferred cage by interaction with the appropriate ligand. Computational studies further supported the fact that coordination interaction of imidazole moiety to Pd-II is enthalpically more preferred compared to pyridine, which drives the selection process. Analysis of crystal packing of both complexes indicated the presence of strong hydrogen bonds between nitrate and water molecules and also H-bonded 3D networks of water. Both complexes exhibit promising proton conductivity (10(-5) to ca. 10(-3) Scm(-1)) at ambient temperature under a relative humidity of circa 98% with low activation energy.

Succinate esters: odd-even effects in melting points

Dialkyl succinates show a pattern of alternating behavior in their melting points, as the number of C atoms in the alkane side chain increases, unlike in the dialkyl oxalates Joseph et al. (2011). Acta Cryst. B67, 525-534]. Dialkyl succinates with odd numbers of C atoms in the alkyl side chain show higher melting points than the immediately adjacent analogues with even numbers. The crystal structures and their molecular packing have been analyzed for a series of dialkyl succinates with 1 - 4 C atoms in the alkyl side chain. The energy difference (Delta E) between the optimized and observed molecular conformations, density, Kitaigorodskii packing index (KPI) and C-H center dot center dot center dot O interactions are considered to rationalize this behavior. In contrast to the dialkyl oxalates where a larger number of moderately strong C-H center dot center dot center dot O interactions were characteristic of oxalates with elevated melting points, here the molecular packing and the density play a major role in raising the melting point. On moving from oxalate to succinate esters the introduction of the C2 spacer adds two activated H atoms to the asymmetric unit, resulting in the formation of stronger C-H center dot center dot center dot O hydrogen bonds in all succinates. As a result the crystallinity of long-chain alkyl substituted esters improves enormously in the presence of hydrogen bonds from activated donors.

Local Structural and Environmental Factors Define the Efficiency of an RNA Pseudoknot Involved in Programmed Ribosomal Frameshift Process

In programmed -1 ribosomal frameshift, an RNA pseudoknot stalls the ribosome at specific sequence and restarts translation in a new reading frame. A precise understanding of structural characteristics of these pseudoknots and their PRF inducing ability has not been clear to date. To investigate this phenomenon, we have studied various structural aspects of a -1 PRF inducing RNA pseudoknot from BWYV using extensive molecular dynamics simulations. A set of functional and poorly functional forms, for which previous mutational data were available, were chosen for analysis. These structures differ from each other by either single base substitutions or base-pair replacements from the native structure. We have rationalized how certain mutations in RNA pseudoknot affect its function; e.g., a specific base substitution in loop 2 stabilizes the junction geometry by forming multiple noncanonical hydrogen bonds, leading to a highly rigid structure that could effectively resist ribosome-induced unfolding, thereby increasing efficiency. While, a CG to AU pair substitution in stem 1 leads to loss of noncanonical hydrogen bonds between stems and loop, resulting in a less stable structure and reduced PRF inducing ability, inversion of a pair in stem 2 alters specific base-pair geometry that might be required in ribosomal recognition of nucleobase groups, negatively affecting pseudoknot functioning. These observations illustrate that the ability of an RNA pseudoknot to induce -1 PRF with an optimal rate depends on several independent factors that contribute to either the local conformational variability or geometry

Depalladation of Neutral Monoalkyne- and Dialkyne-Inserted Palladacycles and Alkyne Insertion/Depalladation Reactions of Cationic Palladacycles Derived from N,N `,N `'-Triarylguanidines as Facile Routes for Guanidine-Containing Heterocycles/Carbocycles: Synthetic, Structural, and Mechanistic Aspects

Depalladation of the monoalkyne-inserted cyclopalldated guanidines (kappa 2(C,N)Pd(2,6-Me2C5H3N)Br] (I and II) in PhCl under reflux conditions and that of the dialkyne-inserted cyclopalladated guanidine kappa(2)(C,N):eta(2)(C=C)PdBr] (III) in pyridine under reflux conditions afforded a guanidine-containing indole (1), imidaziondole (2), and benzazepine (3) in 80%, 67%, and 76%, yields, respectively. trans-L2PdBr2] species (L = 2,6-Me2C5H3N, C5H5N) were also isolated in the aforementioned reactions in 35%, 42%, and 40% yields. Further , the reaction of the cyclopalladated guanidine kappa(2)(C,N)Pd(mu-Br)](2) (IV) with AgBF4 in a CH2Cl2/MeCN mixture afforded the cationic pincer type cyclopalladated guanidine kappa(3)(C,N,O)Pd(MeCN)]BF4] (4) in 85% yield and this palladacycle upon crystallization in MeCN and the reaction of kappa(2)(C,N)Pd(mu-Br)](2) (V) with AgBf(4) in a CH2Cl2/MeCN mixture afforded the cationic palladacycles {kappa(2)(C,N)Pd(MeCN)(2)]BF4](5 and 6) in 89% and 91% yields, respectively. The separate reactions of 4 with 2 equiv of methyl phenylpropiolate (MPP) or diphenylacetylene (DPA) and the reaction of 5 with 2 equiv of MPP in PhCl at 110 degrees C afforded the guanidine-containing quinazolinium tetrafluoroborate 7 in 25-32% yields. The reaction of 6 with 2 equiv of DPA under otherwise identical conditions afforded the unsymmetrically substituted guanidinium tetrafluoroborate 8, containing a highly substituted naphthalene unit, in 82% yield. Compounds 1-8 were characterized by analytical and spectroscopic techniques, and all compounds except 4 were characterized by single-crystal X-ray diffraction. The Molecular structure of 2 and 3 are nove, as the framework in the former arises due to the formation of two C-N bonds upon depalladation while the butadienyl unit in the latter revealed cis,cis stereochemistry, a-feature unprecedented in alkyne insertion chemistry. Plausible pathways for the formation of heterocycles/carbocycles are proposed. the influence of substitutents on the aryl rings fo the cyclopalladated guanidine moiety and those on alkynes upon the nature of the products in addressed. Heterocycles 1 and 7 revealed the presence of two rotamers in about a 1.00:0.43 ratio in CDCl3 and in about a 1.00:0.14 ratio in CD3OD, respectively, as detected by H-1 NMR spectroscopy while in CD3CN and DMSO-d(6) (1) and CD3CN and CDCl3 (7), these heterocycles revealed the presence of a single rotamer. These spectral features are attributed to the restricted C-N single-bond rotation of the CN3 unit of the guanidine moiety, which possibly arises from steric constraint due to the formation of a N-H center dot center dot center dot Cl hydrogen bond with CDCl3 (1) and N-H center dot center dot center dot O and O-D center dot center dot center dot O hydrogen bonds with CD3OD (7).

Sequence and conformational preferences at termini of alpha-helices in membrane proteins: Role of the helix environment

-helices are amongst the most common secondary structural elements seen in membrane proteins and are packed in the form of helix bundles. These -helices encounter varying external environments (hydrophobic, hydrophilic) that may influence the sequence preferences at their N and C-termini. The role of the external environment in stabilization of the helix termini in membrane proteins is still unknown. Here we analyze -helices in a high-resolution dataset of integral -helical membrane proteins and establish that their sequence and conformational preferences differ from those in globular proteins. We specifically examine these preferences at the N and C-termini in helices initiating/terminating inside the membrane core as well as in linkers connecting these transmembrane helices. We find that the sequence preferences and structural motifs at capping (Ncap and Ccap) and near-helical (N' and C') positions are influenced by a combination of features including the membrane environment and the innate helix initiation and termination property of residues forming structural motifs. We also find that a large number of helix termini which do not form any particular capping motif are stabilized by formation of hydrogen bonds and hydrophobic interactions contributed from the neighboring helices in the membrane protein. We further validate the sequence preferences obtained from our analysis with data from an ultradeep sequencing study that identifies evolutionarily conserved amino acids in the rat neurotensin receptor. The results from our analysis provide insights for the secondary structure prediction, modeling and design of membrane proteins. Proteins 2014; 82:3420-3436. (c) 2014 Wiley Periodicals, Inc.

Charge density distribution and electrostatic interactions of ethionamide: an inhibitor of the enoyl acyl carrier protein reductase (inhA) enzyme of Mycobacterium tuberculosis

An experimental charge density analysis of an anti-TB drug ethionamide was carried out from high resolution X-ray diffraction at 100 K to understand its charge density distribution and electrostatic properties. The experimental results were validated from periodic theoretical charge density calculations performed using CRYSTAL09 at the B3LYP/6-31G** level of theory. The electron density rho(bcp)(r) and the Laplacian of electron density del(2)(rho bcp)(r) of the molecule calculated from both the methods display the charge density distribution of the ethionamide molecule in the crystal field. The electrostatic potential map shows a large electropositive region around the pyridine ring and a large electronegative region at the vicinity of the thiol atom. The calculated experimental dipole moment is 10.6D, which is higher than the value calculated from theory (8.2D). The topological properties of C-H center dot center dot center dot S, N-H center dot center dot center dot N and N-H center dot center dot center dot S hydrogen bonds were calculated, revealing their strength. The charge density analysis of the ethionamide molecule determined from both the experiment and theory gives the topological and electrostatic properties of the molecule, which allows to precisely understand the nature of intra and intermolecular interactions.

Directing peptide conformation with centrally positioned pre-organized dipeptide segments: studies of a 12-residue helix and beta-hairpin

Secondary structure formation in oligopeptides can be induced by short nucleating segments with a high propensity to form hydrogen bonded turn conformations. Type I/III turns facilitate helical folding while type II'/I' turns favour hairpin formation. This principle is experimentally verified by studies of two designed dodecapeptides, Boc-Val-Phe-Leu-Phe-Val-Aib-Aib-Val-Phe-Leu-Phe-Val-OMe 1 and Boc-Val-Phe-Leu-Phe-Val- (D) Pro- (L) Pro-Val-Phe-Leu-Phe-Val-OMe 2. The N- and C-terminal flanking pentapeptide sequences in both cases are identical. Peptide 1 adopts a largely alpha-helical conformation in crystals, with a small 3(10) helical segment at the N-terminus. The overall helical fold is maintained in methanol solution as evidenced by NMR studies. Peptide 2 adopts an antiparallel beta-hairpin conformation stabilized by 6 interstrand hydrogen bonds. Key nuclear Overhauser effects (NOEs) provide evidence for the antiparallel beta-hairpin structure. Aromatic proton chemical shifts provide a clear distinction between the conformation of peptides 1 (helical) and 2 (beta-hairpin). The proximity of facing aromatic residues positioned at non-hydrogen bonding positions in the hairpin results in extensively ring current shifted proton resonances in peptide 2.

Hierarchical sampling for metastable conformers determines biomolecular recognition: the case of malectin and diglucosylated N-glycan interactions

Structural information over the entire course of binding interactions based on the analyses of energy landscapes is described, which provides a framework to understand the events involved during biomolecular recognition. Conformational dynamics of malectin's exquisite selectivity for diglucosylated N-glycan (Dig-N-glycan), a highly flexible oligosaccharide comprising of numerous dihedral torsion angles, are described as an example. For this purpose, a novel approach based on hierarchical sampling for acquiring metastable molecular conformations constituting low-energy minima for understanding the structural features involved in a biologic recognition is proposed. For this purpose, four variants of principal component analysis were employed recursively in both Cartesian space and dihedral angles space that are characterized by free energy landscapes to select the most stable conformational substates. Subsequently, k-means clustering algorithm was implemented for geometric separation of the major native state to acquire a final ensemble of metastable conformers. A comparison of malectin complexes was then performed to characterize their conformational properties. Analyses of stereochemical metrics and other concerted binding events revealed surface complementarity, cooperative and bidentate hydrogen bonds, water-mediated hydrogen bonds, carbohydrate-aromatic interactions including CH-pi and stacking interactions involved in this recognition. Additionally, a striking structural transition from loop to beta-strands in malectin CRD upon specific binding to Dig-N-glycan is observed. The interplay of the above-mentioned binding events in malectin and Dig-N-glycan supports an extended conformational selection model as the underlying binding mechanism.