Derivation of a 3D pharmacophore model for the angiotensin-II site one receptor

Synthesis, binding studies and in vivo biological evaluation of novel non-peptide antihypertensive analogues

AT(1) antagonists (SARTANs) constitute the last generation of drugs for the treatment of hypertension, designed and synthesized to mimic the C-terminal segment of the vasoconstrictive hormone angiotensin II (AngII). They exert their action by blocking the binding of AngII on the AT(1) receptor. Up to date eight AT(1) antagonists have been approved for the regulation of high blood pressure. Although these molecules share common structural features and are designed to act under the same mechanism, they have differences in their pharmacological profiles and antihypertensive efficacy. Thus, there is still a need for novel analogues with better pharmacological and financial profiles. An example of a novel synthetic non peptide AT(1) antagonist which devoids the classical template of SARTANs is MM1. In vivo studies showed that MMK molecules, which fall in the same class of MM1, had a significant antihypertensive (40-80% compared to the drug losartan) activity. However, in vitro affinity studies showed that losartan has considerably higher affinity. The theoretical docking studies showed that MM1 acts on the same site of the receptor as losartan. They exert hydrophobic interactions with amino acid Val108 of the third helix of the AT(1) receptor and other hydrophobic amino acids in spatial vicinity. In addition, losartan favours multiple hydrogen bondings between its tetrazole group with Lys199. These additional interactions may in part explain its higher in vitro binding affinity.

Design and synthesis of new tetrazolyl- and carboxy-biphenylylmethyl-quinazolin-4-one derivatives as angiotensin II AT1 receptor antagonists

A series of novel quinazolin-4-ones was designed and their molecular modeling simulation fitting to a new HipHop 3D pharmacophore model using CATALYST was examined. Several compounds showed significant high simulation fit values. The designed compounds were synthesized and eight of them were biologically evaluated in vitro using an AT1 receptor binding assay, where compound XX competed weakly against radiolabeled Sar1Ile8-angiotensin II (Ang II) binding, compounds XIV and XXII showed moderate competition, and compound XXV showed almost equal ability to displace radiolabeled Sar1Ile8-Ang II binding to AT1 receptors as losartan. In vivo biological evaluation study of compounds XIV, XXII, and XXV on both normotensive and hypertensive rats revealed that compound XXV demonstrated higher hypotensive and antihypertensive activity than the reference compound losartan. To obtain a highly active compound from a candidate set of only eight tested compounds illustrates the power and utility of our pharmacophore model.

Losartan's molecular basis of interaction with membranes and AT1 receptor

Physicochemical methods were used to study the thermal and dynamic changes caused by losartan in the membrane bilayers. In addition, molecular modeling was implemented to explore its topography both in membranes and AT(1) receptor. Its incorporation resulted in the modification of thermal profile of dipalmitoyl phosphatidylcholine (DPPC) bilayers in a concentration dependent way up to 20mol% as it is depicted from the combination of differential scanning calorimetry (DSC) and MAS data. In particular, the presence of losartan caused lowering of the phase transition temperature and abolishment of the pretransition. T(1) experiments revealed the location of the drug into the membrane bilayers. The use of a combination of biophysical methods along with docking experiments brought out a possible two-step mechanism which involves incorporation of losartan at the interface of membrane bilayers and diffusion in the upper parts of AT(1) receptor helices IV-VII.

A 3D-QSAR of Angiotensin II (AT1) Receptor Antagonists Based on Receptor Surface Analysis

A hypothetical receptor surface model has been constructed for a set of 38 AT1 antagonists using activity data of each molecule as a weight in the building of the receptor surface. The best model was derived by optimizing various parameters such as atomic partial charges, surface fit, and the manner of representation of electrostatics on the surface. Descriptors such as van der Waals energy, electrostatic energy, and total nonbonded energy were used individually or in combination to derive a family of quantitative structure-activity relationship equations using G/PLS as the statistical method.

Therapeutic effects of angiotensin (AT1) receptor antagonists: potential contribution of mechanisms other than AT1 receptor blockade

Blockade of the renin-angiotensin system improves morbidity and mortality of patients with cardiovascular diseases, e.g. arterial hypertension, renal failure, following myocardial infarction and in congestive heart failure. The angiotensin II type 1 (AT(1)) receptor antagonists (angiotensin receptor blockers; ARBs), i.e. losartan, eprosartan, irbesartan and valsartan were developed by computer-based molecule design. Early observations already indicate that the ARBs elicit pleiotropic effects developing anti-aggregatory, anti-inflammatory and anti-mitogenic effects independent from their actions at the AT(1) receptor. Losartan metabolism indicates a number of known active intermediates and pointed to further interactions of these derivatives with other receptors and cellular signaling systems. Here we discuss a compilation of detailed pharmacokinetic and pharmacodynamic data of active metabolites of ARBs indicating their mode of action and suggest novel therapeutic implications. The clinical observations that ARBs elicit potencies in patients with cardiovascular diseases via the regulation of inflammatory, growth and homeostatic factors lead us to focus on specific, reactive metabolites, which hold potential for future indications and possible drug interactions in cardiovascular diseases.

A comparative molecular field analysis and molecular modelling studies on pyridylimidazole type of angiotensin II antagonists

A large number of compounds known as "AII (Angiotensin II) antagonists" have been developed for the treatment of various heart diseases such as hypertension, congestive heart failure, and chronic renal failure. Most of the currently known AII receptor antagonists share a similar chemical structure, consisting of nitrogen atoms, a lipophilic alkyl side chain and an acidic group. As a new series, we have designed and synthesized various pyridylimidazole derivatives. In this report we would like to discuss the structure-activity relationship of these series of compounds using the comparative molecular field analysis (CoMFA) methods. We could come up with a good CoMFA model (cross-validated and conventional r2 values equal to 0.702 and 0.991, respectively) and the validity of the model was confirmed by synthesizing and measuring their biological activity of additional 6 compounds suggested by the model. This result provides additional information on the structural requirement for structurally diverse group of AII receptor antagonists.

Prospective identification of biologically active structures by topomer shape similarity searching

The principle of bioisosterism-similarly shaped molecules are more likely to share biological properties than are other molecules-has long helped to guide drug discovery. An algorithmic implementation of this principle, based on shape comparisons of a single rule-generated "topomer" conformation per molecule, had been found to be the descriptor most consistently predictive of similar biological properties, in retrospective studies, and also to be well-suited for searching large (>10(12)) "virtual libraries" of potential reaction products. Therefore a prospective trial of this shape similarity searching method was carried out, with synthesis of 425 compounds and testing of them for inhibition of binding of angiotensin II (A-II). The 63 compounds that were identified by shape searching as most similar to any of four query structures included all of the seven compounds found to be highly active, with none of the other 362 structures being highly active (p < 0.001). Additional consistent relations (p < 0.05) were found, among all 425 compounds, between the degree of shape similarity to the nearest query structure and the frequency of various levels of observed activity. Known "SAR" (rules specifying structural features required for A-II antagonism) were also regenerated within the biological data for the 63 shape similar structures.

An effort to understand the molecular basis of hypertension through the study of conformational analysis of losartan and sarmesin using a combination of nuclear magnetic resonance spectroscopy and theoretical calculations

Losartan is the first recently approved drug against hypertension disease that competes with the biological action of angiotensin II (AII) at the AT1 receptor. Its design was based on the mimicry of the C-terminal segment of AII. Due to the biological significance of Losartan, its structure elucidation and conformational properties are reported as determined by NMR spectroscopy and computational analysis. In addition, molecular modeling of the peptide Sarmesin [Sar1Tyr(OMe)4AII], a competitive antagonist of AII, was also developed based on NMR and computational analysis data. Sarmesin's C-terminal was used as a template for superimposition with specific molecular features of interest in the structure of Losartan such as the conformation of biphenyltetrazole, the n-butyl chain, and the orientation of hydroxymethylimidazole relative to the biphenyl template. The major conclusions derived from this study are the following: (a) Sarmesin, like the AII superagonist [Sar1]AII, adopts a conformation which keeps in close proximity the key amino acids Sar1 (or Arg2)-Tyr(OMe)4-His6-Phe8. (b) Losartan favors a low-energy conformation in which imidazole and tetrazole rings are placed in the opposite site relative to the spacer phenyl ring plane; the hydroxymethyl group is placed away from the spacer phenyl ring, the alkyl chain is oriented above the spacer phenyl ring, and the two phenyl rings deviate approximately 60 degrees from being coplanar. (c) Overlay of the C-terminal region of Sarmesin with Losartan using equivalent groups revealed an excellent match. (d) Interestingly, the matching between enantiomeric structures of Losartan was not equivalent, proposing that the chirality of this molecule is significant in order to exert its biological activity. These findings open a new avenue for synthetic chemists to design and synthesize peptidomimetic drugs based on the C-terminal segment of the proposed model of Sarmesin. The new candidate drug molecules are not restricted to structurally resemble Losartan as the design is hitherto focused.

SQ: a program for rapidly producing pharmacophorically relevent molecular superpositions

A new method SQ has been developed to provide fast, automatic, and objective pairwise three-dimensional molecular alignments. SQ uses an atom-based clique-matching step followed by an alignment scoring function that has been parametrized to recognize pharmacologically relevant atomic properties. Molecular alignments from SQ are consistent with known drug-receptor interactions. We demonstrate this with six pairs of receptor-ligand complexes from the Brookhaven Protein Data Bank. The SQ-generated alignment of one isolated ligand onto another is shown to approximate the alignment of the ligands when the receptors are superimposed. SQ appears to be better than its predecessor SEAL (Kearsley and Smith, Tetrahedron Comput. Methodol. 1990, 3, 615-633) in this regard. SQ has been tailored so that, given one molecule as a probe, it can be used to routinely scan large chemical databases for which precomputed conformations have been stored. The SQ score, a measure of 3D similarity of each candidate molecule to the probe, can be used to rank compounds for the purposes of chemical screening. We demonstrate this with three probes (a thrombin inhibitor, an HIV protease inhibitor, and a model for angiotensin II). In each case SQ can preferentially select from the MDDR database other compounds with the same activity as the probe. We further show, using the angiotensin example, how SQ can identify topologically diverse compounds with the same activity.

Metabolites of the angiotensin II antagonist tasosartan: the importance of a second acidic group

Described in this paper is the synthesis and pharmacological activity of five metabolites of the angiotensin II antagonist tasosartan (1). Of particular interest is the effect of the additional acidic group of the enol metabolite (8) on activity. As suggested by the structural-activity relationship of other angiotensin II antagonist series, a second acidic group can improve receptor binding activity but decrease in vivo activity after oral dosing. The metabolic introduction of a second acidic group in tasosartan bypasses this problem and contributes to the excellent profile of the compound. A molecular modeling study provides a rationale for the role of the enol group of 8 in AT1 receptor binding.

Derivation of a three-dimensional pharmacophore model of substance P antagonists bound to the neurokinin-1 receptor

Constrained systematic search was used in an exhaustive conformational analysis of a structurally diverse set of substance P (SP) antagonists to identify a unique hypothesis for their bound conformation at the neurokinin-1 receptor. In this conformation, two aromatic groups essential for high affinity adopt a perpendicular or edge-on arrangement. This pharmacophore hypothesis for the receptor-bound conformation was used in a comparative molecular field analysis (CoMFA) of an expanded set of SP antagonists, and the predictive ability of the resulting three-dimensional quantitative structure-activity relationship (3D-QSAR) was evaluated against a test set of SP antagonists different from those in the training set. This CoMFA model based on the Constrained Search alignment yielded significant cross-validated, conventional, and predictive r2 values equal to 0.70, 0.93, and 0.82, respectively. For comparison, the SP antagonists were forced into an alternative poorer alignment in which the two aromatic rings were parallel and then subjected to a CoMFA analysis. Both the parallel and perpendicular arrangements of the aromatic rings are seen in X-ray structures of SP antagonists and have been proposed as candidates for the receptor-bound conformation. The parallel (or stacked) conformation yielded a poorer correlation with a cross-validated r2 = 0.57, a conventional r2 = 0.90, and a predictive r2 = 0.78. Our results indicate that although both alignments could generate a reasonable CoMFA correlation, the stacked conformation is unlikely to be the receptor-bound conformation, as the covalent structure of the antagonists precludes a common geometry in which the aromatic rings are stacked.

An automated method for predicting the positions of hydrogen-bonding atoms in binding sites

Hydrogen bonds are the most specific, and therefore predictable of the intermolecular interactions involved in ligand-protein binding. Given the structure of a molecule, it is possible to estimate the positions at which complementary hydrogen-bonding atoms could be found. Crystal-survey data are used in the design of a program, HBMAP, that generates a hydrogen-bond map for any given ligand, which contains all the feasible positions at which a complementary atom could be found. On superposition of ligands, the overlapping regions of their maps represent positions of receptor atoms to which each molecule can bind. The certainty of these positions is increased by the incorporation of a larger number and diversity of molecules. In this work, superposition is achieved using the program HBMATCH, which uses simulated annealing to generate the correspondence between points from the hydrogen-bonding maps of the two molecules. Equivalent matches are distinguished on the basis of their steric similarity. The strategy is tested on a number of ligands for which ligand-protein complexes have been solved crystallographically, which allows validation of the techniques. The receptor atom positions of thermolysin are successfully predicted when the correct superposition is obtained.

A 3D QSAR CoMFA study of non-peptide angiotensin II receptor antagonists

A series of non-peptide angiotensin II receptor antagonists was investigated with the aim of developing a 3D QSAR model using comparative molecular field analysis descriptors and approaches. The main goals of the study were dictated by an interest in methodologies and an understanding of the binding requirements to the AT1 receptor. Consistency with the previously derived activity models was always checked to contemporarily test the validity of the various hypotheses. The specific conformations chosen for the study, the procedures invoked to superimpose all structures, the conditions employed to generate steric and electrostatic field values and the various PCA/PLS runs are discussed in detail. The effect of experimental design techniques to select objects (molecules) and variables (descriptors) with respect to the predictive power of the QSAR models derived was especially analysed.

Proposals for the angiotensin II receptor-bound conformation by comparative computer modeling of AII and cyclic analogs

A conformational search using high-temperature molecular dynamics on angiotensin II(AII) and on two cyclic S-S bridged analogs, namely [Hcy3,5]AII and [Cys3,5]AII, in conjunction with a cluster analysis based on the similarities of the three-dimensional patterns of the binding and activation elements, had led to putative AII receptor-bound conformations. These conformations are characterized by a compact folded shape of the peptide backbone, and by particular relative positions of the four pharmacophore groups, namely the aromatic moieties of the Tyr4, His6 and Phe8 residues, and the C-terminal carboxyl group. This compact folded shape, arising from attractive electrostatic interactions between the desolvated N- and C-terminal groups, is similar to the crystallographically determined conformation of AII bound to the antibody Fab receptor.

Tetrazole and carboxylate groups of angiotensin receptor antagonists bind to the same subsite by different mechanisms

To identify specific interactions between either the tetrazole or carboxylate pharmacophores of non-peptide antagonists and the rat AT1 receptor, 6 basic residues were examined by site-directed mutagenesis. Three of the mutants (H183Q, H256Q, and H272Q) appeared to be like wild type. Lys102 and Arg167 mutants displayed reduced binding of the non-peptide antagonist losartan. Examination of their properties employing group-specific angiotensin II analogues indicated that their effects on binding were indirect. Interestingly, the affinity of losartan was not altered by a K199Q mutation, but the same mutation reduced the affinity of angiotensin II, the antagonist [Sar1,Ile8]angiotensin II, and several carboxylate analogues of losartan. An Ala199 substitution reduced the affinity of peptide analogues to a larger extent as compared to the affinity of losartan. Thus, the crucial acidic pharmacophores of angiotensin and losartan appear to occupy the same space within the receptor pocket, but the protonated amino group of Lys199 is not essential for binding the tetrazole anion. The binding of the tetrazole moiety with the AT1 receptor involves multiple contacts with residues such as Lys199 and His256 that constitute the same subsite of the ligand binding pocket. However, this interaction does not involve a conventional salt bridge, but rather an unusual lysine-aromatic interaction.

Structural model of antagonist and agonist binding to the angiotensin II, AT1 subtype, G protein coupled receptor

BACKGROUND

The family of G protein coupled receptors is the largest and perhaps most functionally diverse class of cell-surface receptors. Due to the difficulty of obtaining structural data on membrane proteins there is little information on which to base an understanding of ligand structure-activity relationships, the effects of receptor mutations and the mechanism(s) of signal transduction in this family. We therefore set out to develop a structural model for one such receptor, the human angiotensin II receptor.

RESULTS

An alignment between the human angiotensin II (type 1; hAT1), human beta 2 adrenergic, human neurokinin-1, and human bradykinin receptors, all of which are G protein coupled receptors, was used to generate a three-dimensional model of the hAT1 receptor based on bacteriorhodopsin. We observed a region within the model that was congruent with the biogenic amine binding site of beta 2, and were thus able to dock a model of the hAT1 antagonist L-158,282 (MK-996) into the transmembrane region of the receptor model. The antagonist was oriented within the helical domain by recognising that the essential acid functionality of this antagonist interacts with Lys199. The structural model is consistent with much of the information on structure-activity relationships for both non-peptide and peptide ligands.

CONCLUSIONS

Our model provides an explanation for the conversion of the antagonist L-158,282 (MK-996) to an agonist by the addition of an isobutyl group. It also suggests a model for domain motion during signal transduction. The approach of independently deriving three-dimensional receptor models and pharmacophore models of the ligands, then combining them, is a powerful technique which helps validate both models.