Alternative Conformations Research Articles

Assessing the Conformational Equilibrium of Carboxylic Acid via Quantum Mechanical and Molecular Dynamics Studies on Acetic Acid.

Accurate hydrogen placement in molecular modeling is crucial for studying the interactions and dynamics of biomolecular systems. The carboxyl functional group is a prototypical example of a functional group that requires protonation during structure preparation. To our knowledge, when in their neutral form, carboxylic acids are typically protonated in the syn conformation by default in classical molecular modeling packages, with no consideration of alternative conformations, though we are not aware of any careful examination of this topic. Here, we investigate the general belief that carboxylic acids should always be protonated in the syn conformation. We calculate and compare the relative energetic stabilities of syn and anti acetic acid using ab initio quantum mechanical calculations and atomistic molecular dynamics simulations. We focus on the carboxyl torsional potential and configurations of microhydrated acetic acid from molecular dynamics simulations, probing the effects of solvent, force field (GAFF vs GAFF2), and partial charge assignment of acetic acid. We show that while the syn conformation is the preferred state, the anti state may in some cases also be present under normal NPT conditions in solution.

Fast NMR method to probe solvent accessibility and disordered regions in proteins

Understanding protein structure and dynamics, which govern key cellular processes, is crucial for basic and applied research. Intrinsically disordered protein (IDP) regions display multifunctionality via alternative transient conformations, being key players in disease mechanisms. IDP regions are abundant, namely in small viruses, allowing a large number of functions out of a small proteome. The relation between protein function and structure is thus now seen from a different perspective: as IDP regions enable transient structural arrangements, each conformer can play different roles within the cell. However, as IDP regions are hard and time-consuming to study via classical techniques (optimized for globular proteins with unique conformations), new methods are required. Here, employing the dengue virus (DENV) capsid (C) protein and the immunoglobulin-binding domain of streptococcal protein G, we describe a straightforward NMR method to differentiate the solvent accessibility of single amino acid N-H groups in structured and IDP regions. We also gain insights into DENV C flexible fold region biological activity. The method, based on minimal pH changes, uses the well-established 1H-15N HSQC pulse sequence and is easily implementable in current protein NMR routines. The data generated are simple to interpret, with this rapid approach being an useful first-choice IDPs characterization method.

Alternative conformation induced by substrate binding for Arabidopsis thalianaN6-methyl-AMP deaminase.

Adenosine deaminase is involved in adenosine degradation and salvage pathway, and plays important physiological roles in purine metabolism. Recently, a novel type of adenosine deaminase-like protein has been identified, which displays deamination activity toward N6-methyl-adenosine monophosphate but not adenosine or AMP, and was consequently named N6-methyl-AMP deaminase (MAPDA). The underlying structural basis of MAPDA recognition and catalysis is poorly understood. Here, we present the crystal structures of MAPDA from Arabidopsis thaliana in the free and in the ligand-bound forms. The protein contains a conserved (β/α)8 Tim-barrel domain and a typical zinc-binding site, but it also exhibits idiosyncratic local differences for two flexible helices important for substrate binding. The extensive interactions between the N6-methyl-AMP substrate or the inosine monophosphate product and the enzyme were identified, and subsequently evaluated by the deamination activity assays. Importantly, each structure reported here represents a different stage of the catalytic pathway and their structural differences suggested that the enzyme can exist in two distinct conformational states. The open state switches to the closed one upon the binding of ligands, brought about by the two critical helices. Our structural studies provide the first look of this important metabolic enzyme and shed lights on its catalytic pathway, which holds promise for the structure-based drug design for MAPDA-related diseases.

Abstract B150: STING acquired species-specific motifs to control alternative immune responses

Abstract Stimulator of Interferon Genes (STING) is an adaptor protein critical for downstream signaling during recognition of mislocalized cytosolic DNA. STING activation requires recognition of a unique second messenger 2′–5′, 3′–5′ cyclic GMP–AMP (2′3′ cGAMP) synthesized by the pathogen recognition receptor cyclic GMP–AMP synthase (cGAS). Due to its broad expression in most tissues, and ability to respond to different endogenous and foreign small molecules, STING is emerging as a promising drug target for cancer immunotherapy and treatment of autoimmune diseases. However, it remains unknown how the single adaptor protein STING controls distinct transcriptional outputs leading to production of type I interferon, proinflammatory cytokine responses or autophagy. Here we present surprising evidence for the ancient origins of the cGAS-STING signaling pathway by discovery of bacterial cGAS-like enzymes and complete functional cGAS-STING pathways in lower metazoans. Our previous structural data demonstrate that human and metazoan STING proteins share a common architecture that couples ligand binding with conformational changes and allows signal activation. However, it remained unclear how alternative STING conformations control different downstream signaling outputs, and which motifs are responsible for recruitment of the downstream factors controlling each pathway. We have now combined a phylogenetic and biochemical approach to explain the conserved elements that regulate STING downstream signaling events. Our analysis identifies new immune pathway-specific regulatory motifs acquired by distinct STING species, and we are currently focused on identifying the recruited factors important for controlling human STING signaling. Together, our results explain the molecular basis for distinct STING downstream signaling, and provide new insights for the rational design of STING pathway-specific therapeutics. Citation Format: Carina C. de Oliveira Mann, David S. King, Philip J. Kranzusch. STING acquired species-specific motifs to control alternative immune responses [abstract]. In: Proceedings of the Fourth CRI-CIMT-EATI-AACR International Cancer Immunotherapy Conference: Translating Science into Survival; Sept 30-Oct 3, 2018; New York, NY. Philadelphia (PA): AACR; Cancer Immunol Res 2019;7(2 Suppl):Abstract nr B150.

RIG-I Selectively Discriminates against 5′-Monophosphate RNA

The innate immune sensor RIG-I must sensitively detect and respond to viral RNAs that enter the cytoplasm, while remaining unresponsive to the abundance of structurally similar RNAs that are the products of host metabolism. In the case of RIG-I, these viral and host targets differ by only a few atoms, and a molecular mechanism for such selective differentiation has remained elusive. Using a combination of quantitative biophysical and immunological studies, we show that RIG-I, which is normally activated by duplex RNAs containing a 5'-tri- or diphosphate (5'-ppp or 5'-pp RNAs), is actively antagonized by RNAs containing 5'-monophosphates (5'-p RNAs). This is accomplished by a gating mechanism in which an alternative RIG-I conformation blocks the C-terminal domain (CTD) upon 5'-p RNA binding, thereby short circuiting the activation of signaling.

Surveying the Sequence Space Landscape of Fold-Switching Proteins

A central tenet of biology is that globular proteins have a unique three-dimensional structure under physiological conditions. Over the past 15 years, however, a growing number of proteins have been shown to switch folds, i.e. remodel their secondary and tertiary structures in response to cellular stimuli, leading to biologically relevant functional changes. Previously, we identified 96 of these fold-switching proteins in the Protein Data Bank (PDB). We also found evidence suggesting that their natural abundance likely exceeds their PDB abundance. Thus, we sought to identify other putative fold switchers with only one solved conformation. To do this, we identified two characteristic features of our ∼100 fold-switching proteins—incorrect secondary structure predictions and likely independent folding cooperativity—and searched the PDB for other proteins with similar features. Reassuringly, this method identified the fold-switching regions of 13/15 proteins, each with one solved structure and experimental support for an alternative conformation (p < 2.0∗10−5, hypergeometric test). We are now using a modified version of this method to estimate the size of sequence space for a family of fold-switching proteins. Specifically, we have mapped out sequence space for an E. coli protein known to switch folds and functions (RfaH) and its paralog that maintains a single fold and function (NusG). Preliminary experimental results suggest that our computational method successfully identifies both fold switchers and non-fold switchers in bacterial phyla distantly related to E. coli. Thus, the fold-switching mechanism appears to be conserved among distantly-related bacterial species, highlighting its biological relevance and suggesting that fold switching could be an evolutionarily robust mechanism.

Interaction of Halictine-Related Antimicrobial Peptides with Membrane Models.

We have investigated structural changes of peptides related to antimicrobial peptide Halictine-1 (HAL-1) induced by interaction with various membrane-mimicking models with the aim to identify a mechanism of the peptide mode of action and to find a correlation between changes of primary/secondary structure and biological activity. Modifications in the HAL-1 amino acid sequence at particular positions, causing an increase of amphipathicity (Arg/Lys exchange), restricted mobility (insertion of Pro) and consequent changes in antimicrobial and hemolytic activity, led to different behavior towards model membranes. Secondary structure changes induced by peptide-membrane interaction were studied by circular dichroism, infrared spectroscopy, and fluorescence spectroscopy. The experimental results were complemented by molecular dynamics calculations. An α-helical structure has been found to be necessary but not completely sufficient for the HAL-1 peptides antimicrobial action. The role of alternative conformations (such as β-sheet, PPII or 310-helix) also seems to be important. A mechanism of the peptide mode of action probably involves formation of peptide assemblies (possibly membrane pores), which disrupt bacterial membrane and, consequently, allow membrane penetration.

Journey to the center of the protein: allostery from multitemperature multiconformer X-ray crystallography.

Proteins inherently fluctuate between conformations to perform functions in the cell. For example, they sample product-binding, transition-state-stabilizing and product-release states during catalysis, and they integrate signals from remote regions of the structure for allosteric regulation. However, there is a lack of understanding of how these dynamic processes occur at the basic atomic level. This gap can be at least partially addressed by combining variable-temperature (instead of traditional cryogenic temperature) X-ray crystallography with algorithms for modeling alternative conformations based on electron-density maps, in an approach called multitemperature multiconformer X-ray crystallography (MMX). Here, the use of MMX to reveal alternative conformations at different sites in a protein structure and to estimate the degree of energetic coupling between them is discussed. These insights can suggest testable hypotheses about allosteric mechanisms. Temperature is an easily manipulated experimental parameter, so the MMX approach is widely applicable to any protein that yields well diffracting crystals. Moreover, the general principles of MMX are extensible to other perturbations such as pH, pressure, ligand concentration etc. Future work will explore strategies for leveraging X-ray data across such perturbation series to more quantitatively measure how different parts of a protein structure are coupled to each other, and the consequences thereof for allostery and other aspects of protein function.

The Rich Solid-State Phase Behavior of l-Phenylalanine: Disappearing Polymorphs and High Temperature Forms.

After years of controversy over the solid state structure of the essential amino acid l-phenylalanine, four different polymorphic forms were published recently. The common form I has symmetry P21 with four molecules in the asymmetric unit (Z′ = 4), similar to form III, but with a different arrangement of molecular bilayers. Form II, obtained from the hydrate at very low humidity, is unrelated to forms I and III, as is the high-density form IV. The present investigation demonstrates that this prototype aromatic amino acid has two additional high-temperature phases Ih and IIIh obtained from form I and form III above 458 and 440 K, respectively, when flipping between two alternative side-chain conformations becomes dynamic and causes pairs of molecules, initially crystallographically independent, to become equivalent above a sharp transition temperature. These abrupt and reversible phase changes occur with a reduction of Z′ from 4 (low T) to 2 (high T) and modified crystal symmetry. We furthermore experienced an example of disappearing polymorph for form I which after growing form III in one of our laboratories could no longer be crystallized at room temperature. In contrast, form III crystals may be irreversibly converted to form I crystals as a result of sliding of molecular bilayers in the crystal at elevated temperature. No conversions between the high-temperature forms Ih and IIIh were found. The remarkable crystallographic results are here corroborated by Molecular Dynamics and metadynamics simulations of the form I – form III system.

Morphology Clustering Software for AFM Images, Based on Particle Isolation and Artificial Neural Networks

Advanced microscopy techniques currently allow scientists to visualize biomolecules at high resolution. Among them, atomic force microscopy (AFM) shows the advantage of imaging molecules in their native state, without requiring any staining or coating of the sample. Biopolymers, including proteins and structured nucleic acids, are flexible molecules that can fold into alternative conformations for any given monomer sequence, as exemplified by the different three-dimensional structures adopted by RNA in solution. Therefore, the manual analysis of images visualized by AFM and other microscopy techniques becomes very laborious and time-consuming (and may also be inadvertently biased) when large populations of biomolecules are studied. Here we present a novel morphology clustering software, based on particle isolation and artificial neural networks, which allows the automatic image analysis and classification of biomolecules that can show alternative conformations. It has been tested with a set of AFM images of RNA molecules (a 574 nucleotides-long functional region of the hepatitis C virus genome that contains its internal ribosome entry site element) structured in folding buffers containing 0, 2, 4, 6 or 10 mM Mg 2+ . The developed software shows a broad applicability in the microscopy-based analysis of biopolymers and other complex biomolecules.

The potential of G-quadruplexes as a therapeutic target for neurological diseases

The most common form of DNA is a right-handed helix, the B-form DNA. DNA can also adopt a variety of alternative conformations, termed non-B-form DNA secondary structures, including the G-quadruplex (G4). Furthermore, non-canonical RNA G4 secondary structures are also observed. Recent bioinformatics analysis revealed genomic positions of G4. In addition, G4 formation may be associated with various biological functions, including DNA replication, transcription, epigenetic modification, and RNA metabolism. In this review, we focus on G4 structures in neuronal functions, which may have important roles reveal mechanisms underlying neurological disorders. In addition, we discuss the potential of G4s as a therapeutic target for neurological diseases.

Structural basis for effects of CpA modifications on C/EBPβ binding of DNA.

CCAAT/enhancer binding proteins (C/EBPs) regulate gene expression in a variety of cells/tissues/organs, during a range of developmental stages, under both physiological and pathological conditions. C/EBP-related transcription factors have a consensus binding specificity of 5′-TTG-CG-CAA-3′, with a central CpG/CpG and two outer CpA/TpG dinucleotides. Methylation of the CpG and CpA sites generates a DNA element with every pyrimidine having a methyl group in the 5-carbon position (thymine or 5-methylcytosine (5mC)). To understand the effects of both CpG and CpA modification on a centrally-important transcription factor, we show that C/EBPβ binds the methylated 8-bp element with modestly-increased (2.4-fold) binding affinity relative to the unmodified cognate sequence, while cytosine hydroxymethylation (particularly at the CpA sites) substantially decreased binding affinity (36-fold). The structure of C/EBPβ DNA binding domain in complex with methylated DNA revealed that the methyl groups of the 5mCpA/TpG make van der Waals contacts with Val285 in C/EBPβ. Arg289 recognizes the central 5mCpG by forming a methyl-Arg-G triad, and its conformation is constrained by Val285 and the 5mCpG methyl group. We substituted Val285 with Ala (V285A) in an Ala-Val dipeptide, to mimic the conserved Ala-Ala in many members of the basic leucine-zipper family of transcription factors, important in gene regulation, cell proliferation and oncogenesis. The V285A variant demonstrated a 90-fold binding preference for methylated DNA (particularly 5mCpA methylation) over the unmodified sequence. The smaller side chain of Ala285 permits Arg289 to adopt two alternative conformations, to interact in a similar fashion with either the central 5mCpG or the TpG of the opposite strand. Significantly, the best-studied cis-regulatory elements in RNA polymerase II promoters and enhancers have variable sequences corresponding to the central CpG or reduced to a single G:C base pair, but retain a conserved outer CpA sequence. Our analyses suggest an important modification-dependent CpA recognition by basic leucine-zipper transcription factors.

TDP-43 extracted from frontotemporal lobar degeneration subject brains displays distinct aggregate assemblies and neurotoxic effects reflecting disease progression rates.

Accumulation of abnormally phosphorylated TDP-43 (pTDP-43) is the main pathology in affected neurons of people with amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD). Morphological diversity and neuroanatomical distribution of pTDP-43 accumulations allowed classification of FTLD cases into at least four subtypes, which are correlated with clinical presentations and genetic causes. To understand the molecular basis of this heterogeneity, we developed SarkoSpin, a new method for biochemical isolation of pathological TDP-43. By combining SarkoSpin with mass spectrometry, we revealed proteins beyond TDP-43 that become abnormally insoluble in a disease subtype-specific manner. We show that pTDP-43 extracted from brain forms stable assemblies of distinct densities and morphologies that are associated with disease subtypes. Importantly, biochemically extracted pTDP-43 assemblies showed differential neurotoxicity and seeding that were correlated with disease duration of FTLD subjects. Our data are consistent with the notion that disease heterogeneity could originate from alternate pathological TDP-43 conformations, which are reminiscent of prion strains.

Cysteine Disulfide Traps Reveal Distinct Conformational Ensembles in Dengue Virus NS2B-NS3 Protease

The dengue virus protease (NS2B-NS3pro) plays a critical role in the dengue viral life cycle, making it an attractive drug target for dengue-related pathologies, including dengue hemorrhagic fever. A number of studies indicate that NS2B-NS3pro undergoes a transition between two widely different conformational states: an "open" (inactive) conformation and a "closed" (active) conformation. For the past several years, the equilibrium between these states and the resting conformation of NS2B-NS3pro have been debated, although a strong consensus is emerging. To investigate the importance of such conformational states, we developed versions of NS2B-NS3pro that allow us to trap the enzyme in various distinct conformations. Our data from these variants suggest that the enzymatic activity appears to be dependent on the movement of NS2B and may rely on the flexibility of the protease core. Locking the enzyme into the "closed" conformation dramatically increased activity, strongly suggesting that the "closed" conformation is the active conformation. The observed resting state of the enzyme depends largely on the construct used to express the NS2B-NS3pro complex. In an "unlinked" construct, in which the NS2B and NS3 regions exist as independent, co-expressed polypeptides, the enzyme rests predominantly in a "closed", active conformation. In contrast, in a "linked" construct, in which NS2B and NS3 are attached by a nine-amino acid linker, NS2B-NS3pro adopts a more relaxed, alternative conformation. Nevertheless, even the unlinked construct samples both the "closed" and other alternative conformations. Given our findings, and the more realistic resemblance of NS2B-NS3pro to the native enzyme, these data strongly suggest that studies should focus on the "unlinked" constructs moving forward. Additionally, the results from these studies provide a more detailed understanding of the various poses of the dengue virus NS2B-NS3 protease and should help guide future drug discovery efforts aimed at this enzyme.

QFit-ligand Reveals Widespread Conformational Heterogeneity of Drug-Like Molecules in X-Ray Electron Density Maps.

Proteins and ligands sample a conformational ensemble that governs molecular recognition, activity, and dissociation. In structure-based drug design, access to this conformational ensemble is critical to understand the balance between entropy and enthalpy in lead optimization. However, ligand conformational heterogeneity is currently severely underreported in crystal structures in the Protein Data Bank, owing in part to a lack of automated and unbiased procedures to model an ensemble of protein-ligand states into X-ray data. Here, we designed a computational method, qFit-ligand, to automatically resolve conformationally averaged ligand heterogeneity in crystal structures, and applied it to a large set of protein receptor-ligand complexes. In an analysis of the cancer related BRD4 domain, we found that up to 29% of protein crystal structures bound with drug-like molecules present evidence of unmodeled, averaged, relatively isoenergetic conformations in ligand-receptor interactions. In many retrospective cases, these alternate conformations were adventitiously exploited to guide compound design, resulting in improved potency or selectivity. Combining qFit-ligand with high-throughput screening or multitemperature crystallography could therefore augment the structure-based drug design toolbox.

Hydrolysis of the non-canonical cyclic nucleotide cUMP by PDE9A: kinetics and binding mode.

The non-canonical cyclic nucleotide cUMP and the phosphodiesterase PDE9A both occur in neuronal cells. Using HPLC-coupled tandem mass spectrometry, we characterized the kinetics of PDE9A-mediated cUMP hydrolysis. PDE9A is a low-affinity and high-velocity enzyme for cUMP (Vmax = ~ 6μmol/min/mg; Km = ~ 401μM). The PDE9 inhibitor BAY 73-6691 inhibited PDE9A-catalyzed cUMP hydrolysis (Ki = 590nM). Docking studies indicate two H-bonds between the cUMP uridine moiety and Gln453/Asn405 of PDE9A. By contrast, the guanosine moiety of cGMP forms three H-bonds with Gln453. cCMP is not hydrolyzed at a concentration of 3μM, but inhibits the PDE9A-catalyzed cUMP hydrolysis at concentrations of 100μM or more. The probable main reason is that the cytosine moiety cannot act as H-bond acceptor for Gln453. A comparison of PDE9A with PDE7A suggests that the preference of the former for cGMP and cUMP and of the latter for cAMP and cCMP is due to stabilized alternative conformations of the side chain amide of Gln453 and Gln413, respectively. This so-called glutamine switch is known to be involved in the regulation of cAMP/cGMP selectivity of some PDEs.

Structural Variability of EspG Chaperones from Mycobacterial ESX-1, ESX-3, and ESX-5 Type VII Secretion Systems

Type VII secretion systems (ESX) are responsible for transport of multiple proteins in mycobacteria. How different ESX systems achieve specific secretion of cognate substrates remains elusive. In the ESX systems, the cytoplasmic chaperone EspG forms complexes with heterodimeric PE–PPE substrates that are secreted from the cells or remain associated with the cell surface. Here we report the crystal structure of the EspG1 chaperone from the ESX-1 system determined using a fusion strategy with T4 lysozyme. EspG1 adopts a quasi 2-fold symmetric structure that consists of a central β-sheet and two α-helical bundles. In addition, we describe the structures of EspG3 chaperones from four different crystal forms. Alternate conformations of the putative PE–PPE binding site are revealed by comparison of the available EspG3 structures. Analysis of EspG1, EspG3, and EspG5 chaperones using small-angle X-ray scattering reveals that EspG1 and EspG3 chaperones form dimers in solution, which we observed in several of our crystal forms. Finally, we propose a model of the ESX-3 specific EspG3–PE5–PPE4 complex based on the small-angle X-ray scattering analysis.

208 - Resurrecting an antibody against an alternative conformation of cytochrome c: the recombinant mAb 1D3

Cytochrome c can undergo reversible conformational changes under biologically-relevant conditions including oxidative post-translational modifications, interactions with phospholipids and electric fields. Revealing the existence and relevance of these alternative cytochrome c proteoforms at the cell and tissue level is challenging and has benefitted from immunochemical approaches. Indeed, a monoclonal antibody that recognizes a major alternative conformation of cytochrome c, 1D3, was previously developed and characterized. The alternative conformation involves the destabilization of the Ω-loops encompassing residues 70 to 85 and 40 to 57 and the displacement of the Met80 heme ligand. However, the 1D3 producing hybridoma was lost and several attempts to reproduce it failed. Thus, with the remaining amounts of mAb 1D3 we aimed to “resurrect” 1D3 by the generation of a recombinant version of it. We employed the Valens algorithm for proteomic sequencing of 1D3. Briefly, the antibody was digested with multiple proteases and analyzed via LC-MS/MS using multiple fragmentation modes. By combining de novo peptide sequencing and homology mapping, we were able to obtain full-length sequences for the heavy and light chains of 1D3. Gene fragments representing the predicted heavy and light chain variable domains were synthesized using gBLOCK technology (IDT Inc., http://www.idtdna.com/). These were cloned into the heavy and light chain frameworks, respectively, using standard molecular biology techniques. Expression of light and heavy chains of 1D3 and its molecular assembly in HEK293 and CHO cells resulted in a functional antibody as tested by competitive ELISA with, for example, nitrotyrosine-containing cytochrome c species (NO274Tyr-Cyt c), known to adopt an alternative conformation. The recombinant 1D3 has the same specificity and affinity as the “original” 1D3 and is being successfully applied to various cellular models; enhanced levels of the 1D3-detectable alternative conformation of cytochrome c and their nuclear translocation are observed under conditions of altered cell homeostasis.

Subatomic structure of hyper-sweet thaumatin D21N mutant reveals the importance of flexible conformations for enhanced sweetness

One of the sweetest proteins found in tropical fruits (with a threshold of 50 nM), thaumatin, is also used commercially as a sweetener. Our previous study successfully produced the sweetest thaumatin mutant (D21N), designated hyper-sweet thaumatin, which decreases the sweetness threshold to 31 nM. To investigate why the D21N mutant is sweeter than wild-type thaumatin, we compared the structure of the D21N mutant solved at a subatomic resolution of 0.93 Å with that of wild-type thaumatin determined at 0.90 Å. Although the overall structure of the D21N mutant resembles that of wild-type thaumatin, our subatomic resolution analysis successfully assigned and discriminated the detailed atomic positions of side-chains at position 21. The relative B-factor value of the side-chain at position 21 in the D21N mutant was higher than that of wild-type thaumatin, hinting at a greater flexibility of side-chain at 21 in the hyper-sweet D21N mutant. Furthermore, alternative conformations of Lys19, which is hydrogen-bonded to Asp21 in wild-type, were found only in the D21N mutant. Subatomic resolution analysis revealed that flexible conformations at the sites adjacent to positions 19 and 21 play a crucial role in enhancing sweet potency and may serve to enhance the complementarity of electrostatic potentials for interaction with the sweet taste receptor.

TANGO: A high through-put conformation generation and semiempirical method-based optimization tool for ligand molecules.

Lead optimization is one of the crucial steps in the drug discovery pipeline. After identifying the lead molecule and obtaining its 2D geometry, understanding the best conformation it would attain in 3D still remains one of the most challenging steps in drug discovery. There have been multiple methods and algorithms that are directed toward achieving best conformation for the lead molecules. TANGO focuses on conformation generation and its optimization using semiempirical energy calculations. The conformation generation is based on torsion angle rotation of the exocyclic bonds. The energy calculations are performed using MOPAC. The unique feature of this tool lies in the implementation of Message Passing Interface (MPI) for conformation generation and semiempirical-based optimization. A well-defined architecture handling the input and output generation has been used. The master and slave approach to handle operations involved in torsion angle rotation and energy calculations has helped in load balancing the process of conformation generation. The benchmarking results suggest that TANGO scales significantly well across eight nodes with each node utilizing 16 cores. This tool may prove to very useful in high throughput generation of semiempirically optimized small molecule conformations. The use of semiempirical methods for optimization generates a conformational ensemble thereby helping to obtain stable and alternate stable conformers for a given ligand molecule. © 2018 Wiley Periodicals, Inc.