Enumerative Discovery of Noncanonical Polypeptide Secondary Structures.

Energetically favorable local interactions can overcome the entropic cost of chain ordering and cause otherwise flexible polymers to adopt regularly repeating backbone conformations. A prominent example is the α helix present in many protein structures, which is stabilized by i, i + 4 hydrogen bonds between backbone peptide units. With the increased chemical diversity offered by unnatural amino acids and backbones, it has been possible to identify regularly repeating structures not present in proteins, but to date, there has been no systematic approach for identifying new polymers likely to have such structures despite their considerable potential for molecular engineering. Here we describe a systematic approach to search through dipeptide combinations of 130 chemically diverse amino acids to identify those predicted to populate unique low-energy states. We characterize ten newly identified dipeptide repeating structures using circular dichroism spectroscopy and comparison with calculated spectra. NMR and X-ray crystallographic structures of two of these dipeptide-repeat polymers are similar to the computational models. Our approach is readily generalizable to identify low-energy repeating structures for a wide variety of polymers, and our ordered dipeptide repeats provide new building blocks for molecular engineering.

Design of amyloidogenic peptide traps.

Segments of proteins with high ß-strand propensity can self-associate to form amyloid fibrils implicated in many diseases. We describe a general approach to bind such segments in ß-strand and ß-hairpin conformations using de novo designed scaffolds that contain deep peptide-binding clefts. The designs bind their cognate peptides in vitro with nanomolar affinities. The crystal structure of a designed protein-peptide complex is close to the design model, and NMR characterization reveals how the peptide-binding cleft is protected in the apo state. We use the approach to design binders to the amyloid-forming proteins transthyretin, tau, serum amyloid A1 and amyloid ß1-42 (Aß42). The Aß binders block the assembly of Aß fibrils as effectively as the most potent of the clinically tested antibodies to date and protect cells from toxic Aß42 species.

Improving Protein Expression, Stability, and Function with ProteinMPNN.

Natural proteins are highly optimized for function but are often difficult to produce at a scale suitable for biotechnological applications due to poor expression in heterologous systems, limited solubility, and sensitivity to temperature. Thus, a general method that improves the physical properties of native proteins while maintaining function could have wide utility for protein-based technologies. Here, we show that the deep neural network ProteinMPNN, together with evolutionary and structural information, provides a route to increasing protein expression, stability, and function. For both myoglobin and tobacco etch virus (TEV) protease, we generated designs with improved expression, elevated melting temperatures, and improved function. For TEV protease, we identified multiple designs with improved catalytic activity as compared to the parent sequence and previously reported TEV variants. Our approach should be broadly useful for improving the expression, stability, and function of biotechnologically important proteins.

De novo design of high-affinity binders of bioactive helical peptides.

Many peptide hormones form an α-helix on binding their receptors1-4, and sensitive methods for their detection could contribute to better clinical management of disease5. De novo protein design can now generate binders with high affinity and specificity to structured proteins6,7. However, the design of interactions between proteins and short peptides with helical propensity is an unmet challenge. Here we describe parametric generation and deep learning-based methods for designing proteins to address this challenge. We show that by extending RFdiffusion8 to enable binder design to flexible targets, and to refining input structure models by successive noising and denoising (partial diffusion), picomolar-affinity binders can be generated to helical peptide targets by either refining designs generated with other methods, or completely de novo starting from random noise distributions without any subsequent experimental optimization. The RFdiffusion designs enable the enrichment and subsequent detection of parathyroid hormone and glucagon by mass spectrometry, and the construction of bioluminescence-based protein biosensors. The ability to design binders to conformationally variable targets, and to optimize by partial diffusion both natural and designed proteins, should be broadly useful.

De novo design of highly selective miniprotein inhibitors of integrins αvß6 and αvß8.

The RGD (Arg-Gly-Asp)-binding integrins αvß6 and αvß8 are clinically validated cancer and fibrosis targets of considerable therapeutic importance. Compounds that can discriminate between homologous αvß6 and αvß8 and other RGD integrins, stabilize specific conformational states, and have high thermal stability could have considerable therapeutic utility. Existing small molecule and antibody inhibitors do not have all these properties, and hence new approaches are needed. Here we describe a generalized method for computationally designing RGD-containing miniproteins selective for a single RGD integrin heterodimer and conformational state. We design hyperstable, selective αvß6 and αvß8 inhibitors that bind with picomolar affinity. CryoEM structures of the designed inhibitor-integrin complexes are very close to the computational design models, and show that the inhibitors stabilize specific conformational states of the αvß6 and the αvß8 integrins. In a lung fibrosis mouse model, the αvß6 inhibitor potently reduced fibrotic burden and improved overall lung mechanics, demonstrating the therapeutic potential of de novo designed integrin binding proteins with high selectivity.

De novo design of highly selective miniprotein inhibitors of integrins αvß6 and αvß8.

The RGD (Arg-Gly-Asp)-binding integrins αvß6 and αvß8 are clinically validated cancer and fibrosis targets of considerable therapeutic importance. Compounds that can discriminate between the two closely related integrin proteins and other RGD integrins, stabilize specific conformational states, and have sufficient stability enabling tissue restricted administration could have considerable therapeutic utility. Existing small molecules and antibody inhibitors do not have all of these properties, and hence there is a need for new approaches. Here we describe a method for computationally designing hyperstable RGD-containing miniproteins that are highly selective for a single RGD integrin heterodimer and conformational state, and use this strategy to design inhibitors of αvß6 and αvß8 with high selectivity. The αvß6 and αvß8 inhibitors have picomolar affinities for their targets, and >1000-fold selectivity over other RGD integrins. CryoEM structures are within 0.6-0.7Å root-mean-square deviation (RMSD) to the computational design models; the designed αvß6 inhibitor and native ligand stabilize the open conformation in contrast to the therapeutic anti-αvß6 antibody BG00011 that stabilizes the bent-closed conformation and caused on-target toxicity in patients with lung fibrosis, and the αvß8 inhibitor maintains the constitutively fixed extended-closed αvß8 conformation. In a mouse model of bleomycin-induced lung fibrosis, the αvß6 inhibitor potently reduced fibrotic burden and improved overall lung mechanics when delivered via oropharyngeal administration mimicking inhalation, demonstrating the therapeutic potential of de novo designed integrin binding proteins with high selectivity.

Computational design of constitutively active cGAS.

Cyclic GMP-AMP synthase (cGAS) is a pattern recognition receptor critical for the innate immune response to intracellular pathogens, DNA damage, tumorigenesis and senescence. Binding to double-stranded DNA (dsDNA) induces conformational changes in cGAS that activate the enzyme to produce 2'-3' cyclic GMP-AMP (cGAMP), a second messenger that initiates a potent interferon (IFN) response through its receptor, STING. Here, we combined two-state computational design with informatics-guided design to create constitutively active, dsDNA ligand-independent cGAS (CA-cGAS). We identified CA-cGAS mutants with IFN-stimulating activity approaching that of dsDNA-stimulated wild-type cGAS. DNA-independent adoption of the active conformation was directly confirmed by X-ray crystallography. In vivo expression of CA-cGAS in tumor cells resulted in STING-dependent tumor regression, demonstrating that the designed proteins have therapeutically relevant biological activity. Our work provides a general framework for stabilizing active conformations of enzymes and provides CA-cGAS variants that could be useful as genetically encoded adjuvants and tools for understanding inflammatory diseases.

Accurate de novo design of membrane-traversing macrocycles.

We use computational design coupled with experimental characterization to systematically investigate the design principles for macrocycle membrane permeability and oral bioavailability. We designed 184 6-12 residue macrocycles with a wide range of predicted structures containing noncanonical backbone modifications and experimentally determined structures of 35; 29 are very close to the computational models. With such control, we show that membrane permeability can be systematically achieved by ensuring all amide (NH) groups are engaged in internal hydrogen bonding interactions. 84 designs over the 6-12 residue size range cross membranes with an apparent permeability greater than 1 × 10-6 cm/s. Designs with exposed NH groups can be made membrane permeable through the design of an alternative isoenergetic fully hydrogen-bonded state favored in the lipid membrane. The ability to robustly design membrane-permeable and orally bioavailable peptides with high structural accuracy should contribute to the next generation of designed macrocycle therapeutics.

The nonreceptor tyrosine kinase SRMS inhibits autophagy and promotes tumor growth by phosphorylating the scaffolding protein FKBP51.

Nutrient-responsive protein kinases control the balance between anabolic growth and catabolic processes such as autophagy. Aberrant regulation of these kinases is a major cause of human disease. We report here that the vertebrate nonreceptor tyrosine kinase Src-related kinase lacking C-terminal regulatory tyrosine and N-terminal myristylation sites (SRMS) inhibits autophagy and promotes growth in a nutrient-responsive manner. Under nutrient-replete conditions, SRMS phosphorylates the PHLPP scaffold FK506-binding protein 51 (FKBP51), disrupts the FKBP51-PHLPP complex, and promotes FKBP51 degradation through the ubiquitin-proteasome pathway. This prevents PHLPP-mediated dephosphorylation of AKT, causing sustained AKT activation that promotes growth and inhibits autophagy. SRMS is amplified and overexpressed in human cancers where it drives unrestrained AKT signaling in a kinase-dependent manner. SRMS kinase inhibition activates autophagy, inhibits cancer growth, and can be accomplished using the FDA-approved tyrosine kinase inhibitor ibrutinib. This illuminates SRMS as a targetable vulnerability in human cancers and as a new target for pharmacological induction of autophagy in vertebrates.

KRASQ61H Preferentially Signals through MAPK in a RAF Dimer-Dependent Manner in Non-Small Cell Lung Cancer.

Assembly of RAS molecules into complexes at the cell membrane is critical for RAS signaling. We previously showed that oncogenic KRAS codon 61 mutations increase its affinity for RAF, raising the possibility that KRASQ61H, the most common KRAS mutation at codon 61, upregulates RAS signaling through mechanisms at the level of RAS assemblies. We show here that KRASQ61H exhibits preferential binding to RAF relative to PI3K in cells, leading to enhanced MAPK signaling in in vitro models and human NSCLC tumors. X-ray crystallography of KRASQ61H:GTP revealed that a hyperdynamic switch 2 allows for a more stable interaction with switch 1, suggesting that enhanced RAF activity arises from a combination of absent intrinsic GTP hydrolysis activity and increased affinity for RAF. Disruption of KRASQ61H assemblies by the RAS oligomer-disrupting D154Q mutation impaired RAF dimerization and altered MAPK signaling but had little effect on PI3K signaling. However, KRASQ61H oligomers but not KRASG12D oligomers were disrupted by RAF mutations that disrupt RAF-RAF interactions. KRASQ61H cells show enhanced sensitivity to RAF and MEK inhibitors individually, whereas combined treatment elicited synergistic growth inhibition. Furthermore, KRASQ61H tumors in mice exhibited high vulnerability to MEK inhibitor, consistent with cooperativity between KRASQ61H and RAF oligomerization and dependence on MAPK signaling. These findings support the notion that KRASQ61H and functionally similar mutations may serve as predictive biomarkers for targeted therapies against the MAPK pathway. SIGNIFICANCE: These findings show that oncogenic KRASQ61H forms a cooperative RAS-RAF ternary complex, which renders RAS-driven tumors vulnerable to MEKi and RAFi, thus establishing a framework for evaluating RAS biomarker-driven targeted therapies.

Structure and Characterization of a Covalent Inhibitor of Src Kinase.

Unregulated Src activity promotes malignant processes in cancer, but no Src-directed targeted therapies are used clinically, possibly because early Src inhibitors produce off-target effects leading to toxicity. Improved selective Src inhibitors may enable Src-directed therapies. Previously, we reported an irreversible Src inhibitor, DGY-06-116, based on the hybridization of dasatinib and a promiscuous covalent kinase probe SM1-71. Here, we report biochemical and biophysical characterization of this compound. An x-ray co-crystal structure of DGY-06-116: Src shows a covalent interaction with the kinase p-loop and occupancy of the back hydrophobic kinase pocket, explaining its high potency, and selectivity. However, a reversible analog also shows similar potency. Kinetic analysis shows a slow inactivation rate compared to other clinically approved covalent kinase inhibitors, consistent with a need for p-loop movement prior to covalent bond formation. Overall, these results suggest that a strong reversible interaction is required to allow sufficient time for the covalent reaction to occur. Further optimization of the covalent linker may improve the kinetics of covalent bond formation.

GTP hydrolysis is modulated by Arg34 in the RASopathy-associated KRASP34R.

RAS proteins are commonly mutated in cancerous tumors, but germline RAS mutations are also found in RASopathy syndromes such as Noonan syndrome (NS) and cardiofaciocutaneous (CFC) syndrome. Activating RAS mutations can be subclassified based on their activating mechanisms. Understanding the structural basis for these mechanisms may provide clues for how to manage associated health conditions. We determined high-resolution X-ray structures of the RASopathy mutant KRASP34R seen in NS and CFCS. GTP and GDP-bound KRASP34R crystallized in multiple forms, with each lattice consisting of multiple protein conformations. In all GTP-bound conformations, the switch regions are not compatible with GAP binding, suggesting a structural mechanism for the GAP insensitivity of this RAS mutant. However, GTP-bound conformations are compatible with intrinsic nucleotide hydrolysis, including one that places R34 in a position analogous to the GAP arginine finger or intrinsic arginine finger found in heterotrimeric G proteins, which may support intrinsic GTP hydrolysis. We also note that the affinity between KRASP34R and RAF-RBD is decreased, suggesting another possible mechanism for dampening of RAS signaling. These results may provide a foothold for development of new mutation-specific strategies to address KRASP34R -driven diseases.

Structural basis of the atypical activation mechanism of KRASV14I.

RAS regulation and signaling are largely accomplished by direct protein-protein interactions, making RAS protein dynamics a critical determinant of RAS function. Here, we report a crystal structure of GDP-bound KRASV14I, a mutated KRAS variant associated with the developmental RASopathy disorder Noonan syndrome (NS), at 1.5-1.6 Å resolution. The structure is notable for revealing a marked extension of switch 1 away from the G-domain and nucleotide-binding site of the KRAS protein. We found that this extension is associated with a loss of the magnesium ion and a tilt in the position of the guanine base because of the additional carbon introduced by the isoleucine substitution. Hydrogen-deuterium exchange MS analysis confirmed that this conformation occurs in solution, but also disclosed a difference in kinetics when compared with KRASA146T, another RAS mutant that displays a nearly identical conformation in previously reported crystal structures. This conformational change contributed to a high rate of guanine nucleotide-exchange factor (GEF)-dependent and -independent nucleotide exchange and to an increase in affinity for SOS Ras/Rac GEF 1 (SOS1), which appears to be the major mode of activation for this RAS variant. These results highlight a mechanistic connection between KRASA146T and KRASV14I that may have implications for the regulation of these variants and for the development of therapeutic strategies to manage KRAS variant-associated disorders.

Tissue-Specific Oncogenic Activity of KRASA146T.

KRAS is the most frequently mutated oncogene. The incidence of specific KRAS alleles varies between cancers from different sites, but it is unclear whether allelic selection results from biological selection for specific mutant KRAS proteins. We used a cross-disciplinary approach to compare KRASG12D, a common mutant form, and KRASA146T, a mutant that occurs only in selected cancers. Biochemical and structural studies demonstrated that KRASA146T exhibits a marked extension of switch 1 away from the protein body and nucleotide binding site, which activates KRAS by promoting a high rate of intrinsic and guanine nucleotide exchange factor-induced nucleotide exchange. Using mice genetically engineered to express either allele, we found that KRASG12D and KRASA146T exhibit distinct tissue-specific effects on homeostasis that mirror mutational frequencies in human cancers. These tissue-specific phenotypes result from allele-specific signaling properties, demonstrating that context-dependent variations in signaling downstream of different KRAS mutants drive the KRAS mutational pattern seen in cancer. SIGNIFICANCE: Although epidemiologic and clinical studies have suggested allele-specific behaviors for KRAS, experimental evidence for allele-specific biological properties is limited. We combined structural biology, mass spectrometry, and mouse modeling to demonstrate that the selection for specific KRAS mutants in human cancers from different tissues is due to their distinct signaling properties.See related commentary by Hobbs and Der, p. 696.This article is highlighted in the In This Issue feature, p. 681.

KRAS Switch Mutants D33E and A59G Crystallize in the State 1 Conformation.

KRAS switch loop movements play a crucial role in regulating RAS signaling, and alteration of these sensitive dynamics is a principal mechanism through which disease-associated RAS mutations lead to aberrant RAS activation. Prior studies suggest that despite a high degree of sequence similarity, the switches in KRAS are more dynamic than those in HRAS. We determined X-ray crystal structures of the rare tumorigenic KRAS mutants KRASD33E, in switch 1 (SW1), and KRASA59G, in switch 2 (SW2), bound to GDP and found these adopt nearly identical, open SW1 conformations as well as altered SW2 conformations. KRASA59G bound to a GTP analogue crystallizes in the same conformation. This open conformation is consistent with the inactive "state 1" previously observed for HRAS bound to GTP. For KRASA59G, switch rearrangements may be regulated by increased flexibility in the 57DXXGQ61 motif at codon 59. However, loss of interactions between side chains at codons 33 and 35 in the SW1 33DPT35 motif drives changes for KRASD33E. The 33DPT35 motif is conserved for multiple members of the RAS subfamily but is not found in RAB, RHO, ARF, or Gα families, suggesting that dynamics mediated by this motif may be important for determining the selectivity of RAS-effector interactions. Biochemically, the consequence of altered switch dynamics is the same, showing impaired interaction with the guanine exchange factor SOS and loss of GAP-dependent GTPase activity. However, interactions with the RBD of RAF are preserved. Overall, these observations add to a body of evidence suggesting that HRAS and KRAS show meaningful differences in functionality stemming from differential protein dynamics independent of the hypervariable region.

Takinib, a Selective TAK1 Inhibitor, Broadens the Therapeutic Efficacy of TNF-α Inhibition for Cancer and Autoimmune Disease.

Tumor necrosis factor alpha (TNF-α) has both positive and negative roles in human disease. In certain cancers, TNF-α is infused locally to promote tumor regression, but dose-limiting inflammatory effects limit broader utility. In autoimmune disease, anti-TNF-α antibodies control inflammation in most patients, but these benefits are offset during chronic treatment. TAK1 acts as a key mediator between survival and cell death in TNF-α-mediated signaling. Here, we describe Takinib, a potent and selective TAK1 inhibitor that induces apoptosis following TNF-α stimulation in cell models of rheumatoid arthritis and metastatic breast cancer. We demonstrate that Takinib is an inhibitor of autophosphorylated and non-phosphorylated TAK1 that binds within the ATP-binding pocket and inhibits by slowing down the rate-limiting step of TAK1 activation. Overall, Takinib is an attractive starting point for the development of inhibitors that sensitize cells to TNF-α-induced cell death, with general implications for cancer and autoimmune disease treatment.

In Silico Identification of Bioremediation Potential: Carbamazepine and Other Recalcitrant Personal Care Products.

Emerging contaminants are principally personal care products not readily removed by conventional wastewater treatment and, with an increasing reliance on water recycling, become disseminated in drinking water supplies. Carbamazepine, a widely used neuroactive pharmaceutical, increasingly escapes wastewater treatment and is found in potable water. In this study, a mechanism is proposed by which carbamazepine resists biodegradation, and a previously unknown microbial biodegradation was predicted computationally. The prediction identified biphenyl dioxygenase from Paraburkholderia xenovorans LB400 as the best candidate enzyme for metabolizing carbamazepine. The rate of degradation described here is 40 times greater than the best reported rates. The metabolites cis-10,11-dihydroxy-10,11-dihydrocarbamazepine and cis-2,3-dihydroxy-2,3-dihydrocarbamazepine were demonstrated with the native organism and a recombinant host. The metabolites are considered nonharmful and mitigate the generation of carcinogenic acridine products known to form when advanced oxidation methods are used in water treatment. Other recalcitrant personal care products were subjected to prediction by the Pathway Prediction System and tested experimentally with P. xenovorans LB400. It was shown to biodegrade structurally diverse compounds. Predictions indicated hydrolase or oxygenase enzymes catalyzed the initial reactions. This study highlights the potential for using the growing body of enzyme-structural and genomic information with computational methods to rapidly identify enzymes and microorganisms that biodegrade emerging contaminants.

Insights from Analysis of Binding Sites of Human Meprins: Screening of Inhibitors by Molecular Dynamics Simulation Study.

Human meprin-α and-ß are important regulators of angiogenesis, cancer, inflammation, fibrosis, and neurodegenerative diseases and hence important therapeutic targets. Meprins are the only astacin proteases that are expressed in membrane-bound and secreted form. The cleavage specificity of human meprins is similar in certain cases but differs markedly in others. The inhibitor selectivity of human meprins is controlled by the specific residues involved in binding at the active-site cleft of the proteases. Meprins are inhibited by various small molecular inhibitors as well as macromolecular endogenous inhibitors, making them good drug targets. In the current study, molecular dynamics simulation was performed for 10 ns on ten systems consisting of two apoenzymes of meprin -α/ß and eight complexes of human meprin-α and -ß complexed to four inhibitors with different metal binding moieties and comparable Ki values. These simulation studies helped to elucidate the molecular details of how several parameters influence protein-inhibitor binding affinity. Analysis of the interaction energies of the protein-inhibitor complexes revealed the diverse binding nature of this series of inhibitors. Several structural segments of human meprins exhibited certain conformational changes during the simulation time course. Among the inhibitors studied captopril had a different disposition in the meprin-bound complexes compared to the other three inhibitors, namely Pro- Leu-Gly-hydroxamate, galardin and EDTA. Comparison of the interaction energies for each system helped us to conclude that the hydroxamic acid-based inhibitors are the most potent inhibitors of meprins.

Structure of aminodeoxychorismate synthase from Stenotrophomonas maltophilia.

PabB, aminodeoxychorismate synthase, is the chorismic acid binding component of the heterodimeric PabA-PabB complex that converts chorismic acid to 4-amino-4-deoxychorismate, a precursor of p-aminobenzoate and folic acid in microorganisms. The second component, a glutamine amidotransferase subunit, PabA, generates ammonia that is channeled to the PabB active site where it attacks C4 of a chorismate-derived intermediate that is covalently bound, through C2, to an active site lysine residue. The presence of a PIKGT motif was, until recently, believed to allow discrimination of PabB enzymes from the closely related enzyme anthranilate synthase, which typically contains a PIAGT active site motif and does not form a covalent enzyme-substrate intermediate with chorismate. A subclass of PabB enzymes that employ an alternative mechanism requiring 2 equiv of ammonia from glutamine and that feature a noncovalently bound 2-amino-2-deoxyisochorismate intermediate was recently identified. Here we report the 2.25 Å crystal structure of PabB from the emerging pathogen Stenotrophomonas maltophilia. It is the first reported structure of a PabB that features the PIAGT motif. Surprisingly, no dedicated pabA is evident in the genome of S. maltophilia, suggesting that another cellular amidotransferase is able to fulfill the role of PabA in this organism. Evaluation of the ammonia-dependent aminodeoxychorismate synthase activity of S. maltophilia PabB alone revealed that it is virtually inactive. However, in the presence of a heterologous PabA surrogate, typical levels of activity were observed using either glutamine or ammonia as the nitrogen source. Additionally, the structure suggests that a key segment of the polypeptide can remodel itself to interact with a nonspecialized or shared amidotransferase partner in vivo. The structure and mass spectral analysis further suggest that S. maltophilia PabB, like Escherichia coli PabB, binds tryptophan in a vestigial regulatory site. The observation that the binding site is unoccupied in the crystal structure, however, suggests the affinity may be low relative to that of E. coli PabB.

Conserved water-mediated H-bonding dynamics of catalytic His159 and Asp158: insight into a possible acid-base coupled mechanism in plant thiol protease.

Cysteine protease is ubiquitous in nature. Excess activity of this enzyme causes intercellular proteolysis, muscle tissue degradation, etc. The role of water-mediated interactions in the stabilization of catalytically significant Asp158 and His159 was investigated by performing molecular dynamics simulation studies of 16 three-dimensional structures of plant thiol proteases. In the simulated structures, the hydrophilic W(1), W(2) and WD(1) centers form hydrogen bonds with the OD1 atom of Asp158 and the ND1 atom of His159. In the solvated structures, another water molecule, W(E), forms a hydrogen bond with the NE2 atom of His159. In the absence of the water molecule W(E), Trp177 (NE1) and Gln19 (NE2) directly interact with the NE2 atom of His159. All these hydrophilic centers (the locations of W(1), W(2), WD(1), and W(E)) are conserved, and they play a critical role in the stabilization of His-Asp complexes. In the water dynamics of solvated structures, the water molecules W(1) and W(2) form a water...water hydrogen-bonded network with a few other water molecules. A few dynamical conformations or transition states involving direct (His159 ND1...Asp158 OD1) and water-mediated (His159 ND1...W(2)...Asp158 OD1) hydrogen-bonded complexes are envisaged from these studies.