Identification of Small-Molecule Antagonists Targeting the Growth Hormone Releasing Hormone Receptor (GHRHR).

The growth hormone-releasing hormone receptor (GHRHR) belongs to Class B1 of G protein-coupled receptors (GPCRs). Class B1 GPCR peptides such, as growth hormone-releasing hormone (GHRH), have been proposed to bind in a two-step model, where first the C-terminal region of the peptide interacts with the extracellular domain of the receptor and, subsequently, the N-terminus interacts with the seven transmembrane domain of the receptor, resulting in activation. The GHRHR has recently been highlighted as a promising drug target toward several types of cancer and has been shown to be overexpressed in prostate, breast, pancreatic, and ovarian cancer. Indeed, peptide GHRHR antagonists have displayed promising results in many cancer models. However, no nonpeptide GHRHR-targeting compounds have yet been identified. We have utilized several computational tools to target GHRHR and identify potential small-molecule compounds directed at this receptor. These compounds were validated in vitro using a cyclic adenosine monophosphate (cAMP) ELISA to measure activity at the GHRHR. In vitro results suggest that several of the novel small-molecule compounds could inhibit GHRH-induced cAMP accumulation. Preliminary analysis of the specificity/selectivity of one of the most effective hit compounds indicated that the effect seen was via inhibition of the GHRHR. We therefore report the first nonpeptide antagonists of GHRHR and propose a structural basis for inhibition induced by the compounds, which may assist in the future design of lead GHRHR compounds for treating disorders attributed to dysregulated/aberrant GHRHR signaling.

Density functional theory and enzyme studies support interactions between angiotensin receptor blockers and angiotensin converting enzyme-2: Relevance to coronavirus 2019.

The binding affinities and interactions between eight drug candidates, both commercially available (candesartan; losartan; losartan carboxylic acid; nirmatrelvir; telmisartan) and newly synthesized benzimidazole-N-biphenyltetrazole (ACC519T), benzimidazole bis-N,N'-biphenyltetrazole (ACC519T(2) and 4-butyl-N,N-bis([2-(2H-tetrazol-5-yl)biphenyl-4-yl]) methyl (BV6), and the active site of angiotensin-converting enzyme-2 (ACE2) were evaluated for their potential as inhibitors against SARS-CoV-2 and regulators of ACE2 function through Density Functional Theory methodology and enzyme activity assays, respectively. Notably, telmisartan and ACC519T(2) exhibited pronounced binding affinities, forming strong interactions with ACE2's active center, favorably accepting proton from the guanidinium group of arginine273. The ordering of candidates by binding affinity and reactivity descriptors, emerged as telmisartan > ACC519T(2) > candesartan > ACC519T > losartan carboxylic acid > BV6 > losartan > nirmatrelvir. Proton transfers among the active center amino acids revealed their interconnectedness, highlighting a chain-like proton transfer involving tyrosine, phenylalanine, and histidine. Furthermore, these candidates revealed their potential antiviral abilities by influencing proton transfer within the ACE2 active site. Furthermore, through an in vitro pharmacological assays we determined that candesartan and the BV6 derivative, 4-butyl-N,N0-bis[20-2Htetrazol-5-yl)bipheyl-4-yl]methyl)imidazolium bromide (BV6(K+)2) also contain the capacity to increase ACE2 functional activity. This comprehensive analysis collectively underscores the promise of these compounds as potential therapeutic agents against SARS-CoV-2 by targeting crucial protein interactions.

Natural Compounds for Bone Remodeling: A Computational and Experimental Approach Targeting Bone Metabolism-Related Proteins.

Osteoporosis, characterized by reduced bone density and increased fracture risk, affects over 200 million people worldwide, predominantly older adults and postmenopausal women. The disruption of the balance between bone-forming osteoblasts and bone-resorbing osteoclasts underlies osteoporosis pathophysiology. Standard treatment includes lifestyle modifications, calcium and vitamin D supplementation and specific drugs that either inhibit osteoclasts or stimulate osteoblasts. However, these treatments have limitations, including side effects and compliance issues. Natural products have emerged as potential osteoporosis therapeutics, but their mechanisms of action remain poorly understood. In this study, we investigate the efficacy of natural compounds in modulating molecular targets relevant to osteoporosis, focusing on the Mitogen-Activated Protein Kinase (MAPK) pathway and the gut microbiome's influence on bone homeostasis. Using an in silico and in vitro methodology, we have identified quercetin as a promising candidate in modulating MAPK activity, offering a potential therapeutic perspective for osteoporosis treatment.

Structural Features Influencing the Bioactive Conformation of Angiotensin II and Angiotensin A: Relationship between Receptor Desensitization, Addiction, and the Blood-Brain Barrier.

The N-terminal portion of the octapeptide angiotensin II (DRVYIHPF; AngII), a vasopressor peptide that favorably binds to, and activates, AngII type 1 receptor (AT1R), has an important role in maintaining bioactive conformation. It involves all three charged groups, namely (i) the N-terminal amino group cation, (ii) the Asp sidechain anion and (iii) the Arg guanidino cation. Neutralization of any one of these three charged groups results in a substantial reduction (<5%) in bioactivity, implicating a specialized function for this cluster. In contrast, angiotensin A (ARVYIHPF; AngA) has reduced bioactivity at AT1R; however, replacement of Asp in AngII with sarcosine (N-methyl-glycine) not only restores bioactivity but increases the activity of agonist, antagonist, and inverse agonist analogues. A bend produced at the N-terminus by the introduction of the secondary amino acid sarcosine is thought to realign the functional groups that chaperone the C-terminal portion of AngII, allowing transfer of the negative charge originating at the C-terminus to be transferred to the Tyr hydroxyl-forming tyrosinate anion, which is required to activate the receptor and desensitizes the receptor (tachyphylaxis). Peptide (sarilesin) and nonpeptide (sartans) moieties, which are long-acting inverse agonists, appear to desensitize the receptor by a mechanism analogous to tachyphylaxis. Sartans/bisartans were found to bind to alpha adrenergic receptors resulting in structure-dependent desensitization or resensitization. These considerations have provided information on the mechanisms of receptor desensitization/tolerance and insights into possible avenues for treating addiction. In this regard sartans, which appear to cross the blood-brain barrier more readily than bisartans, are the preferred drug candidates.

Publisher Correction: GPCRmd uncovers the dynamics of the 3D-GPCRome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

GPCRmd uncovers the dynamics of the 3D-GPCRome.

G-protein-coupled receptors (GPCRs) are involved in numerous physiological processes and are the most frequent targets of approved drugs. The explosion in the number of new three-dimensional (3D) molecular structures of GPCRs (3D-GPCRome) over the last decade has greatly advanced the mechanistic understanding and drug design opportunities for this protein family. Molecular dynamics (MD) simulations have become a widely established technique for exploring the conformational landscape of proteins at an atomic level. However, the analysis and visualization of MD simulations require efficient storage resources and specialized software. Here we present GPCRmd (http://gpcrmd.org/), an online platform that incorporates web-based visualization capabilities as well as a comprehensive and user-friendly analysis toolbox that allows scientists from different disciplines to visualize, analyze and share GPCR MD data. GPCRmd originates from a community-driven effort to create an open, interactive and standardized database of GPCR MD simulations.

Identification of compounds that bind the centriolar protein SAS-6 and inhibit its oligomerization.

Centrioles are key eukaryotic organelles that are responsible for the formation of cilia and flagella, and for organizing the microtubule network and the mitotic spindle in animals. Centriole assembly requires oligomerization of the essential protein spindle assembly abnormal 6 (SAS-6), which forms a structural scaffold templating the organization of further organelle components. A dimerization interaction between SAS-6 N-terminal "head" domains was previously shown to be essential for protein oligomerization in vitro and for function in centriole assembly. Here, we developed a pharmacophore model allowing us to assemble a library of low-molecular-weight ligands predicted to bind the SAS-6 head domain and inhibit protein oligomerization. We demonstrate using NMR spectroscopy that a ligand from this family binds at the head domain dimerization site of algae, nematode, and human SAS-6 variants, but also that another ligand specifically recognizes human SAS-6. Atomistic molecular dynamics simulations starting from SAS-6 head domain crystallographic structures, including that of the human head domain which we now resolve, suggest that ligand specificity derives from favorable Van der Waals interactions with a hydrophobic cavity at the dimerization site.

A COVID moonshot: assessment of ligand binding to the SARS-CoV-2 main protease by saturation transfer difference NMR spectroscopy.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the etiological cause of the coronavirus disease 2019, for which no effective antiviral therapeutics are available. The SARS-CoV-2 main protease (Mpro) is essential for viral replication and constitutes a promising therapeutic target. Many efforts aimed at deriving effective Mpro inhibitors are currently underway, including an international open-science discovery project, codenamed COVID Moonshot. As part of COVID Moonshot, we used saturation transfer difference nuclear magnetic resonance (STD-NMR) spectroscopy to assess the binding of putative Mpro ligands to the viral protease, including molecules identified by crystallographic fragment screening and novel compounds designed as Mpro inhibitors. In this manner, we aimed to complement enzymatic activity assays of Mpro performed by other groups with information on ligand affinity. We have made the Mpro STD-NMR data publicly available. Here, we provide detailed information on the NMR protocols used and challenges faced, thereby placing these data into context. Our goal is to assist the interpretation of Mpro STD-NMR data, thereby accelerating ongoing drug design efforts.

Nonconventional glucagon and GLP-1 receptor agonist and antagonist interplay at the GLP-1 receptor revealed in high-throughput FRET assays for cAMP.

G protein-coupled receptors (GPCRs) for glucagon (GluR) and glucagon-like peptide-1 (GLP-1R) are normally considered to be highly selective for glucagon and GLP-1, respectively. However, glucagon secreted from pancreatic α-cells may accumulate at high concentrations to exert promiscuous effects at the ß-cell GLP-1R, as may occur in the volume-restricted microenvironment of the islets of Langerhans. Furthermore, systemic administration of GluR or GLP-1R agonists and antagonists at high doses may lead to off-target effects at other receptors. Here, we used molecular modeling to evaluate data derived from FRET assays that detect cAMP as a read-out for GluR and GLP-1R activation. This analysis established that glucagon is a nonconventional GLP-1R agonist, an effect inhibited by the GLP-1R orthosteric antagonist exendin(9-39) (Ex(9-39)). The GluR allosteric inhibitors LY2409021 and MK 0893 antagonized glucagon and GLP-1 action at the GLP-1R, whereas des-His1-[Glu9]glucagon antagonized glucagon action at the GluR, while having minimal inhibitory action versus glucagon or GLP-1 at the GLP-1R. When testing Ex(9-39) in combination with des-His1-[Glu9]glucagon in INS-1 832/13 cells, we validated a dual agonist action of glucagon at the GluR and GLP-1R. Hybrid peptide GGP817 containing glucagon fused to a fragment of peptide YY (PYY) acted as a triagonist at the GluR, GLP-1R, and neuropeptide Y2 receptor (NPY2R). Collectively, these findings provide a new triagonist strategy with which to target the GluR, GLP-1R, and NPY2R. They also provide an impetus to reevaluate prior studies in which GluR and GLP-1R agonists and antagonists were assumed not to exert promiscuous actions at other GPCRs.

Combined Computational and Structural Approach into Understanding the Role of Peptide Binding and Activation of Melanocortin Receptor 4.

Melanocortin receptor 4 (MC4R) is expressed predominantly in the central nervous system and regulates food intake and sexual function and is also thought to be responsible for effects on mood and cognition. It belongs to the melanocortin receptor subfamily of G protein-coupled receptors (GPCRs). Here, we have synthesized and structurally characterized three peptides that bind to MC4R, producing different signaling events. AgRP is a naturally occurring antagonist, HLWNRS is the minimal sequence of the N-terminal with partial agonist activity, and aMSH is a full agonistic peptide. By implementing molecular dynamics simulations on the different peptide-receptor complexes, we propose their molecular basis of binding to investigate their differential molecular properties regarding the activation states of the receptor. Our analysis shows that the agonist and partial agonist may induce rotation in transmembrane helix 3, which is known to be involved in the key events occurring during GPCR activation, and this movement is impacted by certain aromatic residues and their positioning in the orthosteric binding site of the receptor.

Conformational plasticity of the VEEV macro domain is important for binding of ADP-ribose.

Venezuelan equine encephalitis virus (VEEV) is a new world alphavirus which can be involved in several central nervous system disorders such as encephalitis and meningitis. The VEEV genome codes for 4 non-structural proteins (nsP), of which nsP3 contains a Macro domain. Macro domains (MD) can be found as stand-alone proteins or embedded within larger proteins in viruses, bacteria and eukaryotes. Their most common feature is the binding of ADP-ribose (ADPr), while several macro domains act as ribosylation writers, erasers or readers. Alphavirus MD erase ribosylation but their precise contribution in viral replication is still under investigation. NMR-driven titration experiments of ADPr in solution with the VEEV macro domain (in apo- and complex state) show that it adopts a suitable conformation for ADPr binding. Specific experiments indicate that the flexibility of the loops ß5-α3 and α3-ß6 is critical for formation of the complex and assists a wrapping mechanism for ADPr binding. Furthermore, along with this sequence of events, the VEEV MD undergoes a conformational exchange process between the apo state and a low-populated "dark" conformational state.

Lysophosphatidic acid produced by autotaxin acts as an allosteric modulator of its catalytic efficiency.

Autotaxin (ATX) is a secreted glycoprotein and the only member of the ectonucleotide pyrophosphatase/phosphodiesterase family that converts lysophosphatidylcholine (LPC) into lysophosphatidic acid (LPA). LPA controls key responses, such as cell migration, proliferation, and survival, implicating ATX-LPA signaling in various (patho)physiological processes and establishing it as a drug target. ATX structural and functional studies have revealed an orthosteric and an allosteric site, called the "pocket" and the "tunnel," respectively. However, the mechanisms in allosteric modulation of ATX's activity as a lysophospholipase D are unclear. Here, using the physiological LPC substrate, a new fluorescent substrate, and diverse ATX inhibitors, we revisited the kinetics and allosteric regulation of the ATX catalytic cycle, dissecting the different steps and pathways leading to LPC hydrolysis. We found that ATX activity is stimulated by LPA and that LPA activates ATX lysophospholipase D activity by binding to the ATX tunnel. A consolidation of all experimental kinetics data yielded a comprehensive catalytic model supported by molecular modeling simulations and suggested a positive feedback mechanism that is regulated by the abundance of the LPA products activating hydrolysis of different LPC species. Our results complement and extend the current understanding of ATX hydrolysis in light of the allosteric regulation by ATX-produced LPA species and have implications for the design and application of both orthosteric and allosteric ATX inhibitors.

The size matters? A computational tool to design bivalent ligands.

Motivation: Bivalent ligands are increasingly important such as for targeting G protein-coupled receptor (GPCR) dimers or proteolysis targeting chimeras (PROTACs). They contain two pharmacophoric units that simultaneously bind in their corresponding binding sites, connected with a spacer chain. Here, we report a molecular modelling tool that links the pharmacophore units via the shortest pathway along the receptors van der Waals surface and then scores the solutions providing prioritization for the design of new bivalent ligands. Results: Bivalent ligands of known dimers of GPCRs, PROTACs and a model bivalent antibody/antigen system were analysed. The tool could rapidly assess the preferred linker length for the different systems and recapitulated the best reported results. In the case of GPCR dimers the results suggest that in some cases these ligands might bind to a secondary binding site at the extracellular entrance (vestibule or allosteric site) instead of the orthosteric binding site. Availability and implementation: Freely accessible from the Molecular Operating Environment svl exchange server (https://svl.chemcomp.com/). Supplementary information: Supplementary data are available at Bioinformatics online.

Novel stable analogues of the neurotensin C-terminal hexapeptide containing unnatural amino acids.

Neurotensin (NT) (pGlu-Leu-Tyr-Glu-Asn-Lys-Pro-Arg-Arg-Pro-Tyr-Ile-Leu) exerts a dual function as a neurotransmitter/neuromodulator in the central nervous system and as a hormone/cellular mediator in periphery. This dual function of NT establishes a connection between brain and peripheral tissues that renders this peptide a central player in energy homeostasis. Many biological actions of NT are mediated through its interaction with three types of NT receptors (NTS receptors). Despite its role in energy homeostasis, NT has a short half-life that hampers further determination of the biological actions of this peptide and its receptors in brain and periphery. The short half-life of NT is due to the proteolytic degradation of its C-terminal side by several endopeptidases. Therefore, it is important to synthesize NT analogues with resistant bonds against metabolic deactivation. Based on these findings, we herein report the synthesis of ten linear, two cyclic and two dimeric analogues of NT with modifications in its structure that improve their metabolic stability, while retaining the ability to bind to NTS receptors. Modifications at position 11 (introduction of D-Tyrosine (OEthyl) [D-Tyr(Et)] or D-1-naphtylalanine [D-1-Nal] were combined with introduction of a L-Lysine or a D-Arginine at positions 8 or 9, and 1-[2-(aminophenyl)-2-oxoethyl]-1H-pyrrole-2-carboxylic acid (AOPC) at positions 7 or 8, resulting in compounds NT4-NT21. AOPC is an unnatural amino acid with promise in applications as a building block for the synthesis of peptidomimetic compounds. To biologically evaluate these analogues, we determined their plasma stability and their binding affinities to type 1 NT receptor (NTS1), endogenously expressed in HT-29 cells, Among the fourteen NT analogues, compounds, NT5, NT6, and NT8, which have D-Tyr(Et) at position 11, bound to NTS1 in a dose-response manner and with relatively high affinity but still lower than that of the natural peptide. Despite their lower binding affinities compared to NT, the NT5, NT6, and NT8 exhibited a remarkably higher stability, as a result of their chemistry, which provides protection from enzymatic activity. These results will set the basis for the rational design of novel NT molecules with improved pharmacological properties and enhanced enzymatic stability.

Novel Hit Compounds as Putative Antifungals: The Case of Aspergillus fumigatus.

The prevalence of invasive fungal infections has been dramatically increased as the size of the immunocompromised population worldwide has grown. Aspergillus fumigatus is characterized as one of the most widespread and ubiquitous fungal pathogens. Among antifungal drugs, azoles have been the most widely used category for the treatment of fungal infections. However, increasingly, azole-resistant strains constitute a major problem to be faced. Towards this direction, our study focused on the identification of compounds bearing novel structural motifs which may evolve as a new class of antifungals. To fulfil this scope, a combination of in silico techniques and in vitro assays were implemented. Specifically, a ligand-based pharmacophore model was created and served as a 3D search query to screen the ZINC chemical database. Additionally, molecular docking and molecular dynamics simulations were used to improve the reliability and accuracy of virtual screening results. In total, eight compounds, bearing completely different chemical scaffolds from the commercially available azoles, were proposed and their antifungal activity was evaluated using in vitro assays. Results indicated that all tested compounds exhibit antifungal activity, especially compounds 1, 2, and 4, which presented the most promising minimum inhibitory concentration (MIC) and minimum fungicidal concentration (MFC) values and, therefore, could be subjected to further hit to lead optimization.

Ligand-Triggered Structural Changes in the M2 Muscarinic Acetylcholine Receptor.

The muscarinic M2 acetylcholine receptor, one of the few G-protein coupled receptors that has not only been crystallized in both active and inactive conformations but also in the presence of a positive allosteric modulator, is an interesting system to study the molecular mechanisms of GPCR activation and ligand allosterism. Here, we have employed molecular dynamics (MD) simulations (adding to 14 µs in total) to study conformational changes triggered by the inverse agonist R-(-)-3-quinuclidinyl-benzilate (QNB) in the structure of the active M2 receptor (PBD ID 4MQS ) after replacement of the agonist iperoxo by the inverse agonist QNB. This permitted us to identify the sequence of events in the deactivation mechanism of the M2 acetylcholine receptor, which results first in the rearrangement of the transmission switch, the subsequent opening of the extracellular portion of the receptor and finally, the closure of the intracellular part. We also evaluate the effect of the positive allosteric modulator LY2119620 when bound simultaneously with the orthosteric agonist iperoxo and find that it restricts the conformation of Trp4227.35 in a position that modulates the orientation of the Tyr4267.39 at the orthosteric-binding pocket.

Cyclization of PLP139-151 peptide reduces its encephalitogenic potential in experimental autoimmune encephalomyelitis.

We report the novel synthesis of cyclic PLP139-151 (cPLP) and its application in SJL/J mice to study its encephalitogenic effects. Our results indicate that the cPLP analog is minimally encephalitogenic when administered to induce experimental autoimmune encephalomyelitis (low disease burden, minimal inflammatory, demyelinating and axonopathic pathology compared to its linear counterpart). Proliferation assays confirmed the low stimulatory potential of the cPLP compared to linPLP (2.5-fold lower proliferation) as well as inducing lower antibody responses. Molecular modeling showed a completely different TCR recognition profile of cPLP in regard to linPLP, where H147 replaces W144 and F151-K150 replace H147 as TCR contacts, which may explain the difference on each peptide's response.

Cyclic MOG35-55 ameliorates clinical and neuropathological features of experimental autoimmune encephalomyelitis.

EAE is induced to susceptible mice using linear peptides of myelin proteins of the central nervous system. Specific peptide motifs within the peptide-binding groove of the MHC peptide-complex determines the affinity of the peptide in each animal and the consequent T-cell receptor recognition and activation of the cell. Altered peptide ligand (APL) vaccination is a novel approach based on an effort to induce T-cell tolerance or alter cytokine profile from pro-inflammatory to anti-inflammatory. In the present study we synthesized the MOG35-55 peptide and altered its 3-dimensional conformation to make it a cyclic one (c-MOG35-55). EAE was induced in C57BL/6 mice and pathology was studied on acute and chronic phase of the disease. Our data indicates that c-MOG35-55 peptide alone induces a mild transient acute phase without chronic axonopathy. Administration of the c-MOG35-55 peptide at a 1:1 ratio during disease induction significantly ameliorates clinical disease and underlying pathology, such as demyelination and axonopathy in the acute and chronic phases. Binding and structural studies revealed milder interactions between the c-MOG35-55 and mouse or human MHC class II alleles (H2-IAb and HLA-DR2). Collectively, we provide data supporting for the first time the concept that the cyclic modification of an established encephalitogenic peptide ameliorates the clinical outcomes and underlying pathological processes of EAE. Such a cyclic modification of linear peptides could provide a novel treatment approach for future, patient-selective, immunomodulative treatments of multiple sclerosis.

Cyclic citrullinated MBP87-99 peptide stimulates T cell responses: Implications in triggering disease.

Amino acid mutations to agonist peptide epitopes of myelin proteins have been used to modulate immune responses and experimental autoimmune encephalomyelitis (EAE, animal model of multiple sclerosis). Such amino acid alteration are termed, altered peptide ligands (APL). We have shown that the agonist myelin basic protein (MBP) 87-99 epitope (MBP87-99) with crucial T cell receptor (TCR) substitutions at positions 91 and 96 (K91,P96 (TCR contact residues) to R91,A96; [R91,A96]MBP87-99) results in altered T cell responses and inhibits EAE symptoms. In this study, the role of citrullination of arginines in [R91,A96]MBP87-99 peptide analog was determined using in vivo experiments in combination with computational studies. The immunogenicity of linear [Cit91,A96,Cit97]MBP87-99 and its cyclic analog - cyclo(87-99)[Cit91,A96,Cit97]MBP87-99 when conjugated to the carrier mannan (polysaccharide) were studied in SJL/J mice. It was found that mannosylated cyclo(87-99)[Cit91,A96,Cit97]MBP87-99 peptide induced strong T cell proliferative responses and IFN-gamma cytokine secretion compared with the linear one. Moreover, the interaction of linear and cyclic peptide analogs with the major histocompatibility complex (MHC II, H2-IAs) and TCR was analyzed using molecular dynamics simulations at the receptor level, in order to gain a better understanding of the molecular recognition mechanisms that underly the different immunological profiles of citrullinated peptides compared to its agonist native counterpart MBP87-99 epitope. The results demonstrate that the citrullination of arginine in combination with the backbone conformation of mutated linear and cyclic analogs are significant elements for the immune response triggering the induction of pro-inflammatory cytokines.

Exploring new scaffolds for angiotensin II receptor antagonism.

Nowadays, AT1 receptor (AT1R) antagonists (ARBs) constitute the one of the most prevalent classes of antihypertensive drugs that modulate the renin-angiotensin system (RAS). Their main uses include also treatment of diabetic nephropathy (kidney damage due to diabetes) and congestive heart failure. Towards this direction, our study has been focused on the discovery of novel agents bearing different scaffolds which may evolve as a new class of AT1 receptor antagonists. To fulfill this aim, a combination of computational approaches and biological assays were implemented. Particularly, a pharmacophore model was established and served as a 3D search query to screen the ChEMBL15 database. The reliability and accuracy of virtual screening results were improved by using molecular docking studies. In total, 4 compounds with completely diverse chemical scaffolds from potential ARBs, were picked and tested for their binding affinity to AT1 receptor. Results revealed high nanomolar to micromolar affinity (IC50) for all the compounds. Especially, compound 4 exhibited a binding affinity of 199nM. Molecular dynamics simulations were utilized in an effort to provide a molecular basis of their binding to AT1R in accordance to their biological activities.