Using small, flat aromatic rings as components of fragments or molecules is a common practice in fragment-based drug discovery and lead optimization. With an increasing focus on the exploration of novel biological and chemical space, and their improved synthetic accessibility, 3D fragments are attracting increasing interest. This study presents a detailed analysis of 3D and 2D ring fragments in marketed drugs. Several measures of properties were used, such as the type of ring assemblies and molecular shapes. The study also took into account the relationship between protein classes targeted by each ring fragment, providing target-specific information. The analysis shows the high structural and shape diversity of 3D ring systems and their importance in bioactive compounds. Major differences in 2D and 3D fragments are apparent in ligands that bind to the major drug targets such as GPCRs, ion channels, and enzymes.

Conventionally, most amino acid substitutions at important protein positions are expected to abolish function. However, in several soluble-globular proteins, we identified a class of non-conserved positions for which various substitutions produced progressive functional changes; we consider these evolutionary “rheostats”. Here, we report a strong rheostat position in the integral membrane protein, Na+/taurocholate cotransporting polypeptide (NTCP), at the site of a pharmacologically-relevant polymorphism (S267F). Functional studies were performed for all 20 substitutions (“S267X”) with three substrates (taurocholate, estrone-3-sulfate and rosuvastatin). The S267X set showed strong rheostatic effects on overall transport, and individual substitutions showed varied effects on transport kinetics (Km and Vmax) and substrate specificity. To assess protein stability, we measured surface expression and used the Rosetta software suite to model structure and stability changes of S267X. Although buried near the substrate binding site, S267X substitutions were easily accommodated in the NTCP structure model. Across the modest range of changes, calculated stabilities correlated with surface-expression differences, but neither parameter correlated with altered transport. Thus, substitutions at rheostat position 267 had wide-ranging effects on the phenotype of this integral membrane protein. We further propose that polymorphic positions in other proteins might be locations of rheostat positions.

Traditional hit-to-lead optimization assumes that upon elaboration of chemical structure, the ligand retains its binding mode relative to the receptor. Here, we build a large-scale collection of related ligand pairs solved in complex with the same protein partner: we find that for 41 of 297 pairs (14%), the binding mode changes upon elaboration of the smaller ligand. While certain ligand physiochemical properties predispose changes in binding mode, particularly those properties that define fragments, simple structure-based modeling proves far more effective for identifying substitutions that alter the binding mode. Some ligand pairs change binding mode because the added substituent would irreconcilably conflict with the receptor in the original pose, whereas others change because the added substituent enables new, stronger interactions that are available only in a different pose. Scaffolds that can engage their target using alternate poses may enable productive structure-based optimization along multiple divergent pathways.

AbstractConventionally, most amino acid substitutions at important protein positions are expected to abolish function. However, in several soluble-globular proteins, we identified a class of non-conserved positions for which various substitutions produced progressive functional changes; we consider these evolutionary “rheostats”. Here, we report a strong rheostat position in the integral membrane protein, Na+/taurocholate cotransporting polypeptide (NTCP), at the site of a pharmacologically-relevant polymorphism (S267F). Functional studies were performed for all 20 substitutions (“S267X”) with three substrates (taurocholate, estrone-3-sulfate and rosuvastatin). The S267X set showed strong rheostatic effects on overall transport, and individual substitutions showed varied effects on transport kinetics (Km and Vmax). However, the outcomes were substrate dependent, indicating altered specificity. To assess protein stability, we measured surface expression and used the Rosetta software suite to model structure and stability changes of S267X. Although buried near the substrate binding site, S267X substitutions were easily accommodated in the NTCP structure model. Across the modest range of changes, calculated stabilities correlated with surface-expression differences, but neither parameter correlated with altered transport. Thus, substitutions at rheostat position 267 had wide-ranging effects on the phenotype of this integral membrane protein. We further propose that polymorphic positions in other proteins might be locations of rheostat positions.