QCSPScore: a new scoring function for driving protein-ligand docking with quantitative chemical shifts perturbations

Through the use of information about the biological target structure, the optimization of potential drugs can be improved. In this work I have developed a procedure that uses the quantitative change in the chemical perturbations (CSP) in the protein from NMR experiments for driving protein-ligand docking. The approach is based on a hybrid scoring function (QCSPScore) which combines traditional DrugScore potentials, which describe the interaction between protein and ligand, with Kendall’s rank correlation coefficient, which evaluates docking poses in terms of their agreement with experimental CSP. Prediction of the CSP for a specific ligand pose is done efficiently with an empirical model, taking into account only ring current effects. QCSPScore has been implemented in the AutoDock software package. Compared to previous methods, this approach shows that the use of rank correlation coefficient is robust to outliers. In addition, the prediction of native-like complex geometries improved because the CSP are already being used during the docking process, and not only in a post-filtering setting for generated docking poses. Since the experimental information is guaranteed to be quantitatively used, CSP effectively contribute to align the ligand in the binding pocket. The first step in the development of QCSPScore was the analysis of 70 protein-ligand complexes for which reference CSP were computed. The success rate in the docking increased from 71% without involvement of CSP to 100% if CSP were considered at the highest weighting scheme. In a second step QCSPScore was used in re-docking three test cases, for which reference experimental CSP data was available. Without CSP, i.e. in the use of conventional DrugScore potentials, none of the three test cases could be successfully re-docked. The integration of CSP with the same weighting factor as described above resulted in all three cases successfully re-docked. For two of the three complexes, native-like solutions were only produced if CSP were considered.Conformational changes in the binding pockets of up to 2 Å RMSD did not affect the success of the docking. QCSPScore will be particularly interesting in difficult protein-ligand complexes. They are in particular those cases in which the shape of the binding pocket does not provide sufficient steric restraints such as in flat protein-protein interfaces and in the virtual screening of small chemical fragments.Durch die Verwendung von Information über die biologische Zielstruktur kann die Optimierung potentieller Wirkstoffe verbessert werden. Im Rahmen dieser Arbeit habe ich ein Verfahren entwickelt, das quantitativ die Veränderung der Chemischen Verschieben (CSP) im Protein aus NMR-Experimenten für das Protein-Ligand-Docking verwendet. Der Ansatz basiert auf einer Hybridbewertungsfunktion (QCSPScore) und kombiniert herkömmliche DrugScore-Potentiale, welche die Wechselwirkung zwischen Protein und Ligand beschreiben, mit dem Rangkorrelationskoeffizienten nach Kendall, der die Dockingposen hinsichtlich ihrer Übereinstimmung mit experimentellen CSP. Die Vorhersage der CSP für einen bestimmten Liganden geschieht effizient mit einem empirischen Modell, wobei nur Ringstromeffekte berücksichtigt werden. QCSPScore wurde in das AutoDock Softwarepaket implementiert. Im Vergleich zu früheren Verfahren zeigt dieser Ansatz, dass die Verwendung des Rangkorrelationskoeffizienten robuster ist gegenüber Ausreißern in den vorhergesagten CSP. Außerdem ist die Vorhersage nativ-ähnlicher Komplexgeometrien verbessert, da die CSP bereits während des Docking-Prozesses eingesetzt werden, und nicht erst in einem nachträglichen Filter für generierte Dockingposen. Da die experimentelle Informationen quantitativ benutzt werden wird sichergestellt, dass die CSP effektiv dazu beitragen, den Liganden in der Bindetasche auszurichten. Der erste Schritt bei der Entwicklung des QCSPScore war die Analyse von 70 Protein-Ligand-Komplexen, für die als Referenz CSP vorhergesagt wurden. Die Erfolgsrate im Docking erhöhte sich von 71 %, ohne Einbeziehung von CSP, auf 100 %, wenn CSP mit höchster Gewichtung mit einbezogen wurden. Die globale Optimierung auf der kombinierten Docking-Energiehyperfläche ist also erfolgreich. In einem zweiten Schritt wurde QCSPScore zum Docking dreier Testfälle verwendet, für die als Referenz experimentelle CSP zur Verfügung standen. Ohne CSP, d.h. bei der Verwendung von herkömmlichen DrugScore-Potentialen, konnte keiner der drei Testfälle erfolgreich gedockt werden. Die Einbeziehung von CSP mit dem selben hohen Gewichtungsfaktor wie oben führte in allen drei Fällen zu erfolgreichen Docking-Ergebnissen. Für zwei der drei Komplexe wurden zudem nur bei Einbeziehung der experimentellen Information nativ-ähnliche Geometrien vorhergesagt. Konformationelle Änderungen der Bindetasche bis zu 2 Å RMSD beeinträchtigen den Erfolg des Dockings nicht. Ich bin davon überzeugt, dass mein Verfahren besonders für Protein-Ligand-Komplexe interessant sein wird, für die die Vorhersage nativ-ähnlicher Komplexe bislang schwierig war. Das sind insbesondere solche Fälle, in denen die Form der Bindetasche zur Vorhersage des Komplexes nicht ausreichend, wie das bei flachen Protein-Protein-Wechselwirkungsregionen oder beim virtuellen Screening kleiner Fragmente der Fall ist

Biological Systems Workbook: Data modelling and simulations at molecular level

Nowadays, there are huge quantities of data surrounding the different fields of biology derived from experiments and theoretical simulations, where results are often stored in biological databases that are growing at a vertiginous rate every year. Therefore, there is an increasing research interest in the application of mathematical and physical models able to produce reliable predictions and explanations to understand and rationalize that information. All these investigations are helping to overcome biological questions pushing forward in the solution of problems faced by our society. In this Biological Systems Workbook, we aim to introduce the basic pieces allowing life to take place, from the 3D structural point of view. We will start learning how to look at the 3D structure of molecules from studying small organic molecules used as drugs. Meanwhile, we will learn some methods that help us to generate models of these structures. Then we will move to more complex natural organic molecules as lipid or carbohydrates, learning how to estimate and reproduce their dynamics. Later, we will revise the structure of more complex macromolecules as proteins or DNA. Along this process, we will refer to different computational tools and databases that will help us to search, analyze and model the different molecular systems studied in this course

Theoretical-experimental study on protein-ligand interactions based on thermodynamics methods, molecular docking and perturbation models

The current doctoral thesis focuses on understanding the thermodynamic events of protein-ligand interactions which have been of paramount importance from traditional Medicinal Chemistry to Nanobiotechnology. Particular attention has been made on the application of state-of-the-art methodologies to address thermodynamic studies of the protein-ligand interactions by integrating structure-based molecular docking techniques, classical fractal approaches to solve protein-ligand complementarity problems, perturbation models to study allosteric signal propagation, predictive nano-quantitative structure-toxicity relationship models coupled with powerful experimental validation techniques. The contributions provided by this work could open an unlimited horizon to the fields of Drug-Discovery, Materials Sciences, Molecular Diagnosis, and Environmental Health Sciences

Development and Application of Computational Biology tools for Biomedicine

Biomolecular simulation can be considered as a virtual microscope for molecular biology, allowing to gain insights into the sub-cellular mechanisms of biological relevance at spatial and temporal scales that are difficult to observe experimentally. It provides a powerful tool to link the laws of physics with the complex behavior of biological systems. Dramatic recent advancements in achievable simulation speed and the underlying physical models will increasingly lead to molecular views of large systems. These improvements may largely affect biological sciences. In this thesis, I have applied computational molecular biology approaches to different biological systems using state of the art structural bioinformatics and computational biophysics tools (Chapter 3). My principal focus was on the computational design of molecular imprinted polymers (MIPs), which have recently attracted significant attention as cost effective substitutes for natural antibodies and receptors in chromatography, sensors and assays. I have used molecular modelling in the optimization of polymer compositions to make high affinity synthetic receptors based on Molecular Imprinting. In particular, I developed a new free of charge protocol that can be performed within just few hours that outputs a list of candidate monomers which are capable of strong binding interactions with the template. Furthermore, I have produced a new computational method for the calculation of the ideal monomer: template stoichiometric ratio to be used in the lab for the MIPs synthesis. These protocols have been implemented as a webserver that is available at http://mirate.di.univr.it/. In parallel, I have also investigated the modelling of much more complex MIPs systems by the introduction of some factors e.g. solvent and cross-linker molecules that are also essential in the polymerisation process. A novel algorithm, which mimics a radical polymerization mechanism, has been written for application in the rational design of MIPs (Chapter 4). Moreover, I have been involved in the field of computational molecular biomedicine. Indeed, in Chapters 5 and 6 I describe the work done in collaboration with two labs at the Department of Neuroscience, Biomedicine and Movement Sciences, University of Verona. In Chapter 5, starting from unpublished experimental data I have computationally characterized the interaction of ACOT8 with HIV-1 Nef accessory protein. I have performed a detailed structural and functional characterization of these two proteins in order to infer any possible functional details about their interactions. The bioinformatics predictions were then confirmed by wet-lab experiments. I have also carried out a detailed structural and functional characterization of two pathogenic mutations of AGT-Mi (Chapter 6). In particular, I have used classical molecular dynamics (MD) simulations to study the possible interference with the dimerization process of AGT-Mi exerted by I244T-Mi and F152I-Mi mutants. Those variants are associated with Primary Hyperoxaluria type 1 disease. In Chapter 7, I present the coarse-grained MD simulations of Membrane/Human ileal bile-acid-binding protein Interactions. This study was carried out in collaboration with the NMR group at the University of Verona and it is a part of an extensive research aimed at better understanding of the main biomolecular interactions in crowded cellular environments. MD simulations results were in agreement with experimental findings

The Impact of Dynamics in Protein Assembly

Predicting the assembly of multiple proteins into specific complexes is critical to understanding their biological function in an organism, and thus the design of drugs to address their malfunction. Consequently, a significant body of research and development focuses on methods for elucidating protein quaternary structure. In silico techniques are used to propose models that decode experimental data, and independently as a structure prediction tool. These computational methods often consider proteins as rigid structures, yet proteins are inherently flexible molecules, with both local side-chain motion and larger conformational dynamics governing their behaviour. This treatment is particularly problematic for any protein docking engine, where even a simple rearrangement of the side-chain and backbone atoms at the interface of binding partners complicates the successful determination of the correct docked pose. Herein, we present a means of representing protein surface, electrostatics and local dynamics within a single volumetric descriptor, before applying it to a series of physical and biophysical problems to validate it as representative of a protein. We leverage this representation in a protein-protein docking context and demonstrate that its application bypasses the need to compensate for, and predict, specific side-chain packing at the interface of binding partners for both water-soluble and lipid-soluble protein complexes. We find little detriment in the quality of returned predictions with increased flexibility, placing our protein docking approach as highly competitive versus comparative methods. We then explore the role of larger, conformational dynamics in protein quaternary structure prediction, by exploiting large-scale Molecular Dynamics simulations of the SARS-CoV-2 spike glycoprotein to elucidate possible high-order spike-ACE2 oligomeric states. Our results indicate a possible novel path to therapeutics following the COVID-19 pandemic. Overall, we find that the structure of a protein alone is inadequate in understanding its function through its possible binding modes. Therefore, we must also consider the impact of dynamics in protein assembly

Examination of Molecular Recognition in Protein-Ligand Interactions

This dissertation is a compilation of two main projects that were investigated during my thesis research. The first project was a prospective study which identified and characterized drug-like inhibitors of a prototype of bacterial two-component signal transduction response regulator using computational and experimental methods. The second project was the development and validation of a scoring function, PHOENIX, derived using high-resolution structures and calorimetry measurements to predict binding affinities of protein-ligand interactions. Collectively, my thesis research aimed to better understand the underlying driving forces and principles which govern molecular recognition and molecular design. A prospective study coupled computational predictions with experimental validation resulted in the discovery of first-in-class inhibitors targeting a signal transduction module important for bacterial virulence. Development and validation of the PHOENIX scoring function for binding affinity prediction derived using high-resolution structures and calorimetry measurements should guide future molecular recognition studies and endeavors in computer-aided molecular design. To request for an electronic copy of this dissertation, please email the author: yattang at gmail dot com)

Pharmacogenetic modeling of human cytochrome P450 2D6; On the force of variation in inducing toxicity

Understanding the way in which drugs are metabolized by CYP2D6 and hence the underlying mechanisms that define potential toxicity is crucial to avoid adverse reactions. The high occurrence of CYP2D6 polymorphs enhances the complexity of the toxicity assessment of a drug candidate and should be tackled from early drug discovery phase on. The research described in this PhD thesis has been performed to provide novel fundamental insights regarding the metabolic activity of CYP2D6 wild-type and several polymorphs using various state-of-the-art in silico techniques. The results of the CYP2D6-focused studies enhance our knowledge regarding the enzyme particularities, and can be used to accelerate the development of CYP2D6 modeling tools with more accurate and reliable predictions

NMR studies of transiente protein complexes

Dissertação apresentada para obtenção do Grau de Doutor em Bioquímica,especialidade Bioquímica Física,pela Universidade Nova de Lisboa,Faculdade de Ciências e TecnologiaFundação para a Ciência e Tecnologia - Bolsa de Doutoramento (SFRH/BD/25342/2005)no âmbito do Programa Operacional Potencial Humano, da União Europeia (Fundo Social Europeu

Predicting multibody assembly of proteins

textThis thesis addresses the multi-body assembly (MBA) problem in the context of protein assemblies. [...] In this thesis, we chose the protein assembly domain because accurate and reliable computational modeling, simulation and prediction of such assemblies would clearly accelerate discoveries in understanding of the complexities of metabolic pathways, identifying the molecular basis for normal health and diseases, and in the designing of new drugs and other therapeutics. [...] [We developed] F²Dock (Fast Fourier Docking) which includes a multi-term function which includes both a statistical thermodynamic approximation of molecular free energy as well as several of knowledge-based terms. Parameters of the scoring model were learned based on a large set of positive/negative examples, and when tested on 176 protein complexes of various types, showed excellent accuracy in ranking correct configurations higher (F² Dock ranks the correcti solution as the top ranked one in 22/176 cases, which is better than other unsupervised prediction software on the same benchmark). Most of the protein-protein interaction scoring terms can be expressed as integrals over the occupied volume, boundary, or a set of discrete points (atom locations), of distance dependent decaying kernels. We developed a dynamic adaptive grid (DAG) data structure which computes smooth surface and volumetric representations of a protein complex in O(m log m) time, where m is the number of atoms assuming that the smallest feature size h is [theta](r[subscript max]) where r[subscript max] is the radius of the largest atom; updates in O(log m) time; and uses O(m)memory. We also developed the dynamic packing grids (DPG) data structure which supports quasi-constant time updates (O(log w)) and spherical neighborhood queries (O(log log w)), where w is the word-size in the RAM. DPG and DAG together results in O(k) time approximation of scoring terms where k << m is the size of the contact region between proteins. [...] [W]e consider the symmetric spherical shell assembly case, where multiple copies of identical proteins tile the surface of a sphere. Though this is a restricted subclass of MBA, it is an important one since it would accelerate development of drugs and antibodies to prevent viruses from forming capsids, which have such spherical symmetry in nature. We proved that it is possible to characterize the space of possible symmetric spherical layouts using a small number of representative local arrangements (called tiles), and their global configurations (tiling). We further show that the tilings, and the mapping of proteins to tilings on arbitrary sized shells is parameterized by 3 discrete parameters and 6 continuous degrees of freedom; and the 3 discrete DOF can be restricted to a constant number of cases if the size of the shell is known (in terms of the number of protein n). We also consider the case where a coarse model of the whole complex of proteins are available. We show that even when such coarse models do not show atomic positions, they can be sufficient to identify a general location for each protein and its neighbors, and thereby restricts the configurational space. We developed an iterative refinement search protocol that leverages such multi-resolution structural data to predict accurate high resolution model of protein complexes, and successfully applied the protocol to model gp120, a protein on the spike of HIV and currently the most feasible target for anti-HIV drug design.Computer Science