Publications (2015-present)

71. Harris RC and Shen J*, GPU-accelerated implementation of continuous constant pH molecular dynamics in Amber: pKa predictions with single pH simulations. J Chem Inf Model 59: 4821-4832, 2019.

70. Tsai CC, Yue Z, and Shen J*, How electrostatic coupling enables conformational plasticity in a tyrosine kinase.

J Am Chem Soc 141: 15092-15101, 2019.

We present a computational tool to assess the nucleophilicities of lysines and cysteines to accelerate covalent drug design. The tool enabled us to discover that kinases in a rare inactive conformation can be covalently targeted at catalytic lysines. A new strategy for designing highly selective, lysine-targeted covalent kinase inhibitors was proposed to address the emergent drug resistance issue facing the cysteine-targeted strategy and help usher in the next generation covalent kinase inhibitor design.

69. Kim E, Li J, Kang M, Kelly DL, Chen S, Napolitano A, Panzella L, Shi X, Yan K, Wu S, Shen J, Bentley, WE, Payne GF, Redox is a global biodevice information processing modality. Proc IEEE: Special Issue on "Molecular Communications and Networking" 107: 1402-1424, 2019. 

68. Liu R, Yue Z, Tsai CC, and Shen J*, Assessing lysine and cysteine reactivities for designing targeted covalent kinase inhibitors. J Am Chem Soc 141: 6553-6560, 2019. 

67. Li J, Wu S, Kim E, Yan K, Liu H, Liu C, Dong H, Qu X, Shi X, Shen J, Bentley, WE, Payne GF*, Electrobiofabrication: electrically-based fabrication with biologically-derived materials. Biofabrication 11: 032002, 2019. 

66. Tsai CC, Payne GF, and Shen J*,  Exploring pH-responsive, switchable crosslinking mechanisms for programming reconfigurable hydrogels based on aminopolysaccharides. Chem Mater 30: 8597-8605, 2018

Designing dynamic, stimuli-responsive, reconfigurable materials has growing interest in academia and industry. We report a computational study that provided novel, atomically-detailed physical insights for a pH-responsive, reconfigurable system. pKa gradients are known to be integral to protein structures and functions; our work presents the first theoretical demonstration of a pKa gradient for a highly dynamical polysaccharide system and how it allows a persistent but erasable gradient in the structural and mechanical properties of the hydrogel.

65. Henderson JA, Harris RC, Tsai CC, and Shen J*,  How ligand protonation state controls water in protein-ligand binding. J Phys Chem Lett 9: 5440-5444, 2018.

The role of water in protein-ligand has been intensely studied in recent years; however, how ligand protonation state change perturbs water has not been investigated. Here we use two homologues aspartyl proteases to illustrate how modeling protein-ligand binding while allowing ligand titration can further advance our understanding of the role of water in modulating the inhibitor selectivity for a particular target. 

64. Wu S, Yan K, Zhao Y, Tsai CC, Shen J, Bentley WE, Chen Y, Deng H, Du Y, Payne GF* and Shi X*,

Electrical writing onto a dynamically responsive polysaccharide medium: patterning structure and function into a reconfigurable medium. Adv Funct Mater 28: 1803139, 2018.

Historically, paper is the medium to write information. Here, we expand the concept of writing to the addition of electric energy to a dynamically reconfigurable polysaccharide medium. 

63. Shen J*, Zooming in on a small multidrug transporter reveals details of asymmetric protonation. 

Proc Natl Acad Sci USA 115: 8060-8062, 2018.

Nature is efficient; why would it require two protons if one is enough to do the job? This commentary article discusses a recent computational work in the context of our current understanding of how proton translocation mediates the active transport of drugs and substrates across the biological membrane. 

62. D.A. Case, I.Y. Ben-Shalom, S.R. Brozell, D.S. Cerutti, T.E. Cheatham, III, V.W.D. Cruzeiro, T.A. Darden, R.E. Duke, D. Ghoreishi, M.K. Gilson, H. Gohlke, A.W. Goetz, D. Greene, R Harris, N. Homeyer, S. Izadi, A. Kovalenko, T. Kurtzman, T.S. Lee, S. LeGrand, P. Li, C. Lin, J. Liu, T. Luchko, R. Luo, D.J. Mermelstein, K.M. Merz, Y. Miao, G. Monard, C. Nguyen, H. Nguyen, I. Omelyan, A. Onufriev, F. Pan, R. Qi, D.R. Roe, A. Roitberg, C. Sagui, S. Schott-Verdugo, J. Shen, C.L. Simmerling, J. Smith, R. Salomon-Ferrer, J. Swails, R.C. Walker, J. Wang, H. Wei, R.M. Wolf, X. Wu, L. Xiao, D.M. York and P.A. Kollman (2018), AMBER 2018, University of California, San Francisco.

61. Huang YD, Harris RC, Shen J*, Generalized Born based continuous constant pH molecular dynamics in Amber: implementation, benchmarking, and analysis. J Chem Inf Model 58:1372-1383, 2018.

 

We report the CPU implementation of the CpHMD method based on the GBNeck2 generalized Born (GB) implicit-solvent model in the pmemd engine of the Amber molecular dynamics package. Benchmark data based on 11 proteins are discussed. The related software is disseminated through Amber18. 

60. Huang YD, Yue Z, Tsai CC, Henderson JA, and Shen J*,  Predicting catalytic proton donors and nucleophiles in enzymes: how adding dynamics helps elucidate the structure-function relationships. 

J Phys Chem Lett 9: 1179-1184, 2018.

Despite the importance, prediction of catalytic proton donors and nucleophiles has remained an unsolved problem. Here we tested three types of state-of-the-art computational methods to predict the pKa’s of buried and coupled carboxyl dyads in enzymes. We asked the question what makes a residue proton donor vs. nucleophile. Our study revealed a surprisingly simple trend, which is not apparent from crystal structures, empirical or traditional pKa calculations.

59. Yue Z and Shen J*, pH-dependent cooperativity and existence of a dry molten globule in the folding of a mini-protein BBL. Phys Chem Chem Phys 20, 3523-3530, 2018.

pH plays an important role in protein dynamics, stability, and folding; however, detailed mechanisms remain poorly understood. We carried out the first  explicit-solvent protein folding simulation in a pH-dependent manner. Remarkably, the data revealed an unfolding intermediate bearing the signature of a dry molten globule, which is currently hypothesized as a universal intermediate in protein folding.

58. Yan K, Liu Y, Zhang J, Correa SO, Shang W, Tsai CC, Bentley WE, Shen J, Scarcelli G, Raub CB, Shi XW, Payne GF*,

Electrical programming of soft matter: using temporally varying electrical inputs to spatially control self assembly. Biomacromolecules 19: 364-373, 2018. 

Top-down fabrication methods lack the capabilities to generate complex matrix architectures observed in biology. Here we show that temporally varying electric signals can cue a self-assembling polysaccharide to controllably form a hydrogel with complex internal patterns.

57. Yue Z, Chen W, Zgurskaya, H, and Shen J*, Constant pH molecular dynamics reveals how proton release drives the conformational transition of a transmembrane efflux pump. J Chem Theory Comput 13: 6405-6414, 2017.

AcrB is the inner-membrane transporter of the AcrAB-TolC efflux complex in E. coli which confers resistance to antibiotics. Crystal structures have revealed three distinct states in the drug transport cycle; however, detailed mechanism remains unknown. Our study offered convincing support for the 1:1 drug:proton stoichiometry by depicting how proton release triggers the conformational transition, in agreement with the crystal structure.

56. Harris RC, Tsai CC, Ellis CR, and Shen J*, Proton-coupled conformational allostery modulates the inhibitor selectivity for β-secretase. J Phys Chem Lett 8: 4832-4837, 2017.

Many pharmaceutical targets, such as aspartyl proteases and kinases, exhibit pH-dependent dynamics and functions. Accurate prediction of their binding free energies is challenging, and current computational studies neglect the effects of pH and changes in protonation states.  We demonstrated a new protocol to improve the predictive power of current binding free energy calculations. The allosteric mechanism discovered here could lead to exciting opportunities in drug design.        

55. Tsai CC, Morrow BH, Chen W, Payne GP, and Shen J*, Towards understanding the environmental control of hydrogel film properties: how salt modulates the flexibility of chitosan chains. 

Macromolecules 50: 5946-5952, 2017.

Aminopolysaccharides can form complex structures with the properties highly dependent on environmental conditions; however, molecular details are not understood. We found that, contrary to common belief, chitosan chains have a high intrinsic flexibility and salt increases the flexibility by shifting the population from extended to bent glycosidic backbone conformation. 

54. Deredge DJ, Huang W, Huib C, Matsumura H, Yue Z, Moenne-Loccoz P, Shen J, Wintrode PL, and Wilks A*,

Ligand-induced allostery in the interaction of the Pseudomonas aeruginosa heme binding protein with heme oxygenase. Proc Natl Acad Sci USA 114:3421-3426, 2017.

53. Ellis CR, Tsai CC, Lin FY, and Shen J*, Conformational dynamics of cathepsin D and binding to a small-molecule BACE1 inhibitor. J Comput Chem 38:1260-1269, 2017.

Targeting the aspartyl protease BACE1 with small-molecule inhibitors offers a promising route for treatment of Alzheimer's disease. However, developing inhibitors that can selectively target BACE1 in favor of other proteases, especially cathepsin D, has presented significant challenges. We reported the first computational study that characterizes the conformational dynamics of capthesin D and how it binds a BACE1 small-molecule inhibitor.

52. Chen W, Huang YD and Shen J*, Conformational activation of a transmembrane proton channel from constant pH molecular dynamics. J Phys Chem Lett 7: 3961-3966, 2016.

Proton-coupled transmembrane proteins play important roles in human health and diseases. We performed the first computer simulation that 1) directly describes the proton-coupled conformational activation of a transmembrane channel with fully atomic detail; 2) accurately determines the stepwise acid-base constants of a transmembrane channel; 3) provides the proton-coupled free energy of channel activation. The presented methodologies and major findings are generalizable for studies of proton-coupled channels and transporters.

51. Huang YD, Chen W, and Shen J*, Continuous constant pH molecular dynamics with particle-mesh Ewald and titratabe water. J Chem Theory Comput 12: 5411-5421, 2016.

Development of a pH stat to control solution pH in biomolecular simulations remains a outstanding goal. In recent years, the constant pH molecular dynamics methods have emerged; however, the accuracy and generality have been hampered by the use of implicit-solvent models or truncation-based electrostatic schemes. We reported the first implementation of the all-atom PME-based continuous constant pH method and benchmark results on a set of commonly used proteins.

50. Huang YD, Chen W, Dotson DL, Beckstein O, and Shen J*, Mechanism of pH-dependent activation of the sodium-proton antiporter NhaA. Nat Commun 7: 12940, 2016.

NhaA is the principal sodium-proton antiporter in E. Coli., essential for cellular sodium and pH homeostasis. Abundant experimental and computational studies as well as two high-resolution crystal structures have been published; however, the pH-regulated mechanism at the molecular-level remains elusive. We reported the first computational study that offers an atomically-detailed view of the pH-dependent activation and key conformational events in the sodium-proton exchange cycle of NhaA. 

49. Kim E, Liu Y, Ben-Yoav H, Winkler TE, Yan K, Shi S, Shen J, Kelly DL, Ghodssi R, Bentley WE, and Payne GF*, 

Fusing sensor paradigms to acquire chemical information: an integrative role for smart biopolymeric hydrogels. Adv Healthcare Mater 5: 2595-2616, 2016.

48. Ellis CR, Tsai CC, Hou XJ, and Shen J*, Constant pH molecular dynamics reveals pH-modulated binding of two small-molecule BACE1 inhibitors. J Phys Chem Lett 7: 944-949, 2016.

Targeting the aspartyl protease BACE1 with small-molecule inhibitors offers a promising route for treatment of Alzheimer's disease. However, the intricate pH dependence of BACE1 function and inhibitor efficacy has posed major challenges. We reported the first simulation that uncovers the atomic details of pH-dependent binding of two drug candidates for BACE1. Our work demonstrated the drastic contrast in affinity of the two inhibitors is due to the subtle differences in the protonation behavior of the protein-inhibitor complex.

47.  Ellis CR and Shen J*, pH-dependent population shift regulates BACE1 activity and inhibition. 

J Am Chem Soc 137: 9543-9546, 2015.

The aspartyl protease BACE1 is a major pharmaceutical target for treatment of Alzheimer’s disease. Over the past fifteen years, hundreds of inhibitors have been developed some of which have entered clinical trials. However, the fundamental biology of BACE1 including the pH-dependent regulatory mechanism is not well understood. We reported the first molecular simulations to elucidate the mechanism of the pH-dependent enzymatic activity and inhibitor binding of BACE1.

46. Morrow BH, Payne GF, and Shen J*, Stimuli-responsive self-assembly of polysaccharide through a rugged energy landscape. J Am Chem Soc 137: 13024-13030, 2015.

Stimuli-responsive self-assembling materials based on polysaccharide have found a wide range of applications from regenerative medicine to bioelectronics. However, the lack of detailed understanding limits 

technological efforts to access the broad design space. We presented the first all-atom constant pH simulation that describes the simultaneous titration of 240 ionizable sites, while allowing large-scale conformational transitions.

45. Chen W, Shi C, and Shen J*, Nascent beta-hairpin formation of a natively unfolded fragment peptide reveals the role of hydrophobic contacts. Biophys J 109: 630-638, 2015.

The unfolded states play an important role in protein stability and folding. Here we explored the residual structure in the unfolded NTL9. Interestingly, despite the large number of charged residues, the formation of backbone hydrogen bonds is concomitant with hydrophobic but not electrostatic contacts. The data lends support to a long-standing hypothesis that unfolded states exhibit native-like topology with long-range hydrophobic clusters.

44. Chen W, Shi C, MacKerell Jr. AD, and Shen J*, Conformational dynamics of two natively unfolded fragment peptides: Comparison of the AMBER and CHARMM force fields. J Phys Chem B 119: 7902-7910, 2015.

43. Kim E, Xing Y, Cheng Y, Wu H-C, Liu Y, Morrow BH, Ben-Yoav H, Ghodssi R, Rubloff GW, Shen J, Bentley WE, Shi X, and Payne GF*, Chitosan to connect biology to electronics: Fabricating the bio-device interface and communicating across this interface. Polymers 7: 1-46, 2015.