Many histidine-rich sites in proteins bind transition metal ions such as Zn2+, natively. Such sites can encourage proper protein folding or allow access to enzymatic capabilities such as hydrolysis. Ru(II) and Tc(I) also bind to aromatic amines providing access to unique chemistries not observed in biology. Ru(II) complexes have shown efficacy in fighting cancer and catalysis, while 99mTc complexes are used in radio-imaging. To incorporate such metal-ions' activities into proteins, several mutants have been designed to bind Zn2+, Ru2+, and Tc1+ by introducing three histidines onto their surfaces. The first design, Z0, utilized a chimeric approach by substituting a turn in engrailed homeodomain for the superimposable Zn2+-binding loop of astacin. In the second design, 3HT-C, three histidine residues were incorporated into the N-terminus of the Trp-cage. The final scaffold, ubiquitin, was used to make two mutants: 3HIU with a 3-histidine containing loop inserted between residues 9 and 10, and 3HPU with 3-histidine point mutations near residues 35-38. Z0 proved unstable due to incorporation of a hydrophobic patch onto its surface and was not able to be isolated in sufficient quantities for study. However, the other proteins were stable and soluble. Zn2+-binding by 3HT-C was investigated by intrinsic tryptophan fluorescence quenching, circular dichroism, and RP-HPLC. Binding by 3HIU and 3HPU was studied by CD. All designed proteins bind to Zn2+ with Kd values in the micromolar range. 3HIU and 3HPU were further studied for their ability to bind Ru(tacn)2+ complexes. While addition of Ru-complexes caused oligomerization to various extents depending upon reaction conditions, homogeneous Ru-protein monomers were purified by a combination of size-exclusion, cation-exchange and immobilized-metal affinity chromatographies. Ru-binding was confirmed by ESI-MS, and structural integrity was investigated by CD. Results indicate that Ru(tacn)2+ complexes can be bound to surface binding sites in proteins without disruption of structure, opening the door for the study of catalysis in a protein context. 99mTc(CO)3+-binding studies were performed with 3HIU by RP-HPLC. This protein binds Tc+ and resists substitution by free L-histidine, suggesting that peptidic Tc-binding tags could be designed with this approach to readily incorporate the radionuclide into any protein expression system.

The prediction of energetically favorable side-chain conformations is a fundamental element in homology modeling of proteins and the design of novel protein sequences. The space of side-chain conformations can be approximated by a discrete space of probabilistically representative side-chain conformations (called rotamers). The problem is, then, to find a rotamer selection for each amino acid that minimizes a potential energy function. This is called the Global Minimum Energy Conformation (GMEC) problem. This problem is an NP-hard optimization problem. The Dead-End Elimination theorem together with the A* algorithm (DEE/A*) has been successfully applied to this problem. However, DEE fails to converge for some complex instances. In this paper, we explore two alternatives to DEE/A* in solving the GMEC problem. We use a probabilistic inference method, the max-product (MP) belief-propagation algorithm, to estimate (often exactly) the GMEC. We also investigate integer programming formulations to obtain the exact solution. There are known ILP formulations that can be directly applied to the GMEC problem. We review these formulations and compare their effectiveness using CPLEX optimizers. We also present preliminary work towards applying the branch-and-price approach to the GMEC problem. The preliminary results suggest that the max-product algorithm is very effective for the GMEC problem. Though the max-product algorithm is an approximate method, its speed and accuracy are comparable to those of DEE/A* in large side-chain placement problems and may be superior in sequence design. / Singapore-MIT Alliance (SMA)

Calcium regulates many biological processes through interaction with proteins with different conformational, dynamic, and metal binding properties. Previous studies have shown that the electrostatic environment plays a key role in calcium binding affinity. In this research, we aim to dissect the contribution of the electrostatic environment to calcium binding affinity using protein design. Many natural calcium binding proteins undergo large conformational changes upon calcium binding which hampers the study of these proteins. In addition, cooperativity between multiple calcium binding sites makes it difficult to study site-specific binding affinity. The design of a single calcium binding site into a host system eliminates the difficulties that occur in the study of calcium binding affinity. Using a computer algorithm we have rationally designed several calcium binding sites with a pentagonal bipyramidal geometry in the non-calcium dependent cell adhesion protein CD2 (CD2-D1) to better investigate the key factors that affect calcium binding affinity. The first generation proteins are all in varying electrostatic environments. The conformational and metal binding properties of each of these designed proteins were analyzed. The second generation designed protein, CD2.6D79, was designed based on criteria learned from the first generation proteins. This protein contains a novel calcium binding site with ligands all from the â-strands of the non-calcium dependent cell adhesion protein CD2. The resulting protein maintains native secondary and tertiary packing and folding properties. In addition to its selectivity for calcium over other mono and divalent metal ions, it displays strong metal binding affinities for calcium and its analogues terbium and lanthanum. Furthermore, our designed protein binds CD48, the ligand binding partner of CD2, with an affinity three-fold stronger than CD2. The electrostatic potential of the calcium binding site was modified through mutation to facilitate the study of the effect of electrostatic interactions on calcium binding affinity. Several charge distribution mutants display varying metal binding affinities based on their charge, distance to the calcium binding site, and protein stability. This study will provide insight into the key site factors that control calcium binding affinity and calcium dependent biological function.

<p>Proteins are involved in all of life's processes and are also responsible for many diseases. Thus, engineering proteins to perform new tasks could revolutionize many areas of biomedical research. One promising technique for protein engineering is computational structure-based protein design (CSPD). CSPD algorithms search large protein conformational spaces to approximate biophysical quantities. In this dissertation we present new algorithms to realistically and accurately model how amino acid mutations change protein structure. These algorithms model continuous flexibility, protein ensembles and positive/negative design, while providing guarantees on the output. Using these algorithms and the OSPREY protein design program we design and apply protocols for three biomedically-relevant problems: (i) prediction of new drug resistance mutations in bacteria to a new preclinical antibiotic, (ii) the redesign of llama antibodies to potentially reduce their immunogenicity for use in preclinical monkey studies, and (iii) scaffold-based anti-HIV antibody design. Experimental validation performed by our collaborators confirmed the importance of the algorithms and protocols.</p> / Dissertation

Transient change of cytosolic calcium level leads to physiological actions, which are modulated by the intracellular calcium stores, and gated by membrane calcium channels/pumps. To closely monitor calcium dynamics there is a pressing need to develop calcium sensors that are targeted to high calcium environment such as the ER/SR with relatively low binding affinity and fast kinetic properties to complement the current calcium indicator toolkits. In this dissertation, the development of fast red florescent calcium binding protein using the protein design is reported. The results show the calcium dependent fluorescence increase of mCherry mutant MCD1 (RapidER) and MCD15 (RapidER’) is able to monitor the ER calcium release in several cell lines responding to perturbations of extracellular calcium signaling. The specific targeting to the ER membrane was achieved by fusing the ryanodine receptor 1 transmembrane domains for the spatio-temporal calcium imaging. To understand the underlying mechanism of calcium binding induced fluorescence increase in the designed calcium sensor CatchER, the fluorescence lifetime of CatchER was determined in calcium free and bound forms using time resolved florescence spectroscopy. The results suggest that calcium binding inhibits the geminate quenching, resulting in a longer lifetime when the anionic form is indirectly excited at 395 nm. It is believed that such unique calcium-induced lifetime change can be applied to monitor calcium signaling in cell imaging. NMR spectroscopy was used to investigate the protein-protein/ligand interaction in this dissertation. The residual dipolar coupling and T1, T2, NOE dynamic study were carried out to understand the binding mode of CaM and the N-terminal intracellular loop of connexin 43. The results show that both N and C terminal domains of Ca2+-CaM contact with the peptide, leading to a partially unwound and bending central helix of CaM. The ligand binding induced conformational change was demonstrated by selectively labeled proteins including extracellular domain of calcium sensing receptor and the bacterial membrane protein SecA fragments C34 and N68.

<p>Proteins are among the most complex entities known to science. Composed of just 20 fundamental building blocks arranged in simple linear strings, they nonetheless fold into a dizzying array of architectures that carry out the machinations of life at the molecular level.</p><p>Despite this central role in biology, we cannot reliably predict the structure of a protein from its sequence, and therefore rely on time-consuming and expensive experimental techniques to determine their structures. Although these methods can reveal equilibrium structures with great accuracy, they unfortunately mask much of the inherent molecular flexibility that enables proteins to dynamically perform biochemical tasks. As a result, much of the field of structural biology is mired in a static perspective; indeed, most attempts to naively model increased structural flexibility still end in failure.</p><p>This document details my work to validate alternative protein conformations beyond the primary or equilibrium conformation. The underlying hypothesis is that more realistic modeling of flexibility will enhance our understanding of how natural proteins function, and thereby improve our ability to design new proteins that perform desired novel functions.</p><p>During the course of my work, I used structure validation techniques to validate conformational alternatives in a variety of settings. First, I extended previous work introducing the backrub, a local, sidechain-coupled backbone motion, by demonstrating that backrubs also accompany sequence changes and therefore are useful for modeling conformational changes associated with mutations in protein design. Second, I extensively studied a new local backbone motion, helix shear, by documenting its occurrence in both crystal and NMR structures and showing its suitability for expanding conformational search space in protein design. Third, I integrated many types of local alternate conformations in an ultra-high-resolution crystal structure and discovered the combinatorial complexity that arises when adjacent flexible segments combine into networks. Fourth, I used structural bioinformatics techniques to construct smoothed, multi-dimensional torsional distributions that can be used to validate trial conformations or to propose new ones. Fifth, I participated in judging a structure prediction competition by using validation of geometrical and all-atom contact criteria to help define correctness across thousands of submitted conformations. Sixth, using similar tools plus collation of multiple comparable structures from the public database, I determined that low-energy states identified by the popular structure modeling suite Rosetta sometimes are valid conformations likely to be populated in the cell, but more often are invalid conformations attributable to artifacts in the physical/statistical hybrid energy function.</p><p>Unified by the theme of validating conformational alternatives by reference to high-quality experimental structures, my cumulative work advances our fundamental understanding of protein structural variability, and will benefit future endeavors to design useful proteins for biomedicine or industrial chemistry.</p> / Dissertation

The creation of enzymes displaying desired substrate specificity is an important objective of enzyme engineering. To help achieve this goal, computational protein design (CPD) can be used to identify sequences that can fulfill interactions required to productively bind a desired substrate. Standard CPD protocols find optimal sequences in the context of a single state, for example an enzyme structure with a single substrate bound at its active site. However, many enzymes catalyze reactions requiring them to bind multiple substrates during successive steps of the catalytic cycle. The design of multi-substrate enzyme specificity requires the ability to evaluate sequences in the context of multiple substrate-bound states because mutations designed to enhance activity for one substrate may be detrimental to the binding of a second substrate. Additionally, many enzymes undergo conformational changes throughout their catalytic cycle and the equilibrium between these conformations can have an impact on their substrate specificity. In this thesis, I present the development and implementation of two multistate computational protein design methodologies for the redesign of multi-substrate enzyme specificity and the modulation of enzyme conformational equilibrium. Overall, our approaches open the door to the design of multi-substrate enzymes displaying tailored specificity for any biocatalytic application.

abstract: This thesis explores a wide array of topics related to the protein folding problem, ranging from the folding mechanism, ab initio structure prediction and protein design, to the mechanism of protein functional evolution, using multi-scale approaches. To investigate the role of native topology on folding mechanism, the native topology is dissected into non-local and local contacts. The number of non-local contacts and non-local contact orders are both negatively correlated with folding rates, suggesting that the non-local contacts dominate the barrier-crossing process. However, local contact orders show positive correlation with folding rates, indicating the role of a diffusive search in the denatured basin. Additionally, the folding rate distribution of E. coli and Yeast proteomes are predicted from native topology. The distribution is fitted well by a diffusion-drift population model and also directly compared with experimentally measured half life. The results indicate that proteome folding kinetics is limited by protein half life. The crucial role of local contacts in protein folding is further explored by the simulations of WW domains using Zipping and Assembly Method. The correct formation of N-terminal &beta;-turn turns out important for the folding of WW domains. A classification model based on contact probabilities of five critical local contacts is constructed to predict the foldability of WW domains with 81% accuracy. By introducing mutations to stabilize those critical local contacts, a new protein design approach is developed to re-design the unfoldable WW domains and make them foldable. After folding, proteins exhibit inherent conformational dynamics to be functional. Using molecular dynamics simulations in conjunction with Perturbation Response Scanning, it is demonstrated that the divergence of functions can occur through the modification of conformational dynamics within existing fold for &beta;-lactmases and GFP-like proteins: i) the modern TEM-1 lactamase shows a comparatively rigid active-site region, likely reflecting adaptation for efficient degradation of a specific substrate, while the resurrected ancient lactamases indicate enhanced active-site flexibility, which likely allows for the binding and subsequent degradation of different antibiotic molecules; ii) the chromophore and attached peptides of photocoversion-competent GFP-like protein exhibits higher flexibility than the photocoversion-incompetent one, consistent with the evolution of photocoversion capacity. / Dissertation/Thesis / Ph.D. Physics 2014

abstract: Cyanovirin-N (CVN) is a cyanobacterial lectin with potent anti-HIV activity, mediated by binding to the N-linked oligosaccharide moiety of the envelope protein gp120. CVN offers a scaffold to develop multivalent carbohydrate-binding proteins with tunable specificities and affinities. I present here biophysical calculations completed on a monomeric-stabilized mutant of cyanovirin-N, P51G-m4-CVN, in which domain A binding activity is abolished by four mutations; with comparisons made to CVN<super>mutDB</super>, in which domain B binding activity is abolished. Using Monte Carlo calculations and docking simulations, mutations in CVN<super>mutDB</super> were considered singularly, and the mutations E41A/G and T57A were found to impact the affinity towards dimannose the greatest. <super>15</super>N-labeled proteins were titrated with Man&#945;(1-2)Man&#945;, while following chemical shift perturbations in NMR spectra. The mutants, E41A/G and T57A, had a larger Kd than P51G-m4-CVN, matching the trends predicted by the calculations. We also observed that the N42A mutation affects the local fold of the binding pocket, thus removing all binding to dimannose. Characterization of the mutant N53S showed similar binding affinity to P51G-m4-CVN. Using biophysical calculations allows us to study future iterations of models to explore affinities and specificities. In order to further elucidate the role of multivalency, I report here a designed covalent dimer of CVN, Nested cyanovirin-N (Nested CVN), which has four binding sites. Nested CVN was found to have comparable binding affinity to gp120 and antiviral activity to wt CVN. These results demonstrate the ability to create a multivalent, covalent dimer that has comparable results to that of wt CVN. WW domains are small modules consisting of 32-40 amino acids that recognize proline-rich peptides and are found in many signaling pathways. We use WW domain sequences to explore protein folding by simulations using Zipping and Assembly Method. We identified five crucial contacts that enabled us to predict the folding of WW domain sequences based on those contacts. We then designed a folded WW domain peptide from an unfolded WW domain sequence by introducing native contacts at those critical positions. / Dissertation/Thesis / Doctoral Dissertation Biochemistry 2014

Traditional computational protein design (CPD) calculations model sequence perturbations and evaluate their stabilities using a single fixed protein backbone template in an approach referred to as single‐state design (SSD). However, certain design objectives require the explicit consideration of multiple conformational states. Cases where a multistate framework may be advantageous over the single‐state approach include the computer aided discovery of new enzyme substrates, the prediction of protein stabilities, and the design of protein dynamics. These design objectives can be tackled using multistate design (MSD). However, it is often the case that a design objective requires the consideration of a protein state having no available structure information. For such circumstances the multistate framework cannot be applied. In this thesis I present the development of two template and ensemble preparation methodologies and their application to three projects. The purpose of which is to demonstrate the necessary ensemble modeling strategies to overcome limitations in available structure information. Particular emphasis is placed on the ability to recapitulate experimental data to guide modelling of the design space. Specifically, the use of MSD allowed for the accurate prediction of a methyltransferase recognition motif and new substrates, the prediction of mutant sequence stabilities with quantitative accuracy, and the design of dynamics into the rigid Gβ1 scaffold producing a set of dynamic variants whose tryptophan residue exchanges between two conformations on the millisecond timescale. Implementation of both the ensemble, coordinate perturbation followed by energy minimization (PertMin), and template, rotamer optimization followed by energy minimization (ROM), generation protocols developed here allow for exploration and manipulation of the structure space enabling the success of these applications.