Lecture 17: Prof. Sagar Khare

A New Weapon in the Fight Against Cancer and Viral Infections: Custom-Designed Enzymes

Arising out of natural selection, the structures of proteins (and their complexes with small molecules, nucleic acids, and other proteins) display exquisitely fine-tuned molecular recognition, which is critical for life to operate. Under selection conditions, accurate molecular recognition must be robust to random perturbations such as mutations. Yet, natural proteins are also evolvable — variation in a few amino acids can lead to profound changes in function, e.g. a new enzymatic activity can arise in an “old” enzyme. In other words, these molecular interactions have the fascinating property of being simultaneously functionally robust and plastic.

Fitness scoring function [1].

Characterizing the substrate specificity of protease enzymes is critical for illuminating the molecular basis of their diverse and complex roles in a wide array of biological processes. Rapid and accurate prediction of their extended substrate specificity would also aid in the design of custom proteases capable of selectively and controllably cleaving biotechnologically or therapeutically relevant targets. However, current in silico approaches for protease specificity prediction, rely on, and are therefore limited by, machine learning of sequence patterns in known experimental data. In his talk, Prof. Sagar Khare described a general approach for predicting peptidase substrates de novo using protein structure modeling and biophysical evaluation of enzyme–substrate complexes. His research team constructed atomic resolution models of thousands of candidate substrate–enzyme complexes for each of five model proteases belonging to the four major protease mechanistic classes—serine, cysteine, aspartyl, and metallo-proteases—and develop a discriminatory scoring function using enzyme design modules from Rosetta and AMBER's MMPBSA. They ranked putative substrates based on calculated interaction energy with a modeled near-attack conformation of the enzyme active site. Their results show that the energetic patterns obtained from these simulations can be used to robustly rank and classify known cleaved and uncleaved peptides and that these structural-energetic patterns have greater discriminatory power compared to purely sequence-based statistical inference. Combining sequence and energetic patterns using machine-learning algorithms further improves classification performance, and analysis of structural models provides physical insight into the structural basis for the observed specificities.

Summary of the CPG2 circular permutations [2].

In the second part of his talk, Prof. Khare talked about spatio-temporal design of enzymes. Carboxypeptidase G2 (CPG2) is an Food and Drug Administration (FDA)-approved enzyme drug used to treat methotrexate (MTX) toxicity in cancer patients receiving MTX treatment. It has also been used in directed enzyme-prodrug chemotherapy, but this strategy has been hampered by off-site activation of the prodrug by the circulating enzyme. The development of a tumor protease activatable CPG2, which could be achieved using a circular permutation of CPG2 fused to an inactivating ‘prodomain’, would aid in these applications. The research team reported the development of a protease accessibility-based screen to identify candidate sites for circular permutation in proximity of the CPG2 active site. The resulting six circular permutants showed similar expression, structure, thermal stability, and, in four cases, activity levels compared to the wild-type enzyme. They rationalize these results based on structural models of the permutants obtained using the Rosetta software by developing a cell growth-based selection system, and demonstrated that when fused to periplasm-directing signal peptides, one of the circular permutants confers MTX resistance in Escherichiacoli with equal efficiency as the wild-type enzyme. As the permutants have similar properties to wild-type CPG2, these enzymes are promising starting points for the development of autoinhibited, protease-activatable zymogen forms of CPG2 for use in therapeutic contexts[2]. 

Related References:

[1] Pethe, M. A., Rubenstein, A. B., & Khare, S. D. (2017). Large-scale Structure-based Prediction and Identification of Novel Protease Substrates using Computational Protein Design. Journal of Molecular Biology, 429(2), 220-236

[2] Yachnin, B. J., & Khare, S. D. (2017). Engineering carboxypeptidase G2 circular permutations for the design of an autoinhibited enzyme. Protein engineering, design & selection: PEDS, 1.

____________________________________________________________________________

Sagar Khare is an Assistant Professor at Rutgers University. He teaches Chemistry and Chemical Biology and seeks to understand the structural determinants of enzymatic specificity and reactivity using a combination of computational protein design and experimental characterization. His research team's goal is to develop a quantitative and predictive understanding of specificity at protein-ligand and protein-peptide interfaces; which will inform various therapeutic and synthetic applications. Prior to working at Rutgers, Prof. Sagar Khare was a Postdoctoral Fellow at the University of Washington in Seattle, Washington. He was also a Software Engineer for Affymax Research Institute in Bangalore, India. 

Lecture 16: Dr. Yukun Wang (Martin Ulmschneider Lab)

Experimentally guided simulations reveal the mechanism of action of antimicrobial peptides

Since their discovery over 100 years ago, thousands of pore-forming antimicrobial peptides (AMPs) have been identified, revealing great variations in size, secondary structure and sequence composition. Despite this wealth of data no common pore-forming motif has been discovered, and the molecular basis of antimicrobial activity remains poorly understood. One of the key difficulties has been the transience of pores formed by AMPs, which has proved challenging for experimental structure determination. In the absence of experimental structures a range of theoretical models have been proposed that try to rationalize how soluble AMPs, which carry multiple charged residues, can insert into hydrophobic membranes to form pores and what these pores might look like.

Spontaneous membrane insertion and intra-membrane
oligomerization of maculatin [1]

Scientists from the Johns Hopkins Institue of NanoBio technology reports[1] an experimentally guided unbiased simulation methodology that yields the mechanism of spontaneous pore assembly for the AMP maculatin at atomic resolution. Rather than a single pore, maculatin forms an ensemble of structurally diverse temporarily functional low-oligomeric pores, which mimic integral membrane protein channels in structure. These pores continuously form and dissociate in the membrane. Membrane permeabilization is dominated by hexa-, hepta- and octamers, which conduct water, ions and small dyes. Pores were shown to form by consecutive addition of individual helices to a transmembrane helix or helix bundle, in contrast to current poration models. The diversity of the pore architectures—formed by a single sequence—may be a key feature in preventing bacterial resistance and could explain why sequence–function relationships in AMPs remain elusive.

Related References:

[1] Wang, Y., Chen, C. H., Hu, D., Ulmschneider, M. B., & Ulmschneider, J. P. (2016). Spontaneous formation of structurally diverse membrane channel architectures from a single antimicrobial peptide. Nature Communications, 7.

[2] Wimley, W. C., & Hristova, K. (2011). Antimicrobial peptides: successes, challenges and unanswered questions. The Journal of membrane biology, 239(1-2), 27-34.

[3] Ulmschneider, M. B., Ulmschneider, J. P., Schiller, N., Wallace, B. A., Von Heijne, G., & White, S. H. (2014). Spontaneous transmembrane helix insertion thermodynamically mimics translocon-guided insertion. Nature communications, 5, 4863.

_______________________________________________________________________________

Dr. Yukun Wang did his undergrad in Biology in Nanjing Agricultural University. After graduation, he did his Ph.D. in Biophysics at the Institute of Natural Sciences & School Of Life Sciences And Biotechnology, Shanghai Jiao Tong University. Currently, he is a postdoctoral fellow at the Institute for Nano-Bio-Technology in the Whiting School of Engineering at The Johns Hopkins University working with Prof. Martin Ulmschneider. Dr. Wang's research interests span from the functional mechanism of antimicrobial peptides to developing free energy calculation algorithms for rare event dynamics in the lipid-peptides systems. He enjoys computational designing of membrane-active peptides and proteins.

Lecture 15: Prof. Alex MacKerell

From Lemkul, J. A. et. al., Chemical reviews, 2016 [2].

Recent developments in polarizable force field research
Atomistic molecular dynamics (MD) simulations have become an integral component of the tool set used to examine biomolecular systems. Since the first MD simulation of a protein in 1977 [1], the time scales and sizes of computationally tractable simulations have grown by orders of magnitude. Systems may now be simulated on the microsecond or millisecond time scale and contain over a million atoms due to increased computing power and ever-improving algorithms for parallel and graphical processing unit (GPU)-based computing. Such capabilities also allow for more rigorous testing of the accuracy of the models used for such MD simulations.
Molecular mechanics force fields that explicitly account for induced polarization represent the next generation of physical models for molecular dynamics simulations. Several methods exist for modeling induced polarization, and in this seminar, Prof. MacKerell reviewed the classical Drude oscillator model, in which electronic degrees of freedom are modeled by charged particles attached to the nuclei of their core atoms by harmonic springs. He described the latest developments in Drude force field parametrization and application, primarily in the last 15 years. Emphasis was placed on the Drude-2013 polarizable force field for proteins, DNA, lipids, and carbohydrates. He discussed its parametrization protocol, development history, and recent simulations of biologically interesting systems, highlighting specific studies in which induced polarization plays a critical role in reproducing experimental observables and understanding physical behavior. As the Drude oscillator model is computationally tractable and available in a wide range of simulation packages, it is anticipated that use of these more complex physical models will lead to new and important discoveries of the physical forces driving a range of chemical and biological phenomena.
Related References:
[1] McCammon, J. A., Gelin, B. R., & Karplus, M. (1977). Dynamics of folded proteins. Nature, 267(5612), 585.
[2] Lemkul, J. A., Huang, J., Roux, B., & MacKerell Jr, A. D. (2016). An empirical polarizable force field based on the classical drude oscillator model: development history and recent applications. Chemical reviews, 116(9), 4983-5013.

__________________________________________________________________________

Alex MacKerell received an A.S. in biology in 1979 from Gloucester County College, Sewell, NJ, followed by a B.S. in chemistry in 1981 from the University of Hawaii, Honolulu, HI, and a Ph.D. in biochemistry in 1985 from Rutgers University, New Brunswick, NJ. Subsequent training involved postdoctoral fellowships in the Department of Medical Biophysics, Karolinska Intitutet, Stockholm, Sweden, in the area of experimental and theoretical biophysics and in the Department of Chemistry, Harvard University in theoretical chemistry. Following one year as a visiting professor at Swarthmore College, Swarthmore, PA, he assumed his faculty position in the School of Pharmacy, University of Maryland, Baltimore in 1993. MacKerell is currently the Grollman-Glick Professor of Pharmaceutical Sciences in the School of Pharmacy and the Director of the University of Maryland Computer-Aided Drug Design Center. MacKerell is also cofounder and Chief Scientific Officer of SilcsBio LLC. Research interests include the development of theoretical chemistry methods, with emphasis on empirical force field development; structure–function studies of proteins, carbohydrates, and nucleic acids; and the application of theoretical methods to drug discovery.

Lecture 14: David Holland (Margaret Johnson Lab)

Finding the right partner in a crowded world

Mammalian cells contain ∼21,000 genes encoding upwards to 100,000 protein types. 5-40% of cell volume is occupied by macromolecules, posing challenges for cell proteins to locate functional partners and increasing the risk of nonspecific (nonfunctional) interactions. Overexpressed proteins, in particular, will saturate functional partners, leaving leftovers for nonspecific binding instead. Eukaryotic cells have evolved various methods to help proteins function reliably, including compartmentalization, allostery, and structural and chemical properties of binding sites. Cells may optimize specificity as a function of their binding networks; first through concentration balance, next through network structure. Using the Gillespie algorithm, David Holland and Dr. Margaret Johnson from Johns Hopkins Biophysics simulated specific and nonspecific binding in 500 networks of 90-200 nodes with varying topological properties under equal, random, and stoichiometrically balanced protein concentrations. Binding affinities for all specific and nonspecific interactions were determined using a coarse-grained protein sequence model. The research team found out that the concentration balance significantly reduced the number of nonspecific interactions, as did local topological patterns (motifs) that allowed increased difference between specific and nonspecific binding affinities. 

Schematic illustrating (left) the binding of three proteins, where matching features indicate binding interfaces, (center) the corresponding interface-interaction network with orange and green indicating shared and non-shared binding interfaces, respectively, and (right) the protein-protein interaction network (from [1]).

To test if these motifs are selected for in real networks, they sampled possible interface-interaction networks (IINs) for the Clathrin-mediated endocytosis PPI network in yeast, using a Monte Carlo algorithm and a fitness function to select for features that allow increased specificity. The fitness function reproduced several features of the real IIN, including fragmentation, a scale-free degree distribution, the presence of square motifs and a low number of chain and triangle motifs. These features differed significantly from unbiased sampling. In conclusion, there is selective pressure on the evolution of real IINs to avoid nonfunctional interactions and that non-optimal features of the real IIN (chains, larger network components) likely serve a functional purpose.

Relevant Publications:

[1]Johnson, ME, & G Hummer (2013) “Evolutionary Pressure on the Topology of Protein Interface Interaction Networks.” J Phys Chem B 117:13098-106. 

[2]Johnson, ME, G Hummer (2013). “Interface-resolved network of protein-protein interactions.” PLoS Comput Biol 9(5):e1003065 

[3]Keil, C, E Verschueren, J Yang, & L Serrano (2013) “Integration of protein abundance and structure data reveals competition in the ErbB signaling network.” Science Signaling 6(306): ra109 [4]Johnson, ME, & G Hummer (2011) “Nonspecific binding limits the number of proteins in a cell and shapes their interaction networks.” PNAS 108(2):603-8. 

[5]Zhang, J, S Maslov, & EI Shakhnovich (2008) “Constraints imposed by non-functional protein-protein interactions on gene expression and proteome size.” Mol Sys Biol 4:210 

[6]Vavouri, T, JI Semple, R Garcia-Verdugo, & B Lehner (2009) “Intrinsic protein disorder and interaction promiscuity are widely associated with dosage sensitivity.” Cell 138: 198-208.

_________________________________________________________________________________

David Holland received his B.Sc. in Biomedical Engineering from the University of Virginia in 2011. After graduating, he joined the department of Biomedical Engineering at Johns Hopkins. He is currently earning his Ph.D. under Margaret Johnson where he studies the effects of protein abundance and protein-protein interaction network structure on protein mis-interactions. In his spare time, David practices taekwondo and is also teaching a course on network science. 

Lecture 13: Athena Chen (Margaret Johnson Lab)

Spatial Cell Modeling Methods 

Due to the complexity of cells, it is useful to use computational tools to understand and predict mechanisms of biological processes. Ordinary differential equations (ODEs) and partial differential equations (PDEs) have proven to be successful at modeling various large systems. However, to understand certain biological processes such as bacterial cell division, high spatial resolution is necessary to discern the underlying mechanisms and interactions. ODEs and PDEs, along with the Gillespie algorithm for stochastic modeling, do not provide spatial resolution at a single-particle level; they are all concentration-based methods. MCell, the Free Propagator Reweighting Algorithm (FPR), and Smoldyn are algorithms that provide single-particle resolution, but may be computationally expensive and may not necessarily yield accurate protein dynamics.

Change in concentration of molecule A in an irreversible 3D reaction A+A -> 0. In this parameter set, the initial distances between molecules causes an increased initial reaction rate.

In her recent work, Athena Chen and her mentors Dr. Margaret Johnson and Dr. Osman Yogurtcu in the Johns Hopkins Biophysics department, analyzed the strengths and limitations of ODEs, PDEs, Gillespie, MCell, Smoldyn, and FPR through establishing and performing a set of benchmark tests. Though all modeling methods extracted the correct equilibrium concentrations for most of the tested reactions, the resulting protein dynamics were not necessarily correct. Simulations at high rates and large densities showed that despite providing single-particle resolution, Smoldyn and MCell do not pick up single-particle effects where the distance between two molecules affects the probability of binding. Furthermore, the dynamics given by Smoldyn and MCell are dependent on the time step selected. On the other hand, FPR correctly identified single-particle effects and yielded dynamics independent of the selected timestep.

As an example of the effects of molecular geometry, diffusion, and stochasticity of protein dynamics, we examined a model for bacterial cell division. From oscillations in protein concentrations, the cell can identify the center of the cell to ensure identical offspring and division of genetic information. 

Written by Athena Chen

Relevant Articles:

1-Yogurtcu, Osman N., and Margaret E. Johnson. "Theory of bi-molecular association dynamics in 2D for accurate model and experimental parameterization of binding rates." The Journal of chemical physics 143.8 (2015): 084117.

2-Andrews, Steven S., et al. "Detailed simulations of cell biology with Smoldyn 2.1." PLoS Comput Biol 6.3 (2010): e1000705.

3-Johnson, Margaret E., and Gerhard Hummer. "Free-Propagator Reweighting Integrator for Single-Particle Dynamics in Reaction-Diffusion Models of Heterogeneous Protein-Protein Interaction Systems." Physical Review X 4.3 (2014): 031037.

4-Kerr, Rex A., et al. "Fast Monte Carlo simulation methods for biological reaction-diffusion systems in solution and on surfaces." SIAM journal on scientific computing 30.6 (2008): 3126-3149.

_________________________________________________________________________________

Athena Chen is currently working on her Bachelors of Arts in Biophysics and Bachelors of Science in Applied Mathematics and Statistics from Johns Hopkins University. As part of Dr. Margaret Johnson’s lab in the department of Biophysics, she analyzes the accuracy of methods for modeling the dynamics of protein interactions. In her free time, she enjoys yoga and figure skating. 

Lecture 12: Andrei Kucharavy (Rong Li Lab)

Cellular Adaptation Under Stress

Diagram for a cell population adaptation.

Whether the products of human activity, or naturally occurring social, economic or biological, complex systems share the same properties. Composed of a large number of individual components, their components do not have a straightforward relation to their properties and often interact one with another in unexpected ways. Because of that, different instances of the same complex systems are built from slightly different components. Such differences give rise to heterogeneity within a population, which in turn raises significant difficulties for their study. From the biological perspective, such events have been formalized as Fisher’s geometric model, that has been formalized in the thirties of the last century and has been independently re-discovered in unrelated domains as algorithms for ergodic explorations of multi-dimensional spaces for an optimal value function point. 

Andrei Kucharavy and Dr. Rong Li from Johns Hopkins Medicine propose an enhancement of Fisher’s geometric model, allowing to explain a range of previously unexplained observations in biology. Mathematical analysis of their enhancement provides a set of rules applicable to the optimization of a large class of ergodic exploration algorithms.

Related Journal Articles:
1. H. A. Orr, The genetic theory of adaptation: a brief history. Nat. Rev. Genet. 6, 119–127 (2005) 
2. H. A. Orr, R. L. Unckless, The Population Genetics of Evolutionary Rescue. PLoS Genet. 10, e1004551 (2014). 
3. P. S. Pennings, Standing genetic variation and the evolution of drug resistance in HIV. PLoS Comput. Biol. 8 (2012).
_________________________________________________________________________________

Andrei Kucharavy received his Engineer’s Degree in Physics, Mathematics, Programming and Bioinformatics from Ecole Polytechnique, France in 2011. After graduation, he did masters in computational biology at Ecole Polytechnique Fédérale de Lausanne, Switzerland. Currently, Andrei is a Ph.D. student under the joint direction of Dr. Rong Li from Johns Hopkins and Dr. Gilles Fischer, exploring molecular mechanisms of aneuploidy-enabled stress adaptation and drug resistance it enables in yeast and cancer. He works on developing computational methods enabling systematic analysis of molecular mechanisms underlying complex traits. The interface of biological network analysis and evolution theory is his particular interest.

Lecture 11: Dr. Ana Damjanović

Simulating with the right pH.

Schematic representation of a constant-pH simulation
with the two-dimensional EDS-HREM method.

Solution pH is one of the most important environmental factors that affects the structure and dynamics of proteins. Almost all biologically relevant properties of proteins are affected by pH: stability, folding and assembly, interactions with ligands and other biological molecules, solubility, aggregation properties, and enzymatic activity. A change in pH may induce a change in the protonation state of ionizable groups, which in turn can cause structural changes in proteins. Structural changes triggered by protonation/deprotonation can be exploited for function, such as in the case of ATP synthase, bacteriorhodopsin, cytochrome c oxidase, or the photoactive yellow protein.

Superimposed 1 ns trajectories of snake cardiotoxin
from the (A) 2D and (B) 1D constant-pH
EDS-HREM simulations at pH = 2.

Ana Damjanovic from Johns Hopkins Biophysics and her colleagues from NIH present a new method for enhanced sampling for constant-pH simulations in explicit water based on a two-dimensional (2D) replica exchange scheme. The new method is a significant extension of a previously developed constant-pH simulation method, which is based on enveloping distribution sampling (EDS) coupled with a one-dimensional (1D) Hamiltonian exchange method (HREM). EDS constructs a hybrid Hamiltonian from multiple discrete end state Hamiltonians that, in this case, represent different protonation states of the system. The ruggedness and heights of the hybrid Hamiltonian’s energy barriers can be tuned by the smoothness parameter. Within the context of the 1D EDS-HREM method, exchanges are performed between replicas with different smoothness parameters, allowing frequent protonation-state transitions and sampling of conformations that are favored by the end-state Hamiltonians. In this work, the 1D method is extended to 2D with an additional dimension, external pH. Within the context of the 2D method (2D EDS-HREM), exchanges are performed on a lattice of Hamiltonians with different pH conditions and smoothness parameters. The research team demonstrates that both the 1D and 2D methods exactly reproduce the thermodynamic properties of the semi-grand canonical (SGC) ensemble of a system at a given pH. They have tested the new 2D method on aspartic acid, glutamic acid, lysine, a four-residue peptide (sequence KAAE), and snake cardiotoxin. In all cases, the 2D method converges faster and without loss of precision; the only limitation is a loss of flexibility in how CPU time is employed. The results for snake cardiotoxin demonstrate that the 2D method enhances protonation-state transitions, samples a wider conformational space with the same amount of computational resources, and converges significantly faster overall than the original 1D method.

_________________________________________________________________________________

Ana Damjanovic received her B.Sc. in Physics from Belgrade University (1995), and a Ph.D. in Physics from University of Illinois at Urbana-Champaign (2001). For her Ph.D. she worked with Prof. Klaus Schulten on QM description of energy transfer in light-harvesting complexes in various photosynthetic organisms. She did postdocs at UC Berkeley and Johns Hopkins University. At JHU she worked on understanding hydration and conformational changes in proteins through molecular dynamics simulations. She is presently an Associate Research Scientist and Lecturer in the Dept. of Biophysics at Johns Hopkins University. Her current research interests are in the area of development and applications of molecular dynamics simulations at constant pH.