News & Events


posted Aug 23, 2017, 3:14 PM by Emilio Gallicchio   [ updated Aug 24, 2017, 6:34 AM ]

The Computers in Chemistry Division of the American Chemical Society awarded Prof. Emilio Gallicchio of Brooklyn College the NVIDIA GPU Best Poster Award at the 254th ACS National Meeting & Exposition in Washington DC. Prof. Gallicchio's poster is entitled "Efficient GPU/OpenMM implementation of the AGBNP solvation model for macromolecular binding"


Emilio Gallicchio,, Denise Kilburg, Baofeng Zhang
Brooklyn College, Brooklyn, NY, and Doctoral Programs in Chemistry and Biochemistry, Graduate Center of the City University of New York, New York, NY

The poster describes the GPU implementation of the Analytic Generalized Born plus Non-Polar implicit solvent model (AGBNP) for the OpenMM molecular simulation package. As part of this work, novel recursive computational geometry algorithms have been developed to estimate volumes and surface areas of macromolecules. The Generalized Born pairwise descreening algorithm has been formulated to take advantage of the GPU architecture. The GPU implementation is 50 to 100 times faster than our best CPU implementation. Illustrative applications to binding free energy calculations are presented. OpenMM plugin:

The generous donation of an impressive Titan Xp GPU by NVIDIA is greatly appreciated. We will put it to good use!

Professor Gallicchio is awarded the Henry Wasser Award from CUNY!

posted Jul 6, 2017, 2:28 PM by Rajat Kumar Pal   [ updated Jul 6, 2017, 2:29 PM ]

Professor Gallicchio is the recipient of the 2017 Henry Wasser award from the CUNY Academy. This award is given to outstanding CUNY Assistant Professors by the CUNY Academy for the Humanities and Sciences. Congratulations Dr. Gallicchio!

Presentation on Analytic Model of Binding @CECAM

posted Jun 19, 2017, 2:44 PM by Emilio Gallicchio

I have had the opportunity to present our recent development of an analytic probabilistic model of alchemical binding at the "Beyond Kd's" CECAM conference organized by Sereina Riniker, Lyna Luo, and Chris Chipot at Lausanne in Switzerland. One of the aims of this project is to extract physical parameters such as the volume of the binding site and the local dielectric constant of the binding site from binding free energy profiles. 

Click below to view the slides of the presentation. Feedback is welcome.

A Novel Surface Area Algorithm for GPUs added to OpenMM

posted Jan 8, 2017, 8:08 AM by Emilio Gallicchio   [ updated Jan 8, 2017, 9:30 AM ]

Surface area of a moleculeAfter a year-long effort, we finally developed an efficient implementation of the Gaussian-based volume and surface area model for GPUs. The work, performed in collaboration with Peter Eastman Vijay Pande at Stanford University, is described in a recent publication to appear in the Journal of Computational Chemistry: 

Baofeng Zhang, Denise Kilburg, Peter Eastman, Vijay S. Pande, Emilio Gallicchio. Efficient Gaussian Density Formulation of Volume and Surface Areas of Macromolecules on Graphical Processing Units. J. Comp. Chem. (2017). pdf of submitted manuscript

Volume and surface area of macromolecules are employed in medicinal chemistry to measure structural similarity and complementarity of compounds. They are also the basis of many implicit models of non-polar solvation, which is our primary interest.

With this algorithm we were able to achieve a 50- to 100-fold speed-up on GPU's relative to our best CPU implementation. more ->

The GaussVol code is freely available on github as a plugin of the OpenMM molecular mechanics package.

Support from the National Science Foundation is gratefully acknowledged

Publication on binding free energy study of RNase H inhibitors

posted Sep 18, 2016, 8:16 AM by Emilio Gallicchio   [ updated Sep 18, 2016, 8:25 AM ]

Our collaboration with the group of Ryan Murelli here at Brooklyn College, and John A Beutler and Stuart F Le Grice at NCI, aimed at identifying  α-hydroxytropolones derivatives capable of inhibiting the RNase H enzyme of the HIV virus is now described in two recent publications (below).

The computational work has been made possible by the WEB computational grid at Brooklyn College maintained by the ITS office at Brooklyn College. Additional computing has been performed on the NSF XSEDE SuperMIC cluster at LSU. Thanks!

Good job done! to Baofeng Zhang who completed his postdoctoral residency in the lab.

Baofeng Zhang, Michael P. D'Erasmo, Ryan P. Murelli, Emilio Gallicchio. Free Energy-Based Virtual Screening and Optimization of RNase H Inhibitors of HIV-1 Reverse Transcriptase. ACS Omega. 2016. pdf of submitted manuscript

Ryan P Murelli, Michael P D'Erasmo, Danielle R Hirsch, C Meck, T Masaoka, JA Wilson, Baofeng Zhang, Rajat K Pal, Emilio Gallicchio, John A Beutler, Stuart F Le Grice. Synthetic α-hydroxytropolones as inhibitors of HIV reverse transcriptase ribonuclease H activity. MedChemComm. 2016. doi:10.1039/C6MD00238B

Support from the National Science Foundation is gratefully acknowledged


posted Sep 7, 2016, 6:20 AM by Emilio Gallicchio   [ updated Sep 7, 2016, 6:24 AM ]

We are happy to announce that a team of CUNY investigators led by our laboratory and including Prof. Lauren Wickstrom (Borough of Manhattan Community College), Prof. Tom Kurtzman (Lehman College), and Prof. Wayne Harding (Hunter College) has been awarded a 1-year Interdisciplinary Research Grant award of $39,500 by the CUNY Office of Research to conduct theoretical and experimental studies of the inhibition of the D3 dopamine receptor:

A Combined Treatment of Hydration and Dynamical Effects for the Modeling of Protein-Ligand Binding Thermodynamics

Summary. The project seeks to develop an improved computational protocol for the study of protein-drug interactions marrying the accuracy of detailed explicit solvation models with the versatility of implicit solvation approaches, and to validate it experimentally on the important D3 dopamine receptor drug target. The work will lead to improved descriptions of the displacement of water molecules from the protein surface and of the dynamical response of molecular conformations upon binding. Current computer models of molecular binding lack in one or both of these critical elements. The primary aim of the project is to exploit the unique and complementary expertise of our team of investigators at CUNY to establish a proof-of-concept base and acquire sufficient preliminary data to compete for federal funding at the national level.

The team will work to develop a custom parameterization of the AGBNP solvation model for the D3 dopamine receptor using accurate Hydration Site Analysis data obtained from explicit solvent simulations. The resulting model will be employed to predict binding constants of a virtual library of D3 inhibitors using the BEDAM alchemical methodology. The most promising identified inhibitors will be synthesized and essayed for binding.  Selective inhibition of the dopamine D3 receptor is a promising therapeutic approach for the treatment of drug abuse disorders. Building on the methodology employed in our recent participation to the SAMPL5 blinded challenge, the project will for the first time attempt to apply a model combining enclosed hydration effects and a dynamical description of binding to a complex protein receptor system.

The award will help fund the research of Rajat Pal, a Ph.D. student in our lab, for the next year. Thanks!

Publication on SAMPL5 host-guest blinded challenge

posted Aug 24, 2016, 12:50 PM by Emilio Gallicchio   [ updated Sep 24, 2016, 9:38 AM ]

Our contribution to the SAMPL5 blinded challenge is now described in a paper to appear in the Journal of Computer Aided Molecular Design:

Rajat Kumar Pal, Kamran Haider, Divya Kaur, William Flynn, Junchao Xia,  Ronald M Levy, Tetiana Taran,* Lauren Wickstrom,  Tom Kurtzman, Emilio Gallicchio. A Combined Treatment of Hydration and Dynamical Effects for the Modeling of Host-Guest Binding Thermodynamics: The SAMPL5 Blinded Challenge. J. Comp. Aided. Mol. Des. 2016. pdf of submitted manuscript *Undergraduate student.

The computational work has been performed on WEB computational grid at Brooklyn College maintained by the ITS office at Brooklyn College, and on the NSF XSEDE SuperMIC cluster at LSU. Thanks!

Support from the National Science Foundation is gratefully acknowledged

Review Paper on Computational Modeling of Protein-Peptide Interactions

posted Aug 24, 2016, 12:48 PM by Emilio Gallicchio

Denise Kilburg has recently completed a review paper summarizing the various methodologies in current use for the modeling of protein-peptide binding with examples of their applications:

Denise Kilburg and Emilio Gallicchio. 
Recent Advances in Computational Models for the Study Protein-Peptide Interactions. Adv. Prot. Chem. Struct. Biol. (2016). pdf of submitted manuscript 

Posters at ACS Sci-Mix

posted Aug 24, 2016, 12:24 PM by Emilio Gallicchio   [ updated Aug 24, 2016, 12:31 PM ]

The posters by Baofeng Zhang and Rajat Pal were selected for the prestigious interdisciplinary Sci-Mix event at the Fall ACS meeting in Philadelphia :

Targeting viral receptors using binding free energy-based virtual screening and GPU-accelerated software

Baofeng Zhang et al.

Abstract:We will describe the GPU-accelerated formulation of our single-decoupling free energy perturbation method for the rapid estimation of absolute binding free energies of protein-ligand complexes. The computational system is based on the OpenMM MD engine and the Asynchronous Replica Exchange job coordination system. While computationally intensive, free energy-based virtual screening can more accurately discriminate binders from non-binders because it includes entropic and conformational reorganization effects in addition to energetic complementarity factors. The method is therefore suitable for the secondary screening of focused ligand libraries, for example to improve on docking enrichment. Because it targets absolute binding free energies, this protocol is useful for preliminary free energy scoring of ligand libraries with diverse scaffolds. We will present examples drawn from recent screening campaigns targeting viral targets with alpha-hydroxytropolone derivatives. The combination of binding free energy predictions with structural analysis of dynamical trajectories has pointed out key aspects which we are currently pursuing to optimize leads and identify novel classes of inhibitors.

Alchemical computational methodologies for the estimation of binding free energies of supramolecular complexes

Rajat Kumar Pal et al.

Abstract:We will review the development and current applications of computational binding free energy protocols based on the alchemical Single-Decoupling Method (SDM) and multi-dimensional Asynchronous Replica Exchange sampling framework (AsyncRE). Asynchronous equilibrium techniques of this kind are particularly suited for heterogeneous computational clusters and computational grids. We have started to explore non-equilibrium approaches which, although inherently less statistically efficient than equilibrium ones, have the advantage of not incurring the overhead of communication between computational threads in a dynamic processing environment. We will describe a novel asynchronous methodology for binding free energy estimation which works under non-equilibrium conditions and requires no exchanges between replicas, thereby leading to higher computational throughput than equilibrium approaches on computational grids. We will present illustrative examples of equilibrium and non-equilibrium calculations drawn from recent investigations of Octa-acid and CBClip host-guest supramolecular complexes. The results indicate that non-equilibrium techniques have the potential to accelerate and streamline binding free energy calculations on computational grids.

Presentation at the Vertex Free Energy Workshop

posted Jun 9, 2016, 12:01 PM by Emilio Gallicchio   [ updated Jun 9, 2016, 12:04 PM ]

This Free Energy Workshop has been held at Vertex's headquarters in Boston last May. It attracts the best practitioners in the area of free energy modeling for biological and pharmaceutical applications. I like it especially because people here are not afraid of exposing problems and challenges. I talked about our latest attempts at employing single-decoupling alchemical transformation to compute the absolute binding free energies of protein-peptide complexes. Slides and animations are below:


YouTube Video

YouTube Video

1-10 of 41