Electrostatic components of drug-receptor recognition. II. The DNA-binding antibiotic actinomycin

The semiempirical molecular orbital method PCILO has been used to study the electronic charge distribution and conformation of the antibiotic actinomycin. Molecular electrostatic potential fields have been constructed in three-dimensional space round the molecule and displayed stereoscopically together with contour maps in the orthogonal planes relative to that of the chromophore. It was found that the drug has a considerable dipole moment and that the electrostatic fields can be separated into two large regions; an area of positive potential surrounds the chromophore and one of negative potential is directed away from the peptide rings. In this and the preceding paper these observations are discussed with respect to drug-receptor recognition phenomena. Pattern-matching of complementary electrostatic fields between the drug and polynucleotide receptor can be discerned.     

How does stress affect human being—a molecular dynamic simulation study on cortisol and its glucocorticoid receptor

Stress can be either positive or negative to human beings. Under stressful conditions, the mental and physical conditions of human can be affected. There exists certain relation between stress and illness. The cortisol and other glucocorticoids bind to the same receptor, which is called glucocorticoid receptor. Some evidences indicated that cortisol molecule binding to its glucocorticoid receptor was necessary for the stress response. Up to now, the structure–function relationships between cortisol molecule and its glucocorticoid receptor have not been deliberated from the atomic-level. In order to get a detailed understanding of the structure–function relationships between the cortisol molecule and glucocorticoids receptor, we have carried out molecular dynamic (MD) simulations on glucocorticoid receptor (Apo system) and cortisol with its glucocorticoid receptor complex (HCY system). On the basis of molecular dynamic simulations, a couple of key residues were identified, which were crucial for the binding of cortisol molecule. The results of binding free energy calculations are in good agreement with the experiment data. Our research gives clear insights from atomic-level into the structural–functional aspects of cortisol molecule and its glucocorticoid receptor, and also provides valuable information for the design of drug which can treat stress related illnesses.     

Electrostatic Similarities between Protein and Small Molecule Ligands Facilitate the Design of Protein-Protein Interaction Inhibitors

One of the underlying principles in drug discovery is that a biologically active compound is complimentary in shape and molecular recognition features to its receptor. This principle infers that molecules binding to the same receptor may share some common features. Here, we have investigated whether the electrostatic similarity can be used for the discovery of small molecule protein-protein interaction inhibitors (SMPPIIs). We have developed a method that can be used to evaluate the similarity of electrostatic potentials between small molecules and known protein ligands. This method was implemented in a software called EleKit. Analyses of all available (at the time of research) SMPPII structures indicate that SMPPIIs bear some similarities of electrostatic potential with the ligand proteins of the same receptor. This is especially true for the more polar SMPPIIs. Retrospective analysis of several successful SMPPIIs has shown the applicability of EleKit in the design of new SMPPIIs.     

Probing the effect of electric field in 9,10-dimethoxy-2,6-bis(2-p-tolylethynyl)anthracene molecular nanowire using quantum chemical and charge density analysis

The effect of electric field (EF) in a newly designed molecular nanowire 9,10-dimethoxy-2,6-bis(2-p-tolylethynyl)anthracene has been analysed theoretically from the structural and electronic charge transport properties using quantum chemical and charge density calculations. The applied EF (0–0.36 VÅ^? 1) alters the molecular conformation, charge density distribution, electrostatic properties and the electronic energy levels of the molecule. Furthermore, the applied EF decreases the highest occupied molecular orbital–lowest unoccupied molecular orbital gap significantly from 1.775 to 0.258 eV and it also induces polarisation in the molecule, which leads to increase the dipole moment of the molecule. The electrostatic potential for various levels of applied EF reveals the charge-accumulated regions of the molecule. The I–V characteristics of the molecule have been studied against various applied fields using Landauer formalism.     

Small molecule mimetics of an interferon-α receptor interacting domain

Small molecules that mimic IFN-α epitopes that interact with the cell surface receptor, IFNAR, would be useful therapeutics. One such 8-amino acid region in IFN-α2, designated IRRP-1, was used to derive 11 chemical compounds that belong to 5 distinct chemotypes, containing the molecular features represented by the key residues Leu30, Arg33, and Asp35 in IRRP-1. Three of these compounds exhibited potential mimicry to IRRP-1 and, in cell based assays, as predicted, effectively inhibited IFNAR activation by IFN-α. Of these, compound 3 did not display cell toxicity and reduced IFN-α-inducible STAT1 phosphorylation and STAT-DNA binding. Based on physicochemical properties’ analyses, our data suggest that moieties with acidic pKa on the small molecule may be a necessary element for mimicking the carboxyl group of Asp35 in IRRP-1. Our data confirm the relevance of this strategy of molecular mimicry of ligand–receptor interaction domains of protein partners for small molecule drug discovery.     

Integrin triplets of marine sponges in human D2 receptor heteromers

The evidence for the existence of receptor heteromers opens up a new field for a better understanding of neural transmission. Based on our theory, we have discovered main triplets of amino acid residues in cell-adhesion receptors of marine sponges, which appear also as homologies in several dopamine D2 receptor heteromers of human brain. The obtained results probably mean a general molecular mechanism for receptor–receptor interactions in heteromers originated from the lowest animals (marine sponges).     

血管紧张素转换酶抑制剂的分子设计

目的:以血管紧张素转换酶(ACE)为受体靶蛋白,设计相应的抑制剂分子。方法:采用计算机辅助药物分子设计的方法,在分析已有的关于ACE及其抑制剂的基础上,确定分子设计软件LigBuilder中的相关参数。结果:设计得出了1500个配体分子,并从中筛选出了1个各项得分最高的新型血管紧张素转换酶抑制剂分子,其分子式为C27H36N2O3,分子量为436,logP值为4.95,与受体靶蛋白ACE的结合力得分为11.71。结论:从计算得分来看,所设计出的分子与ACE的亲和性好,且脂溶性高,为发掘高活性的降压药物提供了线索。     

Exploration of the conformational landscape in pregnane X receptor reveals a new binding pocket

Ligand‐regulated pregnane X receptor (PXR), a member of the nuclear receptor superfamily, plays a central role in xenobiotic metabolism. Despite its critical role in drug metabolism, PXR activation can lead to adverse drug‐drug interactions and early stage metabolism of drugs. Activated PXR can induce cancer drug resistance and enhance the onset of malignancy. Since promiscuity in ligand binding makes it difficult to develop competitive inhibitors targeting PXR ligand binding pocket (LBP), it is essential to identify allosteric sites for effective PXR antagonism. Here, molecular dynamics (MD) simulation studies unravelled the existence of two different conformational states, namely “expanded” and “contracted”, in apo PXR ligand binding domain (LBD). Ligand binding events shifted this conformational equilibrium and locked the LBD in a single “ligand‐adaptable” conformational state. Ensemble‐based computational solvent mapping identified a transiently open potential small molecule binding pocket between α5 and α8 helices, named “α8 pocket”, whose opening‐closing mechanism directly correlated with the conformational shift in LBD. A virtual hit identified through structure‐based virtual screening against α8 pocket locks the pocket in its open conformation. MD simulations further revealed that the presence of small molecule at allosteric site disrupts the LBD dynamics and locks the LBD in a “tightly‐contracted” conformation. The molecular details provided here could guide new structural studies to understand PXR activation and antagonism.     

Advancing Aptamers as Molecular Probes for Cancer Theranostic Applications—The Role of Molecular Dynamics Simulation

Theranostics cover emerging technologies for cell biomarking for disease diagnosis and targeted introduction of drug ingredients to specific malignant sites. Theranostics development has become a significant biomedical research endeavor for effective diagnosis and treatment of diseases, especially cancer. An efficient biomarking and targeted delivery strategy for theranostic applications requires effective molecular coupling of binding ligands with high affinities to specific receptors on the cancer cell surface. Bioaffinity offers a unique mechanism to bind specific target and receptor molecules from a range of non‐targets. The binding efficacy depends on the specificity of the affinity ligand toward the target molecule even at low concentrations. Aptamers are fragments of genetic materials, peptides, or oligonucleotides which possess enhanced specificity in targeting desired cell surface receptor molecules. Aptamer–target binding results from several inter‐molecular interactions including hydrogen bond formation, aromatic stacking of flat moieties, hydrophobic interaction, electrostatic, and van der Waals interactions. Advancements in Systematic Evolution of Ligands by Exponential Enrichment (SELEX) assay has created the opportunity to artificially generate aptamers that specifically bind to desired cancer and tumor surface receptors with high affinities. This article discusses the potential application of molecular dynamics (MD) simulation to advance aptamer‐mediated receptor targeting in targeted cancer therapy. MD simulation offers real‐time analysis of the molecular drivers of the aptamer‐receptor binding and generate optimal receptor binding conditions for theranostic applications. The article also provides an overview of different cancer types with focus on receptor biomarking and targeted treatment approaches, conventional molecular probes, and aptamers that have been explored for cancer cells targeting.     

Receptor Activity-modifying Proteins 2 and 3 Generate Adrenomedullin Receptor Subtypes with Distinct Molecular Properties

Adrenomedullin (AM) is a peptide hormone with numerous effects in the vascular systems. AM signals through the AM₁ and AM₂ receptors formed by the obligate heterodimerization of a G protein-coupled receptor, the calcitonin receptor-like receptor (CLR), and receptor activity-modifying proteins 2 and 3 (RAMP2 and RAMP3), respectively. These different CLR-RAMP interactions yield discrete receptor pharmacology and physiological effects. The effective design of therapeutics that target the individual AM receptors is dependent on understanding the molecular details of the effects of RAMPs on CLR. To understand the role of RAMP2 and -3 on the activation and conformation of the CLR subunit of AM receptors, we mutated 68 individual amino acids in the juxtamembrane region of CLR, a key region for activation of AM receptors, and determined the effects on cAMP signaling. Sixteen CLR mutations had differential effects between the AM₁ and AM₂ receptors. Accompanying this, independent molecular modeling of the full-length AM-bound AM₁ and AM₂ receptors predicted differences in the binding pocket and differences in the electrostatic potential of the two AM receptors. Druggability analysis indicated unique features that could be used to develop selective small molecule ligands for each receptor. The interaction of RAMP2 or RAMP3 with CLR induces conformational variation in the juxtamembrane region, yielding distinct binding pockets, probably via an allosteric mechanism. These subtype-specific differences have implications for the design of therapeutics aimed at specific AM receptors and for understanding the mechanisms by which accessory proteins affect G protein-coupled receptor function.     

The changing world of G protein-coupled receptors: from monomers to dimers and receptor mosaics with allosteric receptor–receptor interactions

Based on indications of direct physical interactions between neuropeptide and monoamine receptors in the early 1980s, the term receptor–receptor interactions was introduced and later on the term receptor heteromerization in the early 1990s. Allosteric mechanisms allow an integrative activity to emerge either intramolecularly in G protein-coupled receptor (GPCR) monomers or intermolecularly via receptor–receptor interactions in GPCR homodimers, heterodimers, and receptor mosaics. Stable heteromers of Class A receptors may be formed that involve strong high energy arginine–phosphate electrostatic interactions. These receptor–receptor interactions markedly increase the repertoire of GPCR recognition, signaling and trafficking in which the minimal signaling unit in the GPCR homomers appears to be one receptor and one G protein. GPCR homomers and GPCR assemblies are not isolated but also directly interact with other proteins to form horizontal molecular networks at the plasma membrane.     

Benzodiazepine ligands can act as allosteric modulators of the Type 1 cholecystokinin receptor

The cholecystokinin (CCK₁) receptor is a G protein-coupled receptor important for nutrient homeostasis. The molecular basis of CCK-receptor binding has been debated, with one prominent model suggesting occupation of the same region of the intramembranous helical bundle as benzodiazepines. Here, we used a specific assay of allosteric ligand interaction to probe the mode of binding of devazepide, a prototypic benzodiazepine ligand. Devazepide elicited marked slowing of dissociation of pre-bound CCK, only possible through binding to a topographically distinct allosteric site. This effect was disrupted by chemical modification of a cysteine in the benzodiazepine-binding pocket. Application of an allosteric model to the equilibrium interaction between a series of benzodiazepine ligands and CCK yielded quantitative estimates of each modulator’s affinity for the allosteric site, as well as the degree of negative cooperativity for the interaction between occupied orthosteric and allosteric sites. The allosteric nature of benzodiazepine binding to the CCK₁ receptor provides new opportunities for small molecule drug development.     

Formation of the 67-kDa laminin receptor by acylation of the precursor

Even though the involvement of the 67-kDa laminin receptor (67LR) in tumor invasiveness has been clearly demonstrated, its molecular structure remains an open problem, since only a full-length gene encoding a 37-kDa precursor protein (37LRP) has been isolated so far. A pool of recently obtained monoclonal antibodies directed against the recombinant 37LRP molecule was used to investigate the processing that leads to the formation of the 67-kDa molecule. In soluble extracts of A431 human carcinoma cells, these reagents recognize the precursor molecule as well as the mature 67LR and a 120-kDa molecule. The recovery of these proteins was found to be strikingly dependent upon the cell solubilization conditions: the 67LR is soluble in NP-40-lysis buffer whereas the 37LRP is NP-40-insoluble. Inhibition of 67LR formation by cerulenin indicates that acylation is involved in the processing of the receptor. It is likely a palmitoylation process, as indicated by sensitivity of NP-40-soluble extracts to hydroxylamine treatment. Immunoblotting assays performed with a polyclonal serum directed against galectin3 showed that both the 67- and the 120-kDa proteins carry galectin3 epitopes whereas the 37LRP does not. These data suggest that the 67LR is a heterodimer stabilized by strong intramolecular hydrophobic interactions, carried by fatty acids bound to the 37LRP and to a galectin3 cross-reacting molecule. J. Cell. Biochem. 69:244–251, 1998. © 1998 Wiley-Liss, Inc.     

Modeling the interaction between anti-cancer drug penicillamine and pristine and functionalized carbon nanotubes for medical applications: density functional theory investigation and a molecular dynamics simulation

Abstract

The present study focuses on the prediction and investigation of binding properties of penicillamine with pure (5,5) single-walled carbon nanotube (SWCNT) and functionalized SWCNT (f-SWCNT) through the B3LYP and M06-2X functionals using the 6-31G** basis set. The electronic and structural properties, adsorption energy and frontier molecular orbitals of various configurations are examined. Our theoretical results indicated that the interaction of the nanotubes with penicillamine molecule is weak so that the drug adsorption process is typically physisorption. Also, results of theoretical calculations show that the adsorption of the drug molecule on f-SWCNT is stronger with shorter intermolecular distances in comparison to pure SWCNT. The natural bond orbital (NBO) analysis of studied systems demonstrates that the charge is transferred from penicillamine molecule to the nanotubes. Furthermore, molecular dynamics (MD) simulation is employed to evaluate the dynamic and diffusion behavior of drug molecule on SWCNT and f-SWCNT. Energy results show that drug molecule spontaneously moves toward the carriers, and the van der Waals energy contributions in drug adsorption are more than electrostatic interaction. The obtained results from MD simulation confirm that the functionalization of SWCNT leads to increase in the solubility of the carrier in aqueous solution.

Communicated by Ramaswamy H. Sarma     

A structural study of the synergic envelopment of acetonitrile by a Cd(II) activated molecular receptor formed from cyclen with appended 2-hydroxy-3-phenylpropyl moieties

The newly synthesised metal ion activated molecular receptor [Cd{1,4,7,10-tetrakis((R)-(−)-2-hydroxy-3-phenylpropyl)-1,4,7,10-tetraazacyclododecane}](ClO₄)₂ (hereafter [CdL](ClO₄)₂) acts as a molecular receptor for acetonitrile. The receptor was characterised by X-ray crystallography in its metal free form, as its Cd(II) complex with no included molecule, as its Cd(II) complex with an included acetonitrile molecule and, for comparative purposes, as its Zn(II) complex with a partially included acetonitrile molecule. These structural studies demonstrated that the Cd(II) complex is eight-coordinate, with the potential to form a well defined hydrophobic cavity that can contain one acetonitrile molecule through four hydrogen-bonds, whereas the Zn(II) complex is six-coordinate, with a less rigidly constituted binding cavity, such that when solvated by acetonitrile the solvating molecule is retained by only a single hydrogen-bond.     

Molecular modeling of noncompetitive antagonists of the NMDA receptor: proposal of a pharmacophore and a description of the interaction mode

Since the three-dimensional structure of the NMDA receptor has not been determined experimentally, indirect computer-assisted molecular modeling techniques appear to be of great usefulness in the characterization of the common pharmacophore of all NMDA receptor noncompetitive antagonists, despite their structural differences. Indeed, the conformational analysis of three different chemical families (MK801, PCP, dexoxadrol and their analogues), has allowed us to visualize the different conformations and configurations of each molecule. Superimposition with configurations 1 and 2 of the MK801 molecule has allowed us to propose active conformations and thereafter a geometrical characterization of the pharmacophore, especially the determination of the orientation of the nitrogen lone pair (NLP) related to the phenyl. On the other hand, electrostatic studies, combined with geometrical features, have allowed us to schematize the interaction mode of an active conformation to the binding site. Finally, studies of the molecular lipophilic potential (MLP) have provided us information on the position of lipophilic and hydrophilic zones of the pharmacophore.     

All-Atom Structural Models for Complexes of Insulin-Like Growth Factors IGF1 and IGF2 with Their Cognate Receptor

Type 1 insulin-like growth factor receptor (IGF1R) is a membrane-spanning glycoprotein of the insulin receptor family that has been implicated in a variety of cancers. The key questions related to molecular mechanisms governing ligand recognition by IGF1R remain unanswered, partly due to the lack of testable structural models of apo or ligand-bound receptor complexes. Using a homology model of the IGF1R ectodomain IGF1RΔβ, we present the first experimentally consistent all-atom structural models of IGF1/IGF1RΔβ and IGF2/IGF1RΔβ complexes. Our explicit-solvent molecular dynamics (MD) simulation of apo-IGF1RΔβ shows that it displays asymmetric flexibility mechanisms that result in one of two binding pockets accessible to growth factors IGF1 and IGF2, as demonstrated via an MD-assisted Monte Carlo docking procedure. Our MD-generated ensemble of structures of apo and IGF1-bound IGF1RΔβ agrees reasonably well with published small-angle X-ray scattering data. We observe simultaneous contacts of each growth factor with sites 1 and 2 of IGF1R, suggesting cross-linking of receptor subunits. Our models provide direct evidence in favor of suggested electrostatic complementarity between the C-domain (IGF1) and the cysteine-rich domain (IGF1R). Our IGF1/IGF1RΔβ model provides structural bases for the observation that a single IGF1 molecule binds to IGF1RΔβ at low concentrations in small-angle X-ray scattering studies. We also suggest new possible structural bases for differences in the affinities of insulin, IGF1, and IGF2 for their noncognate receptors.