Transaminase-Catalyzed Synthesis of ß-Branched Noncanonical Amino Acids Driven by a Lysine Amine Donor.

Transaminases are choice biocatalysts for the synthesis of chiral primary amines, including amino acids bearing contiguous stereocenters. In this study, we employ lysine as a "smart" amine donor in transaminase-catalyzed dynamic kinetic resolution reactions to access ß-branched noncanonical arylalanines. Our mechanistic investigation demonstrates that, upon transamination, the lysine-derived ketone byproduct readily cyclizes to a six-membered imine, driving the equilibrium in the desired direction and thus alleviating the need to load superstoichiometric quantities of the amine donor or deploy a multienzyme cascade. Lysine also shows good overall compatibility with a panel of wild-type transaminases, a promising hint of its application as a smart donor more broadly. Indeed, by this approach, we furnished a broad scope of ß-branched arylalanines, including some bearing hitherto intractable cyclopropyl and isopropyl substituents, with high yields and excellent selectivities.

Complex macrocycle exploration: parallel, heuristic, and constraint-based conformer generation using ForceGen.

ForceGen is a template-free, non-stochastic approach for 2D to 3D structure generation and conformational elaboration for small molecules, including both non-macrocycles and macrocycles. For conformational search of non-macrocycles, ForceGen is both faster and more accurate than the best of all tested methods on a very large, independently curated benchmark of 2859 PDB ligands. In this study, the primary results are on macrocycles, including results for 431 unique examples from four separate benchmarks. These include complex peptide and peptide-like cases that can form networks of internal hydrogen bonds. By making use of new physical movements ("flips" of near-linear sub-cycles and explicit formation of hydrogen bonds), ForceGen exhibited statistically significantly better performance for overall RMS deviation from experimental coordinates than all other approaches. The algorithmic approach offers natural parallelization across multiple computing-cores. On a modest multi-core workstation, for all but the most complex macrocycles, median wall-clock times were generally under a minute in fast search mode and under 2 min using thorough search. On the most complex cases (roughly cyclic decapeptides and larger) explicit exploration of likely hydrogen bonding networks yielded marked improvements, but with calculation times increasing to several minutes and in some cases to roughly an hour for fast search. In complex cases, utilization of NMR data to constrain conformational search produces accurate conformational ensembles representative of solution state macrocycle behavior. On macrocycles of typical complexity (up to 21 rotatable macrocyclic and exocyclic bonds), design-focused macrocycle optimization can be practically supported by computational chemistry at interactive time-scales, with conformational ensemble accuracy equaling what is seen with non-macrocyclic ligands. For more complex macrocycles, inclusion of sparse biophysical data is a helpful adjunct to computation.

Conformational Strain of Macrocyclic Peptides in Ligand-Receptor Complexes Based on Advanced Refinement of Bound-State Conformers.

Macrocyclic peptides are an important modality in drug discovery, but molecular design is limited due to the complexity of their conformational landscape. To better understand conformational propensities, global strain energies were estimated for 156 protein-macrocyclic peptide cocrystal structures. Unexpectedly large strain energies were observed when the bound-state conformations were modeled with positional restraints. Instead, low-energy conformer ensembles were generated using xGen that fit experimental X-ray electron density maps and gave reasonable strain energy estimates. The ensembles featured significant conformational adjustments while still fitting the electron density as well or better than the original coordinates. Strain estimates suggest the interaction energy in protein-ligand complexes can offset a greater amount of strain for macrocyclic peptides than for small molecules and non-peptidic macrocycles. Across all molecular classes, the approximate upper bound on global strain energies had the same relationship with molecular size, and bound-state ensembles from xGen yielded favorable binding energy estimates.

XGen: Real-Space Fitting of Complex Ligand Conformational Ensembles to X-ray Electron Density Maps.

We report a new method for X-ray density ligand fitting and refinement that is suitable for a wide variety of small-molecule ligands, including macrocycles. The approach (called "xGen") augments a force field energy calculation with an electron density fitting restraint that yields an energy reward during the restrained conformational search. The resulting conformer pools balance goodness-of-fit with ligand strain. Real-space refinement from pre-existing ligand coordinates of 150 macrocycles resulted in occupancy-weighted conformational ensembles that exhibited low strain energy. The xGen ensembles improved upon electron density fit compared with the PDB reference coordinates without making use of atom-specific B-factors. Similarly, on nonmacrocycles, de novo fitting produced occupancy-weighted ensembles of many conformers that were generally better-quality density fits than the deposited primary/alternate conformational pairs. The results suggest ubiquitous low-energy ligand conformational ensembles in X-ray diffraction data and provide an alternative to using B-factors as model parameters.

Secondary amine selective Petasis (SASP) bioconjugation.

Selective modification of proteins enables synthesis of antibody-drug conjugates, cellular drug delivery and construction of new materials. Many groups have developed methods for selective N-terminal modification without affecting the side chain of lysine by judicious pH control. This is due to lower basicity of the N-terminus relative to lysine side chains. But none of the methods are capable of selective modification of secondary amines or N-terminal proline, which has similar basicity as lysine. Here, we report a secondary amine selective Petasis (SASP) reaction for selective bioconjugation at N-terminal proline. We exploited the ability of secondary amines to form highly electrophilic iminium ions with aldehydes, which rapidly reacted with nucleophilic organoboronates, resulting in robust labeling of N-terminal proline under biocompatible conditions. This is the first time the Petasis reaction has been utilized for selective modification of secondary amines on completely unprotected peptides and proteins under physiological conditions. Peptide screening results showed that the reaction is highly selective for N-terminal proline. There are no other chemical methods reported in literature that are selective for N-terminal proline in both peptides and proteins. This is a multicomponent reaction leading to the synthesis of doubly functionalized bioconjugates in one step that can be difficult to achieve using other methods. The key advantage of the SASP reaction includes its high chemoselective and stereoselective (>99% de) nature, and it affords dual labeled proteins in one pot. The broad utility of this bioconjugation is highlighted for a variety of peptides and proteins, including aldolase and creatine kinase.

Interatomic potentials and solvation parameters from protein engineering data for buried residues.

Van der Waals (vdW) interaction energies between different atom types, energies of hydrogen bonds (H-bonds), and atomic solvation parameters (ASPs) have been derived from the published thermodynamic stabilities of 106 mutants with available crystal structures by use of an originally designed model for the calculation of free-energy differences. The set of mutants included substitutions of uncharged, inflexible, water-inaccessible residues in alpha-helices and beta-sheets of T4, human, and hen lysozymes and HI ribonuclease. The determined energies of vdW interactions and H-bonds were smaller than in molecular mechanics and followed the "like dissolves like" rule, as expected in condensed media but not in vacuum. The depths of modified Lennard-Jones potentials were -0.34, -0.12, and -0.06 kcal/mole for similar atom types (polar-polar, aromatic-aromatic, and aliphatic-aliphatic interactions, respectively) and -0.10, -0.08, -0.06, -0.02, and nearly 0 kcal/mole for different types (sulfur-polar, sulfur-aromatic, sulfur-aliphatic, aliphatic-aromatic, and carbon-polar, respectively), whereas the depths of H-bond potentials were -1.5 to -1.8 kcal/mole. The obtained solvation parameters, that is, transfer energies from water to the protein interior, were 19, 7, -1, -21, and -66 cal/moleA(2) for aliphatic carbon, aromatic carbon, sulfur, nitrogen, and oxygen, respectively, which is close to the cyclohexane scale for aliphatic and aromatic groups but intermediate between octanol and cyclohexane for others. An analysis of additional replacements at the water-protein interface indicates that vdW interactions between protein atoms are reduced when they occur across water.

Small-molecule inhibition of the interaction between the translation initiation factors eIF4E and eIF4G.

Assembly of the eIF4E/eIF4G complex has a central role in the regulation of gene expression at the level of translation initiation. This complex is regulated by the 4E-BPs, which compete with eIF4G for binding to eIF4E and which have tumor-suppressor activity. To pharmacologically mimic 4E-BP function we developed a high-throughput screening assay for identifying small-molecule inhibitors of the eIF4E/eIF4G interaction. The most potent compound identified, 4EGI-1, binds eIF4E, disrupts eIF4E/eIF4G association, and inhibits cap-dependent translation but not initiation factor-independent translation. While 4EGI-1 displaces eIF4G from eIF4E, it effectively enhances 4E-BP1 association both in vitro and in cells. 4EGI-1 inhibits cellular expression of oncogenic proteins encoded by weak mRNAs, exhibits activity against multiple cancer cell lines, and appears to have a preferential effect on transformed versus nontransformed cells. The identification of this compound provides a new tool for studying translational control and establishes a possible new strategy for cancer therapy.