Lecture 4: Dr. Sirish Kaushik Lakkaraju (Alex MacKerell Lab at UMB)

In search for a better drug against asthma (... and many others)

The seven-transmembrane α-helix structure
of a G protein–coupled receptor (wikipedia).

Occluded ligand binding pockets (LBP) in proteins with minimal or no accessibility to the surrounding environment represent a significant, yet challenging opportunity for structure-based and computer-aided drug design approaches. LBPs of more than half of all clinical drug targets, including the G-protein coupled receptors (GPCR) and nuclear receptors (NR) are either partially or fully occluded. As the efficacies of ligands of both GPCRs and NRs are known to be coupled to small conformational changes in their binding sites, accurate modeling of these sites is critical for future development of therapeutic agents.

To determine free energy maps of functional groups of LBPs, Dr. Kaushik Lakkaraju and his colleagues from the University of Maryland Baltimore combined a Grand-Canonical Monte Carlo/Molecular Dynamics (GCMC/MD) strategy [1] with the Site Identification by Ligand Competitive Saturation (SILCS) methodology. In this method, the inclusion of protein flexibility helps identify regions of the binding pockets not accessible in crystal conformations and allows for better quantitative estimates of relative ligand binding affinities in all the proteins tested. The research group’s approach successfully recapitulates locations of functional groups across diverse classes of ligands in the LBPs of many important NR and GPCR proteins. In their study, the differences in functional group requirements of the active and inactive states of the β2AR LBP were used in virtual drug screening to identify high-efficacy agonists targeting β2AR in Airway Smooth Muscle (ASM) cells which are very important drug targets in patients with severe asthma. The experimental tests showed that 7 of the 15 selected ligands effect ASM relaxation, representing a very high hit rate (46%) [2].

FragMaps overlaid on the LBP of AR [2]

The conformational dynamics of a macromolecule can be modulated by a number of factors, including changes in the environment, ligand binding, and interactions with other macromolecules, among others. Therefore, there is always a need for faster extraction and comparison of functionally relevant conformational dynamics information from extended molecular simulations of large and complex macromolecular systems. Dr. Lakkaraju and his colleagues also have developed DIRECT-ID, a method that quantifies the differences in macromolecular conformational dynamics and automatically extracts the structural features responsible for these changes. Given a set of molecular dynamics (MD) simulations of a macromolecule, the method calculates the norms of the differences in covariance matrices for each pair of trajectories. A matrix of these norms thus quantifies the differences in conformational dynamics across the set of simulations. As a demonstration of its applicability to bio-macromolecular systems, DIRECT-ID was used to identify relevant ligand-modulated structural variations in the β2 -adrenergic (β2 AR) G-protein coupled receptor. Micro-second MD simulations of the β2 AR in an explicit lipid bilayer were run in the apo-state and complexed with the ligands: BI-167107 (agonist), epinephrine (agonist), salbutamol (long-acting partial agonist), or carazolol (inverse agonist). Each ligand was found to modulate the conformational dynamics of β2 AR differently and DIRECT-ID analysis of the inverse-agonist vs. agonist-modulated β2 AR identified residues known through previous studies to selectively propagate deactivation/activation information, along with some previously unidentified ligand-specific microswitches across the GPCR [3]. DIRECT-ID is accessible from here.

References:

[1] Lakkaraju, S. K.; Raman, E. P.; Yu, W.; MacKerell, A. D.Sampling of Organic Solutes in Aqueous and Heterogeneous Environments using Oscillating Excess Chemical Potentials in Grand Canonical-Like Monte Carlo-Molecular Dynamics Simulations J. Chem. Theory Comput. 2014,10, 2281– 2290

[2] Lakkaraju, Sirish Kaushik, et al. "Mapping Functional Group Free Energy Patterns at Protein Occluded Sites: Nuclear Receptors and G-Protein Coupled Receptors." Journal of chemical information and modeling 55.3 (2015): 700-708.

[3] Lakkaraju, Sirish Kaushik, et al. "DIRECT‐ID: An automated method to identify and quantify conformational variations—application to β2‐adrenergic GPCR." Journal of computational chemistry 37.4 (2016): 416-425.

_______________________________________________________________________________

Dr. Sirish Kaushik Lakkaraju received his Bachelor of Engineering degree in electronics and instrumentation from the University of Madras at Chennai, India in 2003. He received his Master of Science degree in biomedical engineering from Texas A&M University in 2008 and his Ph.D. in December 2011 with Dr. Wonmuk Hwang. After graduation, he worked as a postdoctoral fellow under the supervision of Dr. Alex MacKerell at University of Maryland, Baltimore on the projects that involve method development and in-silico drug design. Recently, he started his second postdoc at Pfizer. 

Dr. Lakkaraju has published 8 first-author papers and co-authored 2 papers in peer-reviewed international journals. His current research interests lie in studying the role of dynamics of biological macromolecules and macromolecular assemblies on their function and also computational drug design.