Steven Plotkin


Relevant Degree Programs

Affiliations to Research Centres, Institutes & Clusters


Graduate Student Supervision

Doctoral Student Supervision (Jan 2008 - April 2022)
Protein-Solvent Interactions and Classical Density Functional Theory (2016)

We use classical density functional theory to investigate the interactions between solvents and proteins. We examine a diverse experimental literature to establish thermodynamic properties of protein-cosolute interaction, particularly the compensation between transfer entropy and transfer enthalpy. We develop a method of analysing the uncertainties in such measurements and use the method to resolve a long-standing debate over entropy-enthalpy compensation. We develop a classical density functional theory for interactions between proteins and cosolutes. The theory developed here ignores the solvent-solvent interaction but is nonetheless quite accurate. We use this approach to reproduce transfer free energies reported elsewhere, and show that the cDFT model captures the desolvation barrier and the temperature dependence of the transfer free energy. We use experimental values that we have analyzed to define the parameter space of a model density functional theory approach. We then extend the classical density functional theory to capture protein-water interactions, thus developing a new implicit solvent model. Along the way we give a proof that the free energy of a bath of particles in a finite external potential is independent of the external potential in the isothermal-isobaric ensemble. We finally discuss the challenges remaining in implementing our implicit solvent model.

View record

Emergence of Riemannian manifolds from graphs and aspects of Chern-Simons theory (2013)

This thesis is divided into three chapters. Chapter 1 is a study of statistical models of graphs, in order to explore possible realizations of emergent manifolds. Graphs with given numbers of vertices and edges are considered, governed by a Hamiltonian that favors graphs with near-constant valency and local rotational symmetry. The model is simulated numerically in the canonical ensemble. It is found that the model exhibits a first-order phase transition, and that the low energy states are almost triangulations of two dimensional manifolds. The resulting manifold shows topological "handles" and surface intersections in a higher embedding space as well as non-trivial fractal dimension. The model exhibits a phase transition temperature of zero in the bulk limit. We explore the effects of adding long-range interactions to the model, which restore a finite transition temperature in the bulk limit.In Chapter 2, aspects of Chern-Simons theory are studied. The relations between Chern-Simons theory, a model known as BF theory named after the fields that appear in the actions, and 3D gravity, are explored and generalized to the case of non-orientable spacetime manifolds. U(1) Chern-Simons theory is quantized canonically on orientable manifolds, and U(1) BF theory is similarly quantized on non-orientable manifolds. By requiring the quantum states to form a representation of the deformed holonomy group and the deformed large gauge transformation group, we find that the mapping class group of the spacetime manifold can be consistently represented, provided the prefactor k of the Chern-Simon action satisfies quantization conditions which in general are non-trivial. We also find a k 1/k duality for the representations. Motivated by open questions about interpreting the finite size results from Chapter 1, models of finite size scaling for systems with a first-order phase transition are discussed in Chapter 3. Three physics models -- the Potts model, the Go model for protein folding, and the graph model in Chapter 1 -- are simulated. Several finite size scaling models, including three functional forms to fit the energy distributions, and a capillarity model, are compared with simulations of the corresponding physics models.

View record

Generalized distance and applications in protein folding (2013)

The Euclidean distance, D, between two points is generalized to the distance between strings or polymers. The problem is of great mathematical beauty and very rich in structure even for the simplest of cases. The necessary and sufficient conditions for finding minimal distance transformations are presented. Locally minimal solutions for one-link and two-link chains are discussed, and the large N limit of a polymer is studied. Applications of D to protein folding and structural alignment are explored, in particular for finding minimal folding pathways. Non-crossing constraints and the resulting untangling moves in folding pathways are discussed as well. It is observed that, compared to the total distance, these extra untangling moves constitute a small fraction of the total movement. The resulting extra distance from untangling movements (Dnx ) are used to distinguish different protein classes, e.g. knotted proteins from unknotted proteins. By studying the ensembles of untangling moves, dominant folding pathways are constructed for three proteins, in particular a knotted protein. Finally, applications of D, and related metrics to protein folding rate prediction are discussed. It is seen that distance metrics are good at predicting the folding rates of 3-state folders.

View record

Template-directed protein misfolding in neurodegenerative disease (2013)

Protein misfolding diseases represent a large burden to human health for which only symptomatictreatment is generally available. These diseases, such as Creutzfeldt-Jakob disease, amyotrophiclateral sclerosis, and the systemic amyloidoses, are characterized by conversion of globular, nativelyfoldedproteins into pathologic β-sheet rich protein aggregates deposited in affected tissues. Understandingthe thermodynamic and kinetic details of protein misfolding on a molecular level dependson accurately appraising the free energies of the folded, partially unfolded intermediate,and misfolded protein conformers. There are multiple energetic and entropic contributions to thetotal free energy, including nonpolar, electrostatic, solvation, and configurational terms. To accuratelyassess the electrostatic contribution, a method to calculate the spatially-varying dielectricconstant in a protein/water system was developed using a generalization of Kirkwood Frohlich theoryalong with brief all-atom molecular dynamics simulations. This method was combined withpreviously validated models for nonpolar solvation and configurational entropy in an algorithm tocalculate the free energy change on partial unfolding of contiguous protein subsequences. Resultswere compared with those from a minimal, topologically-based Gō model and direct calculationof free energies by steered all-atom molecular dynamics simulations. This algorithm was appliedto understand the early steps in the misfolding mechanism for β₂-microglobulin, prion protein,and superoxide dismutase 1 (SOD1). It was hypothesized that SOD1 misfolding may follow atemplate-directed mechanism like that discovered previously for prion protein, so misfolding ofSOD1 was induced in cell culture by transfection with mutant SOD1 constructs and observed tostably propagate intracellularly and intercellularly much like an infectious prion. A defined minimalassay with recombinant SOD protein demonstrated the sufficiency of mutant SOD1 aloneto trigger wtSOD1 misfolding, reminiscent of the “protein-only” hypothesis of prion spread. Finally,protein misfolding as a feature of disease may extend beyond neurodegeneration and amyloidformation to cancer, in which derangement of protein folding quality control may lead to antibodyrecognizablemisfolded protein present selectively on cancer cell surfaces. The evidence for thishypothesis and possible therapeutic targets are discussed as a future direction.

View record

The effect of osmolytes on protein folding (2011)

Living cells are composed of a variety of biological macromolecules such as nucleic acid, metabolites, proteins and cytoskeletal filaments as well as other particles. The fraction of the cellular interior volume that is taken by these biomolecules is about 30%, leading to a highly crowded environment. Biomolecules present in an extremely dense environment inside a cell have a completely different set of kinetic and thermodynamic behavior than in a test tube. Therefore comprehending the effect of crowding conditions on biological molecules is crucial to broad research fields such as biochemical, medical and pharmaceutical sciences. Experimentally, we are able to mimic such crowded environments; which are of more physiological relevance, by adding high concentrations of synthetic macromolecules into uncrowded buffers. Theoretically, very little attention has been paid to the effects of the dense cellular cytoplasm on biological reactions. The purpose of this work is to investigate analytically the effects of crowding agents on protein folding and stability. We present a new parameter as the measure of the polymer size, which will substitute the traditional measurements of the radius of gyration of the polymer and the end to end distance of a polymeric chain. Using this quantity we derive an expression for the free energy of the polymer which can easily be generalized to provide the free energy of a protein. This mechanism enables us to study the effect of crowding on folding and stability of a protein. The stabilization effects of the crowding particles depend on the concentration and the size of the crowders and also the type of the crowding particles that are present in the system. In our calculations the type of the crowders is controlled by the energetic parameter between the protein and surrounding macromolecules.

View record

Master's Student Supervision (2010 - 2021)
Rational computational design of epitopes and antibodies for neurodegenerative disease and COVID-19 (2020)

Protein computational design uses current knowledge in structural biology and statistical tools to predict amino acid sequences that exhibit targeted properties. In this study, two protein designs and one statistical tool are developed.The first designed protein family in this thesis is chimeric antibodies for Covid-19 and future coronavirus variants. The chimeric antibodies are composed of an IgG1 framework with "ACE2 units" grafted on complementarity-determining regions. ACE2 units were small protein fragments built around the spike-interacting regions of ACE2. Such a chimeric construct is designed to neutralize SARS-COV-2 by binding spike receptor binding domain and is expected to be tolerant to receptor binding domain (RBD) mutations, as long as ACE2 recognition is required for infection. The binding free energy of ACE2 units to the spike RBD was assessed by molecular dynamics simulation. Surprisingly, the computation result showed that some ACE2 units had similar or even stronger RBD binding than full length ACE2. Moreover, it adds validity to the simulations that the calculated binding free energy between the ACE2 and SARS-COV-2 RBD, -52.9 +- 5.0 kJ/mol, is within the range of the experimental results.The second designed protein in this thesis is an immunogen scaffold design for neurodegenerative disease using cyclic peptides. Effectively scaffolding epitopes on immunogens, in order to raise conformationally selective antibodies through active immunization, is a central problem in treating protein misfolding diseases. We seek to selectively target conformations enriched in toxic, oligomeric propagating species while sparing healthy forms of the protein which are often more abundant. To this end, we scaffolded cyclic peptides by varying the number of flanking glycines, to best mimic a misfolding-specific conformation of an epitope of alpha-synuclein enriched in the oligomer ensemble. The cyclic peptide scaffolds of alpha-synuclein are screened in silico based on their ensemble overlap properties with the fibril, oligomer-model, and isolated monomer ensembles.Lastly, a simulation tool, reservoir replica exchange molecular dynamics simulation (R-REMD), was implemented in GROMACS software. The enhanced sampling of the R-REMD was tested on several systems including a cyclic peptide scaffold whose structural ensemble can predict the selectivity of the raised antibody.

View record

Computational study of a quantum mechanical based model of magnetoreception (2012)

In this work we study a quantum mechanical model of magnetoreception in birds based on a pair of electrons undergoing coherent evolution. The magnetoreception arises because there is a different outcome depending on whether the electrons form a singlet or a triplet state, and a magnetic field can influence this configuration. We perform a variety of computational simulations on a simplified model of the above mechanism, with the purpose to study the effect of changing the different variables in which the process depends, ultimately to determine the plausibility of the mechanism.

View record

News Releases

This list shows a selection of news releases by UBC Media Relations over the last 5 years.

If this is your researcher profile you can log in to the Faculty & Staff portal to update your details and provide recruitment preferences.


Sign up for an information session to connect with students, advisors and faculty from across UBC and gain application advice and insight.