Comparative Modelling of Organic Anion Transporting Polypeptides: Structural Insights and Comparison of Binding Modes.

To better understand the functionality of organic anion transporting polypeptides (OATPs) and to design new ligands, reliable structural data of each OATP is needed. In this work, we used a combination of homology model with molecular dynamics simulations to generate a comprehensive structural dataset, that encompasses a diverse set of OATPs but also their relevant conformations. Our OATP models share a conserved transmembrane helix folding harbouring a druggable binding pocket in the shape of an inner pore. Our simulations suggest that the conserved salt bridges at the extracellular region between residues on TM1 and TM7 might influence the entrance of substrates. Interactions between residues on TM1 and TM4 within OATP1 family shown their importance in transport of substrates. Additionally, in transmembrane (TM) 1/2, a known conserved element, interact with two identified motifs in the TM7 and TM11. Our simulations suggest that TM1/2-TM7 interaction influence the inner pocket accessibility, while TM1/2-TM11 salt bridges control the substrate binding stability.

Interruption of Vif/Elongin C interaction: In silico and experimental elucidation of the underlying molecular mechanism of benzimidazole-based APOBEC3G stabilizers.

In 2014, two novel and promising benzimidazole-based APOBEC3G stabilizers MM-1 and MM-2 (MMs) were uncovered with an elusive mechanism of action. Vif-APOBEC3G axis has been recognized as a novel therapeutic target for anti HIV-1 drug development. The unexplored mechanism of MMs hindered their further development into lead compounds. To recognize their underlying mechanism we adopted an exhaustive in silico workflow by which we tested their ability to interrupt Vif complex network formation. The preliminary outcome guided us to a high likelihood of MMs interaction within Elongin C binding site, which in turn, perturbs Vif/Elongin C binding and ultimately undermines Vif action. To validate our estimation, we synthesized only MM-1 as a model to complement our study by in vitro assay for a real-time understanding. An immunoprecipitation experiment confirmed the capacity of MM-1 to interrupt Vif/Elongin C interaction. This is an integral study that lies at the interface between theoretical and experimental approaches showing the potential of molecular modelling to address issues related to drug development.

What makes it difficult to refine protein models further via molecular dynamics simulations?

Protein structure refinement remains a challenging yet important problem as it has the potential to bring already accurate template-based models to near-native resolution. Refinement based on molecular dynamics simulations has been a highly promising approach and the performance of MD-based refinement in the Feig group during CASP12 is described here. During CASP12, sampling was extended well into the microsecond scale, an improved force field was applied, and new protocol variations were tested. Progress over previous rounds of CASP was found to be limited which is analyzed in terms of the quality of the initial models and dependency on the amount of sampling and refinement protocol variations. As current MD-based refinement protocols appear to be reaching a plateau, detailed analysis is presented to provide new insight into the major challenges towards more extensive structure refinement, focusing in particular on sampling with and without restraints.

Structure refinement of membrane proteins via molecular dynamics simulations.

A refinement protocol based on physics-based techniques established for water soluble proteins is tested for membrane protein structures. Initial structures were generated by homology modeling and sampled via molecular dynamics simulations in explicit lipid bilayer and aqueous solvent systems. Snapshots from the simulations were selected based on scoring with either knowledge-based or implicit membrane-based scoring functions and averaged to obtain refined models. The protocol resulted in consistent and significant refinement of the membrane protein structures similar to the performance of refinement methods for soluble proteins. Refinement success was similar between sampling in the presence of lipid bilayers and aqueous solvent but the presence of lipid bilayers may benefit the improvement of lipid-facing residues. Scoring with knowledge-based functions (DFIRE and RWplus) was found to be as good as scoring using implicit membrane-based scoring functions suggesting that differences in internal packing is more important than orientations relative to the membrane during the refinement of membrane protein homology models.

Gating of Connexin Channels by transjunctional-voltage: Conformations and models of open and closed states.

Voltage is an important physiologic regulator of channels formed by the connexin gene family. Connexins are unique among ion channels in that both plasma membrane inserted hemichannels (undocked hemichannels) and intercellular channels (aggregates of which form gap junctions) have important physiological roles. The hemichannel is the fundamental unit of gap junction voltage-gating. Each hemichannel displays two distinct voltage-gating mechanisms that are primarily sensitive to a voltage gradient formed along the length of the channel pore (the transjunctional voltage) rather than sensitivity to the absolute membrane potential (Vm or Vi-o). These transjunctional voltage dependent processes have been termed Vj- or fast-gating and loop- or slow-gating. Understanding the mechanism of voltage-gating, defined as the sequence of voltage-driven transitions that connect open and closed states, first and foremost requires atomic resolution models of the end states. Although ion channels formed by connexins were among the first to be characterized structurally by electron microscopy and x-ray diffraction in the early 1980's, subsequent progress has been slow. Much of the current understanding of the structure-function relations of connexin channels is based on two crystal structures of Cx26 gap junction channels. Refinement of crystal structure by all-atom molecular dynamics and incorporation of charge changing protein modifications has resulted in an atomic model of the open state that arguably corresponds to the physiologic open state. Obtaining validated atomic models of voltage-dependent closed states is more challenging, as there are currently no methods to solve protein structure while a stable voltage gradient is applied across the length of an oriented channel. It is widely believed that the best approach to solve the atomic structure of a voltage-gated closed ion channel is to apply different but complementary experimental and computational methods and to use the resulting information to derive a consensus atomic structure that is then subjected to rigorous validation. In this paper, we summarize our efforts to obtain and validate atomic models of the open and voltage-driven closed states of undocked connexin hemichannels. This article is part of a Special Issue entitled: Gap Junction Proteins edited by Jean Claude Herve.

Structural insights into human microsomal epoxide hydrolase by combined homology modeling, molecular dynamics simulations, and molecular docking calculations.

A new homology model of human microsomal epoxide hydrolase was derived based on multiple templates. The model obtained was fully evaluated, including MD simulations and ensemble-based docking, showing that the quality of the structure is better than that of only previously known model. Particularly, a catalytic triad was clearly identified, in agreement with the experimental information available. Analysis of intermediates in the enzymatic mechanism led to the identification of key residues for substrate binding, stereoselectivity, and intermediate stabilization during the reaction. In particular, we have confirmed the role of the oxyanion hole and the conserved motif (HGXP) in epoxide hydrolases, in excellent agreement with known experimental and computational data on similar systems. The model obtained is the first one that fully agrees with all the experimental observations on the system. Proteins 2017; 85:720-730. © 2016 Wiley Periodicals, Inc.

Escape from adamantane: Scaffold optimization of novel P2X7 antagonists featuring complex polycycles.

The adamantane scaffold, despite being widely used in medicinal chemistry, is not devoid of problems. In recent years we have developed new polycyclic scaffolds as surrogates of the adamantane group with encouraging results in multiple targets. As an adamantane scaffold is a common structural feature in several P2X7 receptor antagonists, herein we report the synthesis and pharmacological evaluation of multiple replacement options of adamantane that maintain a good activity profile. Molecular modeling studies support the binding of the compounds to a site close to the central pore, rather than to the ATP-binding site and shed light on the structural requirements for novel P2X7 antagonists.

New structural insights into the apelin receptor: identification of key residues for apelin binding.

Apelin is the endogenous ligand of the orphan 7-transmembrane domain GPCR APJ, now named the apelin receptor (ApelinR). Apelin plays a prominent role in body fluid and cardiovascular homeostasis. To better understand the structural organization of the ApelinR, we built 3 homology 3-dimensional (3D) models of the human ApelinR using the validated cholecystokinin receptor-1 3D model or the X-ray structures of the ß2-adrenergic and CXCR4 receptors as templates. Docking of the pyroglutamyl form of apelin 13 (pE13F) into these models revealed the conservation at the bottom of the binding site of a hydrophobic cavity in which the C-terminal Phe of pE13F was embedded. In contrast, at the top of the binding site, depending on the model, different interactions were visualized between acidic residues of the ApelinR and the basic residues of pE13F. Using site-directed mutagenesis, we showed that Asp 92, Glu 172, and Asp 282 of rat ApelinR are key residues in apelin binding by interacting with Lys 8, Arg 2, and Arg 4 of pE13F, respectively. These residues are only seen in the CXCR4-based ApelinR 3D model, further validating this model. These findings bring new insights into the structural organization of the ApelinR and the mode of apelin binding.

The hypertrophic cardiomyopathy myosin mutation R453C alters ATP binding and hydrolysis of human cardiac ß-myosin.

The human hypertrophic cardiomyopathy mutation R453C results in one of the more severe forms of the myopathy. Arg-453 is found in a conserved surface loop of the upper 50-kDa domain of the myosin motor domain and lies between the nucleotide binding pocket and the actin binding site. It connects to the cardiomyopathy loop via a long α-helix, helix O, and to Switch-2 via the fifth strand of the central ß-sheet. The mutation is, therefore, in a position to perturb a wide range of myosin molecular activities. We report here the first detailed biochemical kinetic analysis of the motor domain of the human ß-cardiac myosin carrying the R453C mutation. A recent report of the same mutation (Sommese, R. F., Sung, J., Nag, S., Sutton, S., Deacon, J. C., Choe, E., Leinwand, L. A., Ruppel, K., and Spudich, J. A. (2013) Proc. Natl. Acad. Sci. U.S.A. 110, 12607-12612) found reduced ATPase and in vitro motility but increased force production using an optical trap. Surprisingly, our results show that the mutation alters few biochemical kinetic parameters significantly. The exceptions are the rate constants for ATP binding to the motor domain (reduced by 35%) and the ATP hydrolysis step/recovery stroke (slowed 3-fold), which could be the rate-limiting step for the ATPase cycle. Effects of the mutation on the recovery stroke are consistent with a perturbation of Switch-2 closure, which is required for the recovery stroke and the subsequent ATP hydrolysis.

Statin Lactonization by Uridine 5'-Diphospho-glucuronosyltransferases (UGTs).

Statins are cholesterol-lowering drugs that have proven to be effective in lowering the risk of major cardiovascular events. Although well tolerated, statin-induced myopathies are the most common side effects. Compared to their pharmacologically active acid form, statin lactones are more potent inducers of toxicity. They can be formed by glucuronidation mediated by uridine 5'-diphospho-glucuronosyltransferases (UGTs), but a systematic characterization of subtype specificity and kinetics of lactonization is lacking. Here, we demonstrate for six clinically relevant statins that only UGT1A1, 1A3, and 2B7 contribute significantly to their lactonization. UGT1A3 appeared to have the highest lactonization capacity with marked differences in statin conversion rates: pitavastatin ≫ atorvastatin > cerivastatin > lovastatin > rosuvastatin (simvastatin not converted). Using in silico modeling we could identify a probable statin interaction region in the UGT binding pocket. Polymorphisms in these regions of UGT1A1, 1A3, and 2B7 may be a contributing factor in statin-induced myopathies, which could be used in personalization of statin therapy with improved safety.

Synthesis, cholinesterase inhibition and molecular modelling studies of coumarin linked thiourea derivatives.

Alzheimer's disease is among the most widespread neurodegenerative disorder. Cholinesterases (ChEs) play an indispensable role in the control of cholinergic transmission and thus the acetylcholine level in the brain is enhanced by inhibition of ChEs. Coumarin linked thiourea derivatives were designed, synthesized and evaluated biologically in order to determine their inhibitory activity against acetylcholinesterases (AChE) and butyrylcholinesterases (BChE). The synthesized derivatives of coumarin linked thiourea compounds showed potential inhibitory activity against AChE and BChE. Among all the synthesized compounds, 1-(2-Oxo-2H-chromene-3-carbonyl)-3-(3-chlorophenyl)thiourea (2e) was the most potent inhibitor against AChE with an IC50 value of 0.04±0.01µM, while 1-(2-Oxo-2H-chromene-3-carbonyl)-3-(2-methoxyphenyl)thiourea (2b) showed the most potent inhibitory activity with an IC50 value of 0.06±0.02µM against BChE. Molecular docking simulations were performed using the homology models of both cholinesterases in order to explore the probable binding modes of inhibitors. Results showed that the novel synthesized coumarin linked thiourea derivatives are potential candidates to develop for potent and efficacious acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) inhibitors.

Genetic diversity and bioinformatic analysis in the L1 gene of HPV genotypes 31, 33, and 58 circulating in women with normal cervical cytology.

BACKGROUND: HPV-31, -33, and -58, along with HPV-45 and -52, account for almost 11% of HPV-associated cancers. Our previous studies showed that after HPV-16 and -51, HPV-58 was common and HPV-31 was as frequent as HPV-18 among Iranian women with normal cytology. Hence, in this study, we aimed to investigate the intra-type variations in L1 genes of HPV-58, -31, and -33 to find the predominant lineages circulating in women with normal cytology. METHODS: Complete coding sequencing of the L1 gene was amplified and nucleotide and amino acid sequences were compared to those of the references. The selective pressure on L1 protein and whether the variations of the L1 genes embed in L1 loops, or N-glycosylated sites were also investigated. RESULTS: B1, A, and A1 (sub)lineages were common in the HPV-58, -33, and -31 samples, respectively. Ninety nucleotide mutations were observed. Twenty nine nucleotide changes corresponded to nonsynonymous substitutions in which seventeen mutations were located in L1 loops. Only one codon position in HPV-58 sequences was found as the positive selection. No difference was observed in N-glycosylation sites between reference and understudied amino acid sequences. CONCLUSION: In the current study, we reported, for the first time, the (sub) lineages, amino acid, and genetic diversity in the L1 gene of circulating HPV-58, -33, and -31, in women with normal cytology, in Iran. Such studies can not only have epidemiological values, but also aid to set vaccination programs.

ClusPro LigTBM: Automated Template-based Small Molecule Docking.

The template-based approach has been essential for achieving high-quality models in the recent rounds of blind protein-protein docking competition CAPRI (Critical Assessment of Predicted Interactions). However, few such automated methods exist for protein-small molecule docking. In this paper, we present an algorithm for template-based docking of small molecules. It searches for known complexes with ligands that have partial coverage of the target ligand, performs conformational sampling and template-guided energy refinement to produce a variety of possible poses, and then scores the refined poses. The algorithm is available as the automated ClusPro LigTBM server. It allows the user to specify the target protein as a PDB file and the ligand as a SMILES string. The server then searches for templates and uses them for docking, presenting the user with top-scoring poses and their confidence scores. The method is tested on the Astex Diverse benchmark, as well as on the targets from the last round of the D3R (Drug Design Data Resource) Grand Challenge. The server is publicly available as part of the ClusPro docking server suite at https://ligtbm.cluspro.org/.

VGLUT substrates and inhibitors: A computational viewpoint.

The vesicular glutamate transporters (VGLUTs) bind and move glutamate (Glu) from the cytosol into the lumen of synaptic vesicles using a H+-electrochemical gradient (ΔpH and Δψ) generated by the vesicular H+-ATPase. VGLUTs show very low Glu binding and to date, no three-dimensional structure has been elucidated. Prior studies have attempted to identify the key residues involved in binding VGLUT substrates and inhibitors using homology models and docking experiments. Recently, the inward and outward oriented crystal structures of d-galactonate transporter (DgoT) emerged as possible structure templates for VGLUT. In this review, a new homology model for VGLUT2 based on DgoT has been developed and used to conduct docking experiments to identify and differentiate residues and binding orientations involved in ligand interactions. This review describes small molecule-ligand interactions including docking using a VGLUT2 homology model derived from DgoT.

Allosteric Inhibition of Adenylyl Cyclase Type 5 by G-Protein: A Molecular Dynamics Study.

Adenylyl cyclases (ACs) have a crucial role in many signal transduction pathways, in particular in the intricate control of cyclic AMP (cAMP) generation from adenosine triphosphate (ATP). Using homology models developed from existing structural data and docking experiments, we have carried out all-atom, microsecond-scale molecular dynamics simulations on the AC5 isoform of adenylyl cyclase bound to the inhibitory G-protein subunit Gαi in the presence and in the absence of ATP. The results show that Gαi has significant effects on the structure and flexibility of adenylyl cyclase, as observed earlier for the binding of ATP and Gsα. New data on Gαi bound to the C1 domain of AC5 help explain how Gαi inhibits enzyme activity and obtain insight on its regulation. Simulations also suggest a crucial role of ATP in the regulation of the stimulation and inhibition of AC5.

Comparative Structural Analysis of Different Mycobacteriophage-Derived Mycolylarabinogalactan Esterases (Lysin B).

Mycobacteriophage endolysins have emerged as a potential alternative to the current antimycobacterial agents. This study focuses on mycolylarabinogalactan hydrolase (LysB) enzymes of the α/ß-hydrolase family, which disrupt the unique mycolic acid layer of mycobacterium cell wall. Multiple sequence alignment and structural analysis studies showed LysB-D29, the only enzyme with a solved three-dimensional structure, to share several common features with esterases (lacking lid domain) and lipases (acting on long chain lipids). Sequence and structural comparisons of 30 LysB homology models showed great variation in domain organizations and total protein length with major differences in the loop-5 motif harboring the catalytic histidine residue. Docking of different p-nitrophenyl ligands (C4-C18) to LysB-3D models revealed that the differences in length and residues of loop-5 contributed towards wide diversity of active site conformations (long tunnels, deep and superficial funnels, shallow bowls, and a narrow buried cave) resembling that of lipases, cutinases, and esterases. A set of seven LysB enzymes were recombinantly produced; their activity against p-nitrophenyl esters could be related to their active site conformation and acyl binding site. LysB-D29 (long tunnel) showed the highest activity with long chain p-nitrophenyl palmitate followed by LysB-Omega (shallow bowl) and LysB-Saal (deep funnel).

[1]Benzothieno[3,2-d]pyrimidine derivatives as ligands for the serotonergic 5-HT7 receptor.

The present paper describes the synthesis and the binding properties for the serotonergic 5-HT7 and 5-HT1A receptors of three new series (A-C) of (benzo)thienopyrimidinone derivatives. All series exhibit a basic moiety at the 2-position of the heterocyclic scaffold such as N,N-dialkylaminoalkylthio chain in series A and phenylpiperazine, benzylpiperazine, or benzylpiperidine alkyl chain in series B and C. Compounds endowed with the best binding properties at 5-HT7R belong to the B and C types. In particular, derivatives B2 and C1 (RSC4) exhibit notable affinity for 5-HT7R (Ki = 9.08 and 0.85 nM, respectively) and selectivity over the 5-HT1AR (254- and 48-fold, respectively). The structure-affinity relationships for these three new classes of 5-HT7R ligands are discussed and, in order to rationalize and deeply investigate the observed results, molecular modeling studies were performed. In particular, the binding poses of the synthesized compounds were studied by docking them in the two receptor proteins suitably built by homology modeling. The calculated binding energies resulted in an excellent agreement with the experimental measured Ki, further validating the quality of the models.

In Silico Approaches for TRP Channel Modulation.

The implication of several TRP ion channels (e.g., TRPV1) in diverse physiological and pathological processes has signaled them as pivotal drug targets. Consequently, the identification of selective and potent ligands for these channels is of great interest in pharmacology and biomedicine. However, a major challenge in the design of modulators is ensuring the specificity for their intended targets. In recent years, the emergence of high-resolution structures of ion channels facilitates the computer-assisted drug design at molecular levels. Here we describe some computational methods and general protocols to discover channel modulators, including homology modelling, docking and virtual screening, and structure-based peptide design.

GPCR Homology Model Generation for Lead Optimization.

The vast increase of recently solved GPCR X-ray structures forms the basis for GPCR homology modeling to atomistic accuracy. Nowadays, homology models can be employed for GPCR-ligand optimization and have been reported as invaluable tools for drug design in the last few years. Elucidation of the complex GPCR pharmacology and the associated GPCR conformations made clear that different homology models have to be constructed for different activation states of the GPCRs. Therefore, templates have to be chosen accordingly to their sequence homology as well as to their activation state. The subsequent ligand placement is nontrivial, as some recent X-ray structures show very unusual ligand binding sites and solvent involvement, expanding the space of the putative ligand binding site from the generic retinal binding pocket to the whole receptor. In the present study, a workflow is presented starting from the selection of the target sequence, guiding through the GPCR modeling process, and finishing with ligand placement and pose validation.

Natural Polyphenols Selectively Inhibit ß-Carbonic Anhydrase from the Dandruff-Producing Fungus Malassezia globosa: Activity and Modeling Studies.

Around 50 % of the worldwide population is affected by dandruff, which is triggered by a variety of factors. The yeast Malassezia globosa has been labeled as the most probable causative agent for the onset of dandruff. The ß-carbonic anhydrase (CA) of MgCA was recently validated as an anti-dandruff target, with its inhibition being responsible for in vivo growth defects in the fungus. As classical CA inhibitors of the sulfonamide type give rise to permeability problems through biological membranes, finding non-sulfonamide alternatives for MgCA inhibition is of considerable interest in the cosmetic field. We recently screened a large library of human (h) CA inhibitors for MgCA inhibition, including different chemotypes, such as monothiocarbamates, dithiocarbamates, phenols, and benzoxaboroles. Herein, we expanded the research toward new MgCA inhibitors by considering a set of natural polyphenols (including flavones, flavonols, flavanones, flavanols, isoflavones, and depsides) that exhibited MgCA inhibitory activity in the micromolar range, as well as selectivity for the fungal isozyme over off-target human isoforms. The binding mode of representative derivatives within the MgCA catalytic cleft was investigated by docking studies using a homology-built model.