AlphaSpace 2.0: Representing Concave Biomolecular Surfaces Using ß-Clusters.

Modern rational modulator design and structure-function characterization often concentrate on concave regions of biomolecular surfaces, ranging from well-defined small-molecule binding sites to large protein-protein interaction interfaces. Here, we introduce a ß-cluster as a pseudomolecular representation of fragment-centric pockets detected by AlphaSpace [J. Chem. Inf. Model. 2015, 55, 1585], a recently developed computational analysis tool for topographical mapping of biomolecular concavities. By mimicking the shape as well as atomic details of potential molecular binders, this new ß-cluster representation allows direct pocket-to-ligand shape comparison and can be used to guide ligand optimization. Furthermore, we defined the ß-score, the optimal Vina score of the ß-cluster, as an indicator of pocket ligandability and developed an ensemble ß-cluster approach, which allows one-to-one pocket mapping and comparison among aligned protein structures. We demonstrated the utility of ß-cluster representation by applying the approach to a wide variety of problems including binding site detection and comparison, characterization of protein-protein interactions, and fragment-based ligand optimization. These new ß-cluster functionalities have been implemented in AlphaSpace 2.0, which is freely available on the web at http://www.nyu.edu/projects/yzhang/AlphaSpace2.

Computational Strategy for Bound State Structure Prediction in Structure-Based Virtual Screening: A Case Study of Protein Tyrosine Phosphatase Receptor Type O Inhibitors.

Accurate protein structure in the ligand-bound state is a prerequisite for successful structure-based virtual screening (SBVS). Therefore, applications of SBVS against targets for which only an apo structure is available may be severely limited. To address this constraint, we developed a computational strategy to explore the ligand-bound state of a target protein, by combined use of molecular dynamics simulation, MM/GBSA binding energy calculation, and fragment-centric topographical mapping. Our computational strategy is validated against low-molecular weight protein tyrosine phosphatase (LMW-PTP) and then successfully employed in the SBVS against protein tyrosine phosphatase receptor type O (PTPRO), a potential therapeutic target for various diseases. The most potent hit compound GP03 showed an IC50 value of 2.89 µM for PTPRO and possessed a certain degree of selectivity toward other protein phosphatases. Importantly, we also found that neglecting the ligand energy penalty upon binding partially accounts for the false positive SBVS hits. The preliminary structure-activity relationships of GP03 analogs are also reported.

Targeting Unoccupied Surfaces on Protein-Protein Interfaces.

The use of peptidomimetic scaffolds to target protein-protein interfaces is a promising strategy for inhibitor design. The strategy relies on mimicry of protein motifs that exhibit a concentration of native hot spot residues. To address this constraint, we present a pocket-centric computational design strategy guided by AlphaSpace to identify high-quality pockets near the peptidomimetic motif that are both targetable and unoccupied. Alpha-clusters serve as a spatial representation of pocket space and are used to guide the selection of natural and non-natural amino acid mutations to design inhibitors that optimize pocket occupation across the interface. We tested the strategy against a challenging protein-protein interaction target, KIX/MLL, by optimizing a single helical motif within MLL to compete against the full-length wild-type MLL sequence. Molecular dynamics simulation and experimental fluorescence polarization assays are used to verify the efficacy of the optimized peptide sequence.

AlphaSpace: Fragment-Centric Topographical Mapping To Target Protein-Protein Interaction Interfaces.

Inhibition of protein-protein interactions (PPIs) is emerging as a promising therapeutic strategy despite the difficulty in targeting such interfaces with drug-like small molecules. PPIs generally feature large and flat binding surfaces as compared to typical drug targets. These features pose a challenge for structural characterization of the surface using geometry-based pocket-detection methods. An attractive mapping strategy--that builds on the principles of fragment-based drug discovery (FBDD)--is to detect the fragment-centric modularity at the protein surface and then characterize the large PPI interface as a set of localized, fragment-targetable interaction regions. Here, we introduce AlphaSpace, a computational analysis tool designed for fragment-centric topographical mapping (FCTM) of PPI interfaces. Our approach uses the alpha sphere construct, a geometric feature of a protein's Voronoi diagram, to map out concave interaction space at the protein surface. We introduce two new features--alpha-atom and alpha-space--and the concept of the alpha-atom/alpha-space pair to rank pockets for fragment-targetability and to facilitate the evaluation of pocket/fragment complementarity. The resulting high-resolution interfacial map of targetable pocket space can be used to guide the rational design and optimization of small molecule or biomimetic PPI inhibitors.

Revelation of a catalytic calcium-binding site elucidates unusual metal dependence of a human apyrase.

Human soluble calcium-activated nucleotidase 1 (hSCAN-1) represents a new family of apyrase enzymes that catalyze the hydrolysis of nucleotide di- and triphosphates, thereby modulating extracellular purinergic and pyrimidinergic signaling. Among well-characterized phosphoryl transfer enzymes, hSCAN-1 is unique not only in its unusual calcium-dependent activation, but also in its novel phosphate-binding motif. Its catalytic site does not utilize backbone amide groups to bind phosphate, as in the common P-loop, but contains a large cluster of acidic ionizable side chains. By employing a state-of-the-art computational approach, we have revealed a previously uncharacterized catalytic calcium-binding site in hSCAN-1, which elucidates the unusual calcium-dependence of its apyrase activity. In a high-order coordination shell, the newly identified calcium ion organizes the active site residues to mediate nucleotide binding, to orient the nucleophilic water, and to facilitate the phosphoryl transfer reaction. From ab initio QM/MM molecular dynamics simulations with umbrella sampling, we have characterized a reverse protonation catalytic mechanism for hSCAN-1 and determined its free energy reaction profile. Our results are consistent with available experimental studies and provide new detailed insight into the structure-function relationship of this novel calcium-activated phosphoryl transfer enzyme.

Identification of Secondary Binding Sites on Protein Surfaces for Rational Elaboration of Synthetic Protein Mimics.

Minimal mimics of protein conformations provide rationally designed ligands to modulate protein function. The advantage of minimal mimics is that they can be chemically synthesized and coaxed to be proteolytically resistant; a key disadvantage is that minimization of the protein binding epitope may be associated with loss of affinity and specificity. Several approaches to overcome this challenge may be envisioned, including deployment of covalent warheads and use of nonnatural residues to improve contacts with the binding surface. Herein, we describe our computational and experimental efforts to enhance the minimal protein mimics with fragments that can contact undiscovered binding pockets on Mdm2 and MdmX-two well-studied protein partners of p53.

Modulation of virus-induced NF-κB signaling by NEMO coiled coil mimics.

Protein-protein interactions featuring intricate binding epitopes remain challenging targets for synthetic inhibitors. Interactions of NEMO, a scaffolding protein central to NF-κB signaling, exemplify this challenge. Various regulators are known to interact with different coiled coil regions of NEMO, but the topological complexity of this protein has limited inhibitor design. We undertook a comprehensive effort to block the interaction between vFLIP, a Kaposi's sarcoma herpesviral oncoprotein, and NEMO using small molecule screening and rational design. Our efforts reveal that a tertiary protein structure mimic of NEMO is necessary for potent inhibition. The rationally designed mimic engages vFLIP directly causing complex disruption, protein degradation and suppression of NF-κB signaling in primary effusion lymphoma (PEL). NEMO mimic treatment induces cell death and delays tumor growth in a PEL xenograft model. Our studies with this inhibitor reveal the critical nexus of signaling complex stability in the regulation of NF-κB by a viral oncoprotein.

Intuitive Understanding of σ Delocalization in Loose and σ Localization in Tight Helical Conformations of an Oligosilane Chain.

Conformational effects on the σ-electron delocalization in oligosilanes are addressed by Hartree-Fock and time-dependent density functional theory calculations (B3LYP, 6-311G**) at MP2 optimized geometries of permethylated uniformly helical linear oligosilanes (all-ω-Sin R2n+2 ) up to n=16 and for backbone dihedral angles ω=55-180°. The extent of σ delocalization is judged by the partition ratio of the highest occupied molecular orbital and is reflected in the dependence of its shape and energy and of UV absorption spectra on n. The results agree with known spectra of all-transoid loose-helix conformers (all-[±165]-Sin Me2n+2 ) and reveal a transition at ω≈90° from the "σ-delocalized" limit at ω=180° toward and close to the physically non-realizable "σ-localized" tight-helix limit ω=0 with entirely different properties. The distinction is also obtained in the Hückel Ladder H and C models of σ delocalization. An easy intuitive way to understand the origin of the two contrasting limits is to first view the linear chain as two subchains with alternating primary and vicinal interactions (σ hyperconjugation), one consisting of the odd and the other of the even σ(SiSi) bonds, and then allow the two subchains to interact by geminal interactions (σ conjugation).

Resveratrol serves as a protein-substrate interaction stabilizer in human SIRT1 activation.

Resveratrol is a natural compound found in red wine that has been suggested to exert its potential health benefit through the activation of SIRT1, a crucial member of the mammalian NAD+-dependent deacetylases. SIRT1 has emerged as an attractive therapeutic target for many aging related diseases, however, how its activity can only be activated toward some specific substrates by resveratrol has been poorly understood. Herein, by employing extensive molecular dynamics simulations as well as fragment-centric topographical mapping of binding interfaces, we have clarified current controversies in the literature and elucidated that resveratrol plays an important activation role by stabilizing SIRT1/peptide interactions in a substrate-specific manner. This new mechanism highlights the importance of the N-terminal domain in substrate recognition, explains the activity restoration role of resveratrol toward some "loose-binding" substrates of SIRT1, and has significant implications for the rational design of new substrate-specific SIRT1 modulators.

Time-dependent density functional theory treatment of the first UV absorption band in all-transoid permethyloligosilanes SinMe2n + 2 (n = 2-8, 10).

The TD B3LYP/6-311G(d,p) method slightly overestimates the excitation energies of the first UV absorption band of the all-transoid conformers of SinMe2n + 2 (n = 2-8, 10), deduced from temperature-dependent measurements on conformer mixtures in hydrocarbon solvents, by a nearly constant amount (approximately 2000 cm-1). The TD B3LYP/6-31G(d) results are less satisfactory. The first band is calculated to be due to a sigma pi * excitation in Si2Me6 and to a superposition of overlapping sigma sigma * and sigma pi * excitations in the longer oligosilanes. The sigma pi * excitation is calculated to lie a little below the sigma sigma * excitation up to Si4Me10, the two transitions are nearly degenerate in Si5Me12, and the sigma sigma * excitation drops increasingly below the sigma pi * as the chain length is extended. The dipole strength of the sigma sigma * excitation grows by 4.8 D2 (D = debye) per added SiSi bond (more slowly up to n = 5) and the calculation overestimates it by a factor of about three. The sigma pi * excitation is computed to carry no or almost no oscillator strength, but as noted earlier by others, its presence is critical for the interpretation of the observed thermochromism. Upon cooling below room temperature, the first absorption maximum is blue-shifted in short chains and red-shifted in long chains. Unlike the prior investigators, we attribute the blue shift to the disappearance of hot bands built on the sigma pi * origin using intensity borrowing sigma-pi mixing vibrations (b1 in Si3Me8). As usual, the red shift is attributed to the disappearance of twisted conformers, which have higher calculated sigma sigma * excitation energies.