Transition states. Trapping a transition state in a computationally designed protein bottle

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to prevent unphysical rearrangements during geometry relaxation. A simple fixed-atom or "coordinate-lock" approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

Engineering of an in-cell protein crystal for fastening a metastable conformation of a target miniprotein

Protein crystals can be utilized as porous scaffolds to capture exogenous molecules. Immobilization of target proteins using protein crystals is expected to facilitate X-ray structure analysis of proteins that are difficult to be crystallized. One of the advantages of scaffold-assisted structure determination is the analysis of metastable structures that are not observed in solution. However, efforts to fix target proteins within the pores of scaffold protein crystals have been limited due to the lack of strategies to control protein-protein interactions formed in the crystals. In this study, we analyze the metastable structure of the miniprotein, CLN025, which forms a β-hairpin structure in solution, using a polyhedra crystal (PhC), an in-cell protein crystal. CLN025 is successfully fixed within the PhC scaffold by replacing the original loop region. X-ray crystal structure analysis and molecular dynamics (MD) simulation reveal that CLN025 is fixed as a helical structure in a metastable state by non-covalent interactions in the scaffold crystal. These results indicate that modulation of intermolecular interactions can trap various protein conformations in the engineered PhC and provides a new strategy for scaffold-assisted structure determination.

Enzyme redesign and genetic code expansion

Enzyme design is an important application of computational protein design (CPD). It can benefit enormously from the additional chemistries provided by noncanonical amino acids (ncAAs). These can be incorporated into an 'expanded' genetic code, and introduced in vivo into target proteins. The key step for genetic code expansion is to engineer an aminoacyl-transfer RNA (tRNA) synthetase (aaRS) and an associated tRNA that handles the ncAA. Experimental directed evolution has been successfully used to engineer aaRSs and incorporate over 200 ncAAs into expanded codes. But directed evolution has severe limits, and is not yet applicable to noncanonical AA backbones. CPD can help address several of its limitations, and has begun to be applied to this problem. We review efforts to redesign aaRSs, studies that designed new proteins and functionalities with the help of ncAAs, and some of the method developments that have been used, such as adaptive landscape flattening Monte Carlo, which allows an enzyme to be redesigned with substrate or transition state binding as the design target.

Coumarin luciferins and mutant luciferases for robust multi-component bioluminescence imaging

Multi-component bioluminescence imaging requires an expanded collection of luciferase-luciferin pairs that emit far-red or near-infrared light. Toward this end, we prepared a new class of luciferins based on a red-shifted coumarin scaffold. These probes (CouLuc-1s) were accessed in a two-step sequence via direct modification of commercial dyes. The bioluminescent properties of the CouLuc-1 analogs were also characterized, and complementary luciferase enzymes were identified using a two-pronged screening strategy. The optimized enzyme-substrate pairs displayed robust photon outputs and emitted a significant portion of near-infrared light. The CouLuc-1 scaffolds are also structurally distinct from existing probes, enabling rapid multi-component imaging. Collectively, this work provides novel bioluminescent tools along with a blueprint for crafting additional fluorophore-derived probes for multiplexed imaging.

Genetic Incorporation of Selenotyrosine Significantly Improves Enzymatic Activity of Agrobacterium radiobacter Phosphotriesterase

Tyrosine plays important roles in many enzymes. To facilitate enzyme design, mechanistic studies and minimize structural perturbation in the active site, here we report the genetic incorporation of a novel unnatural amino acid selenotyrosine (SeHF), which has single-atom replacement in comparison to tyrosine. The arPTE-(Agrobacterium radiobacter Phosphotriesterase) Tyr309SeHF mutant exhibits a significant 12-fold increase in k_cat and 3.2-fold enhancement in k_cat /K_M at pH 7.0. Molecular dynamics simulations show that the SeHF309 mutation results in a conformational switch which opens up the product release pocket and increases the product release rate, thereby elevating the overall enzyme activity. Significant improvement of the catalytic efficiency at neutral pH by single unnatural amino acid (UAA) mutation broadens the application of this enzyme, and provides valuable insights to the mechanism. Our method represents a new approach for designing enzymes with enhanced activity.

Development and Testing of Force Field Parameters for Phenylalanine and Tyrosine Derivatives

Theoretical analyses are valuable for the exploration of the effects of unnatural amino acids on enzyme functions; however, many necessary parameters for unnatural amino acids remain lacking. In this study, we developed and tested force field parameters compatible with Amber ff14SB for 18 phenylalanine and tyrosine derivatives. The charge parameters were derived from ab initio calculations using the RESP fitting approach and then adjusted to reproduce the benchmark relative energies (at the MP2/TZ level) of the α- and β-backbones for each unnatural amino acid dipeptide. The structures optimized under the proposed force field parameters for the 18 unnatural amino acid dipeptides in both the α- and β-backbone forms were in good agreement with their QM structures, as the average RMSD was as small as 0.1 Å. The force field parameters were then tested in their application to seven proteins containing unnatural amino acids. The RMSDs of the simulated configurations of these unnatural amino acids were approximately 1.0 Å compared with those of the crystal structures. The vital interactions between proteins and unnatural amino acids in five protein–ligand complexes were also predicted using MM/PBSA analysis, and they were largely consistent with experimental observations. This work will provide theoretical aid for drug design involving unnatural amino acids.

Titrating bacterial growth and chemical biosynthesis for efficient N-acetylglucosamine and N-acetylneuraminic acid bioproduction

Metabolic engineering facilitates chemical biosynthesis by rewiring cellular resources to produce target compounds. However, an imbalance between cell growth and bioproduction often reduces production efficiency. Genetic code expansion (GCE)-based orthogonal translation systems incorporating non-canonical amino acids (ncAAs) into proteins by reassigning non-canonical codons to ncAAs qualify for balancing cellular metabolism. Here, GCE-based cell growth and biosynthesis balance engineering (GCE-CGBBE) is developed, which is based on titrating expression of cell growth and metabolic flux determinant genes by constructing ncAA-dependent expression patterns. We demonstrate GCE-CGBBE in genome-recoded Escherichia coli Δ321AM by precisely balancing glycolysis and N-acetylglucosamine production, resulting in a 4.54-fold increase in titer. GCE-CGBBE is further expanded to non-genome-recoded Bacillus subtilis to balance growth and N-acetylneuraminic acid bioproduction by titrating essential gene expression, yielding a 2.34-fold increase in titer. Moreover, the development of ncAA-dependent essential gene expression regulation shows efficient biocontainment of engineered B. subtilis to avoid unintended proliferation in nature.

An imbalance between cell growth and bioproduction of engineered microbes often reduces production efficiency. Here, the authors report genetic code expansion-based cell growth and biosynthesis balance engineering to achieve high levels production of N-acetylglucosamine in E. coli and N-acetylneuraminic acid in B. subtilis.

Structural Insights into How Protein Environments Tune the Spectroscopic Properties of a Noncanonical Amino Acid Fluorophore

Genetically encoded fluorescent noncanonical amino acids (fNCAAs) could be used to develop novel fluorescent sensors of protein function. Previous efforts toward this goal have been limited by the lack of extensive physicochemical and structural characterizations of protein-based sensors containing fNCAAs. Here, we report the steady-state spectroscopic properties and first structural analyses of an fNCAA-containing Fab fragment of the 5c8 antibody, which binds human CD40L. A previously reported 5c8 variant in which the light chain residue Ile^L98 is replaced with the fNCAA l-(7-hydroxycoumarin-4-yl)ethylglycine (7-HCAA) exhibits a 1.7-fold increase in fluorescence upon antigen binding. Determination and comparison of the apparent pK_as of 7-HCAA in the unbound and bound forms indicate that the observed increase in fluorescence is not the result of perturbations in pK_a. Crystal structures of the fNCAA-containing Fab in the apo and bound forms reveal interactions between the 7-HCAA side chain and surrounding residues that are disrupted upon antigen binding. This structural characterization not only provides insight into the manner in which protein environments can modulate the fluorescence properties of 7-HCAA but also could serve as a starting point for the rational design of new fluorescent protein-based reporters of protein function.

Multicomponent Bioluminescence Imaging with a π-Extended Luciferin

Bioluminescence imaging with luciferase-luciferin pairs is commonly used for monitoring biological processes in cells and whole organisms. Traditional bioluminescent probes are limited in scope, though, as they cannot be easily distinguished in biological environments, precluding efforts to visualize multicellular processes. Additionally, many luciferase-luciferin pairs emit light that is poorly tissue penetrant, hindering efforts to visualize targets in deep tissues. To address these issues, we synthesized a set of π-extended luciferins that were predicted to be red-shifted luminophores. The scaffolds were designed to be rotationally labile such that they produced light only when paired with luciferases capable of enforcing planarity. A luciferin comprising an intramolecular "lock" was identified as a viable light-emitting probe. Native luciferases were unable to efficiently process the analog, but a complementary luciferase was identified via Rosetta-guided enzyme design. The unique enzyme-substrate pair is red-shifted compared to well-known bioluminescent tools. The probe set is also orthogonal to other luciferase-luciferin probes and can be used for multicomponent imaging. Four substrate-resolved luciferases were imaged in a single session. Collectively, this work provides the first example of Rosetta-guided design in engineering bioluminescent tools and expands the scope of orthogonal imaging probes.

Protein engineering: the potential of remote mutations

Engineered proteins, especially enzymes, are now commonly used in many industries owing to their catalytic power, specific binding of ligands, and properties as materials and food additives. As the number of potential uses for engineered proteins has increased, the interest in engineering or designing proteins to have greater stability, activity and specificity has increased in turn. With any rational engineering or design pursuit, the success of these endeavours relies on our fundamental understanding of the systems themselves; in the case of proteins, their structure-dynamics-function relationships. Proteins are most commonly rationally engineered by targeting the residues that we understand to be functionally important, such as enzyme active sites or ligand-binding sites. This means that the majority of the protein, i.e. regions remote from the active- or ligand-binding site, is often ignored. However, there is a growing body of literature that reports on, and rationalises, the successful engineering of proteins at remote sites. This minireview will discuss the current state of the art in protein engineering, with a particular focus on engineering regions that are remote from active- or ligand-binding sites. As the use of protein technologies expands, exploiting the potential improvements made possible through modifying remote regions will become vital if we are to realise the full potential of protein engineering and design.

A genetically encoded photosensitizer protein facilitates the rational design of a miniature photocatalytic CO2-reducing enzyme

Photosensitizers, which harness light energy to upgrade weak reductants to strong reductants, are pivotal components of the natural and artificial photosynthesis machineries. However, it has proved difficult to enhance and expand their functions through genetic engineering. Here we report a genetically encoded, 27 kDa photosensitizer protein (PSP), which facilitates the rational design of miniature photocatalytic CO₂-reducing enzymes. Visible light drives PSP efficiently into a long-lived triplet excited state (PSP*), which reacts rapidly with reduced nicotinamide adenine dinucleotide to generate a super-reducing radical (PSP^•), which is strong enough to reduce many CO₂-reducing catalysts. We determined the three-dimensional structure of PSP^• at 1.8 Å resolution by X-ray crystallography. Genetic engineering enabled the site-specific attachment of a nickel-terpyridine complex and the modular optimization of the photochemical properties of PSP, the chromophore/catalytic centre distance and the catalytic centre microenvironment, which culminated in a miniature photocatalytic CO₂-reducing enzyme that has a CO₂/CO conversion quantum efficiency of 2.6%.

An expanded allosteric network in PTP1B by multitemperature crystallography, fragment screening, and covalent tethering

Allostery is an inherent feature of proteins, but it remains challenging to reveal the mechanisms by which allosteric signals propagate. A clearer understanding of this intrinsic circuitry would afford new opportunities to modulate protein function. Here, we have identified allosteric sites in protein tyrosine phosphatase 1B (PTP1B) by combining multiple-temperature X-ray crystallography experiments and structure determination from hundreds of individual small-molecule fragment soaks. New modeling approaches reveal 'hidden' low-occupancy conformational states for protein and ligands. Our results converge on allosteric sites that are conformationally coupled to the active-site WPD loop and are hotspots for fragment binding. Targeting one of these sites with covalently tethered molecules or mutations allosterically inhibits enzyme activity. Overall, this work demonstrates how the ensemble nature of macromolecular structure, revealed here by multitemperature crystallography, can elucidate allosteric mechanisms and open new doors for long-range control of protein function.

Proteins perform many important jobs in each of the cells in our bodies, such as transporting other molecules and helping chemical reactions to occur. The part of the protein directly involved in these tasks is called the active site. Other areas of the protein can communicate with the active site to switch the protein on or off. This method of control is known as allostery.

Switching proteins on and off could help us to develop treatments for certain diseases. For example, a protein called PTP1B reduces how well cells can respond to insulin. Switching this protein off could therefore help to treat diabetes. However, much like it’s hard to guess how a light switch is wired to a light bulb without seeing behind the walls, it is hard to predict which remote areas of a protein are ‘wired’ to the active site.

Keedy, Hill et al. have now used two complementary methods to examine the structure of PTP1B and find new allosteric sites. The first method captured a series of X-ray images from crystallized molecules of the protein held at different temperatures. This revealed areas of PTP1B that can move like windshield wipers to communicate with each other. The second method soaked PTP1B crystals in trays with hundreds of drug-sized molecules and assessed which sites on the protein the molecules bound to. The molecules generally bound to just a few sites of the protein. Further tests on one of these sites showed that it can communicate with the active site to turn the protein on or off.

Further work will be needed to develop drugs that could treat diabetes by binding to the newly identified allosteric sites in PTP1B. More generally, the methods developed by Keedy, Hill et al. could be used to study allostery in other important proteins.

Playing with the Molecules of Life

Our understanding of the complex molecular processes of living organisms at the molecular level is growing exponentially. This knowledge, together with a powerful arsenal of tools for manipulating the structures of macromolecules, is allowing chemists to to harness and reprogram the cellular machinery in ways previously unimaged. Here we review one example in which the genetic code itself has been expanded with new building blocks that allow us to probe and manipulate the structures and functions of proteins with unprecedented precision.

A transition state “trapped”? QM-cluster models of engineered threonyl-tRNA synthetase

QM-cluster models demonstrate how protein bioengineering alters the local energy landscape of p-biphenylalanine to stabilize a transition state analogue.

Expanding and reprogramming the genetic code

Nature uses a limited, conservative set of amino acids to synthesize proteins. The ability to genetically encode an expanded set of building blocks with new chemical and physical properties is transforming the study, manipulation and evolution of proteins, and is enabling diverse applications, including approaches to probe, image and control protein function, and to precisely engineer therapeutics. Underpinning this transformation are strategies to engineer and rewire translation. Emerging strategies aim to reprogram the genetic code so that noncanonical biopolymers can be synthesized and evolved, and to test the limits of our ability to engineer the translational machinery and systematically recode genomes.

Rewriting the Genetic Code

The genetic code-the language used by cells to translate their genomes into proteins that perform many cellular functions-is highly conserved throughout natural life. Rewriting the genetic code could lead to new biological functions such as expanding protein chemistries with noncanonical amino acids (ncAAs) and genetically isolating synthetic organisms from natural organisms and viruses. It has long been possible to transiently produce proteins bearing ncAAs, but stabilizing an expanded genetic code for sustained function in vivo requires an integrated approach: creating recoded genomes and introducing new translation machinery that function together without compromising viability or clashing with endogenous pathways. In this review, we discuss design considerations and technologies for expanding the genetic code. The knowledge obtained by rewriting the genetic code will deepen our understanding of how genomes are designed and how the canonical genetic code evolved.

Mechanism of Formation of the Nonstandard Product in the Prenyltransferase Reaction of the G115T Mutant of FtmPT1: A Case of Reaction Dynamics Calling the Shots?

FtmPT1 is a fungal indole prenyltransferase that affords Tryprostatin B from Brevianamide F and dimethylallyl pyrophosphate; however, when a single residue in the active site is mutated (Gly115Thr), a novel five-membered ring compound is obtained as the major product with Tryprostatin B as the minor product. Herein, we describe detailed studies of the catalysis of the Gly115Thr mutant of FtmPT1 with a focus on the observed regioselectivity of the reaction. We employ one- and two-dimensional potential of mean force simulations to explore the catalytic mechanism, along with molecular dynamics simulations exploring the reaction dynamics of the prenyl transfer reaction. Single-point electronic structure calculations were also used to explore the performance of the self-consistent charge density functional tight-binding method to model specific reaction steps. Importantly, we observe that the two reaction pathways have comparable activation parameters and propose that the origin of the novel product is predicated, at least in part, on the topology of the potential energy surface as revealed by reaction dynamics studies.

Minimum action transition paths connecting minima on an energy surface

Dynamics is essential to the biological functions of many bio-molecules, yet our knowledge of dynamics remains fragmented. Experimental techniques for studying bio-molecules either provide high resolution information on static conformations of the molecule or provide low-resolution, ensemble information that does not shed light on single molecule dynamics. In parallel, bio-molecular dynamics occur at time scale that are not yet attainable through detailed simulation methods. These limitations are especially noticeable when studying transition paths. To address this issue, we report in this paper two methods that derive meaningful trajectories for proteins between two of their conformations. The first method, MinActionPath, uses approximations of the potential energy surface for the molecule to derive an analytical solution of the equations of motion related to the concept of minimum action path. The second method, RelaxPath, follows the same principle of minimum action path but implements a more sophisticated potential, including a mixed elastic potential and a collision term to alleviate steric clashes. Using this new potential, the equations of motion cannot be solved analytically. We have introduced a relaxation method for solving those equations. We describe both the theories behind the two methods and their implementations, focusing on the specific techniques we have used that make those implementations amenable to study large molecular systems. We have illustrated the performance of RelaxPath on simple 2D systems. We have also compared MinActionPath and RelaxPath to other methods for generating transition paths on a well suited test set of large proteins, for which the end points of the trajectories as well as an intermediate conformation between those end points are known. We have shown that RelaxPath outperforms those other methods, including MinActionPath, in its ability to generate trajectories that get close to the known intermediates. We have also shown that the structures along the RelaxPath trajectories remain protein-like. Open source versions of the two programs MinActionPath and RelaxPath are available by request.

Diversity of Secondary Structure in Catalytic Peptides with β-Turn-Biased Sequences

X-ray crystallography has been applied to the structural analysis of a series of tetrapeptides that were previously assessed for catalytic activity in an atroposelective bromination reaction. Common to the series is a central Pro-Xaa sequence, where Pro is either l- or d-proline, which was chosen to favor nucleation of canonical β-turn secondary structures. Crystallographic analysis of 35 different peptide sequences revealed a range of conformational states. The observed differences appear not only in cases where the Pro-Xaa loop-region is altered, but also when seemingly subtle alterations to the flanking residues are introduced. In many instances, distinct conformers of the same sequence were observed, either as symmetry-independent molecules within the same unit cell or as polymorphs. Computational studies using DFT provided additional insight into the analysis of solid-state structural features. Select X-ray crystal structures were compared to the corresponding solution structures derived from measured proton chemical shifts, 3J-values, and 1H-1H-NOESY contacts. These findings imply that the conformational space available to simple peptide-based catalysts is more diverse than precedent might suggest. The direct observation of multiple ground state conformations for peptides of this family, as well as the dynamic processes associated with conformational equilibria, underscore not only the challenge of designing peptide-based catalysts, but also the difficulty in predicting their accessible transition states. These findings implicate the advantages of low-barrier interconversions between conformations of peptide-based catalysts for multistep, enantioselective reactions.

At the Interface of Chemical and Biological Synthesis: An Expanded Genetic Code

The ability to site-specifically incorporate noncanonical amino acids (ncAAs) with novel structures into proteins in living cells affords a powerful tool to investigate and manipulate protein structure and function. More than 200 ncAAs with diverse biological, chemical, and physical properties have been genetically encoded in response to nonsense or frameshift codons in both prokaryotic and eukaryotic organisms with high fidelity and efficiency. In this review, recent advances in the technology and its application to problems in protein biochemistry, cellular biology, and medicine are highlighted.

QM/MM Analysis of Transition States and Transition State Analogues in Metalloenzymes

Enzymology is approaching an era where many problems can benefit from computational studies. While ample challenges remain in quantitatively predicting behavior for many enzyme systems, the insights that often come from computations are an important asset for the enzymology community. Here we provide a primer for enzymologists on the types of calculations that are most useful for mechanistic problems in enzymology. In particular, we emphasize the integration of models that range from small active-site motifs to fully solvated enzyme systems for cross-validation and dissection of specific contributions from the enzyme environment. We then use a case study of the enzyme alkaline phosphatase to illustrate specific application of the methods. The case study involves examination of the binding modes of putative transition state analogues (tungstate and vanadate) to the enzyme. The computations predict covalent binding of these ions to the enzymatic nucleophile and that they adopt the trigonal bipyramidal geometry of the expected transition state. By comparing these structures with transition states found through free energy simulations, we assess the degree to which the transition state analogues mimic the true transition states. Technical issues worth treating with care as well as several remaining challenges to quantitative analysis of metalloenzymes are also highlighted during the discussion.

Molecular Self-Assembly Strategy for Generating Catalytic Hybrid Polypeptides

Recently, catalytic peptides were introduced that mimicked protease activities and showed promising selectivity of products even in organic solvents where protease cannot perform well. However, their catalytic efficiency was extremely low compared to natural enzyme counterparts presumably due to the lack of stable tertiary fold. We hypothesized that assembling these peptides along with simple hydrophobic pockets, mimicking enzyme active sites, could enhance the catalytic activity. Here we fused the sequence of catalytic peptide CP4, capable of protease and esterase-like activities, into a short amyloidogenic peptide fragment of Aβ. When the fused CP4-Aβ construct assembled into antiparallel β-sheets and amyloid fibrils, a 4.0-fold increase in the hydrolysis rate of p-nitrophenyl acetate (p-NPA) compared to neat CP4 peptide was observed. The enhanced catalytic activity of CP4-Aβ assembly could be explained both by pre-organization of a catalytically competent Ser-His-acid triad and hydrophobic stabilization of a bound substrate between the triad and p-NPA, indicating that a design strategy for self-assembled peptides is important to accomplish the desired functionality.

Designed metalloprotein stabilizes a semiquinone radical

Enzymes use binding energy to stabilize their substrates in high-energy states that are otherwise inaccessible at ambient temperature. Here we show that a de novo designed Zn(II) metalloprotein stabilizes a chemically reactive organic radical that is otherwise unstable in aqueous media. The protein binds tightly to and stabilizes the radical semiquinone form of 3,5-di-tert-butylcatechol. Solution NMR spectroscopy in conjunction with molecular dynamics simulations show that the substrate binds in the active site pocket where it is stabilized by metal-ligand interactions as well as by burial of its hydrophobic groups. Spectrochemical redox titrations show that the protein stabilized the semiquinone by reducing the electrochemical midpoint potential for its formation via the one-electron oxidation of the catechol by approximately 400 mV (9 kcal mol(-1)). Therefore, the inherent chemical properties of the radical were changed drastically by harnessing its binding energy to the metalloprotein. This model sets the basis for designed enzymes with radical cofactors to tackle challenging chemistry.

Is it possible to reverse aged acetylcholinesterase inhibited by organophosphorus compounds? Insight from the theoretical study

The main treatment for organophosphorus (OP) compound poisoning in clinics is to restore the activity of acetylcholinesterase (AChE) through oxime-induced reactivation of the phosphorylated OP-AChE adduct. It suffers from a competitive and irreversible aging reaction of the phosphorylated OP-AChE adduct, resulting in permanent inactivity of AChE. However, it was recently reported that N-methyl-2-methoxypyridinium species can act as methylating agents to methylate the methyl methane-phosphonate monoanion, in which the reaction mimics the reverse of the aging reaction of the phosphorylated OP-AChE adduct. If the aging reaction could be really reversed, the efficiency for the OP detoxification should be significantly improved, bringing up the possibility to develop an agent to reverse the aging process of the phosphorylated OP-AChE adduct. However, such a reaction with the N-methyl-2-methoxypyridinium species in the enzyme is still not reported so far. It is of great interest to know whether or not this reaction is observable in the enzyme, and more importantly, if it turns out to be not observable in the enzyme, why such a reaction proceeds quickly in aqueous solution but not in the enzyme. In the present study, we performed DFT calculations and quantum mechanical/molecular mechanical (QM/MM) calculations to reveal the fundamental mechanism for the methylation of both the methyl methane-phosphonate monoanion and the aged sarin-AChE adduct by N-methyl-2-methoxypyridinium species, respectively. The obtained results support the SN2 reaction mechanism, not the stepwise mechanism, for the methylation of the methyl methane-phosphonate monoanion by 9 reported N-methyl-2-methoxypyridinium compounds. The calculated free energy barriers are in good agreement with the experimental data. The methylation of the aged sarin-AChE adduct by one N-methyl-2-methoxypyridinium compound (labeled as compound 2) also employs the SN2 reaction mechanism with an extremely high free energy barrier of 30.4 ± 3.5 (or 26.6) kcal mol(-1), implying that this reaction in the enzyme hardly occurs. Our results clearly show that compound 2 forms a strong π-π stacking interaction with the aromatic ring of the W86 residue of AChE, making itself unable to approach sarin for the reverse of the aging process. On the basis of the structure and mechanism, several possible strategies have been suggested for designing methylating agents with higher activity against the aged sarin-AChE adduct.

Human DNMT1 transition state structure

Human DNA methyltransferase 1 (DNMT1) maintains the epigenetic state of DNA by replicating CpG methylation signatures from parent to daughter strands, producing heritable methylation patterns through cell divisions. The proposed catalytic mechanism of DNMT1 involves nucleophilic attack of Cys(1226) to cytosine (Cyt) C6, methyl transfer from S-adenosyl-l-methionine (SAM) to Cyt C5, and proton abstraction from C5 to form methylated CpG in DNA. Here, we report the subangstrom geometric and electrostatic structure of the major transition state (TS) of the reaction catalyzed by human DNMT1. Experimental kinetic isotope effects were used to guide quantum mechanical calculations to solve the TS structure. Methyl transfer occurs after Cys(1226) attack to Cyt C6, and the methyl transfer step is chemically rate-limiting for DNMT1. Electrostatic potential maps were compared for the TS and ground states, providing the electronic basis for interactions between the protein and reactants at the TS. Understanding the TS of DNMT1 demonstrates the possibility of using similar analysis to gain subangstrom geometric insight into the complex reactions of epigenetic modifications.

Protein sequence design and its applications

Design of proteins has far-reaching potentials in diverse areas that span repurposing of the protein scaffold for reactions and substrates that they were not naturally meant for, to catching a glimpse of the ephemeral proteins that nature might have sampled during evolution. These non-natural proteins, either in synthesized or virtual form have opened the scope for the design of entities that not only rival their natural counterparts but also offer a chance to visualize the protein space continuum that might help to relate proteins and understand their associations. Here, we review the recent advances in protein engineering and design, in multiple areas, with a view to drawing attention to their future potential.

Single-Point Mutation with a Rotamer Library Toolkit: Toward Protein Engineering

Protein engineers have long been hard at work to harness biocatalysts as a natural source of regio-, stereo-, and chemoselectivity in order to carry out chemistry (reactions and/or substrates) not previously achieved with these enzymes. The extreme labor demands and exponential number of mutation combinations have induced computational advances in this domain. The first step in our virtual approach is to predict the correct conformations upon mutation of residues (i.e., rebuilding side chains). For this purpose, we opted for a combination of molecular mechanics and statistical data. In this work, we have developed automated computational tools to extract protein structural information and created conformational libraries for each amino acid dependent on a variable number of parameters (e.g., resolution, flexibility, secondary structure). We have also developed the necessary tool to apply the mutation and optimize the conformation accordingly. For side-chain conformation prediction, we obtained overall average root-mean-square deviations (RMSDs) of 0.91 and 1.01 Å for the 18 flexible natural amino acids within two distinct sets of over 3000 and 1500 side-chain residues, respectively. The commonly used dihedral angle differences were also evaluated and performed worse than the state of the art. These two metrics are also compared. Furthermore, we generated a family-specific library for kinases that produced an average 2% lower RMSD upon side-chain reconstruction and a residue-specific library that yielded a 17% improvement. Ultimately, since our protein engineering outlook involves using our docking software, Fitted/Impacts, we applied our mutation protocol to a benchmarked data set for self- and cross-docking. Our side-chain reconstruction does not hinder our docking software, demonstrating differences in pose prediction accuracy of approximately 2% (RMSD cutoff metric) for a set of over 200 protein/ligand structures. Similarly, when docking to a set of over 100 kinases, side-chain reconstruction (using both general and biased conformation libraries) had minimal detriment to the docking accuracy.

Biomolecular Crowding Arising from Small Molecules, Molecular Constraints, Surface Packing, and Nano-Confinement

The effect of macromolecular crowding on the structure, dynamics, and reactivity of biomolecules is well established and the relevant research has been extensively reviewed. Herein, we focus our discussion on crowding effects arising from small cosolvent molecules and densely packed surface conditions. In addition, we highlight recent efforts that capitalize on the excluded volume effect for various tailored biochemical and biophysical applications. Specifically, we discuss how a targeted increase in local mass density can be exploited to gain insight into the folding dynamics of the protein of interest and how confinement via reverse micelles can be used to study a range of biophysical questions, from protein hydration dynamics to amyloid formation.

Smooth heuristic optimization on a complex chemical subspace

In the present work, several heuristic reordering algorithms for deterministic optimization on a combinatorial chemical compound space are evaluated for performance and efficiency.