Machine Learning to Identify Flexibility Signatures of Class A GPCR Inhibition.

We show that machine learning can pinpoint features distinguishing inactive from active states in proteins, in particular identifying key ligand binding site flexibility transitions in GPCRs that are triggered by biologically active ligands. Our analysis was performed on the helical segments and loops in 18 inactive and 9 active class A G protein-coupled receptors (GPCRs). These three-dimensional (3D) structures were determined in complex with ligands. However, considering the flexible versus rigid state identified by graph-theoretic ProFlex rigidity analysis for each helix and loop segment with the ligand removed, followed by feature selection and k-nearest neighbor classification, was sufficient to identify four segments surrounding the ligand binding site whose flexibility/rigidity accurately predicts whether a GPCR is in an active or inactive state. GPCRs bound to inhibitors were similar in their pattern of flexible versus rigid regions, whereas agonist-bound GPCRs were more flexible and diverse. This new ligand-proximal flexibility signature of GPCR activity was identified without knowledge of the ligand binding mode or previously defined switch regions, while being adjacent to the known transmission switch. Following this proof of concept, the ProFlex flexibility analysis coupled with pattern recognition and activity classification may be useful for predicting whether newly designed ligands behave as activators or inhibitors in protein families in general, based on the pattern of flexibility they induce in the protein.

Protein-ligand interfaces are polarized: discovery of a strong trend for intermolecular hydrogen bonds to favor donors on the protein side with implications for predicting and designing ligand complexes.

Understanding how proteins encode ligand specificity is fascinating and similar in importance to deciphering the genetic code. For protein-ligand recognition, the combination of an almost infinite variety of interfacial shapes and patterns of chemical groups makes the problem especially challenging. Here we analyze data across non-homologous proteins in complex with small biological ligands to address observations made in our inhibitor discovery projects: that proteins favor donating H-bonds to ligands and avoid using groups with both H-bond donor and acceptor capacity. The resulting clear and significant chemical group matching preferences elucidate the code for protein-native ligand binding, similar to the dominant patterns found in nucleic acid base-pairing. On average, 90% of the keto and carboxylate oxygens occurring in the biological ligands formed direct H-bonds to the protein. A two-fold preference was found for protein atoms to act as H-bond donors and ligand atoms to act as acceptors, and 76% of all intermolecular H-bonds involved an amine donor. Together, the tight chemical and geometric constraints associated with satisfying donor groups generate a hydrogen-bonding lock that can be matched only by ligands bearing the right acceptor-rich key. Measuring an index of H-bond preference based on the observed chemical trends proved sufficient to predict other protein-ligand complexes and can be used to guide molecular design. The resulting Hbind and Protein Recognition Index software packages are being made available for rigorously defining intermolecular H-bonds and measuring the extent to which H-bonding patterns in a given complex match the preference key.

Heritable chromatin structure: mapping "memory" in histones H3 and H4.

Telomeric position effect in Saccharomyces cerevisiae is a chromatin-mediated phenomenon in which telomere proximal genes are repressed (silenced) in a heritable, but reversible, fashion. Once a transcriptional state (active or silenced) is established, however, there is a strong tendency for that state to be propagated. Twenty-five years ago, H. Weintraub and colleagues suggested that such heritability could be mediated by posttranslational modification of chromatin [Weintraub, H., Flint, S. J., Leffak, I. M., Groudine, M. & Grainger, R. M. (1977) Cold Spring Harbor Symp. Quant. Biol. 42, 401-407]. To identify potential sites within the chromatin that might act as sources of "memory" for the heritable transmission, we performed a genetic screen to isolate mutant alleles of the histones H3 and H4 genes that would "lock" telomeric marker genes into a silenced state. We identified mutations in the NH(2)-terminal tail and core of both histones; most of the amino acid changes mapped adjacent to lysines that are known sites of acetylation or methylation. We developed a method using MS to quantify the level of acetylation at each lysine within the histone H4 NH(2)-terminal tail in these mutants. We discovered that each of these mutants had a dramatic reduction in the level of acetylation at lysine 12 within the histone H4 tail. We propose that this lysine serves as a "memory mark" for propagating the expression state of a telomeric gene: when it is unacetylated, silent chromatin will be inherited; when it is acetylated an active state will be inherited.