Design, Synthesis and Biological Activities of (Thio)Urea Benzothiazole Derivatives.

(Thio)ureas ((T)Us) and benzothiazoles (BTs) each have demonstrated to have a great variety of biological activities. When these groups come together, the 2-(thio)ureabenzothizoles [(T)UBTs] are formed, improving the physicochemical as well as the biological properties, making these compounds very interesting in medicinal chemistry. Frentizole, bentaluron and methabenzthiazuron are examples of UBTs used for treatment of rheumatoid arthritis and as wood preservatives and herbicides in winter corn crops, respectively. With this antecedent, we recently reported a bibliographic review about the synthesis of this class of compounds, from the reaction of substituted 2-aminobenzothiazoles (ABTs) with iso(thio)cyanates, (thio)phosgenes, (thio)carbamoyl chlorides, 1,1'-(thio)carbonyldiimidazoles, and carbon disulfide. Herein, we prepared a bibliographic review about those features of design, chemical synthesis, and biological activities relating to (T)UBTs as potential therapeutic agents. This review is about synthetic methodologies generated from 1968 to the present day, highlighting the focus to transform (T)UBTs to compounds containing a range substituents, as illustrated with 37 schemes and 11 figures and concluded with 148 references. In this topic, the scientists dedicated to medicinal chemistry and pharmaceutical industry will find useful information for the design and synthesis of this interesting group of compounds with the aim of repurposing these compounds.

Synthetic Procedures to Access 2-Guanidinobenzazoles of Biological Interest.

Benzazoles (Bz) and derivatives are interesting molecules in medicinal chemistry. Several of these compounds display diverse biological activities; some are still used in clinical applications. In this way, synthetic chemists are interested in developing new procedures to access compounds with the guanidine moiety as 2-aminobenzimidazole (2ABI), Astemizole (antihistaminic), Albendazole (anthelmintic) and Carbendazim (fungicide). The guanidine group, considered a super base bonded to a benzoxazole ring, results in the 2-guanidinobenzazoles (2GBZs), which could modify the biological activity of these heterocycles. On these bases, we prepared this review article, which covers chemical aspects of 2-guanidinobenzoazoles as potential therapeutic agents and summarizes the current knowledge on the mechanism of pharmacological activities such as cytotoxic, inhibition of cell proliferation via angiogenesis and apoptosis. Specifically, it highlights the most recent results of synthetic approaches to 2GBZs with variety of modifications and functionalization with aromatic, carbohydrate, and amino-acid moieties as illustrated on 28 schemes and is concluded with 141 references. Additionally, the format of this interesting review is exclusively designed on specifically classified category of chemical reactions with primary precursors such as o-substituted anilines and 2-aminobenzazoles (2ABZs). This will constitute the important goals and novelty of this paper to facilitate synthetic chemists in the investigation about development of new pharmacophores.

N-substitution Reactions of 2-Aminobenzimidazoles to Access Pharmacophores.

Benzimidazole (BI) and its derivatives are interesting molecules in medicinal chemistry because several of these compounds have a diversity of biological activities and some of them are even used in clinical applications. In view of the importance of these compounds, synthetic chemists are still interested in finding new procedures for the synthesis of these classes of compounds. Astemizole (antihistaminic), Omeprazole (antiulcerative), and Rabendazole (fungicide) are important examples of compounds used in medicinal chemistry containing BI nuclei. It is interesting to observe that several of these compounds contain 2-aminobenzimidazole (2ABI) as the base nucleus. The structures of 2ABI derivatives are interesting because they have a planar delocalized structure with a cyclic guanidine group, which have three nitrogen atoms with free lone pairs and labile hydrogen atoms. The 10-π electron system of the aromatic BI ring conjugated with the nitrogen lone pair of the hexocyclic amino group, making these heterocycles to have an amphoteric character. Synthetic chemists have used 2ABI as a building block to produce BI derivatives as medicinally important molecules. In view of the importance of the BIs, and because no review was found in the literature about this topic, we reviewed and summarized the procedures related to the recent methodologies used in the N-substitution reactions of 2ABIs by using aliphatic and aromatic halogenides, dihalogenides, acid chlorides, alkylsulfonic chlorides, carboxylic acids, esters, ethyl chloroformates, anhydrides, SMe-isothioureas, alcohols, alkyl cyanates, thiocyanates, carbon disulfide and aldehydes or ketones to form Schiff bases. The use of diazotized 2ABI as intermediate to obtain 2-diazoBIs was included to produce Nsubstituted 2ABIs of pharmacological interest. Some commentaries about their biological activity were included.

Synthesis and Biological Importance of 2-(thio)ureabenzothiazoles.

The (thio)urea and benzothiazole (BT) derivatives have been shown to have a broad spectrum of biological activities. These groups, when bonded, result in the 2-(thio)ureabenzothizoles (TBT and UBT), which could favor the physicochemical and biological properties. UBTs and TBTs are compounds of great importance in medicinal chemistry. For instance, Frentizole is a UBT derivative used for the treatment of rheumatoid arthritis and systemic lupus erythematosus. The UBTs Bentaluron and Bethabenthiazuron are commercial fungicides used as wood preservatives and herbicides in winter corn crops. On these bases, we prepared this bibliography review, which covers chemical aspects of UBTs and TBTs as potential therapeutic agents as well as their studies on the mechanisms of a variety of pharmacological activities. This work covers synthetic methodologies from 1935 to nowadays, highlighting the most recent approaches to afford UBTs and TBTs with a variety of substituents as illustrated in 42 schemes and 13 figures and concluded with 187 references. In addition, this interesting review is designed on chemical reactions of 2-aminobenzothiazoles (2ABTs) with (thio)phosgenes, iso(thio)cyanates, 1,1'-(thio)carbonyldiimidazoles [(T)CDI]s, (thio)carbamoyl chlorides, and carbon disulfide. This topic will provide information of utility for medicinal chemists dedicated to the design and synthesis of this class of compounds to be tested with respect to their biological activities and be proposed as new pharmacophores.

Current Tools and Methods in Molecular Dynamics (MD) Simulations for Drug Design.

Molecular Dynamics (MD) simulations is a computational method that employs Newton's laws to evaluate the motions of water, ions, small molecules, and macromolecules or more complex systems, for example, whole viruses, to reproduce the behavior of the biological environment, including water molecules and lipid membranes. Specifically, structural motions, such as those that are dependent of the temperature and solute/ solvent are very important to study the recognition pattern of ligandprotein or protein-protein complexes, in that sense, MD simulations are very useful because these motions can be modeled using this methodology. Furthermore, MD simulations for drug design provide insights into the structural cavities required to design novel structures with higher affinity to the target. Also, the employment of MD simulations to drug design can help to refine the three-dimensional (3D) structure of targets in order to obtain a better sampling of the binding poses and more reliable affinity values with better structural advantages, because they incorporate some biological conditions that include structural motions compared to traditional docking procedures. This work analyzes the concepts and applicability of MD simulations for drug design because molecular structural motions are considered, and these help to identify hot spots, decipher structural details in the reported protein sites, as well as to eliminate sites that could be structural artifacts which could be originated from the structural characterization conditions from MD. Moreover, better free energy values for protein ligand recognition can also be obtained, and these can be validated under experimental procedures due to the robustness of the MD simulation methods.

Structural and energetic analysis to provide insight residues of CYP2C9, 2C11 and 2E1 involved in valproic acid dehydrogenation selectivity.

Docking and molecular dynamics (MD) simulation have been two computational techniques used to gain insight about the substrate orientation within protein active sites, allowing to identify potential residues involved in the binding and catalytic mechanisms. In this study, both methods were combined to predict the regioselectivity in the binding mode of valproic acid (VPA) on three cytochrome P-450 (CYP) isoforms CYP2C9, CYP2C11, and CYP2E1, which are involved in the biotransformation of VPA yielding reactive hepatotoxic intermediate 2-n-propyl-4-pentenoic acid (4nVPA). There are experimental data about hydrogen atom abstraction of the C4-position of VPA to yield 4nVPA, however, there are not structural evidence about the binding mode of VPA and 4nVPA on CYPs. Therefore, the complexes between these CYP isoforms and VPA or 4nVPA were studied to explore their differences in binding and energetic stabilization. Docking results showed that VPA and 4nVPA are coupled into CYPs binding site in a similar conformation, but it does not explain the VPA hydrogen atom abstraction. On the other hand, MD simulations showed a set of energetic states that reorient VPA at the first ns, then making it susceptible to a dehydrogenation reaction. For 4nVPA, multiple binding modes were observed in which the different states could favor either undergo other reaction mechanism or ligand expulsion from the binding site. Otherwise, the energetic and entropic contribution point out a similar behavior for the three CYP complexes, showing as expected a more energetically favorable binding free energy for the complexes between CYPs and VPA than with 4nVPA.

Molecular modeling used to evaluate CYP2C9-dependent metabolism: homology modeling, molecular dynamics and docking simulations.

Cytochrome P450 (CYP) 2C9 is the principal isoform of the CYP2C subfamily in the human liver and is involved in the oxidation of several endogenous and xenobiotic compounds, including many therapeutic drugs. The metabolism of drugs by CYP2C9 can yield either safe or toxic products, which may be related to the recognition and binding modes of the substrates to this isoform. These interactions can be studied using in silico methods such as quantum chemistry, molecular dynamics and docking simulations, which can also be useful for predicting the structure of metabolites. In these types of studies, the ligand and the protein must be tridimensional models; thus, the protein can be built by homology modeling or retrieved from the Protein Data Bank. Therefore, the current review emphasizes the importance of using in silico methods to predict the metabolism of CYP2C9 because these computational tools have allowed the description of the principal characteristics of the active site of this isoform at the molecular level and the chemical properties of its ligands.

Homology modeling and molecular dynamics of CYP1A1 and CYP2B1 to explore the metabolism of aryl derivatives by docking and experimental assays.

Since many drugs are metabolized by cytochrome P450 (CYP450), biotransformation studies using these enzymes are valuable in drug development. In this work, the biotransformation by CYP1A1 and CYP2B1 of two acetylcholinesterase (AChE) inhibitors, 4-(4'-hydroxy-phenylamino)-4-oxo propanoic acid (A) and 1H-pyrrolidine-1-(4'-hydroxy-phenyl)-2,5-dione (B), was investigated through docking and molecular dynamics (MD) simulations and by experimental methods using rat liver microsomes pretreated with ß-naphthoflavone and phenobarbital (CYP1A1 and CYP2B1 inducers, respectively). The target proteins were initially built by homology modeling, and the resulting three-dimensional structures were refined by MD to obtain fifteen snapshots of each P450 isoform. These snapshots were used to dock compounds A and B as well as the reference compound acetaminophen (APAP). We confirmed that APAP produces a toxic intermediate (N-acetyl-p-benzoquinone imine) upon interaction of its amide group with the heme iron of CYP1A1. However, neither A nor B presented this kind of interaction within any snapshot with CYP1A1. On the other hand, when APAP, A and B were docked on CYP2B1, their hydroxyl group was located near the heme iron on the snapshot at 3.5 ns. Furthermore, B maintained the same position on all snapshots of this isoform. Therefore, theoretical results suggests that A and B do not generate toxic metabolites. These data were supported by HPLC analysis showing only one metabolite from A and B, which was identified by GC-MS as the hydroxylated product. Altogether, our results suggest that neither test compound is toxic.

Comparing the electronic properties and docking calculations of heme derivatives on CYP2B4.

Cytochrome P-450 is a group of enzymes involved in the biotransformation of many substances, including drugs. These enzymes possess a heme group (1) that when it is properly modified induces several important physicochemical changes that affect their enzymatic activity. In this work, the five structurally modified heme derivatives 2-6 and the native heme 1 were docked on CYP2B4, (an isoform of P450), in order to determine whether such modifications alter their binding form and binding affinity for CYP2B4 apoprotein. In addition, docking calculations were used to evaluate the affinity of CYP2B4 apoprotein-heme complexes for aniline (A) and N-methyl-aniline (NMA). Results showing the CYP2B4 heme 4- and heme 6-apoprotein complexes to be most energetically stable indicate that either hindrance effects or electronic properties are the most important factors with respect to the binding of heme derivatives at the heme-binding site. Furthermore, although all heme-apoprotein complexes demonstrated high affinity for both A and NMA, the CYP2B4 apoprotein-5 complex had higher affinity for A, and the heme 6 complex had higher affinity for NMA. Finally, surface electronic properties (SEP) were calculated in order to explain why certain arginine residues of CYP2B4 apoprotein interact with polarizable functionalities, such as ester groups or sp ( 2 ) carbons, present in some heme derivates. The main physicochemical parameter involved in the recognition process of the heme derivatives, the CYP2B4 apoprotein and A or NMA, are reported.