Per Greisen

Accurate de novo design of hyperstable constrained peptides.

Naturally occurring, pharmacologically active peptides constrained with covalent crosslinks generally have shapes that have evolved to fit precisely into binding pockets on their targets. Such peptides can have excellent pharmaceutical properties, combining the stability and tissue penetration of small-molecule drugs with the specificity of much larger protein therapeutics. The ability to design constrained peptides with precisely specified tertiary structures would enable the design of shape-complementary inhibitors of arbitrary targets. Here we describe the development of computational methods for accurate de novo design of conformationally restricted peptides, and the use of these methods to design 18–47 residue, disulfide-crosslinked peptides, a subset of which are heterochiral and/or N–C backbone-cyclized. Both genetically encodable and non-canonical peptides are exceptionally stable to thermal and chemical denaturation, and 12 experimentally determined X-ray and NMR structures are nearly identical to the computational design models. The computational design methods and stable scaffolds presented here provide the basis for development of a new generation of peptide-based drugs.

Engineering V-type nerve agents detoxifying enzymes using computationally focused libraries.

VX and its Russian (RVX) and Chinese (CVX) analogues rapidly inactivate acetylcholinesterase and are the most toxic stockpile nerve agents. These organophosphates have a thiol leaving group with a choline-like moiety and are hydrolyzed very slowly by natural enzymes. We used an integrated computational and experimental approach to increase Brevundimonas diminuta phosphotriesterases (PTE) detoxification rate of V-agents by 5000-fold. Computational models were built of the complex between PTE and V-agents. On the basis of these models, the active site was redesigned to be complementary in shape to VX and RVX and to include favorable electrostatic interactions with their choline-like leaving group. Small libraries based on designed sequences were constructed. The libraries were screened by a direct assay for V-agent detoxification, as our initial studies showed that colorimetric surrogates fail to report the detoxification rates of the actual agents. The experimental results were fed back to improve the computational models. Overall, five rounds of iterating between experiment and model refinement led to variants that hydrolyze the toxic SP isomers of all three V-agents with kcat/KM values of up to 5 × 10(6) M(-1) min(-1) and also efficiently detoxify G-agents. These new catalysts provide the basis for broad spectrum nerve agent detoxification.

Interaction between sodium dodecyl sulfate and membrane reconstituted aquaporins: A comparative study of spinach SoPIP2;1 and E. coli AqpZ

This study describes the interaction between sodium dodecyl sulfate (SDS) and membrane proteins reconstituted into large unilamellar lipid vesicles and detergent micelles studied by circular dichroism (CD) and polarity sensitive probe labeling. Specifically, we carried out a comparative study of two aquaporins with high structural homology SoPIP2;1 and AqpZ using identical reconstitution conditions. Our CD results indicate that SDS, when added to membrane-reconstituted aquaporins in concentrations below the SDS critical micelle concentration (CMC, ~8mM), causes helical rearrangements of both aquaporins. However, we do not find compelling evidence for unfolding. In contrast when SDS is added to detergent stabilized aquaporins, SoPIP2;1 partly unfolds, while AqpZ secondary structure is unaffected. Using a fluorescent polarity sensitive probe (Badan) we show that SDS action on membrane reconstituted SoPIP2;1 as well as AqpZ is associated with initial increased hydrophobic interactions in protein transmembrane (TM) spanning regions up to a concentration of 0.1× CMC. At higher SDS concentrations TM hydrophobic interactions, as reported by Badan, decrease and reach a plateau from SDS CMC up to 12.5× CMC. Combined, our results show that SDS does not unfold neither SoPIP2;1 nor AqpZ during transition from a membrane reconstituted form to a detergent stabilized state albeit the native folds are changed.

Catalytic efficiencies of directly evolved phosphotriesterase variants with structurally different organophosphorus compounds in vitro.

The nearly 200,000 fatalities following exposure to organophosphorus (OP) pesticides each year and the omnipresent danger of a terroristic attack with OP nerve agents emphasize the demand for the development of effective OP antidotes. Standard treatments for intoxicated patients with a combination of atropine and an oxime are limited in their efficacy. Thus, research focuses on developing catalytic bioscavengers as an alternative approach using OP-hydrolyzing enzymes such as Brevundimonas diminuta phosphotriesterase (PTE). Recently, a PTE mutant dubbed C23 was engineered, exhibiting reversed stereoselectivity and high catalytic efficiency (kcat/KM) for the hydrolysis of the toxic enantiomers of VX, CVX, and VR. Additionally, C23’s ability to prevent systemic toxicity of VX using a low protein dose has been shown in vivo. In this study, the catalytic efficiencies of V-agent hydrolysis by two newly selected PTE variants were determined. Moreover, in order to establish trends in sequence–activity relationships along the pathway of PTE’s laboratory evolution, we examined kcat/KM values of several variants with a number of V-type and G-type nerve agents as well as with different OP pesticides. Although none of the new PTE variants exhibited kcat/KM values >107xa0M−1xa0min−1 with V-type nerve agents, which is required for effective prophylaxis, they were improved with VR relative to previously evolved variants. The new variants detoxify a broad spectrum of OPs and provide insight into OP hydrolysis and sequence–activity relationships.

Computational design of environmental sensors for the potent opioid fentanyl

We describe the computational design of proteins that bind the potent analgesic fentanyl. Our approach employs a fast docking algorithm to find shape complementary ligand placement in protein scaffolds, followed by design of the surrounding residues to optimize binding affinity. Co-crystal structures of the highest affinity binder reveal a highly preorganized binding site, and an overall architecture and ligand placement in close agreement with the design model. We use the designs to generate plant sensors for fentanyl by coupling ligand binding to design stability. The method should be generally useful for detecting toxic hydrophobic compounds in the environment.

Overcoming an optimization plateau in the directed evolution of highly efficient nerve agent bioscavengers

Improving an enzymes initially low catalytic efficiency with a new target substrate by an order of magnitude or two may require only a few rounds of mutagenesis and screening or selection. However, subsequent rounds of optimization tend to yield decreasing degrees of improvement (diminishing returns) eventually leading to an optimization plateau. We aimed to optimize the catalytic efficiency of bacterial phosphotriesterase (PTE) toward V-type nerve agents. Previously, we improved the catalytic efficiency of wild-type PTE toward the nerve agent VX by 500-fold, to a catalytic efficiency (kcat/KM) of 5 × 106 M-1 min-1. However, effective in vivo detoxification demands an enzyme with a catalytic efficiency of >107 M-1 min-1. Here, following eight additional rounds of directed evolution and the computational design of a stabilized variant, we evolved PTE variants that detoxify VX with a kcat/KM ≥ 5 × 107 M-1 min-1 and Russian VX (RVX) with a kcat/KM ≥ 107 M-1 min-1. These final 10-fold improvements were the most time consuming and laborious, as most libraries yielded either minor or no improvements. Stabilizing the evolving enzyme, and avoiding tradeoffs in activity with different substrates, enabled us to obtain further improvements beyond the optimization plateau and evolve PTE variants that were overall improved by >5000-fold with VX and by >17 000-fold with RVX. The resulting variants also hydrolyze G-type nerve agents with high efficiency (GA, GB at kcat/KM > 5 × 107 M-1 min-1) and can thus serve as candidates for broad-spectrum nerve-agent prophylaxis and post-exposure therapy using low enzyme doses.

Sampling and Energy Evaluation Challenges in Ligand Binding Protein Design

The steroid hormone 17α‐hydroxylprogesterone (17‐OHP) is a biomarker for congenital adrenal hyperplasia and hence there is considerable interest in development of sensors for this compound. We used computational protein design to generate protein models with binding sites for 17‐OHP containing an extended, nonpolar, shape‐complementary binding pocket for the four‐ring core of the compound, and hydrogen bonding residues at the base of the pocket to interact with carbonyl and hydroxyl groups at the more polar end of the ligand. Eight of 16 designed proteins experimentally tested bind 17‐OHP with micromolar affinity. A co‐crystal structure of one of the designs revealed that 17‐OHP is rotated 180° around a pseudo‐two‐fold axis in the compound and displays multiple binding modes within the pocket, while still interacting with all of the designed residues in the engineered site. Subsequent rounds of mutagenesis and binding selection improved the ligand affinity to nanomolar range, while appearing to constrain the ligand to a single bound conformation that maintains the same “flipped” orientation relative to the original design. We trace the discrepancy in the design calculations to two sources: first, a failure to model subtle backbone changes which alter the distribution of sidechain rotameric states and second, an underestimation of the energetic cost of desolvating the carbonyl and hydroxyl groups of the ligand. The difference between design model and crystal structure thus arises from both sampling limitations and energy function inaccuracies that are exacerbated by the near two‐fold symmetry of the molecule.

Computational redesign of metalloenzymes for catalyzing new reactions.

The ability to design novel activities in existing metalloenzyme active sites is a stringent test of our understanding of enzyme mechanisms, sheds light on enzyme evolution, and would have many practical applications. Here, we describe a computational method in the context of the macromolecular modeling suite Rosetta to repurpose active sites containing metal ions for reactions of choice. The required inputs for the method are a model of the transition state(s) for the reaction and a set of crystallographic structures of proteins containing metal ions. The coordination geometry associated with the metal ion (Zn(2+), for example) is automatically detected and the transition state model is aligned to the open metal coordination site(s) in the protein. Additional interactions to the transition state model are made using RosettaMatch and the surrounding amino acid side chain identities are optimized for transition state stabilization using RosettaDesign. Validation of the design is performed using docking and molecular dynamics simulations, and candidate designs are generated for experimental validation. Computational metalloenzyme repurposing is complementary to directed evolution approaches for enzyme engineering and allows large jumps in sequence space to make concerted sequence and structural changes for introducing novel enzymatic activities and specificities.

Phosphorus binding sites in proteins: structural preorganization and coordination.

Phosphorus is a ubiquitous element of the cell, which is found throughout numerous key molecules related to cell structure, energy and information storage and transfer, and a diverse array of other cellular functions. In this work, we adopt an approach often used for characterizing metal binding and selectivity of metalloproteins in terms of interactions in a first shell (direct residue interactions with the metal) and a second shell (residue interactions with first shell residues) and use it to characterize binding of phosphorus compounds. Similar analyses of binding have previously been limited to individual structures that bind to phosphate groups; here, we investigate a total of 8307 structures obtained from the RCSB Protein Data Bank (PDB). An analysis of the binding site amino acid propensities reveals very characteristic first shell residue distributions, which are found to be influenced by the characteristics of the phosphorus compound and by the presence of cobound cations. The second shell, which supports the coordinating residues in the first shell, is found to consist mainly of protein backbone groups. Our results show how the second shell residue distribution is dictated mainly by the first shell of the binding site, especially by cobound cations and that the main function of the second shell is to stabilize the first shell residues.

Qm Simulation Of Binding Site In P-type Atpases

P-type ATPases are enzymes that stablish cation gradients across biologcal membranes where ions are transported using the energy from hydrolysis of ATP. The Na,K-ATPase and SERCA1a transport 3 Na+/ 2 K+ and Ca2+ , respectively, against a concentration gradient for each ATP molecule hydrolyzed. The two enzymes have very high sequence homology and show structural similarities in the binding site of the ions. This raises questions concerning the selectivity of the different structural basis and how each protein select its specificion. The present study aims at characterising this selectivety.In order to understand the specificity of the enzymes, we investigate the coordination site in the two crystal structure of SERCA1a(PDB ID: 1SU4) and Na+/K+ ATPase(PDB ID: 3B8) by constructing models of the active site. The models are constructed using analogs of the amino acids from the first and second solvation shell. We optimise the geometry of the models by constraining the C-alpha atoms using semi-empirical methods such as PM6 and AM1. The energies and geometries are investigated further using ab initio, restricted Hartree Fock, and density functional theory with the B3LYP functional using the basis set 6-31G(d,p)∗. We substitute Ca2±ions into the binding pocket of the Na+/K+ ATPase model to investigate how the specificity is related to the coordination environment. Furthermore, the protonation state of side chains coordinating the ions are investigated in order to estimate the influence of pH changes on the binding of ions in the two pumps.