Homology modeling and molecular dynamics simulations of the glycine receptor ligand binding domain

We present a homology based model of the ligand binding domain (LBD) of the homopentameric alpha1 glycine receptor (GlyR). The model is based on multiple sequence alignment with other members of the nicotinicoid ligand gated ion channel superfamily and two homologous acetylcholine binding proteins (AChBP) from the freshwater (Lymnaea stagnalis) and saltwater (Aplysia californica) snails with known high resolution structure. Using two template proteins with known structure to model three dimensional structure of a target protein is especially advantageous for sequences with low homology as in the case presented in this paper. The final model was cross-validated by critical evaluation of experimental and published mutagenesis, functional and other biochemical studies. In addition, a complex structure with strychnine antagonist in the putative binding site is proposed based on docking simulation using Autodock program. Molecular dynamics (MD) simulations with simulated annealing protocol are reported on the proposed LBD of GlyR, which is stable in 5 ns simulation in water, as well as for a deformed LBD structure modeled on the corresponding domain determined in low-resolution cryomicroscopy structure of the alpha subunit of the full-length acetylcholine receptor (AChR). Our simulations demonstrate that the beta-sandwich central core of the protein monomer is fairly rigid in the simulations and resistant to deformations in water.     

Molecular dynamics simulations of ethanol binding to the transmembrane domain of the glycine receptor: implications for the channel potentiation mechanism

The glycine receptor (GlyR) is potentiated by ethanol and other anesthetics. The potentiation mechanism at the molecular level is unknown and remains elusive, but mutagenic studies have shown that ethanol and other volatile anesthetics bind to a pocket between TM1, TM2, and TM3. The present study extends previous studies (Cheng et al., Proteins 2007;68:581-593) wherein we conducted homology modeling and molecular dynamics (MD) simulations to construct models of the homopentameric alpha1 subunits of the GlyR transmembrane domain in open and closed states. To understand the potentiation of GlyR by ethanol we compare the binding of ethanol molecules to the channel in these different states. We observe that ethanol stably resides inside solvent-accessible cavities found in the open state of GlyR that are formed by I229 (of TM1) in one subunit and S267 and A288 (of TM2 and TM3, respectively) in the adjacent subunit. The volume of these putative binding pockets is state-dependent. Selective binding to the open states of receptors has been proposed to explain the potentiating actions of this class of anesthetics. In accordance with this model, our MD simulations suggest that the potentiation of ethanol on GlyR may be effected through preferential binding of ethanol molecules to an inter-subunit binding pocket in the open state.     

Homology modeling and molecular dynamics study on Schwanniomyces occidentalis alpha-amylase

With consumers growing increasingly aware of environmental issues, industries find enzymes as a reasonable alternative over physical conditions and chemical catalysts. Amylases are important hydrolase enzymes, which have been widely used in variety of industrial process such as pharmaceutical, food, and fermentation industries. Among amylases α-Amylase is in maximum demand due to its wide range of applications. The homology modeling study on Schwanniomyces occidentalis amylase (AMY1, UniProt identifier number: P19269) was performed by Modeller using Aspergillus oryzae (6TAA) as the template. The resulting structure was analyzed for validity and subjected to 14 ns of molecular dynamics (MD) simulation trough GROMACS. The validity of obtained model may represent that utilized OPLS force field is suitable for calcium-containing enzymes. DSSP secondary structure and contact map analysis represent the conservation of domain A TIM barrel feature together with calcium ion coordination sphere. Investigating the covariance matrix followed by principle component analyses for the first five eigenvectors of both trajectories indicate a little more flexibility for AMY1 structure. The electrostatic calculation for the final structures shows similar isoelectric point and superimposed buffering zone in the 5–8 pH range.     

Homology modeling and molecular dynamics simulations of lymphotactin

We have modeled the structure of human lymphotactin (hLpnt), by homology modeling and molecular dynamics simulations. This chemokine is unique in having a single disulfide bond and a long C-terminal tail. Because other structural classes of chemokines have two pairs of Cys residues, compared to one in Lpnt, and because it has been shown that both disulfide bonds are required for stability and function, the question arises how the Lpnt maintains its structural integrity. The initial structure of hLpnt was constructed by homology modeling. The first 63 residues in the monomer of hLpnt were modeled using the structure of the human CC chemokine, RANTES, whose sequence appeared most similar. The structure of the long C-terminal tail, missing in RANTES, was taken from the human muscle fatty-acid binding protein. In a Protein Data Bank search, this protein was found to contain a sequence that was most homologous to the long tail. Consequently, the modeled hLpnt C-terminal tail consisted of both alpha-helical and beta-motifs. The complete model of the hLpnt monomer consisted of two alpha-helices located above the five-stranded beta-sheet. Molecular dynamics simulations of the solvated initial model have indicated that the stability of the predicted fold is related to the geometry of Pro78. The five-stranded beta-sheet appeared to be preserved only when Pro78 was modeled in the cis conformation. Simulations were also performed both for the C-terminal truncated forms of the hLpnt that contained one or two (CC chemokine-like) disulfide bonds, and for the chicken Lpnt (cLpnt). Our MD simulations indicated that the turn region (T30-G34) in hLpnt is important for the interactions with the receptor, and that the long C-terminal region stabilizes both the turn (T30-G34) and the five-stranded beta-sheet. The major conclusion from our theoretical studies is that the lack of one disulfide bond and the extension of the C-terminus in hLptn are mutually complementary. It is very likely that removal of two Cys residues sufficiently destabilizes the structure of a chemokine molecule, particularly the core beta-sheet, to abolish its biological function. However, this situation is rectified by the long C-terminal segment. The role of this long region is most likely to stabilize the first beta-turn region and alpha-helix H1, explaining how this chemokine can function with a single disulfide bond.     

Conformational states of the glucocorticoid receptor DNA-binding domain from molecular dynamics simulations

Molecular dynamics simulations (MD) have been performed on variant crystal and NMR-derived structures of the glucocorticoid receptor DNA-binding domain (GR DBD). A loop region five residues long, the so-called D-box, exhibits significant flexibility, and transient perturbations of the tetrahedral geometry of two structurally important Cys4 zinc finger are seen, coupled to conformational changes in the D-box. In some cases, one of the Cys ligands to zinc exchanges with water, although no global distortion of the protein structure is observed. Thus, from MD simulation, dynamics of the D-box could partly be explained by solvent effects in conjunction with structural reformation of the zinc finger.     

Cross-linking of sites involved with alcohol action between transmembrane segments 1 and 3 of the glycine receptor following activation

The glycine receptor is a member of the Cys-loop, ligand-gated ion channel family and is responsible for inhibition in the CNS. We examined the orientation of amino acids I229 in transmembrane 1 (TM1) and A288 in TM3, which are both critical for alcohol and volatile anesthetic action. We mutated these two amino acids to cysteines either singly or in double mutants and expressed the receptors in Xenopus laevis oocytes. We tested whether disulfide bonds could form between A288C in TM3 paired with M227C, Y228C, I229C, or S231C in TM1. Application of cross-linking (mercuric chloride) or oxidizing (iodine) agents had no significant effect on the glycine response of wild-type receptors or the single mutants. In contrast, the glycine response of the I229C/A288C double mutant was diminished after application of either mercuric chloride or iodine only in the presence of glycine, indicating that channel gating causes I229C and A288C to fluctuate to be within 6 Å apart and form a disulfide bond. Molecular modeling was used to thread the glycine receptor sequence onto a nicotinic acetylcholine receptor template, further demonstrating that I229 and A288 are near-neighbors that can cross-link and providing evidence that these residues contribute to a single binding cavity.     

Homology modeling of the Abl-SH3 domain

A tertiary structure model of the Abl-SH3 domain is predicted by using homology modeling techniques coupled to molecular dynamics simulations. Two template proteins were used, Fyn-SH3 and Spc-SH3. The refined model was extensively checked for errors using criteria based on stereochemistry, packing, solvation free-energy, accessible surface areas, and contact analyses. The different checking methods do not totally agree, as each one evaluates a different characteristic of protein structures. Several zones of the protein are more susceptible to incorporating errors. These include residues 13, 15, 35, 39, 45, 46, 50, and 60. An interesting finding is that the measurement of the Cα chirality correlated well with the rest of the criteria, suggesting that this parameter might be a good indicator of correct local conformation. Deviations of more than 4 degrees may be indicative of poor local structure. © 1994 Wiley-Liss, Inc.     

Homology modeling and molecular dynamics simulation of odonthubuthus doriae (Od1) scorpion toxin in comparison to the BmK M1

All of the α-subgroups share similarity in their sequence and structure but different in the toxicity to various voltage-gated sodium channels (VGSCs). We modeled the first 3D structural model of the Od1 based on BmK M1 using homology modeling. The reliability of model for more investigation and compare to BmK M1 has been examined and confirmed. Then the model structure is further refined by energy minimization and molecular dynamics methods. The purpose of this modeling and simulation is comparison toxicity of two mentioned toxins by investigation structural feature of functional regions including core domain, 5-turn and C-terminal which make NC domain. In the one hand, it is intriguing that Od1 in comparison to BmK M1 shows same solvent accessible surface area (SASA) in 5-turn region but a little more exposed and feasibility (more SASA) in C-terminal region and key functional residues of C-terminal such as positive residues Arg58, lys62 and Arg (His)64. These data suggested that Od1 has similarity with BmK M1 but has more toxicity to sodium channel. In the other hand 5-turn proximity of C-terminal to 5-turn in BmK M1with cis peptide bond is less than Od1 without cis peptide bond which is a confirmation with experimental data about BmK M1.A better understanding of the 3-D structure of Od1and comparison to BmK M1 will be helpful for more investigation of functional characters action of natural toxins with a specialized role for VGSCs.     

Homology modeling and molecular dynamics simulations of the active state of the nociceptin receptor reveal new insights into agonist binding and activation

The opioid receptor-like receptor, also known as the nociceptin receptor (NOP), is a class A G protein-coupled receptor (GPCR) in the opioid receptor family. Although NOP shares a significant homology with the other opioid receptors, it does not bind known opioid ligands and has been shown to have a distinct mechanism of activation compared to the closely related opioid receptors mu, delta, and kappa. Previously reported homology models of the NOP receptor, based on the inactive-state GPCR crystal structures, give limited information on the activation and selectivity features of this fourth member of the opioid receptor family. We report here the first active-state homology model of the NOP receptor based on the opsin GPCR crystal structure. An inactive-state homology model of NOP was also built using a multiple template approach. Molecular dynamics simulation of the active-state NOP model and comparison to the inactive-state model suggest that NOP activation involves movements of transmembrane (TM)3 and TM6 and several activation microswitches, consistent with GPCR activation. Docking of the selective nonpeptidic NOP agonist ligand Ro 64-6198 into the active-state model reveals active-site residues in NOP that play a role in the high selectivity of this ligand for NOP over the other opioid receptors. Docking the shortest active fragment of endogenous agonist nociceptin/orphaninFQ (residues 1-13) shows that the NOP extracellular loop 2 (EL2) loop interacts with the positively charged residues (8-13) of N/OFQ. Both agonists show extensive polar interactions with residues at the extracellular end of the TM domain and EL2 loop, suggesting agonist-induced reorganization of polar networks, during receptor activation.     

Pharmacophore-directed Homology Modeling and Molecular Dynamics Simulation of G Protein-coupled Receptor： Study of Possible Binding Modes of 5-HT2c Receptor Agonists

A new pharmacophore-based modeling procedure, including homology modeling, pharmacophore study, flexible molecular docking, and long-time molecular dynamics (MD) simulations, was employed to construct the structure of the human 5-HT_(2C) receptor and determine the characteristics of binding modes of 5-HT_(2C) receptor agonists. An agonist-receptor complex has been constructed based on homology modeling and a pharmacophore hypothesis model based on some high active compounds. Then MD simulations of the ligand-receptor complex in an explicit membrane environment were carried out. The conformation of the 5- HT_(2C) receptor during MD simulation was explored, and the stable binding modes of the studied agonist were determined. Flexible molecular docking of several structurally diverse agonists of the human 5-HT_(2C) receptor was carried out, and the general binding modes of these agonists were investigated. According to the models presented in this work and the results of Flexi-Dock, the involvement of the amino acid residues Asp134, Ser138, Ash210, Asn331, Tyr358, Ile131, Ser132, Val135, Thr139, Ile189, Val202, Val208, Leu209, Phe214, Val215, Gly218, Ser219, Phe223, Trp324, Phe327, and Phe328 in agonist recognition was studied. The obtained binding modes of the human 5-HT_(2C) receptor agonists have good agreement with the site-directed mutagenesis data and other studies.     

Structure prediction and molecular dynamics simulations of a G-protein coupled receptor: human CCR2 receptor

CC chemokine receptor type-2 (CCR2) is a member of G-protein coupled receptors superfamily, expressed on the cell surface of monocytes and macrophages. It binds to the monocyte chemoattractant protein-1, a CC chemokine, produced at the sites of inflammation and infection. A homology model of human CCR2 receptor based on the recently available C-X-C chemokine recepor-4 crystal structure has been reported. Ligand information was used as an essential element in the homology modeling process. Six known CCR2 antagonists were docked into the model using simple and induced fit docking procedure. Docked complexes were then subjected to visual inspection to check their suitability to explain the experimental data obtained from site directed mutagenesis and structure-activity relationship studies. The homology model was refined, validated, and assessed for its performance in docking-based virtual screening on a set of CCR2 antagonists and decoys. The docked complexes of CCR2 with the known antagonists, TAK779, a dual CCR2/CCR5 antagonist, and Teijin-comp1, a CCR2 specific antagonist were subjected to molecular dynamics (MD) simulations, which further validated the binding modes of these antagonists. B-factor analysis of 20?ns MD simulations demonstrated that Cys190 is helpful in providing structural rigidity to the extracellular loop (EL2). Residues important for CCR2 antagonism were recognized using free energy decomposition studies. The acidic residue Glu291 from TM7, a conserved residue in chemokine receptors, is favorable for the binding of Teijin-comp1 with CCR2 by ΔG of ?11.4?kcal/mol. Its contribution arises more from the side chains than the backbone atoms. In addition, Tyr193 from EL2 contributes ?0.9?kcal/mol towards the binding of the CCR2 specific antagonist with the receptor. Here, the homology modeling and subsequent molecular modeling studies proved successful in probing the structure of human CCR2 chemokine receptor for the structure-based virtual screening and predicting the binding modes of CCR2 antagonists.     

Molecular dynamics simulations on SDF-1alpha: binding with CXCR4 receptor

Insights into the interacting mode of CXCR4 with SDF-1alpha are crucial in understanding the structural and functional characteristics of CXCR4 receptor. In this paper a computational pipeline, integrating protein structure prediction, molecular dynamics simulations, automated molecular docking, and Brownian dynamics simulations were employed to investigate the dynamic and energetic aspects of CXCR4 associating with SDF-1alpha. The entire simulation revealed the surface distribution feature of electrostatic potentials and conformational "open-close" process of the receptor. The possible binding conformation of CXCR4 was identified, and the CXCR4-SDF-1alpha binding complex was generated. Arg188-Glu277 salt bridge plays an important role for both the extracellular domain conformational change and SDF-1alpha binding. Two binding sites were mapped at the extracellular domain (Site 1) and inside the transmembrane domain (Site 2), which are composed of conserved residues. Sites 1 and 2 contribute approximately 60% and 40% to the binding affinity with SDF-1alpha, respectively. The binding model is in agreement with most of the experimental data. Transmembrane VI has more significant motion in the harmonious conformational transition of CXCR4 during SDF-1alpha binding, which may be possibly associated with signal transduction. Based on the modeling and simulation, a binding mechanism hypothesis between CXCR4 and SDF-1alpha and its relationship to the signal transduction has been proposed.     

Ligand unbinding pathways from the vitamin D receptor studied by molecular dynamics simulations

Molecular dynamics simulation techniques have been used to study the unbinding pathways of 1α,25-dihydroxyvitamin D₃ from the ligand-binding pocket of the vitamin D receptor (VDR). The pathways observed in a large number of relatively short (<200 ps) random acceleration molecular dynamics (RAMD) trajectories were found to be in fair agreement, both in terms of pathway locations and deduced relative preferences, compared to targeted molecular dynamics (TMD) and streered molecular dynamics simulations (SMD). However, the high-velocity ligand expulsions of RAMD tend to favor straight expulsion trajectories and the observed relative frequencies of different pathways were biased towards the probability of entering a particular exit channel. Simulations indicated that for VDR the unbinding pathway between the H1–H2 loop and the β-sheet between H5 and H6 is more favorable than the pathway located between the H1–H2 loop and H3. The latter pathway has been suggested to be the most likely unbinding path for thyroid hormone receptors (TRs) and a likely path for retinoic acid receptor. Ligand entry/exit through these two pathways would not require displacement of H12 from its agonistic position. Differences in the packing of the H1, H2, H3 and β-sheet region explain the changed relative preference of the two unbinding pathways in VDR and TRs. Based on the crystal structures of the ligand binding domains of class 2 nuclear receptors, whose members are VDR and TRs, this receptor class can be divided in two groups according to the packing of the H1, H2, H3 and β-sheet region. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A coupling of homology modeling with multiple molecular dynamics simulation for identifying representative conformation of GPCR structures: a case study on human bombesin receptor subtype-3

In this study, a computational pipeline was therefore devised to overcome homology modeling (HM) bottlenecks. The coupling of HM with molecular dynamics (MD) simulation is useful in that it tackles the sampling deficiency of dynamics simulations by providing good-quality initial guesses for the native structure. Indeed, HM also relaxes the severe requirement of force fields to explore the huge conformational space of protein structures. In this study, the interaction between the human bombesin receptor subtype-3 and MK-5046 was investigated integrating HM, molecular docking, and MD simulations. To improve conformational sampling in typical MD simulations of GPCRs, as in other biomolecules, multiple trajectories with different initial conditions can be employed rather than a single long trajectory. Multiple MD simulations of human bombesin receptor subtype-3 with different initial atomic velocities are applied to sample conformations in the vicinity of the structure generated by HM. The backbone atom conformational space distribution of replicates is analyzed employing principal components analysis. As a result, the averages of structural and dynamic properties over the twenty-one trajectories differ significantly from those obtained from individual trajectories.     

Constitutive activity and ligand-dependent activation of the nuclear receptor CAR-insights from molecular dynamics simulations

The constitutive androstane receptor (CAR) possesses, unlike most other nuclear receptors, a pronounced basal activity in vitro whose structural basis is still not fully understood. Using comparative molecular dynamics simulations of CAR X-ray crystal structures, we evaluated the molecular basis for constitutive activity and ligand-dependent receptor activation. Our results suggest that a combination of van der Waals interactions and hydrogen bonds is required to maintain the activation helix in the active conformation also in absence of a ligand. Furthermore, we identified conformational rearrangements within the ligand-binding pocket upon agonist binding and an influence of CAR inducers pregnanedione and CITCO on the helical conformation of the activation helix. Based on the results a model for ligand-dependent CAR activation is suggested.     

The gating mechanism of the bacterial mechanosensitive channel MscL revealed by molecular dynamics simulations

One of the ultimate goals of the study on mechanosensitive (MS) channels is to understand the biophysical mechanisms of how the MS channel protein senses forces and how the sensed force induces channel gating. The bacterial MS channel MscL is an ideal subject to reach this goal owing to its resolved 3D protein structure in the closed state on the atomic scale and large amounts of electrophysiological data on its gating kinetics. However, the structural basis of the dynamic process from the closed to open states in MscL is not fully understood. In this study, we performed molecular dynamics (MD) simulations on the initial process of MscL opening in response to a tension increase in the lipid bilayer. To identify the tension-sensing site(s) in the channel protein, we calculated interaction energy between membrane lipids and candidate amino acids (AAs) facing the lipids. We found that Phe78 has a conspicuous interaction with the lipids, suggesting that Phe78 is the primary tension sensor of MscL. Increased membrane tension by membrane stretch dragged radially the inner (TM1) and outer (TM2) helices of MscL at Phe78, and the force was transmitted to the pentagon-shaped gate that is formed by the crossing of the neighboring TM1 helices in the inner leaflet of the bilayer. The radial dragging force induced radial sliding of the crossing portions, leading to a gate expansion. Calculated energy for this expansion is comparable to an experimentally estimated energy difference between the closed and the first subconductance state, suggesting that our model simulates the initial step toward the full opening of MscL. The model also successfully mimicked the behaviors of a gain of function mutant (G22N) and a loss of function mutant (F78N), strongly supporting that our MD model did simulate some essential biophysical aspects of the mechano-gating in MscL.     

米曲霉木糖醇脱氢酶的分子模建与对接

[目的]研究米曲霉木糖醇脱氢酶基因的结构与功能.[方法]克隆测序来源于米曲霉的木糖醇脱氢酶(XDH)基因,利用Swiss-MODEL和Modeller对XDH进行三级结构模建,通过PROCHECK和Prosa2003对得到的4个目标模型进行评价,从中得到一个最佳模型.在同源建模的基础上,通过分子对接软件MolsoftICM-Pro,对辅因子进行对接,预测了XDH与NAD+、Zn2+作用的相关残基.寻找底物木糖醇与XDH结合的可能活性口袋,用Molsoft模拟XDH与木糖醇的对接,预测了酶与底物作用的关键氨基酸残基.[结果]结构分析显示,米曲霉XDH含有醇脱氢酶家族锌指纹结构和典型醇脱氢酶Rossmann折叠的辅酶结合域,属于Medium-chain脱氢酶(MDR)家族.通过对接研究,预测了XDH与NAD+之间形成氢键的氨基酸有Asp206、Arg211、Ser255、Ser301和Arg303,这些氨基酸位于结合域,与Zn2+形成氢键的氨基酸有His72和Glu73,位于催化域,与天然底物木糖醇形成氢键的氨基酸有Ile46、Ile349、Lys350和Thr351,位于催化域.[结论]所得信息对XDH分子定向改造、拓展米曲霉工业应用范围有重要意义.     

Using molecular dynamics for the refinement of atomistic models of GPCRs by homology modeling

Despite GPCRs sharing a common seven helix bundle, analysis of the diverse crystallographic structures available reveal specific features that might be relevant for ligand design. Despite the number of crystallographic structures of GPCRs steadily increasing, there are still challenges that hamper the availability of new structures. In the absence of a crystallographic structure, homology modeling remains one of the important techniques for constructing 3D models of proteins. In the present study we investigated the use of molecular dynamics simulations for the refinement of GPCRs models constructed by homology modeling. Specifically, we investigated the relevance of template selection, ligand inclusion as well as the length of the simulation on the quality of the GPCRs models constructed. For this purpose we chose the crystallographic structure of the rat muscarinic M3 receptor as reference and constructed diverse atomistic models by homology modeling, using different templates. Specifically, templates used in the present work include the human muscarinic M2; the more distant human histamine H1 and the even more distant bovine rhodopsin as shown in the GPCRs phylogenetic tree. We also investigated the use or not of a ligand in the refinement process. Hence, we conducted the refinement process of the M3 model using the M2 muscarinic as template with tiotropium or NMS docked in the orthosteric site and compared with the results obtained with a model refined without any ligand bound.     

Homology modeling of human serum paraoxonase1 and its molecular interaction studies with aspirin and cefazolin

Human serum paraoxonase1 (HuPON1) belongs to the family of A-esterases (EC.3.1.8.1). It is associated with HDL particle and prevents atherosclerosis by cleaving lipid hydroperoxides and other proatherogenic molecules of oxidized low density lipoproteins (LDL). Since the precise structure of HuPON1 is not yet available, the structure-function relationship between HuPON1 and activators/inhibitors is still unknown. Therefore, a theoretical model of HuPON1 was generated using homology modelling and precise molecular interactions of an activator aspirin and an inhibitor cefazolin with PON1 were studied using Autodock software. The ligand binding residues were found to be similar to the predicted active site residues. Both cefazolin and aspirin were found to dock in the vicinity of the predicted active sites of PON1; cefazolin bound at residues N166, S193 and Y71, while aspirin at residues N309, I310 and L311. Binding region in the PON1 by prediction (3D2GO server) and docking studies provide useful insight into mechanism of substrate and inhibitor binding to the enzyme active site.     

Ligand-supported homology modeling of the human angiotensin II type 1 (AT(1)) receptor: insights into the molecular determinants of telmisartan binding

Angiotensin II type 1 (AT(1)) receptor belongs to the super-family of G-protein-coupled receptors, and antagonists of the AT(1) receptor are effectively used in the treatment of hypertension. To understand the molecular interactions of these antagonists, such as losartan and telmisartan, with the AT(1) receptor, a homology model of the human AT(1) (hAT(1)) receptor with all connecting loops was constructed from the 2.6 A resolution crystal structure (PDB i.d., 1L9H) of bovine rhodopsin. The initial model generated by MODELLER was subjected to a stepwise ligand-supported model refinement. This protocol involved initial docking of non-peptide AT(1) antagonists in the putative binding site, followed by several rounds of iterative energy minimizations and molecular dynamics simulations. The final model was validated based on its correlation with several structure-activity relationships and site-directed mutagenesis data. The final model was also found to be in agreement with a previously reported AT(1) antagonist pharmacophore model. Docking studies were performed for a series of non-peptide AT(1) receptor antagonists in the active site of the final hAT(1) receptor model. The docking was able to identify key molecular interactions for all the AT(1) antagonists studied. Reasonable correlation was observed between the interaction energy values and the corresponding binding affinities of these ligands, providing further validation for the model. In addition, an extensive unrestrained molecular dynamics simulation showed that the docking-derived bound pose of telmisartan is energetically stable. Knowledge gained from the present studies can be used in structure-based drug design for developing novel ligands for the AT(1) receptor.