We present the performances of our mathematical deep learning (MathDL) models for D3R Grand Challenge 4 (GC4). This challenge involves pose prediction, affinity ranking, and free energy estimation for beta secretase 1 (BACE) as well as affinity ranking and free energy estimation for Cathepsin S (CatS). We have developed advanced mathematics, namely differential geometry, algebraic graph, and/or algebraic topology, to accurately and efficiently encode high dimensional physical/chemical interactions into scalable low-dimensional rotational and translational invariant representations. These representations are integrated with deep learning models, such as generative adversarial networks (GAN) and convolutional neural networks (CNN) for pose prediction and energy evaluation, respectively. Overall, our MathDL models achieved the top place in pose prediction for BACE ligands in Stage 1a. Moreover, our submissions obtained the highest Spearman correlation coefficient on the affinity ranking of 460 CatS compounds, and the smallest centered root mean square error on the free energy set of 39 CatS molecules. It is worthy to mention that our method for docking pose predictions has significantly improved from our previous ones. * Corresponding to Guo-Wei Wei. The Drug Design Data Resource (D3R) offers blind communitywide challenges of ligand pose and binding affinity ranking predictions. 1-3 Benchmarks in D3R contests contain high quality structures and reliable binding energies supplied by experimental groups before the publication. These challenges provide computer-aided drug design (CADD) community a great opportunity to validate, calibrate, and develop drug virtual screening (VS) models. The latest D3R Grand Challenge 4 (GC4), took place from September 4th 2018 to December 4th, 2018. GC4 presented two different protein targets, Cathepsin S (CatS) and beta secretase 1 (BACE), which were generously supplied by Janssen Pharmaceuticals and Novartis, respectively. There were two stages in GC4. The first one has two subchallenges, namely Stage 1a and Stage 1b. In Stage 1a, participants were asked to predict the pose, rank the affinity, and estimate the free energy of BACE ligands. Following Stage 1a, Stage 1b revealed the receptor structures and participants were asked again to predict the crystallographic poses of 20 BACE ligands. There was no affinity calculation in this stage 1b. The second part of GC4 was called Stage 2 which contained the affinity rankings and free energy challenges for both BACE and CatS compounds. In this last stage, participants were able to take advantage of experimental structures of BACE complexes released right after stage 1b.A successful VS model requires a reliable ligand conformation generation and highly accurate scoring function to predict binding affinities. There are several state-of-the-art software packages to take care of the first component of VS, for example, Autodock Vina, 4 GOLD, 5 GLIDE, 6 ICM, 7 etc. Unfortunately, one may fail dramatically to achieve decent poses if blindly using t...