Accurate positioning of functional residues with robotics-inspired computational protein design.

SignificanceComputational protein design promises to advance applications in medicine and biotechnology by creating proteins with many new and useful functions. However, new functions require the design of specific and often irregular atom-level geometries, which remains a major challenge. Here, we develop computational methods that design and predict local protein geometries with greater accuracy than existing methods. Then, as a proof of concept, we leverage these methods to design new protein conformations in the enzyme ketosteroid isomerase that change the protein's preference for a key functional residue. Our computational methods are openly accessible and can be applied to the design of other intricate geometries customized for new user-defined protein functions.

Computational design of a modular protein sense-response system.

Sensing and responding to signals is a fundamental ability of living systems, but despite substantial progress in the computational design of new protein structures, there is no general approach for engineering arbitrary new protein sensors. Here, we describe a generalizable computational strategy for designing sensor-actuator proteins by building binding sites de novo into heterodimeric protein-protein interfaces and coupling ligand sensing to modular actuation through split reporters. Using this approach, we designed protein sensors that respond to farnesyl pyrophosphate, a metabolic intermediate in the production of valuable compounds. The sensors are functional in vitro and in cells, and the crystal structure of the engineered binding site closely matches the design model. Our computational design strategy opens broad avenues to link biological outputs to new signals.

Optimization of Orally Bioavailable Enhancer of Zeste Homolog 2 (EZH2) Inhibitors Using Ligand and Property-Based Design Strategies: Identification of Development Candidate (R)-5,8-Dichloro-7-(methoxy(oxetan-3-yl)methyl)-2-((4-methoxy-6-methyl-2-oxo-1,2-dihydropyridin-3-yl)methyl)-3,4-dihydroisoquinolin-1(2H)-one (PF-06821497).

A new series of lactam-derived EZH2 inhibitors was designed via ligand-based and physicochemical-property-based strategies to address metabolic stability and thermodynamic solubility issues associated with previous lead compound 1. The new inhibitors incorporated an sp3 hybridized carbon atom at the 7-position of the lactam moiety present in lead compound 1 as a replacement for a dimethylisoxazole group. This transformation enabled optimization of the physicochemical properties and potency compared to compound 1. Analysis of relationships between calculated log D (clogD) values and in vitro metabolic stability and permeability parameters identified a clogD range that afforded an increased probability of achieving favorable ADME data in a single molecule. Compound 23a exhibited the best overlap of potency and pharmaceutical properties as well as robust tumor growth inhibition in vivo and was therefore advanced as a development candidate (PF-06821497). A crystal structure of 23a in complex with the three-protein PRC2 complex enabled understanding of the key structural features required for optimal binding.

Design and Synthesis of Pyridone-Containing 3,4-Dihydroisoquinoline-1(2H)-ones as a Novel Class of Enhancer of Zeste Homolog 2 (EZH2) Inhibitors.

A new enhancer of zeste homolog 2 (EZH2) inhibitor series comprising a substituted phenyl ring joined to a dimethylpyridone moiety via an amide linkage has been designed. A preferential amide torsion that improved the binding properties of the compounds was identified for this series via computational analysis. Cyclization of the amide linker resulted in a six-membered lactam analogue, compound 18. This transformation significantly improved the ligand efficiency/potency of the cyclized compound relative to its acyclic analogue. Additional optimization of the lactam-containing EZH2 inhibitors focused on lipophilic efficiency (LipE) improvement, which provided compound 31. Compound 31 displayed improved LipE and on-target potency in both biochemical and cellular readouts relative to compound 18. Inhibitor 31 also displayed robust in vivo antitumor growth activity and dose-dependent de-repression of EZH2 target genes.

Comparison of a High-Throughput Mass Spectrometry Method and Radioactive Filter Binding to Assay the Protein Methyltransferase PRMT5.

Conformational remodeling of chromatin in cells is known to alter gene expression. The histone code hypothesis postulates that multiple modifications present on histone tails can regulate gene expression both through direct effects on chromatin compaction as well as through recruitment of unique complexes that signal specific downstream functions. Histone methylation is an important component of the histone code, and the dysregulation of histone methylation in disease makes methyltransferases and demethylases viable targets for drug discovery. We developed a biochemical assay platform, which takes advantage of the fact that protein methyltransferases (PMTs) all utilize the cofactor S-Adenosyl-L-methionine (SAM) as the methyl donor. The platform utilizes the High-throughput Mass Spectrometry (MS) technology to measure SAM and the S-Adenosyl-L-homocysteine product in a label-free manner. The platform has all the advantages of a label-free system coupled with the benefit of substrate agnostic measurements making it an ideal setup for PMT biochemical studies and drug discovery. In addition, MS is ideally suited for detecting multiple modification events within the same substrate. The ability to adjust the detection to monitor the methyl acceptor product allows for real-time measurements of multiple product species simultaneously, a distinct advantage over other commonly used assay formats.

SAH derived potent and selective EZH2 inhibitors.

A series of novel enhancer of zeste homolog 2 (EZH2) inhibitors was designed based on the chemical structure of the histone methyltransferase (HMT) inhibitor SAH (S-adenosyl-l-homocysteine). These nucleoside-based EZH2 inhibitors blocked the methylation of nucleosomes at H3K27 in biochemical assays employing both WT PRC2 complex as well as a Y641N mutant PRC2 complex. The most potent compound, 27, displayed IC50's against both complexes of 270 nM and 70 nM, respectively. To our knowledge, compound 27 is the most potent SAH-derived inhibitor of the EZH2 PRC2 complex yet identified. This compound also displayed improved potency, lipophilic efficiency (LipE), and selectivity profile against other lysine methyltransferases compared with SAH.

Distal structural elements coordinate a conserved base flipping network.

One of the most dramatic illustrations of enzymatic promotion of a high-energy intermediate is observed in DNA modification and repair enzymes where an individual base is rotated (flipped) 180° around the deoxyribose-phosphate backbone and into the active site. While the end states have been extensively characterized, experimental techniques have yet to yield a full description of the base flipping process and the role played by the enzyme. The C5 cytosine methyltransferase M.HhaI coordinates an ensemble of reciprocal DNA and enzyme rearrangements to efficiently flip the target cytosine from the DNA helix. We sought to understand the role of individual amino acids during base flipping. Our results demonstrate that M.HhaI initiates base flipping before closure of the catalytic loop and utilizes the conserved serine 85 in the catalytic loop to accelerate flipping and maintain distortion of the DNA backbone. Serine 87, which forms specific contacts within the DNA helix after base flipping, is not involved in the flipping process or in maintaining the catalytically competent complex. At the base of the catalytic loop, glycine 98 acts as a hinge to allow conformational dynamism of the loop and mutation to alanine inhibits stabilization of the closed loop. Our results illustrate how an enzyme utilizes numerous, distal residues in concert to transform substrate recognition into catalysis.