Protein-ligand interactions: probing the energetics of a putative cation-π interaction

Discerning Tyrosine Phosphorylation from Multiple Phosphorylations Using a Nanofluidic Logic Platform

Discerning tyrosine phosphorylation (pTyr) catalyzed by Tyr kinase is central to the revelation of oncogenic mechanisms and the development of targeted anticancer drugs. Despite some techniques, this goal remains challenging, especially when faced with the interference of multiple phosphorylation events, including serine (pSer) and threonine phosphorylation (pThr). We describe here a functional polymer-modified artificial ion nanochannel, which enables the sensitive and selective recognition of phosphotyrosine (pY) peptide by the distinct ionic current change. Such a recognition effect allows for the nanochannel to work in a complex protein digest condition. Further, the implementation of nanofluidic logic functions with the addition of Ca²⁺ dramatically improves the selectivity of the nanochannel to pY peptide and thus can discern pTyr by the Tyr kinase from pSer by the Ser/Thr kinase through simultaneously monitoring multisite phosphorylation at the same or different peptide substrates in one-pot. This logic sensing platform displays the potential in differentiating Tyr kinase and Ser/Thr kinase and assessing multi-kinase activities in multi-targeted drug design.

Synthesis and structural characterization of a monocarboxylic inhibitor for GRB2 SH2 domain

A monocarboxylic inhibitor was designed and synthesized to disrupt the protein-protein interaction (PPI) between GRB2 and phosphotyrosine-containing proteins. Biochemical characterizations show compound 7 binds with the Src homology 2 (SH2) domain of GRB2 and is more potent than EGFR₁₀₆₈ phosphopeptide 14-mer. X-ray crystallographic studies demonstrate compound 7 occupies the GRB2 binding site for phosphotyrosine-containing sequences and reveal key structural features for GRB2-inhibitor binding. This compound with a -1 formal charge offers a new direction for structural optimization to generate cell-permeable inhibitors for this key protein target of the aberrant Ras-MAPK signaling cascade.

Potential of ATP toward Prevention of hIAPP Oligomerization and Destabilization of hIAPP Protofibrils: An In Silico Perspective

The aggregation of an intrinsically disordered protein, human islet amyloid polypeptide (hIAPP), leads to one of the most prevalent endocrine disorders, type II diabetes mellitus (T2DM). Hence inhibition of hIAPP aggregation provides a possible therapeutic approach for the treatment of T2DM. In this regard, a new aspect of adenosine triphosphate (ATP), which is widely known as the energy source for biological reactions, has recently been discovered, where it can inhibit the formation of protein aggregates and simultaneously dissolve preformed aggregates at a millimolar concentration scale. In this work, we investigate the effect of ATP on the aggregation of an amyloidogenic segment of hIAPP, hIAPP_22-28, and also of the full length sequence. Using all-atom classical molecular dynamics simulations, we observe that the tendency of hIAPP to oligomerize into β-sheet conformers is inhibited by ATP, due to which the peptides remain distant, loosely packed random monomers. Moreover, it can also disassemble preformed hIAPP protofibrils. ATP preferentially interacts with the hydrophobic residues of hIAPP_22-28 fragment and the terminal and turn residues of the full length peptide. The hydrogen bonding, hydrophobic, π-π, and N-H-π stacking interactions are the driving forces for the ATP induced inhibition of hIAPP aggregation. Interestingly, the hydrophobic adenosine of ATP is found to be more in contact with the peptide residues than the hydrophilic triphosphate moiety. The insight into the inhibitory mechanism of ATP on hIAPP aggregation can prove to be beneficial for the design of novel amyloid inhibitors in the future.

Some thermodynamic effects of varying nonpolar surfaces in protein-ligand interactions

Understanding how making structural changes in small molecules affects their binding affinities for targeted proteins is central to improving strategies for rational drug design. To assess the effects of varying the nature of nonpolar groups upon binding entropies and enthalpies, we designed and prepared a set of Grb2-SH2 domain ligands, Ac-pTyr-Ac6c-Asn-(CH2)n-R, in which the size and electrostatic nature of R groups at the pTyr+3 site were varied. The complexes of these ligands with the Grb2-SH2 domain were evaluated in a series of studies in which the binding thermodynamics were determined using isothermal titration calorimetry, and binding interactions were examined in crystallographic studies of two different complexes. Notably, adding nonpolar groups to the pTyr+3 site leads to higher binding affinities, but the magnitude and energetic origins of these effects vary with the nature of the R substituent. For example, enhancements to binding affinities using aliphatic R groups are driven by more favorable changes in binding entropies, whereas aryl R groups improve binding free energies through a combination of more favorable changes in binding enthalpies and entropies. However, enthalpy/entropy compensation plays a significant role in these associations and mitigates against any significant variation in binding free energies, which vary by only 0.8 kcal•mol-1, with changes in the electrostatic nature and size of the R group. Crystallographic studies show that differences in ΔG° or ΔH° correlate with buried nonpolar surface area, but they do not correlate with the total number of polar or van der Waals contacts. The relative number of ordered water molecules and relative order in the side chains at pTyr+3 correlate with differences in -TΔS°. Overall, these studies show that burial of nonpolar surface can lead to enhanced binding affinities arising from dominating entropy- or enthalpy-driven hydrophobic effects, depending upon the electrostatic nature of the apolar R group.

Functional Nanochannels for Sensing Tyrosine Phosphorylation

Tyrosine phosphorylation (pTyr), much of which occurred on localized multiple sites, initiates cellular signaling, governs cellular functions, and its dysregulation is implicated in many diseases, especially cancers. pTyr-specific sensing is of great significance for understanding disease states and developing targeted anticancer drugs, however, it is very challenging due to the slight difference from serine (pSer) or threonine phosphorylation (pThr). Here we present polyethylenimine-g-phenylguanidine (PEI-PG)-modified nanochannels that can address the challenge. Rich guanidinium groups enabled PEI-PG to form multiple interactions with phosphorylated residues, especially pTyr residue, which triggered the conformational change of PEI-PG. By taking advantage of the "OFF-ON" change of the ion flux arising from the conformational shrinkage of the grafted PEI-PG, the nanochannels could distinguish phosphorylated peptide (PP) from nonmodified peptide, recognize PPs with pSer, pThr, or pTyr residue and PPs with different numbers of identical residues, and importantly could sense pTyr peptides in a biosample. Benefiting from the strong interaction between the guanidinium group and the pTyr side-chain, the specific sensing of pTyr peptide was achieved by performing a simple logic operation based on PEI-PG-modified nanochannels when Ca²⁺ was introduced as an interferent. The excellent pTyr sensing capacity makes the nanochannels available for real-time monitoring of the pTyr process by c-Abl kinase on a peptide substrate, even under complicated conditions, and the proof-of-concept study of monitoring the kinase activity demonstrates its potential in kinase inhibitor screening.

Structure-thermodynamics-relationships of hepatitis C viral NS3 protease inhibitors

Thermodynamic parameters were determined for structurally-related inhibitors of HCV NS3 protease to assess how binding entropies and enthalpies vary with incremental changes at the P2 and P3 inhibitor subsites. Changing the heterocyclic substituent at P2 from a pyridyl to a 7-methoxy-2-phenyl-4-quinolyl group leads to a 710-fold increase in affinity. Annelating a benzene ring onto a pyridine ring leads to quinoline-derived inhibitors having higher affinities, but the individual enthalpy and entropy contributions are markedly different for each ligand pair. Introducing a phenyl group at C2 of the heterocyclic ring at P2 uniformly leads to higher affinity analogs with more favorable binding entropies, while adding a methoxy group at C7 of the quinoline ring at P2 provides derivatives with more favorable binding enthalpies. Significant enthalpy/entropy compensation is observed for structural changes made to inhibitors lacking a 2-phenyl substituent, whereas favorable changes in both binding enthalpies and entropies accompany structural modifications when a 2-phenyl group is present. Overall, binding energetics of inhibitors having a 2-phenyl-4-quinolyl group at P2 are dominated by entropic effects, whereas binding of the corresponding norphenyl analogs are primarily enthalpy driven. Notably, the reversal from an entropy driven association to an enthalpy driven one for this set of inhibitors also correlates with alternate binding modes. When the steric bulk of the side chain at P3 is increased from a hydrogen atom to a tert-butyl group, there is a 770-fold improvement in affinity. The 30-fold increase resulting from the first methyl group is solely the consequence of a more favorable change in entropy, whereas subsequent additions of methyl groups leads to modest increases in affinity that arise primarily from incremental improvements in binding enthalpies accompanied with smaller favorable entropic contributions.

Thermodynamic consequences of Tyr to Trp mutations in the cation-π-mediated binding of trimethyllysine by the HP1 chromodomain

Evolution has converged on cation–π interactions for recognition of quaternary alkyl ammonium groups such as trimethyllysine (Kme3). While computational modelling indicates that Trp provides the strongest cation–π interaction of the native aromatic amino acids, there is limited corroborative data from measurements within proteins. Herein we investigate a Tyr to Trp mutation in the binding pocket of the HP1 chromodomain, a reader protein that recognizes Kme3. Binding studies demonstrate that the Trp-mediated cation–π interaction is about −5 kcal mol⁻¹ stronger, and the Y24W crystal structure shows that the mutation is not perturbing. Quantum mechanical calculations indicate that greater enthalpic binding is predominantly due to increased cation–π interactions. NMR studies indicate that differences in the unbound state of the Y24W mutation lead to enthalpy–entropy compensation. These results provide direct experimental quantification of Trp versus Tyr in a cation–π interaction and afford insight into the conservation of aromatic cage residues in Kme3 reader domains.

In this work, we experimentally validate that tryptophan provides the strongest cation–π binding interaction among aromatic amino acids and also lend insight into the importance of residue identity in trimethyllysine recognition by reader proteins.

Norepinephrine Inhibits Alzheimer's Amyloid-β Peptide Aggregation and Destabilizes Amyloid-β Protofibrils: A Molecular Dynamics Simulation Study

The abnormal self-assembly of amyloid-β (Aβ) peptides into toxic fibrillar aggregates is associated with the pathogenesis of Alzheimer's disease (AD). The inhibition of β-sheet-rich oligomer formation is considered as the primary therapeutic strategy for AD. Previous experimental studies reported that norepinephrine (NE), one of the neurotransmitters, is able to inhibit Aβ aggregation and disaggregate the preformed fibrils. Moreover, exercise can markedly increase the level of NE. However, the underlying inhibitory and disruptive mechanisms remain elusive. In this work, we performed extensive replica-exchange molecular dynamic (REMD) simulations to investigate the conformational ensemble of Aβ_1-42 dimer with and without NE molecules. Our results show that without NE molecules, Aβ_1-42 dimer transiently adopts a β-hairpin-containing structure, and the β-strand regions of this β-hairpin (residues 15QKLVFFA21 and 33GLMVGGVV40) strongly resemble those of the Aβ fibril structure (residues 15QKLVFFA21 and 30AIIGLMVG37) reported in an electron paramagnetic resonance spectroscopy study. NE molecules greatly reduce the interpeptide β-sheet content and suppress the formation of the above-mentioned β-hairpin, leading to a more disordered coil-rich Aβ dimer. Five dominant binding sites are identified, and the central hydrophobic core 16KLVFFA21 site and C-terminal 31IIGLMV36 hydrophobic site are the two most favorable ones. Our data reveal that hydrophobic, aromatic stacking, hydrogen-bonding and cation-π interactions synergistically contribute to the binding of NE molecules to Aβ peptides. MD simulations of Aβ_1-42 protofibril show that NE molecules destabilize Aβ protofibril by forming H-bonds with residues D1, A2, D23, and A42. This work reveals the molecular mechanism by which NE molecules inhibit Aβ_1-42 aggregation and disaggregate Aβ protofibrils, providing valuable information for developing new drug candidates and exercise therapy against AD.

Structural and functional properties of Grb2 SH2 dimer in CD28 binding

Growth factor receptor-bound protein 2 (Grb2) is an adaptor protein that plays a critical role in cellular signal transduction. It contains a central Src homology 2 (SH2) domain flanked by two Src homology 3 (SH3) domains. Binding of Grb2 SH2 to the cytoplasmic region of CD28, phosphorylated Tyr (pY) containing the peptide motif pY-X-N-X, is required for costimulatory signaling in T cells. In this study, we purified the dimer and monomer forms of Grb2 SH2, respectively, and analyzed their structural and functional properties. Size exclusion chromatography analysis showed that both dimer and monomer exist as stable states. Thermal stability analysis using circular dichroism showed that the dimer mostly dissociates into the monomer around 50°C. CD28 binding experiments showed that the affinity of the dimer to the phosphopeptide was about three fold higher than that of the monomer, possibly due to the avidity effect. The present crystal structure analysis of Grb2 SH2 showed two forms; one is monomer at 1.15 Å resolution, which is currently the highest resolution analysis, and another is dimer at 2.00 Å resolution. In the dimer structure, the C-terminal region, comprising residues 123–152, was extended towards the adjacent molecule, in which Trp121 was the hinge residue. The stable dimer purified using size exclusion chromatography would be due to the C-terminal helix “swapping”. In cases where a mutation caused Trp121 to be replaced by Ser in Grb2 SH2, this protein still formed dimers, but lost the ability to bind CD28.

Amino acid modifications for conformationally constraining naturally occurring and engineered peptide backbones: Insights from the Protein Data Bank

Extensive efforts invested in understanding the rules of protein folding are now being applied, with good effect, in de novo design of proteins/peptides. For proteins containing standard α-amino acids alone, knowledge derived from experimentally determined three-dimensional (3D) structures of proteins and biologically active peptides are available from the Protein Data Bank (PDB), and the Cambridge Structural Database (CSD). These help predict and design protein structures, with reasonable confidence. However, our knowledge of 3D structures of biomolecules containing backbone modified amino acids is still evolving. A major challenge in de novo protein/peptide design concerns the engineering of conformationally constrained molecules with specific structural elements and chemical groups appropriately positioned for biological activity. This review explores four classes of amino acid modifications that constrain protein/peptide backbone structure. Systematic analysis of peptidic molecule structures (eg, bioactive peptides, inhibitors, antibiotics, and designed molecules), containing these backbone-modified amino acids, found in the PDB and CSD are discussed. The review aims to provide structure-function insights that will guide future design of proteins/peptides.

Natural Products and Their Mimics as Targets of Opportunity for Discovery

Diverse structural types of natural products and their mimics have served as targets of opportunity in our laboratory to inspire the discovery and development of new methods and strategies to assemble polyfunctional and polycyclic molecular architectures. Furthermore, our efforts toward identifying novel compounds having useful biological properties led to the creation of new targets, many of which posed synthetic challenges that required the invention of new methodology. In this Perspective, selected examples of how we have exploited a diverse range of natural products and their mimics to create, explore, and solve a variety of problems in chemistry and biology will be discussed. The journey was not without its twists and turns, but the unexpected often led to new revelations and insights. Indeed, in our recent excursion into applications of synthetic organic chemistry to neuroscience, avoiding the more-traveled paths was richly rewarding.

Probing the role of cation-π interaction in the thermotolerance and catalytic performance of endo-polygalacturonases

Understanding the dynamics of the key pectinase, polygalacturonase, and improving its thermotolerance and catalytic efficiency are of importance for the cost-competitive bioconversion of pectic materials. By combining structure analysis and molecular dynamics (MD) simulations, eight mutagenesis sites having the potential to form cation-π interactions were identified in the widely used fungal endo-polygalacturonase PG63. In comparison to the wild-type, three single mutants H58Y, T71Y and T304Y showed improved thermostability (the apparent Tms increased by 0.6-3.9 °C) and catalytic efficiency (by up to 32-fold). Chromatogram analysis of the hydrolysis products indicated that a larger amount of shorter sugars were released from the polygalacturonic acid by these three mutants than by the wild-type. MD analysis of the enzyme-substrate complexes illustrated that the mutants with introduced cation-π interaction have modified conformations of catalytic crevice, which provide an enviable environment for the catalytic process. Moreover, the lower plasticity of T3 loop 2 at the edge of the subsite tunnel appears to recruit the reducing ends of oligogalacturonide into the active site tunnel and initiates new hydrolysis reactions. This study demonstrates the importance of cation-π interaction in protein conformation and provides a realistic strategy to enhance the thermotolerance and catalytic performance of endo-polygalacturonases.

Conformational Ensemble of hIAPP Dimer: Insight into the Molecular Mechanism by which a Green Tea Extract inhibits hIAPP Aggregation

Small oligomers formed early along human islet amyloid polypeptide (hIAPP) aggregation is responsible for the cell death in Type II diabetes. The epigallocatechin gallate (EGCG), a green tea extract, was found to inhibit hIAPP fibrillation. However, the inhibition mechanism and the conformational distribution of the smallest hIAPP oligomer – dimer are mostly unknown. Herein, we performed extensive replica exchange molecular dynamic simulations on hIAPP dimer with and without EGCG molecules. Extended hIAPP dimer conformations, with a collision cross section value similar to that observed by ion mobility-mass spectrometry, were observed in our simulations. Notably, these dimers adopt a three-stranded antiparallel β-sheet and contain the previously reported β-hairpin amyloidogenic precursor. We find that EGCG binding strongly blocks both the inter-peptide hydrophobic and aromatic-stacking interactions responsible for inter-peptide β-sheet formation and intra-peptide interaction crucial for β-hairpin formation, thus abolishes the three-stranded β-sheet structures and leads to the formation of coil-rich conformations. Hydrophobic, aromatic-stacking, cation-π and hydrogen-bonding interactions jointly contribute to the EGCG-induced conformational shift. This study provides, on atomic level, the conformational ensemble of hIAPP dimer and the molecular mechanism by which EGCG inhibits hIAPP aggregation.