Role of computational methods in pharmaceutical sciences

Experimental and Computational Approaches to Improve Binding Affinity in Chemical Biology and Drug Discovery

Drug discovery is one of the most complicated processes and establishment of a single drug may require multidisciplinary attempts to design efficient and commercially viable drugs. The main purpose of drug design is to identify a chemical compound or inhibitor that can bind to an active site of a specific cavity on a target protein. The traditional drug design methods involved various experimental based approaches including random screening of chemicals found in nature or can be synthesized directly in chemical laboratories. Except for the long cycle design and time, high cost is also the major issue of concern. Modernized computer-based algorithm including structure-based drug design has accelerated the drug design and discovery process adequately. Surprisingly from the past decade remarkable progress has been made concerned with all area of drug design and discovery. CADD (Computer Aided Drug Designing) based tools shorten the conventional cycle size and also generate chemically more stable and worthy compounds and hence reduce the drug discovery cost. This special edition of editorial comprises the combination of seven research and review articles set emphasis especially on the computational approaches along with the experimental approaches using a chemical synthesizing for the binding affinity in chemical biology and discovery as a salient used in de-novo drug designing. This set of articles exfoliates the role that systems biology and the evaluation of ligand affinity in drug design and discovery for the future.

Ordering of Trotterization: Impact on Errors in Quantum Simulation of Electronic Structure

Trotter–Suzuki decompositions are frequently used in the quantum simulation of quantum chemistry. They transform the evolution operator into a form implementable on a quantum device, while incurring an error—the Trotter error. The Trotter error can be made arbitrarily small by increasing the Trotter number. However, this increases the length of the quantum circuits required, which may be impractical. It is therefore desirable to find methods of reducing the Trotter error through alternate means. The Trotter error is dependent on the order in which individual term unitaries are applied. Due to the factorial growth in the number of possible orderings with respect to the number of terms, finding an optimal strategy for ordering Trotter sequences is difficult. In this paper, we propose three ordering strategies, and assess their impact on the Trotter error incurred. Initially, we exhaustively examine the possible orderings for molecular hydrogen in a STO-3G basis. We demonstrate how the optimal ordering scheme depends on the compatibility graph of the Hamiltonian, and show how it varies with increasing bond length. We then use 44 molecular Hamiltonians to evaluate two strategies based on coloring their incompatibility graphs, while considering the properties of the obtained colorings. We find that the Trotter error for most systems involving heavy atoms, using a reference magnitude ordering, is less than 1 kcal/mol. Relative to this, the difference between ordering schemes can be substantial, being approximately on the order of millihartrees. The coloring-based ordering schemes are reasonably promising—particularly for systems involving heavy atoms—however further work is required to increase dependence on the magnitude of terms. Finally, we consider ordering strategies based on the norm of the Trotter error operator, including an iterative method for generating the new error operator terms added upon insertion of a term into an ordered Hamiltonian.

Designing modulators of monoamine transporters using virtual screening techniques

The plasma-membrane monoamine transporters (MATs), including the serotonin (SERT), norepinephrine (NET) and dopamine (DAT) transporters, serve a pivotal role in limiting monoamine-mediated neurotransmission through the reuptake of their respective monoamine neurotransmitters. The transporters are the main target of clinically used psychostimulants and antidepressants. Despite the availability of several potent and selective MAT substrates and inhibitors the continuing need for therapeutic drugs to treat brain disorders involving aberrant monoamine signaling provides a compelling reason to identify novel ways of targeting and modulating the MATs. Designing novel modulators of MAT function have been limited by the lack of three dimensional structure information of the individual MATs. However, crystal structures of LeuT, a bacterial homolog of MATs, in a substrate-bound occluded, substrate-free outward-open, and an apo inward-open state and also with competitive and non-competitive inhibitors have been determined. In addition, several structures of the Drosophila DAT have also been resolved. Together with computational modeling and experimental data gathered over the past decade, these structures have dramatically advanced our understanding of several aspects of SERT, NET, and DAT transporter function, including some of the molecular determinants of ligand interaction at orthosteric substrate and inhibitor binding pockets. In addition progress has been made in the understanding of how allosteric modulation of MAT function can be achieved. Here we will review all the efforts up to date that has been made through computational approaches employing structural models of MATs to design small molecule modulators to the orthosteric and allosteric sites using virtual screening techniques.

Methods to enable the design of bioactive small molecules targeting RNA

RNA is an immensely important target for small molecule therapeutics or chemical probes of function. However, methods that identify, annotate, and optimize RNA-small molecule interactions that could enable the design of compounds that modulate RNA function are in their infancies. This review describes recent approaches that have been developed to understand and optimize RNA motif-small molecule interactions, including structure-activity relationships through sequencing (StARTS), quantitative structure-activity relationships (QSAR), chemical similarity searching, structure-based design and docking, and molecular dynamics (MD) simulations. Case studies described include the design of small molecules targeting RNA expansions, the bacterial A-site, viral RNAs, and telomerase RNA. These approaches can be combined to afford a synergistic method to exploit the myriad of RNA targets in the transcriptome.

The future of crystallography in drug discovery

INTRODUCTION

X-ray crystallography plays an important role in structure-based drug design (SBDD), and accurate analysis of crystal structures of target macromolecules and macromolecule-ligand complexes is critical at all stages. However, whereas there has been significant progress in improving methods of structural biology, particularly in X-ray crystallography, corresponding progress in the development of computational methods (such as in silico high-throughput screening) is still on the horizon. Crystal structures can be overinterpreted and thus bias hypotheses and follow-up experiments. As in any experimental science, the models of macromolecular structures derived from X-ray diffraction data have their limitations, which need to be critically evaluated and well understood for structure-based drug discovery.

AREAS COVERED

This review describes how the validity, accuracy and precision of a protein or nucleic acid structure determined by X-ray crystallography can be evaluated from three different perspectives: i) the nature of the diffraction experiment; ii) the interpretation of an electron density map; and iii) the interpretation of the structural model in terms of function and mechanism. The strategies to optimally exploit a macromolecular structure are also discussed in the context of 'Big Data' analysis, biochemical experimental design and structure-based drug discovery.

EXPERT OPINION

Although X-ray crystallography is one of the most detailed 'microscopes' available today for examining macromolecular structures, the authors would like to re-emphasize that such structures are only simplified models of the target macromolecules. The authors also wish to reinforce the idea that a structure should not be thought of as a set of precise coordinates but rather as a framework for generating hypotheses to be explored. Numerous biochemical and biophysical experiments, including new diffraction experiments, can and should be performed to verify or falsify these hypotheses. X-ray crystallography will find its future application in drug discovery by the development of specific tools that would allow realistic interpretation of the outcome coordinates and/or support testing of these hypotheses.

Rational design of chemical genetic probes of RNA function and lead therapeutics targeting repeating transcripts

RNA is an important yet vastly underexploited target for small molecule chemical probes or lead therapeutics. Small molecules have been used successfully to modulate the function of the bacterial ribosome, viral RNAs and riboswitches. These RNAs are either highly expressed or can be targeted using substrate mimicry, a mainstay in the design of enzyme inhibitors. However, most cellular RNAs are neither highly expressed nor have a lead small molecule inhibitor, a significant challenge for drug discovery efforts. Herein, I describe the design of small molecules targeting expanded repeating transcripts that cause myotonic muscular dystrophy (DM). These test cases illustrate the challenges of designing small molecules that target RNA and the advantages of targeting repeating transcripts. Lastly, I discuss how small molecules might be more advantageous than oligonucleotides for targeting RNA.

Clarifying off-target effects for torcetrapib using network pharmacology and reverse docking approach

Background

Torcetrapib, a cholesteryl ester transfer protein (CETP) inhibitor which raises high-density lipoprotein (HDL) cholesterol and reduces low-density lipoprotein (LDL) cholesterol level, has been documented to increase mortality and cardiac events associated with adverse effects. However, it is still unclear the underlying mechanisms of the off-target effects of torcetrapib.

Results

In the present study, we developed a systems biology approach by combining a human reassembled signaling network with the publicly available microarray gene expression data to provide unique insights into the off-target adverse effects for torcetrapib. Cytoscape with three plugins including BisoGenet, NetworkAnalyzer and ClusterONE was utilized to establish a context-specific drug-gene interaction network. The DAVID functional annotation tool was applied for gene ontology (GO) analysis, while pathway enrichment analysis was clustered by ToppFun. Furthermore, potential off-targets of torcetrapib were predicted by a reverse docking approach. In general, 10503 nodes were retrieved from the integrative signaling network and 47660 inter-connected relations were obtained from the BisoGenet plugin. In addition, 388 significantly up-regulated genes were detected by Significance Analysis of Microarray (SAM) in adrenal carcinoma cells treated with torcetrapib. After constructing the human signaling network, the over-expressed microarray genes were mapped to illustrate the context-specific network. Subsequently, three conspicuous gene regulatory networks (GRNs) modules were unearthed, which contributed to the off-target effects of torcetrapib. GO analysis reflected dramatically over-represented biological processes associated with torcetrapib including activation of cell death, apoptosis and regulation of RNA metabolic process. Enriched signaling pathways uncovered that IL-2 Receptor Beta Chain in T cell Activation, Platelet-Derived Growth Factor Receptor (PDGFR) beta signaling pathway, IL2-mediated signaling events, ErbB signaling pathway and signaling events mediated by Hepatocyte Growth Factor Receptor (HGFR, c-Met) might play decisive characters in the adverse cardiovascular effects associated with torcetrapib. Finally, a reverse docking algorithm in silico between torcetrapib and transmembrane receptors was conducted to identify the potential off-targets. This screening was carried out based on the enriched signaling network analysis.

Conclusions

Our study provided unique insights into the biological processes of torcetrapib-associated off-target adverse effects in a systems biology visual angle. In particular, we highlighted the importance of PDGFR, HGFR, IL-2 Receptor and ErbB1tyrosine kinase might be direct off-targets, which were highly related to the unfavorable adverse effects of torcetrapib and worthy of further experimental validation.