Hydrophobic Complementarity Research Articles

Structure-Based Stabilization of Non-native Protein-Protein Interactions of Coronavirus Nucleocapsid Proteins in Antiviral Drug Design.

Structure-based stabilization of protein–protein interactions (PPIs) is a promising strategy for drug discovery. However, this approach has mainly focused on the stabilization of native PPIs, and non-native PPIs have received little consideration. Here, we identified a non-native interaction interface on the three-dimensional dimeric structure of the N-terminal domain of the MERS-CoV nucleocapsid protein (MERS-CoV N-NTD). The interface formed a conserved hydrophobic cavity suitable for targeted drug screening. By considering the hydrophobic complementarity during the virtual screening step, we identified 5-benzyloxygramine as a new N protein PPI orthosteric stabilizer that exhibits both antiviral and N-NTD protein-stabilizing activities. X-ray crystallography and small-angle X-ray scattering showed that 5-benzyloxygramine stabilizes the N-NTD dimers through simultaneous hydrophobic interactions with both partners, resulting in abnormal N protein oligomerization that was further confirmed in the cell. This unique approach based on the identification and stabilization of non-native PPIs of N protein could be applied toward drug discovery against CoV diseases.

Docking Troponin T onto the Tropomyosin Overlapping Domain of Thin Filaments

Complete description of thin filament conformational transitions accompanying muscle regulation requires ready access to atomic structures of actin-bound tropomyosin-troponin. To date, several molecular-docking protocols have been employed to identify troponin interactions on actin-tropomyosin because high-resolution experimentally determined structures of filament-associated troponin are not available. However, previously published all-atom models of the thin filament show chain separation and corruption of components during our molecular dynamics simulations of the models, implying artifactual subunit organization, possibly due to incorporation of unorthodox tropomyosin-TnT crystal structures and complex FRET measurements during model construction. For example, the recent Williams et al. (2016) atomistic model of the thin filament displays a paucity of salt bridges and hydrophobic complementarity between the TnT tail (TnT1) and tropomyosin, which is difficult to reconcile with the high, 20 nM Kd binding of TnT onto tropomyosin. Indeed, our molecular dynamics simulations show the TnT1 component in their model partially dissociates from tropomyosin in under 100 ns, whereas actin-tropomyosin and TnT1 models themselves remain intact. We therefore revisited computational work aiming to improve TnT1-thin filament models by employing unbiased docking methodologies, which test billions of trial rotations and translations of TnT1 over three-dimensional grids covering end-to-end bonded tropomyosin alone or tropomyosin on F-actin. We limited conformational searches to the association of well-characterized TnT1 helical domains and either isolated tropomyosin or actin-tropomyosin yet avoided docking TnT domains that lack known or predicted structure. The docking programs PIPER and ClusPro were used, followed by interaction energy optimization and extensive molecular dynamics. TnT1 docked to either side of isolated tropomyosin but uniquely onto one location of actin-bound tropomyosin. The antiparallel interaction with tropomyosin contained abundant salt bridges and intimately integrated hydrophobic networks joining TnT1 and the tropomyosin N-/C-terminal overlapping domain. The TnT1-tropomyosin linkage yields well-defined molecular crevices. Interaction energy measurements strongly favor this TnT1-tropomyosin design over previously proposed models.

Abstract 2174: Inhibition of CDK8 kinase with SEL120-34A allows for a personalized approach in AML

Abstract Inhibition of oncogenic transcriptional programs is recognized to be a promising therapeutic strategy. SEL120-34A is a novel inhibitor of Cyclin- dependent kinase 8 (CDK8), which regulates transcription by associating with the Mediator complex. SEL120-34A interacts with the ATP binding site of CDK8 in type I inhibitor fashion and forms several types of interactions, including halogen bonds with the protein’s hinge region and hydrophobic complementarities within its front pocket. Although the compound was only modestly active in solid tumor cell lines, it repressed phosphorylation of STAT5 Ser726 and could differentially inhibit viability of AML and ALL cell lines in vitro and in vivo, along with other type I CDK8 inhibitors. Transcriptomic analysis identified major transcriptional programs altered in responder cell lines, which strongly indicated that apart from repression of survival pathways, CDK8 inhibitors could induce differentiation in cell lines with leukemia stem cells characteristics. Further studies on a large panels of responder and non-responder cell lines identified robust biomarkers which could be used with high confidence for stratification and personalized approach in CDK8-dependent AML cases. Favorable pharmacokinetics, confirmed safety and in vivo efficacy in leukemia models provide the rationale for further clinical development of SEL120-34A. Citation Format: Tomasz Rzymski, Michal Mikula, Eliza Zylkiewicz, Agnieszka Dreas, Katarzyna Wiklik, Aniela Golas, Katarzyna Wojcik, Magdalena Masiejczyk, Iga Dudzicz, Katarzyna Kucwaj, Malgorzata Statkiewicz, Krzysztof Goryca, Aleksandra Grochowska, Aleksandra Cabaj, Jerzy Ostrowski, Urszula Kukliniska, Krzysztof Brzozka. Inhibition of CDK8 kinase with SEL120-34A allows for a personalized approach in AML [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 2174. doi:10.1158/1538-7445.AM2017-2174

SEL120-34A is a novel CDK8 inhibitor active in AML cells with high levels of serine phosphorylation of STAT1 and STAT5 transactivation domains.

Inhibition of oncogenic transcriptional programs is a promising therapeutic strategy. A substituted tricyclic benzimidazole, SEL120-34A, is a novel inhibitor of Cyclin-dependent kinase 8 (CDK8), which regulates transcription by associating with the Mediator complex. X-ray crystallography has shown SEL120-34A to be a type I inhibitor forming halogen bonds with the protein's hinge region and hydrophobic complementarities within its front pocket. SEL120-34A inhibits phosphorylation of STAT1 S727 and STAT5 S726 in cancer cells in vitro. Consistently, regulation of STATs- and NUP98-HOXA9- dependent transcription has been observed as a dominant mechanism of action in vivo. Treatment with the compound resulted in a differential efficacy on AML cells with elevated STAT5 S726 levels and stem cell characteristics. In contrast, resistant cells were negative for activated STAT5 and revealed lineage commitment. In vivo efficacy in xenotransplanted AML models correlated with significant repression of STAT5 S726. Favorable pharmacokinetics, confirmed safety and in vivo efficacy provide a rationale for the further clinical development of SEL120-34A as a personalized therapeutic approach in AML.

Delineation of autoantibody repertoire through differential proteogenomics in hepatitis C virus-induced cryoglobulinemia.

Antibodies cross-reactive to pathogens and autoantigens are considered pivotal in both infection control and accompanying autoimmunity. However, the pathogenic roles of autoantibodies largely remain elusive without a priori knowledge of disease-specific autoantigens. Here, through a novel quantitative proteogenomics approach, we demonstrated a successful identification of immunoglobulin variable heavy chain (VH) sequences highly enriched in pathological immune complex from clinical specimens obtained from a patient with hepatitis C virus-induced cryoglobulinemia (HCV-CG). Reconstructed single-domain antibodies were reactive to both HCV antigens and potentially liver-derived human proteins. Moreover, over the course of antiviral therapy, a substantial “de-evolution” of a distinct sub-repertoire was discovered, to which proteomically identified cryoprecipitation-prone autoantibodies belonged. This sub-repertoire was characterized by IGHJ6*03-derived, long, hydrophobic complementarity determining region (CDR-H3). This study provides a proof-of-concept of de novo mining of autoantibodies and corresponding autoantigen candidates in a disease-specific context in human, thus facilitating future reverse-translational research for the discovery of novel biomarkers and the development of antigen-specific immunotherapy against various autoantibody-related disorders.

The role of shape complementarity in the protein-protein interactions.

We use a dissipative particle dynamic simulation to investigate the effects of shape complementarity on the protein-protein interactions. By monitoring different kinds of protein shape-complementarity modes, we gave a clear mechanism to reveal the role of the shape complementarity in the protein-protein interactions, i.e., when the two proteins with shape complementarity approach each other, the conformation of lipid chains between two proteins would be restricted significantly. The lipid molecules tend to leave the gap formed by two proteins to maximize the configuration entropy, and therefore yield an effective entropy-induced protein-protein attraction, which enhances the protein aggregation. In short, this work provides an insight into understanding the importance of the shape complementarity in the protein-protein interactions especially for protein aggregation and antibody–antigen complexes. Definitely, the shape complementarity is the third key factor affecting protein aggregation and complex, besides the electrostatic-complementarity and hydrophobic complementarity.

Structure–kinetic relationship study of CDK8/CycC specific compounds

In contrast with the very well explored concept of structure-activity relationship, similar studies are missing for the dependency between binding kinetics and compound structure of a protein ligand complex, the structure-kinetic relationship. Here, we present a structure-kinetic relationship study of the cyclin-dependent kinase 8 (CDK8)/cyclin C (CycC) complex. The scaffold moiety of the compounds is anchored in the kinase deep pocket and extended with diverse functional groups toward the hinge region and the front pocket. These variations can cause the compounds to change from fast to slow binding kinetics, resulting in an improved residence time. The flip of the DFG motif ("DMG" in CDK8) to the inactive DFG-out conformation appears to have relatively little influence on the velocity of binding. Hydrogen bonding with the kinase hinge region contributes to the residence time but has less impact than hydrophobic complementarities within the kinase front pocket.

Semi‐Empirical Predictions on the Structure and Properties of ent‐Kaurenoic Acid and Derivatives

The physicochemical properties of ent‐ kaurenoic acid model derivatives, which possibly influence its therapeutic application, were calculated. Results revealed that the molecule possess favourable attributes which renders it possible to be considered as a drug lead only that its very hydrophobic nature can result to poor bioavailabilty, low absorption and poor systemic circulation. In silico simulations revealed that this setback can be overcome by introduction of hydroxyl group to the tertiary carbon of ent‐kaurenoic acid employing m‐CPBA catalyzed hydroxylation, thus, unleashing its full drug potency. Moreover, molecular similarity analyses derived from semi‐empirical calculations between ent‐kaurenoic acid and a set of kaurane diterpenoids showed differences in hydrophobic complementarity, size and electronic properties despite possessing nearly identical molecular frameworks, thus, arriving in a generalization for their observed mechanistic differences on acting on different targets.

Analysis of hydrophobic interactions of antagonists with the beta2-adrenergic receptor

The adrenergic receptors mediate a wide variety of physiological responses, including vasodilatation and vasoconstriction, heart rate modulation, and others. Beta-adrenergic antagonists (‘beta-blockers’) thus constitute a widely used class of drugs in cardiovascular medicine as well as in management of anxiety, migraine, and glaucoma. The importance of the hydrophobic effect has been evidenced for a wide range of beta-blocker properties. To better understand the role of the hydrophobic effect in recognition of beta-blockers by their receptor, we carried out a molecular docking study combined with an original approach to estimate receptor–ligand hydrophobic interactions. The proposed method is based on automatic detection of molecular fragments in ligands and the analysis of their interactions with receptors separately. A series of beta-blockers, based on phenylethanolamines and phenoxypropanolamines, were docked to the beta2-adrenoceptor binding site in the crystal structure. Hydrophobic complementarity between the ligand and the receptor was calculated using the PLATINUM web-server (http://model.nmr.ru/platinum). Based on the analysis of the hydrophobic match for molecular fragments of beta-blockers, we have developed a new scoring function which efficiently predicts dissociation constant (pKd) with strong correlations (r 2∼ 0.8) with experimental data.

Electrostatics in computer-aided drug design

Selective binding of a drug by a biopolymer is influenced by three major factors: steric, Coulombic, and hydrophobic complementarity. Although the steric fit plays a determining role in binding, the other two factors may be of primary importance, too. We treat Coulombic complementarity as the matching between positive and negative potential regions of the biopolymer and its ligand. On the other hand, hydrophobic complementarity, related to hydration ability, is denned using electrostatic fields. Highland low-field regions of the interacting species should match in order to maximize the net free energy of association in the biophase. According to the above definitions, Coulombic and hydrophobic similarity, or the binding of different molecules to the same receptor by about the same strength, correspond to the more-or-less approximate congruence of the electrostatic potential and field patterns of the molecules compared. The above principles are discussed for a variety of cases studied in European laboratories.

Protein Folding in Membranes: Insights from Neutron Diffraction Studies of a Membrane β-Sheet Oligomer

Studies of the assembly of the hexapeptide Acetyl-Trp-Leu 5 (AcWL 5) into β-sheets in membranes have provided insights into membrane protein folding. Yet, the exact structure of the oligomer in the lipid bilayer is unknown. Here we use neutron diffraction to study the disposition of the peptides in bilayers. We find that pairs of adjacent deuterium-labeled leucines have no well-defined peak or dip in the transmembrane distribution profiles, indicative of heterogeneity in the depth of membrane insertion. At the same time, the monomeric homolog AcWL 4 exhibits a homogeneous, well-defined, interfacial location in neutron diffraction experiments. Thus, although the bilayer location of monomeric AcWL 4 is determined by hydrophobicity matching or complementarity within the bilayer, the AcWL 5 molecules in the oligomer are positioned at different depths within the bilayer because they assemble into a staggered transmembrane β-sheet. The AcWL 5 assembly is dominated by protein-protein interactions rather than hydrophobic complementarity. These results have implications for the structure and folding of proteins in their native membrane environment and highlight the importance of the interplay between hydrophobic complementarity and protein-protein interactions in determining the structure of membrane proteins.

Ligand-specific scoring functions: improved ranking of docking solutions1

Molecular docking is a powerful computational method that has been widely used in many biomolecular studies to predict geometry of a protein-ligand complex. However, while its conformational search algorithms are usually able to generate correct conformation of a ligand in the binding site, the scoring methods often fail to discriminate it among many false variants. We propose to treat this problem by applying more precise ligand-specific scoring filters to re-rank docking solutions. In this way specific features of interactions between protein and different types of compounds can be implicitly taken into account. New scoring functions were constructed including hydrogen bonds, hydrophobic and hydrophilic complementarity terms. These scoring functions also discriminate ligands by the size of the molecule, the total hydrophobicity, and the number of peptide bonds for peptide ligands. Weighting coefficients of the scoring functions were adjusted using a training set of 60 protein–ligand complexes. The proposed method was then tested on the results of docking obtained for an additional 70 complexes. In both cases the success rate was 5–8% better compared to the standard functions implemented in popular docking software. 1Presented at CMTPI 2007: Computational Methods in Toxicology and Pharmacology Integrating Internet Resources (Moscow, Russia, September 1–5, 2007).

Exploring blocker binding to a homology model of the open hERG K + channel using docking and molecular dynamics methods

Binding of blockers to the human voltage-gated hERG potassium channel is studied using a combination of homology modelling, automated docking calculations and molecular dynamics simulations, where binding affinities are evaluated using the linear interaction energy method. A homology model was constructed based on the available crystal structure of the bacterial KvAP channel and the affinities of a series of sertindole analogues predicted using this model. The calculations reproduce the relative binding affinities of these compounds very well and indicate that both polar interactions near the intracellular opening of the selectivity filter as well as hydrophobic complementarity in the region around F656 are important for blocker binding. These results are consistent with recent alanine scanning mutation experiments on the blocking of the hERG channel by other compounds.

Structural Requirements for Factor Xa Inhibition by 3-Oxybenzamides with Neutral P1 Substituents: Combining X-ray Crystallography, 3D-QSAR, and Tailored Scoring Functions

The design, synthesis, and structure-activity relationship of 3-oxybenzamides as potent inhibitors of the coagulation protease factor Xa are described on the basis of X-ray structures, privileged structure motifs, and SAR information. A total of six X-ray structures of fXa/inhibitor complexes led us to identify the major protein-ligand interactions. The binding mode is characterized by a lipophilic dichlorophenyl substituent interacting with Tyr228 in the protease S1 pocket, while polar parts are accommodated in S4. This alignment in combination with docking allowed derivation of 3D-QSAR models and tailored scoring functions to rationalize biological affinity and provide guidelines for optimization. The resulting models showed good correlation coefficients and predictions of external test sets. Furthermore, they correspond to binding site topologies in terms of steric, electrostatic, and hydrophobic complementarity. Two approaches to derive tailored scoring functions combining binding site and ligand information led to predictive models with acceptable predictions of the external set. Good correlations to experimental affinities were obtained for both AFMoC (adaptation of fields for molecular comparison) and the novel TScore function. The SAR information from 3D-QSAR and tailored scoring functions agrees with all experimental data and provides guidelines and reasonable activity estimations for novel fXa inhibitors.

Hydrophobic complementarity in protein-protein docking.

Formation of hydrophobic contacts across a newly formed interface is energetically favorable. Based on this observation we developed a geometric-hydrophobic docking algorithm that estimates quantitatively the hydrophobic complementarity at protein-protein interfaces. Each molecule to be docked is represented as a grid of complex numbers, storing information regarding the shape of the molecule in the real part and information regarding the hydropathy of the surface in the imaginary part. The grid representations are correlated using fast Fourier transformations. The algorithm is used to compare the extent of hydrophobic complementarity in oligomers (represented by D2 tetramers) and in hetero-dimers of soluble proteins (complexes). We also test the implication of hydrophobic complementarity in distinguishing correct from false docking solutions. We find that hydrophobic complementarity at the interface exists in oligomers and in complexes, and in both groups the extent of such complementarity depends on the size of the interface. Thus, the non-polar portions of large interfaces are more often juxtaposed than non-polar portions of small interfaces. Next we find that hydrophobic complementarity helps to point out correct docking solutions. In oligomers it significantly improves the ranks of nearly correct reassembled and modeled tetramers. Combining geometric, electrostatic and hydrophobic complementarity for complexes gives excellent results, ranking a nearly correct solution < 10 for 5 of 23 tested systems, < 100 for 8 systems and < 1000 for 19 systems.

On proteins, grids, correlations, and docking

The activity of a living cell can be portrayed as a network of interactions involving proteins and nucleic acids that transfer biological information. Intervention in cellular processes requires thorough understanding of the interactions between the molecules, which can be provided by docking techniques. Docking methods attempt to predict the structures of complexes given the structures of the component molecules. We focus hereby on protein–protein docking procedures that employ grid representations of the molecules, and use correlation for searching the solution space and evaluating putative complexes. Geometric surface complementarity is the dominant descriptor in docking. Inclusion of electrostatics often improves the results of geometric docking for soluble proteins, whereas hydrophobic complementarity is more important in construction of oligomers. Using binding-site information in the scan or as a filter helps to identify and up-rank nearly correct solutions. To cite this article: M. Eisenstein, E. Katchalski-Katzir, C. R. Biologies 327 (2004).

Structural requirements for inhibition of the neuronal nitric oxide synthase (NOS-I): 3D-QSAR analysis of 4-oxo- and 4-amino-pteridine-based inhibitors.

The family of homodimeric nitric oxide synthases (NOS I-III) catalyzes the generation of the cellular messenger nitric oxide (NO) by oxidation of the substrate L-arginine. The rational design of specific NOS inhibitors is of therapeutic interest in regulating pathological NO levels associated with sepsis, inflammatory, and neurodegenerative diseases. The cofactor (6R)-5,6,7,8-tetrahydrobiopterin (H(4)Bip) maximally activates all NOSs and stabilizes enzyme quaternary structure by promoting and stabilizing dimerization. Here, we describe the synthesis and three-dimensional (3D) quantitative structure-activity relationship (QSAR) analysis of 65 novel 4-amino- and 4-oxo-pteridines (antipterins) as inhibitors targeting the H(4)Bip binding site of the neuronal NOS isoform (NOS-I). The experimental binding modes for two inhibitors complexed with the related endothelial NO synthase (NOS-III) reveal requirements of biological affinity and form the basis for ligand alignment. Different alignment rules were derived by building other compounds accordingly using manual superposition or a genetic algorithm for flexible superposition. Those alignments led to 3D-QSAR models (comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA)), which were validated using leave-one-out cross-validation, multiple analyses with two and five randomly chosen cross-validation groups, perturbation of biological activities by randomization or progressive scrambling, and external prediction. An iterative realignment procedure based on rigid field fit was used to improve the consistency of the resulting partial least squares models. This led to consistent and highly predictive 3D-QSAR models with good correlation coefficients for both CoMFA and CoMSIA, which correspond to experimentally determined NOS-II and -III H(4)Bip binding site topologies as well as to the NOS-I homology model binding site in terms of steric, electrostatic, and hydrophobic complementarity. These models provide clear guidelines and accurate activity predictions for novel NOS-I inhibitors.

Design and quantitative structure-activity relationship of 3-amidinobenzyl-1H-indole-2-carboxamides as potent, nonchiral, and selective inhibitors of blood coagulation factor Xa.

A series of 138 nonchiral 3-amidinobenzyl-1H-indole-2-carboxamides and analogues as inhibitors of the blood coagulation enzyme factor Xa (fXa) were designed, synthesized, and investigated by X-ray structure analysis and 3D quantitative structure-activity relationship (QSAR) studies (CoMFA, CoMSIA) in order to identify important protein-ligand interactions responsible for biological affinity and selectivity. Several compounds from this series are highly potent and selective inhibitors of this important enzyme linking extrinsic and intrinsic coagulation pathways. To rationalize biological affinity and to provide guidelines for further design, all compounds were docked into the factor Xa binding site. Those docking studies were based on X-ray structures of factor Xa in complex with literature-known inhibitors. It was possible to validate those binding modes by four X-ray crystal structures of representative ligands in factor Xa, while one ligand was additionally crystallized in trypsin to rationalize requirements for selective factor Xa inhibition. The 3D-QSAR models based on a superposition rule derived from these docking studies were validated using conventional and cross-validated r(2) values using the leave-one-out method and repeated analyses using two randomly chosen cross-validation groups plus randomization of biological activities. This led to consistent and highly predictive 3D-QSAR models with good correlation coefficients for both CoMFA and CoMSIA, which were found to correspond to experimentally determined factor Xa binding site topology in terms of steric, electrostatic, and hydrophobic complementarity. Subsets selected as smaller training sets using 2D fingerprints and maximum dissimilarity methods resulted in 3D-QSAR models with remarkable correlation coefficients and a high predictive power. The final quantitative SAR information agrees with all experimental data for the binding topology and thus provides reasonable activity predictions for novel factor Xa inhibitors.

Design of three-dimensional, millimeter-scale models for molecular folding.

This communication describes the fabrication of three-dimensional structures of organic polymers using principles of design inspired by protein folding. The structures consist of rigid polyhedral components with dimensions of a few millimeters ("microdomains"), representing alpha-helical and beta-sheet secondary structures, connected with flexible linkers representing loops or turns. These structures were fabricated from polyurethane using photolithographic and soft lithographic techniques. The surfaces of the microdomains were patterned into hydrophobic and hydrophilic regions, and a hydrophobic photocurable liquid (serving both as lubricant and adhesive) was selectively precipitated onto the hydrophobic areas. The unfolded structures were suspended in water and agitated by tumbling. Self-assembly occurred through coalescence of the thin films of hydrophobic liquid, and was caused by minimization of the free energy of the interface between the liquid adhesive and the water. The self-assembled structures were locked in place by curing the adhesive with UV light. These results demonstrate the use of concepts abstracted from the study of proteins-including attractive hydrophobic interactions, shape complementarity, and conformational constraint-in the self-assembly of complex, three-dimensional structures on the millimeter scale.

Scoring Protein-Substrate Complex Decoys Based on Electrostatic and Hydrophobic Complementarities: Evaluation of Hormone-Receptor Interactions

Frequently, hitherto automatic systems for rigid body docking of a substrate into a receptor output series of protein-substrate decoys differing from each other in the orientation of the ligand molecule within the receptor cleft [3, 4, 5]. Then, using several schemes based mainly on the electrostatics of the atoms composing the cleft and the ligand, a scoring of the decoys is carried out, the most plausible structure being that with the highest score. Commonly used properties to score the decoys are the analysis of putative hydrogen bonds, ion-ion interactions, hydrophobic interaction and other electrophysical properties of the atoms, the evaluations being carried on an atomic bases which most of the time is costly in terms of processing times besides relying heavily on a specialist knowledge. Here we have developed a completely automatic method to perform the scoring of protein-substrate decoys, which is a combination of a well known pattern recognition method and a filtering process in order to detect complementary patterns from the electrostatic point of view of regions on the surfaces of both the ligand molecule and the protein cleft. The result is a swift and concise methodology to score the decoys by degree of electrostatic and hydrophobic complementarity.