Tetrahydroisoquinolines acting as dopaminergic ligands. A molecular modeling study using MD simulations and QM calculations

A molecular modeling study on 16 1-benzyl tetrahydroisoquinolines (BTHIQs) acting as dopaminergic ligands was carried out. By combining molecular dynamics simulations with ab initio and density functional theory (DFT) calculations, a simple and generally applicable procedure to evaluate the binding energies of BTHIQs interacting with the human dopamine D2 receptor (D2 DR) is reported here, providing a clear picture of the binding interactions of BTHIQs from both structural and energetic viewpoints. Molecular aspects of the binding interactions between BTHIQs and the D2 DR are discussed in detail. A significant correlation between binding energies obtained from DFT calculations and experimental pKi values was obtained, predicting the potential dopaminergic effect of non-synthesized BTHIQs.     

Synergistic effects in the designs of neuraminidase ligands: Analysis from docking and molecular dynamics studies

Docking and molecular dynamics were used to study the nine ligands (see Scheme 1) at the neuraminidase (NA) active sites. Their binding modes are structurally and energetically different, with details given in the text. Compared with 1A (oseltamivir carboxylate), the changes of core template or/and functional groups in the other ligands cause the reductions of interaction energies and numbers of H-bonds with the NA proteins. Nonetheless, all these ligands occupy the proximity space at the NA active sites and share some commonness in their binding modes. The fragment approach was then used to analyze and understand the binding specificities of the nine ligands. The contributions of each core template and functional group were evaluated. It was found that the core templates rather than functional groups play a larger role during the binding processes; in addition, the binding qualities are determined by the synergistic effects of the core templates and functional groups. Among the nine ligands, 1A (oseltamivir carboxylate) has the largest synergistic energy and its functional groups fit perfectly with the NA active site, consistent with the largest interaction energy, numerous H-bonds with the NA active-site residues as well as experimentally lowest IC₅₀ value. Owing to the poorer metabolizability than oseltamivir, large contribution of the benzene core template and fine synergistic effects of the functional groups, the 4-(N-acetylamino)-5-guanidino-3-(3-pentyloxy)benzoic acid should be an ideal lead compound for optimizing NA drugs.     

Human aryl-hydrocarbon receptor and its interaction with dioxin and physiological ligands investigated by molecular modelling and docking simulations

Molecular structure of the ligand binding domain of hAhR has been modelled by homology modelling techniques and used for docking simulations with dioxin and nine more xenobiotics and endogenous ligands. The study evidences that different sites may bind these ligands, whereas only one binding site has been previously indicated by past studies on the mouse homologous receptor. The differences in the sequence of mouse and human AhR ligand binding domain may explain this observation, being most of them in the additional sites observed. Preferences of the evaluated ligands for the different sites are reported and discussed in view of their functional role.     

Identification of nonplanar small molecule for G-quadruplex grooves: molecular docking and molecular dynamic study

DNA G-quadruplex is an attractive drug target for anticancer therapy. Most G-quadruplex ligands have little selectivity, due to π-stacking interaction with common G-tetrads surface. Thanks to the varieties of G-quadruplex grooves, the groove-binding ligand is expected to create high selectivity. Therefore, developing novel molecular geometries that target G-quadruplex groove has been paid growing attention. In this work, steroid FG, a special nonplanar and nonaromatic small molecule, interacting with different conformations of G-quadruplexes has been studied by molecular docking and molecular dynamics simulations. The results showed the selectivity of the hydrophobic group of steroid FG for the wide groove of antiparallel G-quadruplex. The methyl groups on the tetracyclic ring of steroid represent the specific binding ability for the small hydrophobic cavity formed by reversed stacking of G-tetrads in antiparallel G-quadruplex groove. This work provides new insight for developing new classes of G-quadruplex groove-binding ligands.     

Modeling molecular mechanisms of binding of the anaphylatoxin C5a to the C5a receptor

This study presents the 3D model of the complex between the anaphylatoxin C5a and its specific receptor, C5aR. This is the first 3D model of a G-protein-coupled receptor (GPCR) complex with a peptide ligand deduced by a molecular modeling procedure analyzing various conformational possibilities of the extracellular loops and the N-terminal segment of the GPCR. The modeling results indicated two very different ways of interacting between C5a and C5aR at the two interaction sites suggested earlier based on the data of site-directed mutagenesis. Specifically, C5a and C5aR can be involved in "mutual-induced fit", where the interface between the molecules is determined by both the receptor and the ligand. The rigid core of the C5a ligand selects the proper conformations of the highly flexible N-terminal segment of C5aR (the first interaction site). At the same time, the binding conformation of the flexible C-terminal fragment of C5a is selected by well-defined interactions with the TM region of the C5aR receptor (the second interaction site). The proposed 3D model of C5a/C5aR complex was built without direct use of structural constraints derived from site-directed mutagenesis reserving those data for validation of the model. The available data of site-directed mutagenesis of C5a and C5aR were successfully rationalized with the help of the model. Also, the modeling results predicted that the full-length C5a and C5a-des74 metabolite would have different binding modes with C5aR. Modeling approaches employed in this study are readily applicable for studies of molecular mechanisms of binding of other polypeptide ligands to their specific GPCRs.     

Functional interpretation of APN receptor from M.sexta using a molecular model

Insect pests are the major cause of damage to commercially important agricultural crops. The continuous application of synthetic pesticides resulted in severe insect resistance by plants. This causes irreversible damage to the environment. Bacillus thuringiensis (Bt) emerged as a valuable biological alternative in pest control. However, insect resistance against Bt has been reported in many cases. Insects develop resistance to insecticides through mechanisms that reduce the binding of toxins to gut receptors. Nonetheless, the molecular mechanism of insect resistance is not fully understood. Therefore, it is important to study the mechanism of toxin resistance by analyzing amino‐peptidase‐N (APN) receptor of the insect M. sexta. A homology model of APN was constructed using Insight II molecular modeling software and the model was further evaluated using the PROCHECK program. Oligosaccharides participating in post translational modification were constructed and docked onto specific APN functional sites. Post analyses of the APN model provide insights on the functional properties of APN towards the understanding of receptor and toxin interactions. We also discuss the predicted binding sites for ligands, metals and Bt toxins in M. sexta APN receptor. These data help in the development of a roadmap for the design and synthesis of novel insect resistant Cry toxins.     

Structure and localisation of drug binding sites on neurotransmitter transporters

The dopamine (DAT), serotontin (SERT) and noradrenalin (NET) transporters are molecular targets for different classes of psychotropic drugs. The crystal structure of Aquifex aeolicus LeuT_Aa was used as a template for molecular modeling of DAT, SERT and NET, and two putative drug binding sites (pocket 1 and 2) in each transporter were identified. Cocaine was docked into binding pocket 1 of DAT, corresponding to the leucine binding site in LeuT_Aa, which involved transmembrane helices (TMHs) 1, 3, 6 and 8. Clomipramine was docked into binding pocket 2 of DAT, involving TMHs 1, 3, 6, 10 and 11, and extracellular loops 4 and 6, corresponding to the clomipramine binding site in a crystal structure of a LeuT_Aa–clomipramine complex. The structures of the proposed cocaine- and tricyclic antidepressant-binding sites may be of particular interest for the design of novel DAT interacting ligands.     

Regulation of the black bullhead hepatic beta-adrenoceptors

Black bullhead catfish (Ameiurus melas) were exposed to air for 1 h to examine the effect of an acute stress on the distribution and function of the hepatic beta-adrenoceptors (beta-ARs). Air exposure significantly reduced both adrenaline (ADR)- and noradrenaline (NADR)-stimulated glucose production in isolated hepatocytes with no effect on either receptor affinity (K(d)) or number of binding sites (B(max)). A 24 h exposure of isolated hepatocytes to the beta-agonist isoproterenol also had no significant impact on either binding parameter. Competition studies using selective agonists and antagonists suggest that the hepatic beta-AR in this species is pharmacologically beta(2)-like. However in addition to the beta(2)-AR, molecular evidence provides support for the existence of hepatic beta-ARs that phylogenetically group with the beta(3)-ARs and the beta(1)-ARs. Despite the presence of several potential phosphorylation sites in the third intracellular loop and cytoplasmic tail of the bullhead beta(2)-AR, no significant changes were observed in the binding parameters. While physiological data supports the presence of only a single subtype, molecular data supports the existence of multiple beta-AR subtypes in this species. The mechanisms thought to regulate mammalian beta-ARs exist in the bullhead ARs reported here but these mechanisms are not as effective in this fish system as in mammals.     

Characterization of the modes of binding between human sweet taste receptor and low-molecular-weight sweet compounds

One of the most distinctive features of human sweet taste perception is its broad tuning to chemically diverse compounds ranging from low-molecular-weight sweeteners to sweet-tasting proteins. Many reports suggest that the human sweet taste receptor (hT1R2-hT1R3), a heteromeric complex composed of T1R2 and T1R3 subunits belonging to the class C G protein-coupled receptor family, has multiple binding sites for these sweeteners. However, it remains unclear how the same receptor recognizes such diverse structures. Here we aim to characterize the modes of binding between hT1R2-hT1R3 and low-molecular-weight sweet compounds by functional analysis of a series of site-directed mutants and by molecular modeling-based docking simulation at the binding pocket formed on the large extracellular amino-terminal domain (ATD) of hT1R2. We successfully determined the amino acid residues responsible for binding to sweeteners in the cleft of hT1R2 ATD. Our results suggest that individual ligands have sets of specific residues for binding in correspondence with the chemical structures and other residues responsible for interacting with multiple ligands.     

Theoretical aspects of DNA-protein interactions: co-operative and non-co-operative binding of large ligands to a one-dimensional homogeneous lattice

The interaction of proteins binding non-specifically to DNA, as well as the properties of many other interacting ligand-lattice systems important in molecular biology, requires a fundamentally different type of theoretical analysis than that provided by the classical Scatchard independent-binding-site treatment. Exact and relatively simple equations describing the binding of both non-interacting and interacting (co-operative) ligands to a homogeneous one-dimensional lattice are derived in terms of ligand site size, intrinsic binding constant and ligand-ligand co-operativity (equations (10) and (15) in the text). The mathematical approach is based on simple conditional probabilities, and reveals some largely unrecognized characteristics of such lattice binding systems. The results indicate that the binding of any non-interacting ligand covering more than one lattice residue results in non-linear (convex downward) Scatchard plots. The introduction of positive ligand-ligand co-operativity antagonizes this non-linearity, and eventually leads to plots of the opposite curvature. The maxima, limiting slopes, and intercepts of such plots can be used to estimate the required binding parameters. The method can be extended to systems involving heterogeneous ligands, and some types of heterogeneous lattices. Procedures for applying the method to a variety of interacting systems are presented, and a preliminary analysis is carried out for some selected sets of data from the literature.     

Expanding GPCR homology model binding sites via a balloon potential: A molecular dynamics refinement approach

Homology modeling of G protein-coupled receptors is becoming a widely used tool in drug discovery. However, unrefined models built using the bovine rhodopsin crystal structure as the template, often have binding sites that are too small to accommodate known ligands. Here, we present a novel systematic method to refine model active sites based on a pressure-guided molecular dynamics simulation. A distinct advantage of this approach is the ability to introduce systematic perturbations in model backbone atoms in addition to side chain adjustments. The method is validated on two test cases: (1) docking of retinal into an MD-relaxed structure of opsin and (2) docking of known ligands into a homology model of the CCR2 receptor. In both cases, we show that the MD expansion algorithm makes it possible to dock the ligands in poses that agree with the crystal structure or mutagenesis data.     

Identification of the ligand binding sites on the molecular surface of proteins

Identification of protein biochemical functions based on their three-dimensional structures is now required in the post-genome-sequencing era. Ligand binding is one of the major biochemical functions of proteins, and thus the identification of ligands and their binding sites is the starting point for the function identification. Previously we reported our first trial on structure-based function prediction, based on the similarity searches of molecular surfaces against the functional site database. Here we describe the extension of our first trial by expanding the search database to whole heteroatom binding sites appearing within the Protein Data Bank (PDB) with the new analysis protocol. In addition, we have determined the similarity threshold line, by using 10 structure pairs with solved free and complex structures. Finally, we extensively applied our method to newly determined hypothetical proteins, including some without annotations, and evaluated the performance of our methods.     

Sequence analysis of serum albumins reveals the molecular evolution of ligand recognition properties

Serum albumin (SA) is a circulating protein providing a depot and carrier for many endogenous and exogenous compounds. At least seven major binding sites have been identified by structural and functional investigations mainly in human SA. SA is conserved in vertebrates, with at least 49 entries in protein sequence databases. The multiple sequence analysis of this set of entries leads to the definition of a cladistic tree for the molecular evolution of SA orthologs in vertebrates, thus showing the clustering of the considered species, with lamprey SAs (Lethenteron japonicum and Petromyzon marinus) in a separate outgroup. Sequence analysis aimed at searching conserved domains revealed that most SA sequences are made up by three repeated domains (about 600 residues), as extensively characterized for human SA. On the contrary, lamprey SAs are giant proteins (about 1400 residues) comprising seven repeated domains. The phylogenetic analysis of the SA family reveals a stringent correlation with the taxonomic classification of the species available in sequence databases. A focused inspection of the sequences of ligand binding sites in SA revealed that in all sites most residues involved in ligand binding are conserved, although the versatility towards different ligands could be peculiar of higher organisms. Moreover, the analysis of molecular links between the different sites suggests that allosteric modulation mechanisms could be restricted to higher vertebrates.     

Binding of ibuprofen to human hemoglobin: elucidation of their molecular recognition by spectroscopy,calorimetry, and molecular modeling techniques

Ibuprofen, used for the treatment of acute and chronic pain, osteoarthritis, rheumatoid arthritis, and related conditions has ample affinity to globular proteins. Here we have explored this fundamental study pertaining to the interaction of ibuprofen with human hemoglobin (HHb), using multispectroscopic, calorimetric, and molecular modeling techniques to gain insights into molecular aspects of binding mechanism. Ibuprofen-induced graded decrease in absorption spectra indicates protein disruption along with sedimentation of HHb particle. Red shifting of absorption peak at 195 nm indicates alteration in the secondary structure of HHb upon interaction with ibuprofen. Flouremetric and isothermal titration calorimetric (ITC) studies suggested one binding site in HHb for ibuprofen at 298.15 K. However, with increase in temperature, ITC revealed increasing number of binding sites. The negative values of Gibbs energy change (ΔG⁰) and enthalpy change (ΔH⁰) along with positive value of entropy change (ΔS⁰) strongly suggest that it is entropy-driven spontaneous exothermic reaction. Moreover, hydrophobic interaction, hydrogen bonding, and π–π interaction play major role in this binding process as evidenced from ANS (8-anilino-1-napthalenesulphonic acid), sucrose binding, and molecular modeling studies. The interaction impacts on structural integrity and functional aspects of HHb as confirmed by CD spectroscopy, increased free iron release, increased rate of co-oxidation and decreased rate of esterase activity. These findings suggest us to conclude that ibuprofen upon interaction perturbs both structural and functional aspects of HHb.     

Functional molecular mass of binding sites for [3H]dihydrotetrabenazine and [3H]reserpine and of dopamine beta-hydroxylase and cytochrome b561 from chromaffin granule membrane as determined by radiation inactivation

The monoamine transporter of chromaffin granule membrane has two distinct high-affinity binding sites for tetrabenazine and reserpine, which can be assayed by [3H]dihydrotetrabenazine and [3H]reserpine binding, respectively. The functional molecular mass of the components bearing these sites has been investigated by the radiation inactivation technique. The decline of [3H]dihydrotetrabenazine binding activity with increasing radiation doses followed a single exponential, from which a functional molecular mass of 68 kDa was derived for tetrabenazine binding sites. [3H]Reserpine binding activity declined in a more complex way; however, under conditions where high-affinity reserpine binding sites were specifically assayed, the decline was also exponential, corresponding to a functional molecular mass of 37 kDa for these sites. The figures obtained for high-affinity tetrabenazine and reserpine binding sites are consistent with previous values obtained by photoaffinity of tetrabenazine and serotonin binding sites, respectively. It is thus concluded that the monoamine transporter has an oligomeric structure. By the radiation inactivation technique, cytochrome b561 and dopamine beta-hydroxylase have functional molecular masses of 25 and 123 kDa, respectively. The latter value might be attributed to the dimeric form of the enzyme.     

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.     

A new protein binding pocket similarity measure based on comparison of clouds of atoms in 3D: application to ligand prediction

Background  

Predicting which molecules can bind to a given binding site of a protein with known 3D structure is important to decipher the protein function, and useful in drug design. A classical assumption in structural biology is that proteins with similar 3D structures have related molecular functions, and therefore may bind similar ligands. However, proteins that do not display any overall sequence or structure similarity may also bind similar ligands if they contain similar binding sites. Quantitatively assessing the similarity between binding sites may therefore be useful to propose new ligands for a given pocket, based on those known for similar pockets.