Accelerated molecular dynamics simulations of the octopamine receptor using GPUs: discovery of an alternate agonist‐binding position

Octopamine receptors (OARs) perform key biological functions in invertebrates, making this class of G‐protein coupled receptors (GPCRs) worth considering for insecticide development. However, no crystal structures and very little research exists for OARs. Furthermore, GPCRs are large proteins, are suspended in a lipid bilayer, and are activated on the millisecond timescale, all of which make conventional molecular dynamics (MD) simulations infeasible, even if run on large supercomputers. However, accelerated Molecular Dynamics (aMD) simulations can reduce this timescale to even hundreds of nanoseconds, while running the simulations on graphics processing units (GPUs) would enable even small clusters of GPUs to have processing power equivalent to hundreds of CPUs. Our results show that aMD simulations run on GPUs can successfully obtain the active and inactive state conformations of a GPCR on this reduced timescale. Furthermore, we discovered a potential alternate active‐state agonist‐binding position in the octopamine receptor which has yet to be observed and may be a novel GPCR agonist‐binding position. These results demonstrate that a complex biological system with an activation process on the millisecond timescale can be successfully simulated on the nanosecond timescale using a simple computing system consisting of a small number of GPUs. Proteins 2016; 84:1480–1489. © 2016 Wiley Periodicals, Inc.     

Predicted structures of agonist and antagonist bound complexes of adenosine A3 receptor

We used the GEnSeMBLE Monte Carlo method to predict ensemble of the 20 best packings (helix rotations and tilts) based on the neutral total energy (E) from a vast number (10 trillion) of potential packings for each of the four subtypes of the adenosine G protein-coupled receptors (GPCRs), which are involved in many cytoprotective functions. We then used the DarwinDock Monte Carlo methods to predict the binding pose for the human A(3) adenosine receptor (hAA(3)R) for subtype selective agonists and antagonists. We found that all four A(3) agonists stabilize the 15th lowest conformation of apo-hAA(3)R while also binding strongly to the 1st and 3rd. In contrast the four A(3) antagonists stabilize the 2nd or 3rd lowest conformation. These results show that different ligands can stabilize different GPCR conformations, which will likely affect function, complicating the design of functionally unique ligands. Interestingly all agonists lead to a trans χ1 angle for W6.48 that experiments on other GPCRs associate with G-protein activation while all 20 apo-AA(3)R conformations have a W6.48 gauche+ χ1 angle associated experimentally with inactive GPCRs for other systems. Thus docking calculations have identified critical ligand-GPCR structures involved with activation. We found that the predicted binding site for selective agonist Cl-IB-MECA to the predicted structure of hAA(3)R shows favorable interactions to three subtype variable residues, I253(6.58), V169(EL2), and Q167(EL2), while the predicted structure for hAA(2A)R shows weakened to the corresponding amino acids: T256(6.58), E169(EL2), and L167(EL2), explaining the observed subtype selectivity.     

In silico approaches towards understanding CALB using molecular dynamics simulation and docking

Candida antarctica lipase B (CALB), a serine protease, is involved in the hydrolysis of substrates at the aqueous lipid interface. There is a significant role played by the helices in serine proteases including acting as a flap covering the active site region. The α5 and α10 helices in the path to the active site of CALB appear to play an important role in the region. This study investigates these helices by mutational studies, docking and molecular dynamics simulations. The mutations were selected based on their proximity to the active site and their presence at the α10-helix in the path of the active site. Molecular dynamics studies reveal the flexibility, stability and hydrogen bonding ability of the α5 helix. The radius of gyration (R _g) clearly showed the compactness of the structure. Docking studies show the changes occurring at the protein's binding site before and after 15 ns of simulation. Results from the study demonstrate the importance of the two helices α5 and α10 in the stability of CALB.     

Modeling and simulation of the human delta opioid receptor

A model for the human delta opioid receptor has been generated via sequence alignment, structure building using the crystal structure of bovine rhodopsin as a template, and refinement by molecular dynamics simulation. The model building suggested that, in addition to the previously postulated interaction between D128 and Y308, an internal salt bridge also exists between residues D128 and R192, both of which are conserved in all the opioid receptors. The model and salt bridge were then shown to be stable during a 20-nsec simulation in a lipid bilayer. It is therefore proposed that both of these interactions play a role in stabilizing the inactive state of the receptor. The model is also used in an effort to rationalize many of the mutational studies performed on delta opioid receptors, and to suggest a plausible explanation for the differences between known delta opioid agonists and antagonists.     

The mechanism of high affinity pentasaccharide binding to antithrombin,insights from Gaussian accelerated molecular dynamics simulations

Abstract

The activity of antithrombin (AT), a serpin protease inhibitor, is enhanced by heparin and heparin analogs against its target proteases, mainly thrombin, factors Xa and IXa. Considerable amount of information is available on the multistep mechanism of the heparin pentasaccharide binding and conformational activation. However, much of the details were inferred from ‘static’ structures obtained by X-ray diffraction. Moreover, limited information is available for the early steps of binding mechanism other than kinetic studies with various ligands. To gain insights into these processes, we performed enhanced sampling molecular dynamics (MD) simulations using the Gaussian Accelerated Molecular Dynamics (GAMD) method, applied previously in drug binding studies. We were able to observe the binding of the pentasaccharide idraparinux to a ‘non-activated’ AT conformation in two separate trajectories with low root mean square deviation (RMSD) values compared to X-ray structures of the bound state. These trajectories along with further simulations of the AT-pentasaccharide complex provided insights into the mechanisms of multiple conformational transitions, including the expulsion of the hinge region, the extension of helix D and the conformational behavior of the reactive center loop (RCL). We could also confirm the high stability of helix P in non-activated AT conformations, such states might play an important role in heparin binding. ‘Generalized correlation’ matrices revealed possible paths of allosteric signal propagation to the binding sites for the target proteases, factors Xa and IXa. Enhanced MD simulations of ligand binding to AT may assist the design of new anticoagulant drugs.

Communicated by Ramaswamy H. Sarma     

Molecular dynamics of NPY Y1 receptor activation

A three-dimensional model of the human neuropeptide Y(NPY)Y₁ receptor (hY₁) was constructed, energy refined and used to simulate molecular receptor interactions of the peptide ligands NPY, [L31, P34]NPY, peptide YY (PYY) and pancreatic polypeptide (PP), and of the nonpeptide antagonist R-N²-(diphenylacetyl)-N-(4-hydroxyphenyl)methyl-argininamide (BIBP3226) and its S-enantiomer BIBP3435. The best complementarity in charges between the receptor and the peptides, and the best structural accordance with experimental studies, was obtained with amino acid 1–4 of the peptides interacting with Asp194, Asp200, Gln201, Phe202 and Trp288 in the receptor. Arg33 and Arg35 of the peptides formed salt bridges with Asp104 and Asp287, respectively, while Tyr36 interacted in a binding pocket formed by Phe41, Thr42, Tyr100, Asn297, His298 and Phe302. Calculated electrostatic potentials around NPY and hY₁ molecules indicated that ligand binding is initiated by electrostatic interactions between a highly positive region in the N- and C-terminal parts of the peptides, and a negative region in the extracellular receptor domains. Molecular dynamics simulations of NPY and BIBP3226 interactions with the receptor indicated rigid body motions of TMH5 and TMH6 upon NPY binding as mechanisms of receptor activation, and that BIBP3226 may act as an antagonist by constraining these motions.     

Molecular dynamics simulation of human neurohypophyseal hormone receptors complexed with oxytocin-modeling of an activated state.

The neurohypophyseal hormone oxytocin (CYIQNCPLG-NH(2), OT) is involved in the control of labor, secretion of milk and many social and behavioral functions via interaction with its receptors (OTR) located in the uterus, mammary glands and peripheral tissues, respectively. In this paper we propose the interactions responsible for OT binding and selectivity to OTR versus vasopressin ([F3,R8]OT, AVP) receptors: V1aR and V2R, all three belonging to the Class A G protein-coupled receptors (GPCRs). Three-dimensional models of the activated receptors were constructed using a multiple sequence alignment and the activated rhodopsin-transducin (MII-Gt) prototype [Slusarz and Ciarkowski, 2004] as a template. The 1 ns unconstrained molecular dynamics (MD) of three pairs of receptor-OT complexes (two complexes per each receptor) immersed in the fully hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC) lipid bilayer was conducted in the AMBER 7.0 force field. The relaxed models of ligand-receptor complexes were used to identify the putative binding sites of OT. The stabilizing interactions with conserved Gln residues in all complexes were identified. The nonconserved hydrophobic residues were proposed as responsible for OTR-OT selectivity and ligand recognition. These results provide guidelines for experimental site-directed mutagenesis and if confirmed, they may be helpful in designing new selective OT analogs with both agonistic or antagonistic properties.     

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

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

In silico studies of conformational dynamics of Mu opioid receptor performed using gaussian accelerated molecular dynamics

G-protein-coupled receptors play a crucial role in various signaling pathways and function as targets for treating a wide spectrum of diseases. Since the twentieth century, extensive research has been conducted on the Mu opioid receptor (MOR) as a drug target. We examined the MOR inactivation and activation processes using an enhanced sampling method (Gaussian accelerated molecular dynamics), the binding pocket site area method, the root mean square deviation method, and the free energy (potential of mean force) method. This study revealed two important intermediate MOR structures (intermediate and intermediate inactive), and the results suggest that the intermediate MOR structure is responsible for the selectivity of opioids.     

Molecular dynamics characterization of the active cavity of carboxypeptidase A and some of its inhibitor adducts.

Molecular dynamics (MD) calculations have been performed on carboxypeptidase A and on its adducts with inhibitors, such as d-phenylalanine (dPhe) and acetate. The catalytically essential zinc ion present in the protein was explicitly included in all the simulations. The simulation was carried out over a sphere of 15 A centered on the zinc ion. The crystallographic water molecules were explicitly taken into account; then the protein was solvated with a 18 A sphere of water molecules. MD calculations were carried out for 45-60 ps. There is no large deviation from the available X-ray structures of native and the dPhe adduct for the MD structures. Average MD structures were calculated starting from the X-ray structure of the dPhe adduct, and, from a structure obtained by docking the inhibitor in the native structure. Comparison between these two structures and with that of the native protein shows that some of the key variations produced by inhibitor binding are reproduced by MD calculations. Addition of acetate induces structural changes relevant for the understanding of the interaction network in the active cavity. The structural variations induced by different inhibitors are examined. The effects of these interactions on the catalytic mechanism and on the binding of substrate are discussed.     

Deciphering the GPER/GPR30-agonist and antagonists interactions using molecular modeling studies,molecular dynamics,and docking simulations

The G-protein coupled estrogen receptor 1 GPER/GPR30 is a transmembrane seven-helix (7TM) receptor involved in the growth and proliferation of breast cancer. Due to the absence of a crystal structure of GPER/GPR30, in this work, molecular modeling studies have been carried out to build a three-dimensional structure, which was subsequently refined by molecular dynamics (MD) simulations (up to 120 ns). Furthermore, we explored GPER/GPR30’s molecular recognition properties by using reported agonist ligands (G1, estradiol (E2), tamoxifen, and fulvestrant) and the antagonist ligands (G15 and G36) in subsequent docking studies. Our results identified the E2 binding site on GPER/GPR30, as well as other receptor cavities for accepting large volume ligands, through GPER/GPR30 π–π, hydrophobic, and hydrogen bond interactions. Snapshots of the MD trajectory at 14 and 70 ns showed almost identical binding motifs for G1 and G15. It was also observed that C107 interacts with the acetyl oxygen of G1 (at 14 ns) and that at 70 ns the residue E275 interacts with the acetyl group and with the oxygen from the other agonist whereas the isopropyl group of G36 is oriented toward Met141, suggesting that both C107 and E275 could be involved in the protein activation. This contribution suggest that GPER1 has great structural changes which explain its great capacity to accept diverse ligands, and also, the same ligand could be recognized in different binding pose according to GPER structural conformations.     

Molecular modeling,dynamics studies and density functional theory approaches to identify potential inhibitors of SIRT4 protein from Homo sapiens : a novel target for the treatment of type 2 diabetes

Type 2 diabetes is one of the biggest health challenges in the world and WHO projects it to be the 7th leading cause of death in 2030. It is a chronic condition affecting the way our body metabolizes sugar. Insulin resistance is high risk factor marked by expression of Lipoprotein Lipases and Peroxisome Proliferator-Activated Receptor that predisposes to type 2 diabetes. AMP-dependent protein kinase in AMPK signaling pathway is a central sensor of energy status. Deregulation of AMPK signaling leads to inflammation, oxidative stress, and deactivation of autophagy which are implicated in pathogenesis of insulin resistance. SIRT4 protein deactivates AMPK as well as directly inhibits insulin secretion. SIRT4 overexpression leads to dyslipidimeia, decreased fatty acid oxidation, and lipogenesis which are the characteristic features of insulin resistance promoting type 2 diabetes. This makes SIRT4 a novel therapeutic target to control type 2 diabetes. Virtual screening and molecular docking studies were performed to obtain potential ligands. To further optimize the geometry of protein–ligand complexes Quantum Polarized Ligand Docking was performed. Binding Free Energy was calculated for the top three ligand molecules. In view of exploring the stereoelectronic features of the ligand, density functional theory approach was implemented at B3LYP/6-31G* level. 30 ns MD simulation studies of the protein–ligand complexes were done. The present research work proposes ZINC12421989 as potential inhibitor of SIRT4 with docking score (?7.54 kcal/mol), docking energy (?51.34 kcal/mol), binding free energy (?70.21 kcal/mol), and comparatively low energy gap (?0.1786 eV) for HOMO and LUMO indicating reactivity of the lead molecule.     

Computational studies of LXR molecular interactions reveal an allosteric communication pathway

The liver X receptor, LXRα, is an important regulator of genes involved in metabolism and inflammation. The mechanism of communication between the cofactor peptide and the ligand in the ligand-binding pocket is a crucial and often discussed issue for the nuclear receptors (NRs), but such allosteric signaling pathways are difficult to detect and the transmission mechanism remains elusive. Here, we apply the anisotropic thermal diffusion method to the LXRα with bound coactivator and ligand. We detected a possible communication pathway between the coactivator peptide and the ligand. The signal is transmitted both through the receptor backbone and side chains. A key signaling residue is the first leucine in the cofactor peptide recognition motif LXXLL, which is conserved within the NR cofactors, suggesting a general mechanism for allosteric signaling. Furthermore, we studied the LXR receptor and cofactor molecular interactions in detail using molecular dynamics simulations. The protein-protein interaction patterns in the complexes of nine different cofactor peptides and holo-LXRα were characterized, revealing the importance of the receptor-cofactor charge clamp interaction. Specific, but infrequently occurring interactions were observed toward the cofactor peptide C-terminal residues. Thus, additional specificity between LXRα and its cofactors is likely to be found in molecular interactions outside the cofactor peptide or in other biological factors.     

Analysis of interactions responsible for vasopressin binding to human neurohypophyseal hormone receptors-molecular dynamics study of the activated receptor-vasopressin-G(alpha) systems.

Vasopressin (CYFQNCPRG-NH(2), AVP) is a semicyclic endogenous peptide, which exerts a variety of biological effects in mammals. The main physiological roles of AVP are the regulation of water balance and the control of blood pressure and adrenocorticotropin hormone (ACTH) secretion, mediated via three different subtypes of vasopressin receptors: V1a, V1b and V2 receptors (V1aR, V1bR and V2R, respectively). They are the members of the class A, G-protein-coupled receptors (GPCRs). AVP also modulates several behavioral and social functions. In this study, the interactions responsible for AVP binding to vasopressin V1a and V2 receptors versus the closely related oxytocin ([I3,L8]AVP, OT) receptor (OTR) have been investigated. Three-dimensional models of the activated receptors were constructed using multiple sequence alignment, followed by homology modeling using the complex of activated rhodopsin with Gt(alpha) C-terminal peptide of transducin MII-Gt(338-350) prototype as a template. AVP was docked into the receptor-G(alpha) systems. The three lowest-energy pairs of receptor-AVP-G(alpha) (two complexes per each receptor) were selected. The 1-ns unconstrained molecular dynamics (MD) of complexes embedded into the fully hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC) lipid bilayer was conducted in the AMBER 7.0 force field. Six relaxed receptor-AVP-G(alpha) models were obtained. The residues responsible for AVP binding to vasopressin receptors have been identified and a different mechanism of AVP binding to V2R than to V1aR has been proposed.     

Structural dynamics of an ionotropic glutamate receptor

Ionotropic glutamate receptors (iGluRs) are postsynaptic ion channels involved in excitatory neurotransmission. iGluRs play important roles in development and in forms of synaptic plasticity that underlie higher order processes such as learning and memory. Neurobiological and biochemical studies have long characterized iGluRs in detail. However, the structural basis for the function of iGluRs has not yet been investigated, because there is insufficient information about their three-dimensional structures. In 1998, a crystal structure called S1S2 lobes was first solved for the extracellular bilobed ligand-binding domain of the GluR2 subunit. Since then, the crystal structures for the S1S2 lobes both in the apo and in various liganded states have been reported, and recent biophysical studies have further elucidated the dynamic aspects of the structure of the S1S2 lobes. In this review, the dynamic structures of the S1S2 lobes and their ligands are summarized, and the importance of their structural flexibility and fluctuation is discussed in light of the mechanisms of ligand recognition, activation, and desensitization of the receptor.     

Glutamate and Glycine Binding to the NMDA Receptor

1. Download high-res image (273KB)
2. Download full-size image

     

Location and nature of the residues important for ligand recognition in G-protein coupled receptors

The overall structure of the biogenic amine subclass of the G-protein-coupled receptors, and of their ligand binding sites, is discussed with the aim of highlighting the major structural features of these receptors that are responsible for ligand recognition. A comparison is made between biogenic amine receptors, peptide receptors of the rhodopsin class, and the secretin receptors which all have peptide ligands. The question of where the peptide ligands bind, whether at extracellular sites or within the transmembrane helix bundle, is discussed. The suitability of the rhodopsin crystal structure as a template for construction of homology models is discussed and it is concluded that there are many reasons why a caution should be issued against using it uncritically.