Association of Beta 2 adrenergic receptor (Thr164Ile) polymorphism with Salbutamol refractoriness in severe asthmatics from Indian population

The impact of genomic variants on patient response to inhaled bronchodilators: a comprehensive update

INTRODUCTION

The bronchodilator response (BDR) depends on many factors, including genetic ones. Numerous single nucleotide polymorphisms (SNPs) influencing BDR have been identified. However, despite several studies in this field, genetic variations are not currently being utilized to support the use of bronchodilators.

AREAS COVERED

In this narrative review, the possible impact of genetic variants on BDR is discussed.

EXPERT OPINION

Pharmacogenetic studies of β₂-agonists have mainly focused on ADRB2 gene. Three SNPs, A46G, C79G, and C491T, have functional significance. However, other uncommon variants may contribute to individual variability in salbutamol response. SNPs haplotypes in ADRB2 may have a role. Many variants in genes coding for muscarinic ACh receptor (mAChR) have been reported, particularly in the M₂ and, to a lesser degree, M₃ mAChRs, but no consistent evidence for a pharmacological relevance of these SNPs has been reported. Moreover, there is a link between SNPs and ethnic and/or age profiles regarding BDR. Nevertheless, replication of pharmacogenetic results is limited and often, BDR is dissociated from what is expected based on SNP identification. Pharmacogenetic studies on bronchodilators must continue. However, they must integrate data derived from a multi-omics approach with epigenetic factors that may modify BDR.

Molecular determinants of GPCR pharmacogenetics: Deconstructing the population variants in β2-adrenergic receptor

G protein-coupled receptors (GPCRs) are membrane proteins that play a central role in cell signaling and constitute one of the largest classes of drug targets. The molecular mechanisms underlying GPCR function have been characterized by several experimental and computational methods and provide an understanding of their role in physiology and disease. Population variants arising from nsSNPs affect the native function of GPCRs and have been implicated in differential drug response. In this chapter, we provide an overview on GPCR structure and activation, with a special focus on the β₂-adrenergic receptor (β₂-AR). First, we discuss the current understanding of the structural and dynamic features of the wildtype receptor. Subsequently, the population variants identified in this receptor from clinical and large-scale genomic studies are described. We show how computational approaches such as bioinformatics tools and molecular dynamics simulations can be used to characterize the variant receptors in comparison to the wildtype receptor. In particular, we discuss three examples of clinically important variants and discuss how the structure and function of these variants differ from the wildtype receptor at a molecular level. Overall, the chapter provides an overview of structure and function of GPCR variants and is a step towards the study of inter-individual differences and personalized medicine.

Ligand-specific pharmacogenetic effects of nonsynonymous mutations

In pharmacogenomics, variable receptor phenotypes, resulting from genetic polymorphisms, are often described as a change in protein function or regulation observed upon exposure to a drug. However, in some instances, phenotypes are defined using a class of medications rather than individual drugs. This paradigm assumes that a variation associated with a drug response phenotype will retain the magnitude and direction of the effect for other drugs with the same mechanism of action. However, nonsynonymous polymorphisms may have ligand-specific effects. The purpose of this study was to investigate the potential for point mutations to asymmetrically affect the binding of different drugs to a common target. Ligand binding data from site-directed mutagenesis studies on five G-protein coupled receptors (beta-1 and -2 adrenergic, dopamine D2, angiotensin II and mu-opioid receptor) were collected and analyzed. Binding data from 81 studies for 253 ligands with 447 mutant proteins, including 10 naturally occurring human variants, were analyzed, yielding 1989 mutation-ligand pairs. Fold change in binding affinity for mutant proteins, relative to the wild-type, for different drugs was examined for ligand-specific effects, with a fold-change difference of one or more orders of magnitude between agents considered significant. Of the mutations examined, 49% were associated with ligand-specific effects. One human variant (T164I, beta-2 adrenergic receptor) showed ligand-specific effects for antiasthmatic agents. These results indicate that ligand-specific changes in binding are a possible consequence of missense mutations. This implies that caution needs to be exercised when grouping drugs together during design or interpretation of genotype-phenotype association studies.

A Multi-target Drug Designing for BTK, MMP9, Proteasome and TAK1 for the Clinical Treatment of Mantle Cell Lymphoma

BACKGROUND

Mantle cell lymphoma (MCL) is a type of non-Hodgkin lymphoma characterized by the mutation and overexpression of the cyclin D1 protein by the reciprocal chromosomal translocation t(11;14)(q13:q32).

AIM

The present study aims to identify potential inhibition of MMP9, Proteasome, BTK, and TAK1 and determine the most suitable and effective protein target for the MCL.

METHODOLOGY

Nine known inhibitors for MMP9, 24 for proteasome, 15 for BTK and 14 for TAK1 were screened. SB-3CT (PubChem ID: 9883002), oprozomib (PubChem ID: 25067547), zanubrutinib (PubChem ID: 135565884) and TAK1 inhibitor (PubChem ID: 66760355) were recognized as drugs with high binding capacity with their respective protein receptors. 41, 72, 102 and 3 virtual screened compounds were obtained after the similarity search with compound (PubChem ID:102173753), PubChem compound SCHEMBL15569297 (PubChem ID:72374403), PubChem compound SCHEMBL17075298 (PubChem ID:136970120) and compound CID: 71814473 with best virtual screened compounds.

RESULT

MMP9 inhibitors show commendable affinity and good interaction profile of compound holding PubChem ID:102173753 over the most effective established inhibitor SB-3CT. The pharmacophore study of the best virtual screened compound reveals its high efficacy based on various interactions. The virtual screened compound's better affinity with the target MMP9 protein was deduced using toxicity and integration profile studies.

CONCLUSION

Based on the ADMET profile, the compound (PubChem ID: 102173753) could be a potent drug for MCL treatment. Similar to the established SB-3CT, the compound was non-toxic with LD50 values for both the compounds lying in the same range.

Loss of a water-mediated network results in reduced agonist affinity in a β2-adrenergic receptor clinical variant

The β₂-adrenergic receptor (β₂AR) is a member of the G protein-coupled receptor (GPCR) family that is an important drug target for asthma and COPD. Clinical studies coupled with biochemical data have identified a critical receptor variant, Thr164Ile, to have a reduced response to agonist-based therapy, although the molecular mechanism underlying this seemingly "non-deleterious" substitution is not clear. Here, we couple molecular dynamics simulations with network analysis and free-energy calculations to identify the molecular determinants underlying the differential drug response. We are able to identify hydration sites in the transmembrane domain that are essential to maintain the integrity of the binding site but are absent in the variant. The loss of these hydration sites in the variant correlates with perturbations in the intra-protein interaction network and rearrangements in the orthosteric ligand binding site. In conjunction, we observe an altered binding and reduced free energy of a series of agonists, in line with experimental trends. Our work identifies a functional allosteric pathway connected by specific hydration sites in β₂AR that has not been reported before and provides insight into water-mediated networks in GPCRs in general. Overall, the work is one of the first step towards developing variant-specific potent and selective agonists.

Most Promising Compounds for Treating COVID-19 and Recent Trends in Antimicrobial & Antifungal Agents

Multidrug resistance in microbes poses a major health crisis and demands for the discovery of novel antimicrobial agents. The recent pandemic of SARS-CoV-2 has raised a public health emergency in almost all the countries of the world. Unlike viruses, a bacterium plays a significant role in various environmental issues such as bioremediation. Furthermore, biosurfactants produced by various bacterial species have an edge over traditionally produced chemical surfactants for its biodegradability, low toxicity and better interfacial activity with various applications in agriculture and industry. This special issue focuses on the global perspective of drug discovery for various antimicrobial, antiviral, and antifungal agents for infectious diseases. The issue also emphasizes the ongoing developments and the role of microbes in environmental remediation. We wish the articles published in this issue will enhance the current understanding in microbiology among the readers, and serve as the "seed of an idea" for drug development for ongoing COVID-19 pandemic.

Shape-based Machine Learning Models for the Potential Novel COVID-19 Protease Inhibitors Assisted by Molecular Dynamics Simulation

BACKGROUND

The vast geographical expansion of novel coronavirus and an increasing number of COVID-19 affected cases have overwhelmed health and public health services. Artificial Intelligence (AI) and Machine Learning (ML) algorithms have extended their major role in tracking disease patterns, and in identifying possible treatments.

OBJECTIVE

This study aims to identify potential COVID-19 protease inhibitors through shape-based Machine Learning assisted by Molecular Docking and Molecular Dynamics simulations.

METHODS

31 Repurposed compounds have been selected targeting the main coronavirus protease (6LU7) and a machine learning approach was employed to generate shape-based molecules starting from the 3D shape to the pharmacophoric features of their seed compound. Ligand-Receptor Docking was performed with Optimized Potential for Liquid Simulations (OPLS) algorithms to identify highaffinity compounds from the list of selected candidates for 6LU7, which were subjected to Molecular Dynamic Simulations followed by ADMET studies and other analyses.

RESULTS

Shape-based Machine learning reported remdesivir, valrubicin, aprepitant, and fulvestrant as the best therapeutic agents with the highest affinity for the target protein. Among the best shape-based compounds, a novel compound identified was not indexed in any chemical databases (PubChem, Zinc, or ChEMBL). Hence, the novel compound was named 'nCorv-EMBS'. Further, toxicity analysis showed nCorv-EMBS to be suitable for further consideration as the main protease inhibitor in COVID-19.

CONCLUSION

Effective ACE-II, GAK, AAK1, and protease 3C blockers can serve as a novel therapeutic approach to block the binding and attachment of the main COVID-19 protease (PDB ID: 6LU7) to the host cell and thus inhibit the infection at AT2 receptors in the lung. The novel compound nCorv- EMBS herein proposed stands as a promising inhibitor to be evaluated further for COVID-19 treatment.

Computer-Aided Structure Prediction of Bluetongue Virus Coat Protein VP2 Assisted by Optimized Potential for Liquid Simulations (OPLS)

Background:

The capsid coated protein of Bluetongue virus (BTV) VP2 is responsible for BTV transmission by the Culicoides vector to vertebrate hosts. Besides, VP2 is responsible for BTV entry into permissive cells and hence plays a major role in disease progression. However, its mechanism of action is still unknown.

Objective:

The present investigation aimed to predict the 3D structure of Viral Protein 2 of the bluetongue virus assisted by Optimized Potential for Liquid Simulations (OPLS), structure validation, and an active site prediction.

Methods:

The 3D structure of the VP2 protein was built using a Python-based Computational algorithm. The templates were identified using Smith waterman’s Local alignment. The VP2 protein structure validated using PROCHECK. Molecular Dynamics Simulation (MDS) studies were performed using an academic software Desmond, Schrodinger dynamics, for determining the stability of a model protein. The Ligand-Binding site was predicted by structure comparison using homology search and proteinprotein network analysis to reveal their stability and inhibition mechanism, followed by the active site identification.

Results:

The secondary structure of the VP2 reveals that the protein contains 220 alpha helix atoms, 40 310 helix, 151 beta sheets, 134 coils and 424 turns, whereas the 3D structure of Viral Protein 2 of BTV has been found to have 15774 total atoms in the structure. However, 961 amino acids were found in the final model. The dynamical cross-correlation matrix (DCCM) analysis tool identifies putative protein domains and also confirms the stability of the predicted model and their dynamical behavior difference with the correlative fluctuations in motion.

Conclusion:

The biological interpretation of the Viral Protein 2 was carried out. DCCM maps were calculated, using a different coordinate reference frame, through which, protein domain boundaries and protein domain residue constituents were identified. The obtained model shows good reliability. Moreover, we anticipated that this research should play a promising role in the identification of novel candidates with the target protein to inhibit their functional significance.

In Silico Insights on GD2 : A Potential Target for Pediatric Neuroblastoma

Background:

Originating from the abnormal growth of neuroblasts, pediatric neuroblastoma affects the age group below 15 years. It is an aggressive heterogenous cancer with a high morbidity rate. Biological marker GD2 synthesised by the GD2 gene acts as a powerful predictor of neuroblastoma cells. GD2 gangliosides are sialic acid-containing glycosphingolipids. Differential expression during brain development governs the function of the GD2. The present study explains the interaction of the GD2 with its established inhibitors and discovers the compound having a high binding affinity against the target protein. Technically, during the development of new compounds through docking studies, the best drug among all pre-exist inhibitors was filtered. Hence in reference to the best docked compound, the study proceeded further.

Methodology:

The In silico approach provides a platform to determine and establish potential inhibitor against GD2 in Pediatric neuroblastoma. The 3D structure of GD2 protein was modelled by homology base fold methods using Smith-Watermans’ Local alignment. A total of 18 established potent compounds were subjected to molecular docking and Etoposide (CID: 36462) manifested the highest affinity. The similarity search presented 336 compounds similar to Etoposide.

Results:

Through virtual screening, the compound having PubChem ID 10254934 showed a better affinity towards GD2 than the established inhibitor. The comparative profiling of the two compounds based on various interactions such as H-bond interaction, aromatic interactions, electrostatic interactions and ADMET profiling and toxicity studies were performed using various computational tools.

Conclusion:

The docking separated the virtual screened drug (PubChemID: 10254934) from the established inhibitor with a better re-rank score of -136.33. The toxicity profile of the virtual screened drug was also lesser (less lethal) than the established drug. The virtual screened drug was observed to be bioavailable as it does not cross the blood-brain barrier. Conclusively, the virtual screened compound obtained in the present investigation is better than the established inhibitor and can be further augmented by In vitro analysis, pharmacodynamics and pharmacokinetic studies.

Identification of Potent VEGF Inhibitors for the Clinical Treatment of Glioblastoma, A Virtual Screening Approach

Vascular endothelial growth factor (VEGF) expression could be found in all glioblastomas. VEGF takes part in numerous changes including the endothelial cell proliferation, the vasculature of solid tumor: its survival invasion, and migration, chemotaxis of bone marrow-derived progenitor cells, vasodilation and vascular permeability. VEGF inhibition can be a smart therapeutic strategy because it is extremely specific and less toxic than cytotoxic therapy. To establish better inhibition of VEGF than the current inhibitors, present study approach is by molecular docking, virtual screening to illustrate the inhibitor with superior affinity against VEGF to have a cautious pharma profile. To retrieve the best established and high-affinity high affinity molecule, Molegro Virtual Docker software was executed. The high-affinity scoring compounds were subjected to further similarity search to retrieve the drugs with similar properties from pubchem database. The completion of virtual screening reveals that PubChem compound SCHEMBL1250485 (PubChem CID: 66965667) has the highest affinity. The study of the drug-likeness was verified using OSIRIS Property Explorer software which supported the virtual screened result. Further ADMET study and drug comparative study strongly prove the superiority of the new established inhibitor with lesser rerank score and toxicity. Overall, the new inhibitor has higher potential to stop the expression of VEGF in glioblastoma and positively can be further analysed through In vitro studies.

Identification of High-affinity Small Molecules Targeting Gamma Secretase for the Treatment of Alzheimer’s Disease

Background:

Alzheimers Disease (AD) is a neurodegenerative disease which is characterized by the deposition of amyloid plaques in the brain- a concept supported by most of the researchers worldwide. The main component of the plaques being amyloid-beta (Aβ42) results from the sequential cleavage of Amyloid precursor protein (APP) by beta and gamma secretase. This present study intends to inhibit the formation of amyloid plaques by blocking the action of gamma secretase protein with Inhibitors (GSI).

Methods:

A number of Gamma Secretase Inhibitors (GSI) were targeted to the protein by molecular docking. The inhibitor having the best affinity was used as a subject for further virtual screening methods to obtain similar compounds. The generated compounds were docked again at the same docking site on the protein to find a compound with higher affinity to inhibit the protein. The highlights of virtually screened compound consisted of Pharmacophore Mapping of the docking site. These steps were followed by comparative assessments for both the compounds, obtained from the two aforesaid docking studies, which included interaction energy descriptors, ADMET profiling and PreADMET evaluations.

Results:

111 GSI classified as azepines, sulfonamides and peptide isosteres were used in the study. By molecular docking an amorpholino-amide, compound (22), was identified to be the high affinity compound GSI along with its better interaction profiles.The virtually screened pubchem compound AKOS001083915 (CID:24462213) shows the best affinity with gamma secretase. Collective Pharmacophore mapping (H bonds, electrostatic profile, binding pattern and solvent accesibility) shows a stable interaction. The resulting ADMETand Descriptor values were nearly equivalent.

Conclusion:

These compounds identified herein hold a potential as Gamma Secretase inhibitors.According to PreADMET values the compound AKOS001083915 is effective and specific to the target protein. Its BOILED-egg plot analysis infers the compound permeable to blood brain barrier.Comparative study for both the compounds resulted in having nearly equivalent properties. These compounds have the capacity to inhibit the protein which is indirectly responsible for the formation of amyloid plaques and can be further put to in vitro pharmacokinetic and dynamic studies.

A Computer - Aided Drug Designing for Pharmacological Inhibition of Mutant ALK for the Treatment of Non-small Cell Lung Cancer

Background:

Lung cancer is the most common among all the types of cancer worldwide with 1.8 million people diagnosed every year, leading to 1.6 million deaths every year according to the American cancer society. The involvement of mutated Anaplasic Lymphoma Kinase (ALK) positive fusion protein in the progression of NSCLC has made a propitious target to inhibit and treat NSCLC. In the present study, the main motif is to screen the most effective inhibitor against ALK protein with the potential pharmacological profile. The ligands selected were docked with Molegro Virtual Docker (MVD) and CEP-37440 (PubChem CID- 71721648) was the best docked pre-established compound with a permissible pharmacological profile.

Methods:

The selected ligands were docked with Molegro Virtual Docker (MVD). With reference to the obtained compound with the lowest re-rank score, PubChem database was virtually screened to retrieve a large set of similar compounds which were docked to find the compound with higher affinity. Further comparative studies and in silico prediction included pharmacophore studies, proximity energy parameters, ADMET and BOILED-egg plot analysis.

Results:

CEP-37440 (PubChem CID- 71721648) was the best docked pre-established compound with preferable pharmacological profile and PubChem compound CID-123449015 came out as the most efficient virtually screened inhibitor. Interestingly, the contours of the virtual screened compound PubChem CID- 123449015 fall within our desired high volume cavity of protein having admirable property to control the ALK regulation to prevent carcinogenesis in NSCLC. BOILED-Egg plot analysis depicts that both the compounds have analogous characteristics in the divergent aspects. Moreover, in the evaluations of Blood Brain Barrier, Human Intestinal Absorption, AMES toxicity, and LD50, the virtually screened compound (PubChem CID-123449015) was found within high optimization.

Conclusion:

These investigations denote that the virtually screened compound (PubChem CID- 123449015) is more efficient to be a better prospective candidate for NSCLC treatment having good pharmacological profile than the pre-established compound CEP-37440 (PubChem CID- 71721648) with low re-rank score. The identified virtually screened compound has high potential to act as an ALK inhibitor and can show promising results in the research of non-small cell lung cancer (NSCLC).

Pharmacogenetic Factors Affecting Asthma Treatment Response. Potential Implications for Drug Therapy

Asthma is a frequent disease, mainly characterized by airway inflammation, in which drug therapy is crucial in its management. The potential of pharmacogenomics testing in asthma therapy has been, to date, little explored. In this review, we discuss pharmacogenetic factors affecting asthma treatment, both related to drugs used as controller medications for regular maintenance, such as inhaled corticosteroids, anti-leukotriene agents, long-acting beta-agonists, and the new biologic agents used to treat severe persistent asthma. In addition, we discuss current pharmacogenomics knowledge for rescue medications provided to all patients for as-needed relief, such as short-acting beta-agonists. Evidence for genetic variations as a factor related to drugs response has been provided for the following genes and groups of drugs: Inhaled corticosteroids: FCER2; anti-leukotriene agents: ABCC1, and LTC4S; beta-agonists: ADRB2. However, the following genes require further studies confirming or rejecting association with the response to asthma therapy: ADCY9, ALOX5, ARG1, ARG2, CRHR1, CRHR2, CYP3A4, CYP3A5, CYSLTR1, CYSLTR2, GLCCI1, IL4RA, LTA4H, ORMDL3, SLCO2B1, SPATS2L, STIP1, T, TBX21, THRA, THRB, and VEGFA. Although only a minority of these genes are, at present, listed as associated with drugs used in asthma therapy, in the Clinical Pharmacogenomics Implementation Consortium gene-drug pair list, this review reveals that sufficient evidence to start testing the potential of clinical pharmacogenomics in asthma therapy already exists. This evidence supports the inclusion in pilot pharmacogenetics tests of at least four genes. Hopefully these tests, if proven useful, will increase the efficiency and the safety of asthma therapy.

Design of PD-L1 inhibitors for lung cancer

The progression of lung cancer is associated with inactivation of programmed cell death protein 1, abbreviated as PD- 1 which regulates the suppression of the body's immune system by suppressing T- cell inflammatory activity and is responsible for preventing cancer cell growth. It is of interest to identify inhibitors for PD-L1 dimeric structure through molecular docking and virtual screening. The virtual screened compound XGIQBUNWFCCMAS-UHFFFAOYSA-N (PubChem CID: 127263272) displays a high affinity with the target protein. ADMET analysis and cytotoxicity studies further add weight to this compound as a potential inhibitor of PD-L1. The established compound BMS-202 still shows the high re-rank score, but the virtual screened drug possesses a better ADMET profile with a higher intestinal absorption value and lower toxicity.

FLT3 inhibitor design using molecular docking based virtual screening for acute myeloid leukemia

Acute Myeloid Leukaemia (AML) is a blood cancer, which affects the red blood cells in the bone marrow. Of the possible proteins that are affected in AML, fms-like tyrosine kinase 3 (FLT3) has long been recognized as a potential therapeutic target as it affects the other signaling pathways and leads to a cascade of events. First-generation inhibitors sorafenib and midostaurin, as well as secondgeneration agents such as quizartinib and crenolanib are known. It is of interest to identify new compounds against FLT3 with improved activity using molecular docking and virtual screening. Molecular docking of existing inhibitors selected a top scoring bestestablished candidate Quizartinib having PubChem CID: 24889392. Similarity searching resulted in compound XGIQBUNWFCCMASUHFFFAOYSA-NPubChemCID: 44598530 which shows higher affinity scores. A comparative study of both the compounds using a drug-drug comparison, ADMET studies, boiled egg plot and pharmacophore parameters and properties confirmed the result and predicted the ligand to be an efficient inhibitor of FLT3.

Design of novel JAK3 Inhibitors towards Rheumatoid Arthritis using molecular docking analysis

Multiple cytokines play a pivotal role in the pathogenesis of Rheumatoid Arthritis by inducing intracellular signaling and it is known that the members of the Janus kinase (JAK) family are essential for such signal transduction. Janus kinase 3 is a tyrosine kinase that belongs to the Janus family of kinases. Drugs targeting JAK3 in the treatment of Rheumatoid arthritis is relevant. Therefore, it is of interest to design suitable inhibitors for JAK3 dimer using molecular docking with Molegro Virtual Docker. The compound possessing the highest affinity score is subjected to virtual screening to retrieve inhibitors. The compound SCHEMBL19100243 (PubChem CID- 76749591) displays a high affinity with the target protein. The affinity scores of this compound are more than known drugs. ADMET analysis and BOILED Egg plot provide insights into this compound as a potent inhibitor of JAK3.

Virtual Screening of IL-6 Inhibitors for Idiopathic Arthritis

Juvenile idiopathic arthritis (JIA) is a heterogeneous disease characterized by the arthritis of unknown origin and IL6 is a known target for JIA. 20 known inhibitors towards IL-6 were screened and Methotrexate (MTX) having PubChem ID: 126941 showed high binding capacity with the receptor IL-6. The similarity searching with this compound gave 269 virtual screened compounds. The said screening presented 269 possible drugs having structural similarity to Methotrexate. The docking studies of the screened drugs separated the compound having PubChem CID: 122677576 (re-rank value of -140.262). Toxicity and interaction profile validated this compound for having a better affinity with the target protein. Conclusively, this study shows that according to ADMET profile and BOILED-Egg plot, the compound (PubChem CID: 122677576) obtained from Virtual Screen could be the best drug in future during the prevention of juvenile idiopathic arthritis. In the current study, the drug CID: 122677576 is a potent candidate for treating JIA. The pharmacophore study revealed that the drug CID: 122677576 is a non-inhibitor of CYP450 microsomal enzymes and was found to be non-toxic, similar to the established drug Methotrexate (CID: 126941). It has a lower LD50 value of 2.6698mol/kg as compared to the established compound having LD50 value as 23.4955mol/kg. Moreover, the compound was found to be non-carcinogenic.

Molecular dynamic simulations reveal suboptimal binding of salbutamol in T164I variant of β2 adrenergic receptor

The natural variant C491T (rs1800088) in ADRB2 gene substitutes Threonine to Isoleucine at 164^th position in β2AR and results in receptor sequestration and altered binding of agonists. Present investigation pursues to identify the effect of T164I variation on function and structure of β2AR through systematic computational approaches. The study, in addition, addresses altered binding of salbutamol in T164I variant through molecular dynamic simulations. Methods involving changes in free energy, solvent accessibility surface area, root mean square deviations and analysis of binding cavity revealed structural perturbations in receptor to incur upon T164I substitution. For comprehensive understanding of receptor upon substitution, OPLS force field aided molecular dynamic simulations were performed for 10 ns. Simulations revealed massive structural departure for T164I β2AR variant from the native state along with considerably higher root mean square fluctuations of residues near the cavity. Affinity prediction by molecular docking showed two folds reduced affinity of salbutamol in T164I variant. To validate the credibility docking results, simulations for ligand-receptor complex were performed which demonstrated unstable salbutamol-T164I β2AR complex formation. Further, analysis of interactions in course of simulations revealed reduced ligand-receptor interactions of salbutamol in T164I variant. Taken together, studies herein provide structural rationales for suboptimal binding of salbutamol in T164I variant through integrated molecular modeling approaches.