AI is a viable alternative to high throughput screening: a 318-target study

High throughput screening (HTS) is routinely used to identify bioactive small molecules. This requires physical compounds, which limits coverage of accessible chemical space. Computational approaches combined with vast on-demand chemical libraries can access far greater chemical space, provided that the predictive accuracy is sufficient to identify useful molecules. Through the largest and most diverse virtual HTS campaign reported to date, comprising 318 individual projects, we demonstrate that our AtomNet® convolutional neural network successfully finds novel hits across every major therapeutic area and protein class. We address historical limitations of computational screening by demonstrating success for target proteins without known binders, high-quality X-ray crystal structures, or manual cherry-picking of compounds. We show that the molecules selected by the AtomNet® model are novel drug-like scaffolds rather than minor modifications to known bioactive compounds. Our empirical results suggest that computational methods can substantially replace HTS as the first step of small-molecule drug discovery.

Structural polymorphisms within a common powdery mildew effector scaffold as a driver of coevolution with cereal immune receptors

In plants, host-pathogen coevolution often manifests in reciprocal, adaptive genetic changes through variations in host nucleotide-binding leucine-rich repeat immune receptors (NLRs) and virulence-promoting pathogen effectors. In grass powdery mildew (PM) fungi, an extreme expansion of a RNase-like effector family, termed RALPH, dominates the effector repertoire, with some members recognized as avirulence (AVR) effectors by cereal NLR receptors. We report the structures of the sequence-unrelated barley PM effectors AVRA6, AVRA7, and allelic AVRA10/AVRA22 variants, which are detected by highly sequence-related barley NLRs MLA6, MLA7, MLA10, and MLA22 and of wheat PM AVRPM2 detected by the unrelated wheat NLR PM2. The AVR effectors adopt a common scaffold, which is shared with the RNase T1/F1 family. We found striking variations in the number, position, and length of individual structural elements between RALPH AVRs, which is associated with a differentiation of RALPH effector subfamilies. We show that all RALPH AVRs tested have lost nuclease and synthetase activities of the RNase T1/F1 family and lack significant binding to RNA, implying that their virulence activities are associated with neo-functionalization events. Structure-guided mutagenesis identified six AVRA6 residues that are sufficient to turn a sequence-diverged member of the same RALPH subfamily into an effector specifically detected by MLA6. Similar structure-guided information for AVRA10 and AVRA22 indicates that MLA receptors detect largely distinct effector surface patches. Thus, coupling of sequence and structural polymorphisms within the RALPH scaffold of PMs facilitated escape from NLR recognition and potential acquisition of diverse virulence functions.

Targeting of bone morphogenetic protein complexes to heparin/heparan sulfate glycosaminoglycans in bioactive conformation

Bone morphogenetic proteins (BMP) are powerful regulators of cellular processes such as proliferation, differentiation, and apoptosis. However, the specific molecular requirements controlling the bioavailability of BMPs in the extracellular matrix (ECM) are not yet fully understood. Our previous work showed that BMPs are targeted to the ECM as growth factor-prodomain (GF-PD) complexes (CPLXs) via specific interactions of their PDs. We showed that BMP-7 PD binding to the extracellular microfibril component fibrillin-1 renders the CPLXs from an open, bioactive V-shape into a closed, latent ring shape. Here, we show that specific PD interactions with heparin/heparan sulfate glycosaminoglycans (GAGs) allow to target and spatially concentrate BMP-7 and BMP-9 CPLXs in bioactive V-shape conformation. However, targeting to GAGs may be BMP specific, since BMP-10 GF and CPLX do not interact with heparin. Bioactivity assays on solid phase in combination with interaction studies showed that the BMP-7 PD protects the BMP-7 GF from inactivation by heparin. By using transmission electron microscopy, molecular docking, and site-directed mutagenesis, we determined the BMP-7 PD-binding site for heparin. Further, fine-mapping of the fibrillin-1-binding site within the BMP-7 PD and molecular modeling showed that both binding sites are mutually exclusive in the open V- versus closed ring-shape conformation. Together, our data suggest that targeting exquisite BMP PD-binding sites by extracellular protein and GAG scaffolds integrates BMP GF bioavailability in a contextual manner in development, postnatal life, and connective tissue disease.

Collagen's primary structure determines collagen:HSP47 complex stoichiometry

Collagens play important roles in development and homeostasis in most higher organisms. In order to function, collagens require the specific chaperone HSP47 for proper folding and secretion. HSP47 is known to bind to the collagen triple helix, but the exact positions and numbers of binding sites are not clear. Here, we employed a collagen II peptide library to characterize high-affinity binding sites for HSP47. We show that many previously predicted binding sites have very low affinities due to the presence of a negatively charged amino acid in the binding motif. In contrast, large hydrophobic amino acids such as phenylalanine at certain positions in the collagen sequence increase binding strength. For further characterization, we determined two crystal structures of HSP47 bound to peptides containing phenylalanine or leucine. These structures deviate significantly from previously published ones in which different collagen sequences were used. They reveal local conformational rearrangements of HSP47 at the binding site to accommodate the large hydrophobic side chain from the middle strand of the collagen triple helix and, most surprisingly, possess an altered binding stoichiometry in the form of a 1:1 complex. This altered stoichiometry is explained by steric collisions with the second HSP47 molecule present in all structures determined thus far caused by the newly introduced large hydrophobic residue placed on the trailing strand. This exemplifies the importance of considering all three sites of homotrimeric collagen as independent interaction surfaces and may provide insight into the formation of higher oligomeric complexes at promiscuous collagen-binding sites.

Structure of a collagen VI α3 chain VWA domain array: adaptability and functional implications of myopathy causing mutations

Collagen VI is a ubiquitous heterotrimeric protein of the extracellular matrix (ECM) that plays an essential role in the proper maintenance of skeletal muscle. Mutations in collagen VI lead to a spectrum of congenital myopathies, from the mild Bethlem myopathy to the severe Ullrich congenital muscular dystrophy. Collagen VI contains only a short triple helix and consists primarily of von Willebrand factor type A (VWA) domains, protein-protein interaction modules found in a range of ECM proteins. Disease-causing mutations occur commonly in the VWA domains, and the second VWA domain of the α3 chain, the N2 domain, harbors several such mutations. Here, we investigate structure-function relationships of the N2 mutations to shed light on their possible myopathy mechanisms. We determined the X-ray crystal structure of N2, combined with monitoring secretion efficiency in cell culture of selected N2 single-domain mutants, finding that mutations located within the central core of the domain severely affect secretion efficiency. In longer α3 chain constructs, spanning N6-N3, small-angle X-ray scattering demonstrates that the tandem VWA array has a modular architecture and samples multiple conformations in solution. Single-particle EM confirmed the presence of multiple conformations. Structural adaptability appears intrinsic to the VWA domain region of collagen VI α3 and has implications for binding interactions and modulating stiffness within the ECM.

Triple-Helix-Stabilizing Effects in Collagen Model Peptides Containing PPII-Helix-Preorganized Diproline Modules

Collagen model peptides (CMPs) serve as tools for understanding stability and function of the collagen triple helix and have a potential for biomedical applications. In the past, interstrand cross-linking or conformational preconditioning of proline units through stereoelectronic effects have been utilized in the design of stabilized CMPs. To further study the effects determining collagen triple helix stability we investigated a series of CMPs containing synthetic diproline-mimicking modules (ProMs), which were preorganized in a PPII-helix-type conformation by a functionalizable intrastrand C2 bridge. Results of CD-based denaturation studies were correlated with calculated (DFT) conformational preferences of the ProM units, revealing that the relative helix stability is mainly governed by an interplay of main-chain preorganization, ring-flip preference, adaptability, and steric effects. Triple helix integrity was proven by crystal structure analysis and binding to HSP47.

Structural and biophysical characterization of the type VII collagen vWFA2 subdomain leads to identification of two binding sites

Type VII collagen is an extracellular matrix protein, which is important for skin stability; however, detailed information at the molecular level is scarce. The second vWFA (von Willebrand factor type A) domain of type VII collagen mediates important interactions, and immunization of mice induces skin blistering in certain strains. To understand vWFA2 function and the pathophysiological mechanisms leading to skin blistering, we structurally characterized this domain by X-ray crystallography and NMR spectroscopy. Cell adhesion assays identified two new interactions: one with β1 integrin via its RGD motif and one with laminin-332. The latter interaction was confirmed by surface plasmon resonance with a K_D of about 1 mm. These data show that vWFA2 has additional functions in the extracellular matrix besides interacting with type I collagen.

COMP and TSP-4 interact specifically with the novel GXKGHR motif only found in fibrillar collagens

COMP (cartilage oligomeric matrix protein) is a member of the thrombospondin family and forms homopentamers as well as mixed heterooligomers with its closely related family member TSP-4. COMP is long known to bind to collagens and to influence collagen fibril formation. Recent work indicates that already intracellular interaction with collagen is important for collagen secretion. However, the exact binding site of COMP on the collagen triple helix has not been described up to now. In this study we have identified a GXKGHR motif on the collagen II helix to bind to COMP, using a recombinantly expressed collagen II peptide library. This binding sequence is conserved throughout evolution and we demonstrate that TSP-4 binds to the same sequence. The identified binding motif overlaps with the recognition sites of many other collagen-binding partners (e.g. PEDF, Heparin) and also spans the lysine residues, which form collagen cross-links. COMP might thereby protect collagen helices from premature modification and cross-linking. Interestingly, this motif is only found in classical fibrillar collagens, although COMP is known to also bind other types. This might indicate that COMP has a unique interface for fibrillar collagens, thus making it an interesting target for the development of antifibrotic drugs.

The pH-dependent Client Release from the Collagen-specific Chaperone HSP47 Is Triggered by a Tandem Histidine Pair

Heat shock protein 47 (HSP47) is an endoplasmic reticulum (ER)-resident collagen-specific chaperone and essential for proper formation of the characteristic collagen triple helix. It preferentially binds to the folded conformation of its clients and accompanies them from the ER to the Golgi compartment, where it releases them and is recycled back to the ER. Unlike other chaperones, the binding and release cycles are not governed by nucleotide exchange and hydrolysis, but presumably the dissociation of the HSP47-procollagen complex is triggered by the lower pH in the Golgi (pH 6.3) compared with the ER (pH 7.4). Histidine residues have been suggested as triggers due to their approximate textbook pKa value of 6.1 for their side chains. We present here an extensive theoretical and experimental study of the 14 histidine residues present in canine HSP47, where we have mutated all histidine residues in the collagen binding interface and additionally all of those that were predicted to undergo a significant change in protonation state between pH 7 and 6. These mutants were characterized by biolayer interferometry for their pH-dependent binding to a collagen model. One mutant (H238N) loses binding, which can be explained by a rearrangement of the Arg(222) and Asp(385) residues, which are crucial for specific collagen recognition. Most of the other mutants were remarkably silent, but a double mutant with His(273) and His(274) exchanged for asparagines exhibits a much less pronounced pH dependence of collagen binding. This effect is mainly caused by a lower koff at the low pH values.

AMACO is a component of the basement membrane-associated Fraser complex

Fraser syndrome (FS) is a phenotypically variable, autosomal recessive disorder characterized by cryptophthalmus, cutaneous syndactyly, and other malformations resulting from mutations in FRAS1, FREM2, and GRIP1. Transient embryonic epidermal blistering causes the characteristic defects of the disorder. Fras1, Frem1, and Frem2 form the extracellular Fraser complex, which is believed to stabilize the basement membrane. However, several cases of FS could not be attributed to mutations in FRAS1, FREM2, or GRIP1, and FS displays high clinical variability, suggesting that there is an additional genetic, possibly modifying contribution to this disorder. An extracellular matrix protein containing VWA-like domains related to those in matrilins and collagens (AMACO), encoded by the VWA2 gene, has a very similar tissue distribution to the Fraser complex proteins in both mouse and zebrafish. Here, we show that AMACO deposition is lost in Fras1-deficient zebrafish and mice and that Fras1 and AMACO interact directly via their chondroitin sulfate proteoglycan (CSPG) and P2 domains. Knockdown of vwa2, which alone causes no phenotype, enhances the phenotype of hypomorphic Fras1 mutant zebrafish. Together, our data suggest that AMACO represents a member of the Fraser complex.

Genetic analysis of fin development in zebrafish identifies furin and hemicentin1 as potential novel fraser syndrome disease genes

Using forward genetics, we have identified the genes mutated in two classes of zebrafish fin mutants. The mutants of the first class are characterized by defects in embryonic fin morphogenesis, which are due to mutations in a Laminin subunit or an Integrin alpha receptor, respectively. The mutants of the second class display characteristic blistering underneath the basement membrane of the fin epidermis. Three of them are due to mutations in zebrafish orthologues of FRAS1, FREM1, or FREM2, large basement membrane protein encoding genes that are mutated in mouse bleb mutants and in human patients suffering from Fraser Syndrome, a rare congenital condition characterized by syndactyly and cryptophthalmos. Fin blistering in a fourth group of zebrafish mutants is caused by mutations in Hemicentin1 (Hmcn1), another large extracellular matrix protein the function of which in vertebrates was hitherto unknown. Our mutant and dose-dependent interaction data suggest a potential involvement of Hmcn1 in Fraser complex-dependent basement membrane anchorage. Furthermore, we present biochemical and genetic data suggesting a role for the proprotein convertase FurinA in zebrafish fin development and cell surface shedding of Fras1 and Frem2, thereby allowing proper localization of the proteins within the basement membrane of forming fins. Finally, we identify the extracellular matrix protein Fibrillin2 as an indispensable interaction partner of Hmcn1. Thus we have defined a series of zebrafish mutants modelling Fraser Syndrome and have identified several implicated novel genes that might help to further elucidate the mechanisms of basement membrane anchorage and of the disease's aetiology. In addition, the novel genes might prove helpful to unravel the molecular nature of thus far unresolved cases of the human disease.

There are a large number of human genetic syndromes with limb and digit deformities. It has been shown that the genes underlying these syndromes are well conserved in evolution, and most perform the same role even in the fins of fish. One such human syndrome is Fraser Syndrome, characterized by a number of defects including fusion of the fingers (syndactyly). Data obtained with corresponding mouse mutants suggest that all of these defects are due to transient basement membrane disruptions and epithelial blistering during development. Whilst some of the Fraser Syndrome genes have been identified, others are unknown. We show that mutation of the known Fraser Syndrome genes in zebrafish generate comparable blistering defects in the fins. Importantly, we have also identified additional genes and mechanisms required for the same processes. Included in this are hemicentin1, a gene whose function had thus far only been studied in nematodes, and furinA, encoding a proprotein convertase, for which we reveal a novel role in ectodomain shedding of Fras/Frem proteins. This work thus expands our understanding, not only of Fraser Syndrome, but also of the common processes of basement membrane formation and function during fin and limb development.

Expression of the AMACO (VWA2 protein) ortholog in zebrafish

AMACO is a basement membrane associated protein that belongs to the VWA domain-containing protein superfamily. In addition to three VWA domains it contains two EGF-like domains, a cysteine-rich domain and a unique domain. Mouse AMACO has been partially characterized, but its function remains unknown. The zebrafish genome contains a single AMACO ortholog gene on chromosome 12. The domain structure is completely conserved between zebrafish and mouse and the first EGF-like domain, carrying a rare O-glucosylation and O-fucosylation consensus sequence, has the highest identity at the protein level. RT-PCR shows strongest AMACO expression during development, starting at the 5 somite stage. An antibody specific for zebrafish AMACO detected expression mainly in myosepta but also in skin, pronephros, pituitary gland, otic capsule and gills. In situ hybridization revealed that the muscle precursor cells of the somites express the protein that is laid down in the myosepta.

O-glucosylation and O-fucosylation occur together in close proximity on the first epidermal growth factor repeat of AMACO (VWA2 protein)

AMACO (VWA2 protein) is an extracellular matrix protein of unknown function associated with certain basement membranes in skin, lung, and kidney. AMACO is a member of the von Willebrand factor A-like (VWA) domain containing protein superfamily and in addition to three VWA domains it also contains two epidermal growth factor-like domains. One of these contains the rare, overlapping consensus sequences for both O-glucosylation and O-fucosylation. In earlier studies of other proteins the attachment of either core glucose and fucose moieties or of the respective elongated glycans starting with these monosaccharides has been described. By a detailed mass spectrometric analysis we show that both elongated O-glucosylated (Xyl1-3Xyl1-3Glc) and elongated O-fucosylated glycan chains (NeuAc2-3Gal1-4GlcNAc1-3Fuc) can be attached to AMACO in close proximity on the same epidermal growth factor-like domain. It has been reported that the lack of O-fucosylation can markedly decrease secretion of proteins. However, the secretion of AMACO is not significantly affected when the glycosylation sites are mutated. The number of extracellular matrix proteins carrying the overlapping consensus sequence is very limited and it could be that these modifications have a new, yet unknown function.

Matrilins mediate weak cell attachment without promoting focal adhesion formation

The matrilins form a family of non-collagenous adaptor proteins in the extracellular matrix. The extracellular ligand interactions of matrilins have been studied in some detail, while the potential interplay between matrilins and cells has been largely neglected. Except for matrilin-4, all matrilins mediate cell attachment, but only for matrilin-1 and -3 the binding is clearly dose dependent and seen already at moderate coating concentrations. Even so, much higher concentrations of matrilin-1 or -3 than of fibronectin are required for cell attachment to reach plateau values. Integrins contribute to the matrilin-mediated cell attachment, but the binding does not lead to formation of focal contacts and reorganisation of the actin cytoskeleton. Cells deficient in beta1 integrins are able to adhere, although weaker, and matrilins do not bind the soluble integrin alpha1beta1 and alpha2beta1 ectodomains. Cell surface proteoglycans may promote the attachment, as cells deficient in glycosaminoglycan biosynthesis adhere less well to matrilin-3. Even so, exogenous glycosaminoglycans are not able to compete for the attachment of HaCaT cells to matrilins.