Expression, purification, and circular dichroism analysis of human CDK9

Application of Machine Learning Techniques to Predict Binding Affinity for Drug Targets: A Study of Cyclin-Dependent Kinase 2

BACKGROUND

The elucidation of the structure of cyclin-dependent kinase 2 (CDK2) made it possible to develop targeted scoring functions for virtual screening aimed to identify new inhibitors for this enzyme. CDK2 is a protein target for the development of drugs intended to modulate cellcycle progression and control. Such drugs have potential anticancer activities.

OBJECTIVE

Our goal here is to review recent applications of machine learning methods to predict ligand- binding affinity for protein targets. To assess the predictive performance of classical scoring functions and targeted scoring functions, we focused our analysis on CDK2 structures.

METHODS

We have experimental structural data for hundreds of binary complexes of CDK2 with different ligands, many of them with inhibition constant information. We investigate here computational methods to calculate the binding affinity of CDK2 through classical scoring functions and machine- learning models.

RESULTS

Analysis of the predictive performance of classical scoring functions available in docking programs such as Molegro Virtual Docker, AutoDock4, and Autodock Vina indicated that these methods failed to predict binding affinity with significant correlation with experimental data. Targeted scoring functions developed through supervised machine learning techniques showed a significant correlation with experimental data.

CONCLUSION

Here, we described the application of supervised machine learning techniques to generate a scoring function to predict binding affinity. Machine learning models showed superior predictive performance when compared with classical scoring functions. Analysis of the computational models obtained through machine learning could capture essential structural features responsible for binding affinity against CDK2.

How Docking Programs Work

Protein-ligand docking simulations are of central interest for computer-aided drug design. Docking is also of pivotal importance to understand the structural basis for protein-ligand binding affinity. In the last decades, we have seen an explosion in the number of three-dimensional structures of protein-ligand complexes available at the Protein Data Bank. These structures gave further support for the development and validation of in silico approaches to address the binding of small molecules to proteins. As a result, we have now dozens of open source programs and web servers to carry out molecular docking simulations. The development of the docking programs and the success of such simulations called the attention of a broad spectrum of researchers not necessarily familiar with computer simulations. In this scenario, it is essential for those involved in experimental studies of protein-ligand interactions and biophysical techniques to have a glimpse of the basics of the protein-ligand docking simulations. Applications of protein-ligand docking simulations to drug development and discovery were able to identify hits, inhibitors, and even drugs. In the present chapter, we cover the fundamental ideas behind protein-ligand docking programs for non-specialists, which may benefit from such knowledge when studying molecular recognition mechanism.

Van der Waals Potential in Protein Complexes

Van der Waals forces are determinants of the formation of protein-ligand complexes. Physical models based on the Lennard-Jones potential can estimate van der Waals interactions with considerable accuracy and with a computational complexity that allows its application to molecular docking simulations and virtual screening of large databases of small organic molecules. Several empirical scoring functions used to evaluate protein-ligand interactions approximate van der Waals interactions with the Lennard-Jones potential. In this chapter, we present the main concepts necessary to understand van der Waals interactions relevant to molecular recognition of a ligand by the binding pocket of a protein target. We describe the Lennard-Jones potential and its application to calculate potential energy for an ensemble of structures to highlight the main features related to the importance of this interaction for binding affinity.

Homology Modeling of Protein Targets with MODELLER

Homology modeling is a computational approach to generate three-dimensional structures of protein targets when experimental data about similar proteins are available. Although experimental methods such as X-ray crystallography and nuclear magnetic resonance spectroscopy successfully solved the structures of nearly 150,000 macromolecules, there is still a gap in our structural knowledge. We can fulfill this gap with computational methodologies. Our goal in this chapter is to explain how to perform homology modeling of protein targets for drug development. We choose as a homology modeling tool the program MODELLER. To illustrate its use, we describe how to model the structure of human cyclin-dependent kinase 3 using MODELLER. We explain the modeling procedure of CDK3 apoenzyme and the structure of this enzyme in complex with roscovitine.

Hydrogen Bonds in Protein-Ligand Complexes

Fast and reliable evaluation of the hydrogen bond potential energy has a significant impact in the drug design and development since it allows the assessment of large databases of organic molecules in virtual screening projects focused on a protein of interest. Semi-empirical force fields implemented in molecular docking programs make it possible the evaluation of protein-ligand binding affinity where the hydrogen bond potential is a common term used in the calculation. In this chapter, we describe the concepts behind the programs used to predict hydrogen bond potential energy employing semi-empirical force fields as the ones available in the programs AMBER, AutoDock4, TreeDock, and ReplicOpter. We described here the 12-10 potential and applied it to evaluate the binding affinity for an ensemble of crystallographic structures for which experimental data about binding affinity are available.

Human G-CSF synthesis using stress-responsive bacterial proteins

We previously reported that under the stress condition caused by the addition of 2-hydroxyethyl disulfide, a thiol-specific oxidant, to growing cultures of Escherichia coli BL21(DE3), a population of stress-responsive proteins [peptidyl-prolyl cis–trans isomerase B (PpiB), bacterioferritin (Bfr), putative HTH-type transcriptional regulator yjdC (YjdC), dihydrofolate reductase (FolA), chemotaxis protein cheZ (CheZ), and glutathione synthetase (GshB)] were significantly upregulated when compared with the nonstress condition. When those stress-responsive proteins were used as fusion partners for the expression of human granulocyte colony-stimulating factor (hG-CSF), the solubility of hG-CSF was dramatically enhanced in E. coli cytoplasm, whereas almost all of the directly expressed hG-CSF were aggregated to inclusion bodies. In addition, the spectra of circular dichroism measured with the purified hG-CSF were identical to that of standard hG-CSF, implying that the synthesized hG-CSF has native conformation. These results indicate that the bacterial stress-responsive proteins could be potent fusion expression partners for aggregation-prone heterologous proteins in E. coli cytoplasm.