Investigation of Inhibition Mechanism of Chemokine Receptor CCR5 by Micro-second Molecular Dynamics Simulations

Targeting CCR3 with antagonist SB 328437 sensitizes 5‑fluorouracil‑resistant gastric cancer cells: Experimental evidence and computational insights

Gastric cancer (GC) ranks fifth globally in cancer diagnoses and third for cancer-related deaths. Chemotherapy with 5-fluorouracil (5-FU), a primary treatment, faces challenges due to the development of chemoresistance. Tumor microenvironment factors, including C-C motif chemokine receptor 3 (CCR3), can contribute to chemoresistance. The present study evaluated the effect of CCR3 receptor inhibition using the antagonist SB 328437 and the molecular dynamics of this interaction on resistance to 5-FU in gastric cancer cells. The 5-FU-resistant AGS cell line (AGS R-5FU) demonstrated notable tolerance to higher concentrations of 5-FU, with a 2.6-fold increase compared with the parental AGS cell line. Furthermore, the mRNA expression levels of thymidylate synthase (TS), a molecular marker for 5-FU resistance, were significantly elevated in AGS R-5FU cells. CCR3 was shown to be expressed at significantly higher levels in these resistant cells. Combining SB 328437 with 5-FU resulted in a significant decrease in cell viability, particularly at higher concentrations of 5-FU. Furthermore, when SB 328437 was combined with 5-FU at a high concentration, the relative mRNA expression levels of CCR3 and TS decreased significantly. Computational analysis of CCR3 demonstrated dynamic conformational changes, especially in extracellular loop 2 region, which indicated potential alterations in ligand recognition. Docking simulations demonstrated that SB 328437 bound to the allosteric site of CCR3, inducing a conformational change in ECL2 and hindering ligand recognition. The present study provides comprehensive information on the molecular and structural aspects of 5-FU resistance and CCR3 modulation, highlighting the potential for therapeutic application of SB 328437 in GC treatment.

Insights of ligand binding in modeled h5-HT1A receptor: homology modeling, docking, MM-GBSA, screening and molecular dynamics

The pharmacologically characterized receptor subtype of the serotonin family, the 5HT1A receptor is implicated in the pathophysiology and treatment of depression and anxiety-related disorders. Being the most extensively targeted receptor for developing novel antidepressants and anxiolytics, a near-ideal theoretical model can aid in high-throughput screening of promising drug candidates. However, the design of potential drug candidates suffers owing to a lack of complete structural information. In this work, homology models of 5-HT1A receptor are generated using two distinct alignments (CW and PSTA) and model building methods (KB and EB). The developed models are validated for virtual screening using a ligand dataset of agonists and antagonists. The best-suited model was efficient in discriminating agonist/antagonist binding. Correlation plots between pKi and docking (R²_agonist≥ 0.6, R²_antagonist≥ 0.7) and MM-GBSA dG bind values (R²_agonist≥ 0.5, R²_antagonist≥ 0.7) revealed optimum corroboration between in vitro and in silico outcomes, which further suggested the usefulness of the developed model for the design of high-affinity probes for the neurological disorders.Communicated by Ramaswamy H. Sarma.

Pharmacophore-Oriented Identification of Potential Leads as CCR5 Inhibitors to Block HIV Cellular Entry

Cysteine-cysteine chemokine receptor 5 (CCR5) has been discovered as a co-receptor for cellular entry of human immunodeficiency virus (HIV). Moreover, the role of CCR5 in a variety of cancers and various inflammatory responses was also discovered. Despite the fact that several CCR5 antagonists have been investigated in clinical trials, only Maraviroc has been licensed for use in the treatment of HIV patients. This indicates that there is a need for novel CCR5 antagonists. Keeping this in mind, the present study was designed. The active CCR5 inhibitors with known IC₅₀ value were selected from the literature and utilized to develop a ligand-based common feature pharmacophore model. The validated pharmacophore model was further used for virtual screening of drug-like databases obtained from the Asinex, Specs, InterBioScreen, and Eximed chemical libraries. Utilizing computational methods such as molecular docking studies, molecular dynamics simulations, and binding free energy calculation, the binding mechanism of selected inhibitors was established. The identified Hits not only showed better binding energy when compared to Maraviroc, but also formed stable interactions with the key residues and showed stable behavior throughout the 100 ns MD simulation. Our findings suggest that Hit1 and Hit2 may be potential candidates for CCR5 inhibition, and, therefore, can be considered for further CCR5 inhibition programs.

Computational study of the structural ensemble of CC chemokine receptor type 5 (CCR5) and its interactions with different ligands

CC Chemokine receptor 5 (CCR5), a member of the Superfamily of G Protein-Coupled Receptors (GPCRs), is an important effector in multiple physiopathological processes such as inflammatory and infectious entities, including central nervous system neuroinflammatory diseases such as Alzheimer's disease, recovery from nervous injuries, and in the HIV-AIDS infective processes. Thus, CCR5 is an attractive target for pharmacological modulation. Since maraviroc was described as a CCR5 ligand that modifies the HIV-AIDS progression, multiple efforts have been developed to describe the functionality of the receptor. In this work, we characterized key structural features of the CCR5 receptor employing extensive atomistic molecular dynamics (MD) in its apo form and in complex with an endogenous agonist, the chemokine CCL5/RANTES, an HIV entry inhibitor, the partial inverse agonist maraviroc, and the experimental antagonists Compound 21 and 34, aiming to elucidate the structural features and mechanistic processes that constitute its functional states, contributing with structural details and a general understanding of this relevant system.

Identification of 1H-purine-2,6-dione derivative as a potential SARS-CoV-2 main protease inhibitor: molecular docking, dynamic simulations, and energy calculations

The rapid spread of the coronavirus since its first appearance in 2019 has taken the world by surprise, challenging the global economy, and putting pressure on healthcare systems across the world. The introduction of preventive vaccines only managed to slow the rising death rates worldwide, illuminating the pressing need for developing effective antiviral therapeutics. The traditional route of drug discovery has been known to require years which the world does not currently have. In silico approaches in drug design have shown promising results over the last decade, helping to decrease the required time for drug development. One of the vital non-structural proteins that are essential to viral replication and transcription is the SARS-CoV-2 main protease (Mpro). Herein, using a test set of recently identified COVID-19 inhibitors, a pharmacophore was developed to screen 20 million drug-like compounds obtained from a freely accessible Zinc database. The generated hits were ranked using a structure based virtual screening technique (SBVS), and the top hits were subjected to in-depth molecular docking studies and MM-GBSA calculations over SARS-COV-2 Mpro. Finally, the most promising hit, compound (1), and the potent standard (III) were subjected to 100 ns molecular dynamics (MD) simulations and in silico ADME study. The result of the MD analysis as well as the in silico pharmacokinetic study reveal compound 1 to be a promising SARS-Cov-2 MPro inhibitor suitable for further development.

Dynamic insights into the effects of nonsynonymous polymorphisms (nsSNPs) on loss of TREM2 function

Single nucleotide variations in Triggering Receptor Expressed on Myeloid Cells 2 (TREM2) are associated with many neurodegenerative diseases, including Nasu-Hakola disease (NHD), frontotemporal dementia (FTD), and late-onset Alzheimer's disease because they disrupt ligand binding to the extracellular domain of TREM2. However, the effects of nonsynonymous single nucleotide polymorphisms (nsSNPs) in TREM2 on disease progression remain unknown. In this study, we identified several high-risk nsSNPs in the TREM2 gene using various deleterious SNP predicting algorithms and analyzed their destabilizing effects on the ligand recognizing region of the TREM2 immunoglobulin (Ig) domain by molecular dynamics (MD) simulation. Cumulative prediction by all tools employed suggested the three most deleterious nsSNPs involved in loss of TREM2 function are rs549402254 (W50S), rs749358844 (R52C), and rs1409131974 (D104G). MD simulation showed that these three variants cause substantial structural alterations and conformational remodeling of the apical loops of the TREM2 Ig domain, which is responsible for ligand recognition. Detailed analysis revealed that these variants substantially increased distances between apical loops and induced conformation remodeling by changing inter-loop nonbonded contacts. Moreover, all nsSNPs changed the electrostatic potentials near the putative ligand-interacting region (PLIR), which suggested they might reduce specificity or loss of binding affinity for TREM2 ligands. Overall, this study identifies three potential high-risk nsSNPs in the TREM2 gene. We propose further studies on the molecular mechanisms responsible for loss of TREM2 function and the associations between TREM2 nsSNPs and neurodegenerative diseases.

Structural Consequence of Non-Synonymous Single-Nucleotide Variants in the N-Terminal Domain of LIS1

Disruptive neuronal migration during early brain development causes severe brain malformation. Characterized by mislocalization of cortical neurons, this condition is a result of the loss of function of migration regulating genes. One known neuronal migration disorder is lissencephaly (LIS), which is caused by deletions or mutations of the LIS1 (PAFAH1B1) gene that has been implicated in regulating the microtubule motor protein cytoplasmic dynein. Although this class of diseases has recently received considerable attention, the roles of non-synonymous polymorphisms (nsSNPs) in LIS1 on lissencephaly progression remain elusive. Therefore, the present study employed combined bioinformatics and molecular modeling approach to identify potential damaging nsSNPs in the LIS1 gene and provide atomic insight into their roles in LIS1 loss of function. Using this approach, we identified three high-risk nsSNPs, including rs121434486 (F31S), rs587784254 (W55R), and rs757993270 (W55L) in the LIS1 gene, which are located on the N-terminal domain of LIS1. Molecular dynamics simulation highlighted that all variants decreased helical conformation, increased the intermonomeric distance, and thus disrupted intermonomeric contacts in the LIS1 dimer. Furthermore, the presence of variants also caused a loss of positive electrostatic potential and reduced dimer binding potential. Since self-dimerization is an essential aspect of LIS1 to recruit interacting partners, thus these variants are associated with the loss of LIS1 functions. As a corollary, these findings may further provide critical insights on the roles of LIS1 variants in brain malformation.

Computational strategies to identify new drug candidates against neuroinflammation

The even more increasing application of computational approaches in these last decades has deeply modified the process of discovery and commercialization of new therapeutic entities. This is especially true in the field of neuroinflammation, in which both the peculiar anatomical localization and the presence of the blood-brain barrier makeit mandatory to finely tune the candidates' physicochemical properties from the early stages of the discovery pipeline. The aim of this review is therefore to provide a general overview to the readers about the topic of neuroinflammation, together with the most common computational strategies that can be exploited to discover and design small molecules controlling neuroinflammation, especially those based on the knowledge of the three-dimensional structure of the biological targets of therapeutic interest. The techniques used to describe the molecular recognition mechanisms, such as molecular docking and molecular dynamics, will therefore be eviscerated, highlighting their advantages and their limitations. Finally, we report several case studies in which computational methods have been applied in drug discovery on neuroinflammation, focusing on the last decade's research.

A two-stage computational approach to predict novel ligands for a chemosensory receptor

Olfactory receptor (OR) 1A2 is the member of largest superfamily of G protein-coupled receptors (GPCRs). OR1A2 is an ectopically expressed receptor with only 13 known ligands, implicated in reducing hepatocellular carcinoma progression, with enormous therapeutic potential. We have developed a two-stage screening approach to identify novel putative ligands of OR1A2. We first used a pharmacophore model based on atomic property field (APF) to virtually screen a library of 5942 human metabolites. We then carried out structure-based virtual screening (SBVS) for predicting the potential agonists, based on a 3D homology model of OR1A2. This model was developed using a biophysical approach for template selection, based on multiple parameters including hydrophobicity correspondence, applied to the complete set of available GPCR structures to pick the most appropriate template. Finally, the membrane-embedded 3D model was refined by molecular dynamics (MD) simulations in both the apo and holo forms. The refined model in the apo form was selected for SBVS. Four novel small molecules were identified as strong binders to this olfactory receptor on the basis of computed binding energies.

Assessment of the PETase conformational changes induced by poly(ethylene terephthalate) binding

Recently, a bacterium strain of Ideonella sakaiensis was identified with the uncommon ability to degrade the poly(ethylene terephthalate) (PET). The PETase from I. sakaiensis strain 201-F6 (IsPETase) catalyzes the hydrolysis of PET converting it to mono(2-hydroxyethyl) terephthalic acid (MHET), bis(2-hydroxyethyl)-TPA (BHET), and terephthalic acid (TPA). Despite the potential of this enzyme for mitigation or elimination of environmental contaminants, one of the limitations of the use of IsPETase for PET degradation is the fact that it acts only at moderate temperature due to its low thermal stability. Besides, molecular details of the main interactions of PET in the active site of IsPETase remain unclear. Herein, molecular docking and molecular dynamics (MD) simulations were applied to analyze structural changes of IsPETase induced by PET binding. Results from the essential dynamics revealed that the β1-β2 connecting loop is very flexible. This loop is located far from the active site of IsPETase and we suggest that it can be considered for mutagenesis to increase the thermal stability of IsPETase. The free energy landscape (FEL) demonstrates that the main change in the transition between the unbound to the bound state is associated with the β7-α5 connecting loop, where the catalytic residue Asp206 is located. Overall, the present study provides insights into the molecular binding mechanism of PET into the IsPETase structure and a computational strategy for mapping flexible regions of this enzyme, which can be useful for the engineering of more efficient enzymes for recycling plastic polymers using biological systems.

Structural and Functional Characterization of Allatostatin Receptor Type-C of Thaumetopoea pityocampa, a Potential Target for Next-Generation Pest Control Agents

Insect neuropeptide receptors, including allatostatin receptor type C (AstR-C), a G protein-coupled receptor, are among the potential targets for designing next-generation pesticides that despite their importance in offering a new mode-of-action have been overlooked. Focusing on AstR-C of Thaumetopoea pityocampa, a common pest in Mediterranean countries, by employing resonance energy transfer-based methods, we showed Gαi/o coupling and β-arrestin recruitment of the receptor at sub-nanomolar and nanomolar ranges of the endogenous ligand, AST-C, respectively. Molecular docking and molecular dynamics simulation studies revealed the importance of extracellular loop 2 in AstRC/AST-C interaction, and a combination of in silico and in vitro approaches showed the substantial role of Q271^6.55 in G protein-dependent activation of AstR-C possibly via contributing to the flexibility of the receptor's structure. The functional and structural insights obtained on T. pit AstR-C positively assist future efforts in developing environmentally friendly pest control agents that are needed urgently.

Mechanistic insights into the deleterious roles of Nasu-Hakola disease associated TREM2 variants

Recently, the critical roles played by genetic variants of TREM2 (Triggering Receptor Expressed on Myeloid cells 2) in Alzheimer's disease have been aggressively highlighted. However, few studies have focused on the deleterious roles of Nasu-Hakola disease (NHD) associated TREM2 variants. In order to get insights into the contributions made by these variants to neurodegeneration, we investigated the influences of four NHD associated TREM2 mutations (Y38C, W50C, T66M, and V126G) on loss-of-function, and followed this with in silico prediction and conventional molecular dynamics simulation. NHD mutations were predicted to be highly deleterious by eight different in silico bioinformatics tools and found to induce conformational changes by molecular dynamics simulation. As compared with the wild-type, the four variants produced substantial differences in the collective motions of loop regions, which not only promoted structural remodeling in the CDR2 (complementarity-determining region 2) loop but also in the CDR1 loop, by changing inter- and intra-loop hydrogen bonding networks. In addition, structural studies in a free energy landscape analysis showed that Y38, T66, and V126 are crucial for maintaining the structural features of CDR1 and CDR2 loops, and that mutations in these positions produced steric clashes and loss of ligand binding. These results showed the presence of mutations in the TREM2 ectodomain induced flexibility and caused structural alterations. Dynamical scenarios, as provided by the present study, may be critical to our understanding of the roles of these TREM2 mutations in neurodegenerative diseases.

Unveiling the Structural Insights into the Selective Inhibition of Protein Kinase D1

BACKGROUND

Although protein kinase D1 (PKD1) has been proved to be an efficient target for anticancer drug development, lack of structural details and substrate binding mechanisms are the main obstacles for the development of selective inhibitors with therapeutic benefits.

OBJECTIVE

The present study described the in silico dynamics behaviors of PKD1 in binding with selective and non-selective inhibitors and revealed the critical binding site residues for the selective kinase inhibition.

METHODS

Here, the three dimensional model of PKD1 was initially constructed by homology modeling along with binding site characterization to explore the non-conserved residues. Subsequently, two known inhibitors were docked to the catalytic site and the detailed ligand binding mechanisms and post binding dyanmics were investigated by molecular dynamics simulation and binding free energy calculations.

RESULTS

According to the binding site analysis, PKD1 serves several non-conserved residues in the G-loop, hinge and catalytic subunits. Among them, the residues including Leu662, His663, and Asp665 from hinge region made polar interactions with selective PKD1 inhibitor in docking simulation, which were further validated by the molecular dynamics simulation. Both inhibitors strongly influenced the structural dynamics of PKD1 and their computed binding free energies were in accordance with experimental bioactivity data.

CONCLUSION

The identified non-conserved residues likely to play critical role on molecular reorganization and inhibitor selectivity. Taken together, this study explained the molecular basis of PKD1 specific inhibition, which may help to design new selective inhibitors for better therapies to overcome cancer and PKD1 dysregulated disorders.

Structural and Dynamic Characterizations Highlight the Deleterious Role of SULT1A1 R213H Polymorphism in Substrate Binding

Sulfotransferase 1A1 (SULT1A1) is responsible for catalyzing various types of endogenous and exogenous compounds. Accumulating data indicates that the polymorphism rs9282861 (R213H) is responsible for inefficient enzymatic activity and associated with cancer progression. To characterize the detailed functional consequences of this mutation behind the loss-of-function of SULT1A1, the present study deployed molecular dynamics simulation to get insights into changes in the conformation and binding energy. The dynamics scenario of SULT1A1 in both wild and mutated types as well as with and without ligand showed that R213H induced local conformational changes, especially in the substrate-binding loop rather than impairing overall stability of the protein structure. The higher conformational changes were observed in the loop3 (residues, 235-263), turning loop conformation to A-helix and B-bridge, which ultimately disrupted the plasticity of the active site. This alteration reduced the binding site volume and hydrophobicity to decrease the binding affinity of the enzyme to substrates, which was highlighted by the MM-PBSA binding energy analysis. These findings highlight the key insights of structural consequences caused by R213H mutation, which would enrich the understanding regarding the role of SULT1A1 mutation in cancer development and also xenobiotics management to individuals in the different treatment stages.

Membrane partitioning of peptide aggregates: coarse-grained molecular dynamics simulations

Coarse-grained molecular dynamics (CGMD) simulation technique (MARTINI force field) is applied to monitor the aggregation of helical peptides representing the transmembrane sequence and its extension of bone marrow stromal cell antigen 2 (BST-2). One of the peptides is coupled with a protein transducing domain (PTD) of nine arginine residues (R9) at its N-terminal side as well as a peptide, pep11**, which has been shown to bind to human papilloma virus 16 (HPV16) E6 oncoprotein. A short hydrophobic stretch of the transmembrane domain (TMD) of BST-2 aggregates the fastest and inserts into a lipid membrane. An aggregate of R9-pep11** attaches to the membrane via simultaneous contact of many arginine residues. Monomers from the aggregates of the shortest of the hydrophobic TMDs dissolve into the opposing leaflet when the aggregate spans the bilayer. A 'flipping' of the individual monomeric peptides is not observed.Communicated by Ramaswamy H. Sarma.

Identification and structural characterization of deleterious non-synonymous single nucleotide polymorphisms in the human SKP2 gene

In SCF (Skp, Cullin, F-box) ubiquitin-protein ligase complexes, S-phase kinase 2 (SKP2) is one of the major players of F-box family, that is responsible for the degradation of several important cell regulators and tumor suppressor proteins. Despite of having significant evidence for the role of SKP2 on tumorgenesis, there is a lack of available data regarding the effect of non-synonymous polymorphisms. In this communication, the structural and functional consequences of non-synonymous single nucleotide polymorphisms (nsSNPs) of SKP2 have been reported by employing various computational approaches and molecular dynamics simulation. Initially, several computational tools like SIFT, PolyPhen-2, PredictSNP, I-Mutant 2.0 and ConSurf have been implicated in this study to explore the damaging SNPs. In total of 172 nsSNPs, 5 nsSNPs were identified as deleterious and 3 of them were predicted to be decreased the stability of protein. Guided from ConSurf analysis, P101L (rs761253702) and Y346C (rs755010517) were categorized as the highly conserved and functional disrupting mutations. Therefore, these mutations were subjected to three dimensional model building and molecular dynamics simulation study for the detailed structural consequences upon the mutations. The study revealed that P101L and Y346C mutations increased the flexibility and changed the structural dynamics. As both these mutations are located in the most functional regions of SKP2 protein, these computational insights might be helpful to consider these nsSNPs for wet-lab confirmatory analysis as well as in rationalizing future population based studies and structure based drug design against SKP2.

CCR5 Revisited: How Mechanisms of HIV Entry Govern AIDS Pathogenesis

The chemokine receptor CCR5 has been the focus of intensive studies since its role as a coreceptor for HIV entry was discovered in 1996. These studies lead to the development of small molecular drugs targeting CCR5, with maraviroc becoming in 2007 the first clinically approved chemokine receptor inhibitor. More recently, the apparent HIV cure in a patient transplanted with hematopoietic stem cells devoid of functional CCR5 rekindled the interest for inactivating CCR5 through gene therapy and pharmacological approaches. Fundamental research on CCR5 has also been boosted by key advances in the field of G-protein coupled receptor research, with the realization that CCR5 adopts a variety of conformations, and that only a subset of these conformations may be targeted by chemokine ligands. In addition, recent genetic and pathogenesis studies have emphasized the central role of CCR5 expression levels in determining the risk of HIV and SIV acquisition and disease progression. In this article, we propose to review the key properties of CCR5 that account for its central role in HIV pathogenesis, with a focus on mechanisms that regulate CCR5 expression, conformation, and interaction with HIV envelope glycoproteins.

BIM Binding Remotely Regulates BAX Activation: Insights from the Free Energy Landscapes

Activation of the pro-apoptotic BAX protein, a BCL-2 family member, is known to trigger apoptosis by forming pores in the mitochondrial outer membrane (MOM). While in the cytosol, release of its transmembrane C-terminal helix (called α9 helix) from a well-characterized binding pocket (BC groove) and subsequent permeabilization of the MOM are understood to be the initiating events of the activation. Concerning what initiates BAX activation, so far one plausible suggestion has been that the transient attachment of BH3-only peptide at a distal site from the BC groove triggers the activation process. Yet how this pivotal step displaces α9 from the BC groove has remained poorly understood. Using a combination of standard molecular dynamics and enhanced sampling methods, the energy landscape of BIM (BH3-only peptide) induced BAX activation has been computed, and the molecular origin of those events is hereby reported in atomistic detail. The simulated transition pathway of α9 release reveals that BIM subdues the energetic cost of the process by reducing the activation energy barrier to some extent but mostly by minimizing the free energy difference between the active (α9-released) and inactive (α9-bound) states. Interestingly, the flexibility of the α9 helix itself plays a decisive role in this mechanism. The impact of BIM encounter at the distal site is found to propagate to the α9 (BC groove bound) mostly through conserved pathways of residue level interactions. Overall, the thermodynamic basis of the "hit-and-run" mechanism for activation of the BCL-2 family is presented reconciling the available biochemical observations.

Using accelerated molecular dynamics simulation to shed light on the mechanism of activation/deactivation upon mutations for CCR5

In this work, accelerated molecular dynamics (aMD) simulations were used to study different effects of G286F and R126 mutations on the activity of CCR5. Potential of Mean Force (PMF) results indicate that there are stable inactive-like states and active-like ones existing in the conformation space of the wild type (WT), confirming that CCR5 could possess to some extent constitutive activity. But the R126N mutation could constrain CCR5 in the inactive state through influencing the TXP motif and limiting the movements of TM5 and TM6. In contrast, the G286F mutation promotes the activity of the receptor by increasing the distance of TM2–TM6 and the flexibility of the intracellular part of TM5 and changing the H-bonding in the TXP motif. The observations from the cross correlation analysis further show that the R126N mutation dramatically reduces the motion correlations between TMs, which should partly contribute to the deactivation of CCR5. Compared with the WT system, TM6 and TM7 in the G286F mutant are loosely correlated with other regions, which should be conducive to drive the movement of TM6 and TM7 toward the active conformation. In addition, the result from the protein structure network (PSN) analysis reveals that the shortest pathways connecting the extracellular and the intracellular domains are highly conserved in the three systems despite the different mutations, in which the hydrogen bond plays a pivotal role. However, the G286F mutation shortens the lifetime of the pathway with respect to the R126N mutation, which may be associated with the different activities of the two mutants. The pathway connecting the ligand-binding site and the G-protein region reveals that the allosteric communication between TM6 and TM7 is enhanced by the R126N mutation while the G286F mutation induces the activation of the G-protein pocket by arousing more residues in the NPxxY region to participate in the pathway.

In this work, accelerated molecular dynamics (aMD) simulations were used to study different effects of G286F and R126 mutations on the activity of CCR5.

Molecular mechanisms of allosteric probe dependence in μ opioid receptor

Allostery is one of the most important features of proteins. It greatly contributes to the complexity of life, since it enables possibility of precise tuning of protein function, as well as performing more than one function per protein. Probe dependence is one of the unique features of allostery. It allows a protein to respond differently to the same allosteric modulator when different drugs or transmitters are bound. Unfortunately, allosteric mechanisms are difficult to investigate experimentally. Instead, they can be reproduced artificially in simulations. We simulated in silico a native-like cell membrane fragment with an active-state human μ opioid receptor (MOR) in order to investigate diverse effects of a receptor's positive allosteric modulator on various agonists. Particular emphasis on native-likeness of the environment was put. We managed to reproduce the experimentally observed effects, which allowed us to take deeper insight into their underlying mechanisms. We found an allosteric pathway in the receptor, leading from the ligand binding site to the intracellular, effector site. We observed that the modulator affected the pathway, inducing different resultant responses for full and partial agonists.

Improvising 5-HT7R homology model for design of high affinity ligands: model validation with docking, embrace minimization, MM-GBSA, and molecular dynamic simulations

The subtype, 5-HT₇R has been implicated in neurological disorders and presents itself as a promising target for antidepressant drugs. Design of targeted selective ligands, require a sound knowledge of 3D-receptor structure. In absence of receptor structure, structure-based design of targeted ligands relies on generation of 5-HT₇R homology model. In this study, the impact of template choice, alignment, and model building methods on the homology model of 5-HT₇R is addressed. The compactness and model quality due to the presence of cholesterol (lipidic receptor) have also been observed. The results suggest good stereochemical quality of the final model. Ramachandran Plot Analysis indicated more than 97.5% amino acid residues in the favorable region. The overall quality factor was 91.8% using ERRAT. The Z-score for backbone confirmation and packing quality were -1.248 and -1.427, respectively, using WHATCHECK. The RMS Z-score for side chain planarity was .711. Other validation results for the final model include binding site analysis in which Asp162, Val163, Phe343, Phe344, Arg350, Arg367, and Leu370 conserved residues were found in the active site, correlation coefficient of .82 in ligand-based screening and .85 in embrace minimization. Further, the model showed good correlation for agonist and antagonist in docking ([Formula: see text] ≈ .76, [Formula: see text] ≈ .82), embrace minimization ([Formula: see text] ≈ .73, [Formula: see text] ≈ .72), and MM-GBSA ([Formula: see text] ≈ .69, [Formula: see text] ≈ .75) studies. The model was subjected to Molecular Dynamics (MD) simulation of 20 ns both in ligand-free and ligand-bound receptor (agonist and antagonist) system in order to assess its stability.

Exploration of binding and inhibition mechanism of a small molecule inhibitor of influenza virus H1N1 hemagglutinin by molecular dynamics simulation

Influenza viruses are a major public health threat worldwide. The influenza hemagglutinin (HA) plays an essential role in the virus life cycle. Due to the high conservation of the HA stem region, it has become an especially attractive target for inhibitors for therapeutics. In this study, molecular simulation was applied to study the mechanism of a small molecule inhibitor (MBX2329) of influenza HA. Behaviors of the small molecule under neutral and acidic conditions were investigated, and an interesting dynamic binding mechanism was found. The results suggested that the binding of the inhibitor with HA under neutral conditions facilitates only its intake, while it interacts with HA under acidic conditions using a different mechanism at a new binding site. After a series of experiments, we believe that binding of the inhibitor can prevent the release of HA1 from HA2, further maintaining the rigidity of the HA2 loop and stabilizing the distance between the long helix and short helices. The investigated residues in the new binding site show high conservation, implying that the new binding pocket has the potential to be an effective drug target. The results of this study will provide a theoretical basis for the mechanism of new influenza virus inhibitors.

Understanding the molecular mechanism for the differential inhibitory activities of compounds against MTH1

MTH1 can hydrolyze oxidized nucleotides and is required for cancer survival. The IC₅₀ values were 0.8 nM for TH287 with a methyl substitution, 5.0 nM for TH588 with a cyclopropyl substitution, and 2.1 μM for TH650 with an oxetanyl substitution. Thus, it is very significant to understand inhibitory mechanisms of these structurally similar compounds against MTH1 and influences of the substituent on the bioactivities. Our MD researches indicate that TH287 maintains significant hydrogen bonds with Asn33 and Asp119, stabilizes the binding site, and induces MTH1 adopt a closed motion, leading to a high inhibitory activity. When bound with TH588, the binding site can be partially stabilized and take a semi-closed state, which is because the cyclopropyl group in TH588 has larger steric hindrance than a methyl group in TH287. So TH588 has a slightly reduced inhibitory activity compared to TH287. TH650 induces greater conformation fluctuations than TH588 and the binding site adopts an opening state, which is caused by the large bulk of oxetanyl group and the interference of solvent on the oxetanyl substituent, leading to the lowest inhibitory activity. Thus, the inhibitory activity follows a TH287 > TH588 > TH650 trend, which well matches with the experimental finding.

Virtual screening of eighteen million compounds against dengue virus: Combined molecular docking and molecular dynamics simulations study

Dengue virus is a major issue of tropical and sub-tropical regions. Dengue virus has been the cause behind the major alarming epidemics in the history with mass causalities from the decades. Unavailability of on-shelf drugs for the prevention of further proliferation of virus inside the human body results in immense number of deaths each year. This issue necessitates the design of novel anti-dengue drug. The protease enzyme pathway is the critical target for drug design due to its significance in the replication, survival and other cellular activities of dengue virus. Therefore, approximately eighteen million compounds from the ZINC database have been virtually screened against nonstructural protein 3 (NS3). The incremental construction algorithm of Glide docking program has been used with its features high throughput virtual screening (HTVS), standard precision (SP), extra precision (XP) and in combination of Prime module, induced fit docking (IFD) approach has also been applied. Five top-ranked compounds were then selected from the IFD results with better predicted binding energies with the catalytic triad residues (His51, Asp75, and Ser135) that may act as potential inhibitors for the underlying target protease enzyme. The top-ranked compounds ZINC95518765, ZINC44921800, ZINC71917414, ZINC39500661, ZINC36681949 have shown the predicted binding energies of -7.55, -7.36, -8.04, -8.41, -9.18kcal/mol, respectively, forming binding interactions with three catalytically important amino acids. Top-docking poses of compounds are then used in molecular dynamics (MD) simulations. In computational studies, our proposed compounds confirm promising results against all the four serotypes of dengue virus, strengthening the opportunity of these compounds to work as potential on-shelf drugs against dengue virus. Further experimentation on the proposed compounds can result in development of strong inhibitors.

The Structural Basis of Oncogenic Mutations G12, G13 and Q61 in Small GTPase K-Ras4B

Ras mediates cell proliferation, survival and differentiation. Mutations in K-Ras4B are predominant at residues G12, G13 and Q61. Even though all impair GAP-assisted GTP → GDP hydrolysis, the mutation frequencies of K-Ras4B in human cancers vary. Here we aim to figure out their mechanisms and differential oncogenicity. In total, we performed 6.4 μs molecular dynamics simulations on the wild-type K-Ras4B (K-Ras4B^WT-GTP/GDP) catalytic domain, the K-Ras4B^WT-GTP–GAP complex, and the mutants (K-Ras4B^{G12C/G12D/G12V}-GTP/GDP, K-Ras4B^G13D-GTP/GDP, K-Ras4B^Q61H-GTP/GDP) and their complexes with GAP. In addition, we simulated ‘exchanged’ nucleotide states. These comprehensive simulations reveal that in solution K-Ras4B^WT-GTP exists in two, active and inactive, conformations. Oncogenic mutations differentially elicit an inactive-to-active conformational transition in K-Ras4B-GTP; in K-Ras4B^G12C/G12D-GDP they expose the bound nucleotide which facilitates the GDP-to-GTP exchange. These mechanisms may help elucidate the differential mutational statistics in K-Ras4B-driven cancers. Exchanged nucleotide simulations reveal that the conformational transition is more accessible in the GTP-to-GDP than in the GDP-to-GTP exchange. Importantly, GAP not only donates its R789 arginine finger, but stabilizes the catalytically-competent conformation and pre-organizes catalytic residue Q61; mutations disturb the R789/Q61 organization, impairing GAP-mediated GTP hydrolysis. Together, our simulations help provide a mechanistic explanation of key mutational events in one of the most oncogenic proteins in cancer.

Investigation of allosteric modulation mechanism of metabotropic glutamate receptor 1 by molecular dynamics simulations, free energy and weak interaction analysis

Metabotropic glutamate receptor 1 (mGlu₁), which belongs to class C G protein-coupled receptors (GPCRs), can be coupled with G protein to transfer extracellular signal by dimerization and allosteric regulation. Unraveling the dimer packing and allosteric mechanism can be of great help for understanding specific regulatory mechanism and designing more potential negative allosteric modulator (NAM). Here, we report molecular dynamics simulation studies of the modulation mechanism of FITM on the wild type, T815M and Y805A mutants of mGlu₁ through weak interaction analysis and free energy calculation. The weak interaction analysis demonstrates that van der Waals (vdW) and hydrogen bonding play an important role on the dimer packing between six cholesterol molecules and mGlu₁ as well as the interaction between allosteric sites T815, Y805 and FITM in wild type, T815M and Y805A mutants of mGlu₁. Besides, the results of free energy calculations indicate that secondary binding pocket is mainly formed by the residues Thr748, Cys746, Lys811 and Ser735 except for FITM-bound pocket in crystal structure. Our results can not only reveal the dimer packing and allosteric regulation mechanism, but also can supply useful information for the design of potential NAM of mGlu₁.

Binding Interactions of Dopamine and Apomorphine in D2High and D2Low States of Human Dopamine D2 Receptor Using Computational and Experimental Techniques

We have recently reported G-protein coupled receptor (GPCR) model structures for the active and inactive states of the human dopamine D2 receptor (D2R) using adrenergic crystal structures as templates. Since the therapeutic concentrations of dopamine agonists that suppress the release of prolactin are the same as those that act at the high-affinity state of the D2 receptor (D2High), D2High in the anterior pituitary gland is considered to be the functional state of the receptor. In addition, the therapeutic concentrations of anti-Parkinson drugs are also related to the dissociation constants in the D2High form of the receptor. The discrimination between the high- and low-affinity (D2Low) components of the D2R is not obvious and requires advanced computer-assisted structural biology investigations. Therefore, in this work, the derived D2High and D2Low receptor models (GPCR monomer and dimer three-dimensional structures) are used as drug-binding targets to investigate binding interactions of dopamine and apomorphine. The study reveals a match between the experimental dissociation constants of dopamine and apomorphine at their high- and low-affinity sites of the D2 receptor in monomer and dimer and their calculated dissociation constants. The allosteric receptor-receptor interaction for dopamine D2R dimer is associated with the accessibility of adjacent residues of transmembrane region 4. The measured negative cooperativity between agonist ligand at dopamine D2 receptor is also correctly predicted using the D2R homodimerization model.

Mutated form (G52E) of inactive diphtheria toxin CRM197: molecular simulations clearly display effect of the mutation to NAD binding

Mutated form (G52E) of diphtheria toxin (DT) CRM197 is an inactive and nontoxic enzyme. Here, we provided a molecular insight using comparative molecular dynamics (MD) simulations to clarify the influence of a single point mutation on overall protein and active-site loop. Post-processing MD analysis (i.e. stability, principal component analysis, hydrogen-bond occupancy, etc.) is carried out on both wild and mutated targets to investigate and to better understand the mechanistic differences of structural and dynamical properties on an atomic scale especially at nicotinamide adenine dinucleotide (NAD) binding site when a single mutation (G52E) happens at the DT. In addition, a docking simulation is performed for wild and mutated forms. The docking scoring analysis and docking poses results revealed that mutant form is not able to properly accommodate the NAD molecule.