Binding and Dynamics Demonstrate the Destabilization of Ligand Binding for the S688Y Mutation in the NMDA Receptor GluN1 Subunit

Encephalopathies are brain dysfunctions that lead to cognitive, sensory, and motor development impairments. Recently, the identification of several mutations within the N-methyl-D-aspartate receptor (NMDAR) have been identified as significant in the etiology of this group of conditions. However, a complete understanding of the underlying molecular mechanism and changes to the receptor due to these mutations has been elusive. We studied the molecular mechanisms by which one of the first mutations within the NMDAR GluN1 ligand binding domain, Ser688Tyr, causes encephalopathies. We performed molecular docking, randomly seeded molecular dynamics simulations, and binding free energy calculations to determine the behavior of the two major co-agonists: glycine and D-serine, in both the wild-type and S688Y receptors. We observed that the Ser688Tyr mutation leads to the instability of both ligands within the ligand binding site due to structural changes associated with the mutation. The binding free energy for both ligands was significantly more unfavorable in the mutated receptor. These results explain previously observed in vitro electrophysiological data and provide detailed aspects of ligand association and its effects on receptor activity. Our study provides valuable insight into the consequences of mutations within the NMDAR GluN1 ligand binding domain.

A Case Report of Poor Response to Selpercatinib in the Presence of a 632_633 RET Deletion

Background: Genomic deletions in medullary thyroid cancer (MTC) are rare. Selpercatinib is a highly selective RET inhibitor for treatment of metastatic RET-altered MTC. We report a 35-year-old male with an aggressive metastatic MTC harboring p.632_633del RET that was poorly responsive to RET kinase inhibitor selpercatinib. Objective: Our objective was to understand the clinical phenotype of p.632_633del RET in MTC in the context of novel RET kinase inhibitor treatment. Methods: Wild-type and p.632_633del RET sequences were modeled using a lighter version of the AlphaFold2 (AF2) software. Functional studies were performed on transfected HEK 293 cells (pCMV6-Entry, pCMV6-RET, or pCMV6-RET(p.632_633del) treated with inhibitors for 24 hours and analyzed on luciferase assays. Results: Structural modeling revealed a paucity of disulfide bridge between Cys630-Cys634 in p.632_633del RET sequences, apparent in wild-type, while forming an intermolecular disulfide bridge between two Cys656. Proximity of juxtamembrane segments of each dimer may impede Tyr687 phosphorylation and stable conformation of intracellular RET that hosts selpercatinib. In vitro experiments confirmed a reduction in efficacy of selpercatinib upon p.632_633del RET compared with wild-type RET control. Conclusion: Clinical presentation together with structural modeling and functional studies suggests that p.632_633del RET results in poor response to selpercatinib.

Functional significance of germline EPAS1 variants

Mosaic or somatic EPAS1 mutations are associated with a range of phenotypes including pheochromocytoma and/or paraganglioma (PPGL), polycythemia and somatostatinoma. The pathogenic potential of germline EPAS1 variants however is not well understood. We report a number of germline EPAS1 variants occurring in patients with PPGL, including a novel variant c.739C>A (p.Arg247Ser); a previously described variant c.1121T>A (p.Phe374Tyr); several rare variants, c.581A>G (p.His194Arg), c.2353C>A (p.Pro785Thr) and c.2365A>G (p.Ile789Val); a common variant c.2296A>C (p.Thr766Pro). We performed detailed functional studies to understand their pathogenic role in PPGL. In transient transfection studies, EPAS1/HIF-2α p.Arg247Ser, p.Phe374Tyr and p.Pro785Thr were all stable in normoxia. In co-immunoprecipitation assays, only the novel variant p.Arg247Ser showed diminished interaction with pVHL. A direct interaction between HIF-2α Arg247 and pVHL was confirmed in structural models. Transactivation was assessed by means of a HRE-containing reporter gene in transiently transfected cells, and significantly higher reporter activity was only observed with EPAS1/HIF-2α p.Phe374Tyr and p.Pro785Thr. In conclusion, three germline EPAS1 variants (c.739C>A (p.Arg247Ser), c.1121T>A (p.Phe374Tyr) and c.2353C>A (p.Pro785Thr)) all have some functional features in common with somatic activating mutations. Our findings suggest that these three germline variants are hypermorphic alleles that may act as modifiers to the expression of PPGLs.

Structural and Functional Aspects of Targeting the Secreted Human Group IIA Phospholipase A2

Human group IIA secretory phospholipase A₂ (hGIIA) promotes the proliferation of cancer cells, making it a compelling therapeutic target, but it is also significant in other inflammatory conditions. Consequently, suitable inhibitors of hGIIA have always been sought. The activation of phospholipases A₂ and the catalysis of glycerophospholipid substrates generally leads to the release of fatty acids such as arachidonic acid (AA) and lysophospholipid, which are then converted to mediator compounds, including prostaglandins, leukotrienes, and the platelet-activating factor. However, this ability of hGIIA to provide AA is not a complete explanation of its biological role in inflammation, as it has now been shown that it also exerts proinflammatory effects by a catalysis-independent mechanism. This mechanism is likely to be highly dependent on key specific molecular interactions, and the full mechanistic descriptions of this remain elusive. The current candidates for the protein partners that may mediate this catalysis-independent mechanism are also introduced in this review. A key discovery has been that selective inhibition of the catalysis-independent activity of hGIIA is achieved with cyclised derivatives of a pentapeptide, FLSYK, derived from the primary sequence of hGIIA. The effects of hGIIA on cell function appear to vary depending on the pathology studied, and so its mechanism of action is complex and context-dependent. This review is comprehensive and covers the most recent developments in the understanding of the many facets of hGIIA function and inhibition and the insight they provide into their clinical application for disease treatment. A cyclic analogue of FLSYK, c2, the most potent analogue known, has now been taken into clinical trials targeting advanced prostate cancer.

Structural Studies based on two Lysine Dioxygenases with Distinct Regioselectivity Brings Insights Into Enzyme Specificity within the Clavaminate Synthase-Like Family

Iron(II)/α-ketoacid-dependent oxygenases (αKAOs) are enzymes that catalyze the oxidation of unactivated C-H bonds, mainly through hydroxylation. Among these, those that are active towards amino-acids and their derivatives are grouped in the Clavaminate Synthase Like (CSL) family. CSL enzymes exhibit high regio- and stereoselectivities with strict substrate specificity. This study reports the structural elucidation of two new regiodivergent members, KDO1 and KDO5, active towards lysine, and the structural and computational analysis of the whole family through modelling and classification of active sites. The structures of KDO1 and KDO5 in complex with their ligands show that one exact position in the active site controls the regioselectivity of the reaction. Our results suggest that the substrate specificity and high stereoselectivity typical of this family is linked to a lid that closes up in order to form a sub-pocket around the side chain of the substrate. This dynamic lid is found throughout the family with varying sequence and length and is associated with a conserved stable dimeric interface. Results from this study could be a starting-point for exploring the functional diversity of the CSL family and direct in vitro screening in the search for new enzymatic activities.

Parallel evolution of non-homologous isofunctional enzymes in methionine biosynthesis

Experimental validation of enzyme function is crucial for genome interpretation, but it remains challenging because it cannot be scaled up to accommodate the constant accumulation of genome sequences. We tackled this issue for the MetA and MetX enzyme families, phylogenetically unrelated families of acyl-L-homoserine transferases involved in L-methionine biosynthesis. Members of these families are prone to incorrect annotation because MetX and MetA enzymes are assumed to always use acetyl-CoA and succinyl-CoA, respectively. We determined the enzymatic activities of 100 enzymes from diverse species, and interpreted the results by structural classification of active sites based on protein structure modeling. We predict that >60% of the 10,000 sequences from these families currently present in databases are incorrectly annotated, and suggest that acetyl-CoA was originally the sole substrate of these isofunctional enzymes, which evolved to use exclusively succinyl-CoA in the most recent bacteria. We also uncovered a divergent subgroup of MetX enzymes in fungi that participate only in L-cysteine biosynthesis as O-succinyl-L-serine transferases.

3-Keto-5-aminohexanoate cleavage enzyme: a common fold for an uncommon Claisen-type condensation

The exponential increase in genome sequencing output has led to the accumulation of thousands of predicted genes lacking a proper functional annotation. Among this mass of hypothetical proteins, enzymes catalyzing new reactions or using novel ways to catalyze already known reactions might still wait to be identified. Here, we provide a structural and biochemical characterization of the 3-keto-5-aminohexanoate cleavage enzyme (Kce), an enzymatic activity long known as being involved in the anaerobic fermentation of lysine but whose catalytic mechanism has remained elusive so far. Although the enzyme shows the ubiquitous triose phosphate isomerase (TIM) barrel fold and a Zn(2+) cation reminiscent of metal-dependent class II aldolases, our results based on a combination of x-ray snapshots and molecular modeling point to an unprecedented mechanism that proceeds through deprotonation of the 3-keto-5-aminohexanoate substrate, nucleophilic addition onto an incoming acetyl-CoA, intramolecular transfer of the CoA moiety, and final retro-Claisen reaction leading to acetoacetate and 3-aminobutyryl-CoA. This model also accounts for earlier observations showing the origin of carbon atoms in the products, as well as the absence of detection of any covalent acyl-enzyme intermediate. Kce is the first representative of a large family of prokaryotic hypothetical proteins, currently annotated as the "domain of unknown function" DUF849.

A novel acyl-CoA beta-transaminase characterized from a metagenome

BACKGROUND

Bacteria are key components in all ecosystems. However, our knowledge of bacterial metabolism is based solely on the study of cultivated organisms which represent just a tiny fraction of microbial diversity. To access new enzymatic reactions and new or alternative pathways, we investigated bacterial metabolism through analyses of uncultivated bacterial consortia.

METHODOLOGY/PRINCIPAL FINDINGS

We applied the gene context approach to assembled sequences of the metagenome of the anaerobic digester of a municipal wastewater treatment plant, and identified a new gene which may participate in an alternative pathway of lysine fermentation.

CONCLUSIONS

We characterized a novel, unique aminotransferase that acts exclusively on Coenzyme A (CoA) esters, and proposed a variant route for lysine fermentation. Results suggest that most of the lysine fermenting organisms use this new pathway in the digester. Its presence in organisms representative of two distinct bacterial divisions indicate that it may also be present in other organisms.

On the characterization and selection of diverse conformational ensembles with applications to flexible docking

To address challenging flexible docking problems, a number of docking algorithms pregenerate large collections of candidate conformers. To remove the redundancy from such ensembles, a central problem in this context is to report a selection of conformers maximizing some geometric diversity criterion. We make three contributions to this problem. First, we resort to geometric optimization so as to report selections maximizing the molecular volume or molecular surface area (MSA) of the selection. Greedy strategies are developed, together with approximation bounds. Second, to assess the efficacy of our algorithms, we investigate two conformer ensembles corresponding to a flexible loop of four protein complexes. By focusing on the MSA of the selection, we show that our strategy matches the MSA of standard selection methods, but resorting to a number of conformers between one and two orders of magnitude smaller. This observation is qualitatively explained using the Betti numbers of the union of balls of the selection. Finally, we replace the conformer selection problem in the context of multiple-copy flexible docking. On the aforementioned systems, we show that using the loops selected by our strategy can improve the result of the docking process.

Accounting for large amplitude protein deformation during in silico macromolecular docking

Rapid progress of theoretical methods and computer calculation resources has turned in silico methods into a conceivable tool to predict the 3D structure of macromolecular assemblages, starting from the structure of their separate elements. Still, some classes of complexes represent a real challenge for macromolecular docking methods. In these complexes, protein parts like loops or domains undergo large amplitude deformations upon association, thus remodeling the surface accessible to the partner protein or DNA. We discuss the problems linked with managing such rearrangements in docking methods and we review strategies that are presently being explored, as well as their limitations and success.

Identification of subfamily-specific sites based on active sites modeling and clustering

MOTIVATION

Current computational approaches to function prediction are mostly based on protein sequence classification and transfer of annotation from known proteins to their closest homologous sequences relying on the orthology concept of function conservation. This approach suffers a major weakness: annotation reliability depends on global sequence similarity to known proteins and is poorly efficient for enzyme superfamilies that catalyze different reactions. Structural biology offers a different strategy to overcome the problem of annotation by adding information about protein 3D structures. This information can be used to identify amino acids located in active sites, focusing on detection of functional polymorphisms residues in an enzyme superfamily. Structural genomics programs are providing more and more novel protein structures at a high-throughput rate. However, there is still a huge gap between the number of sequences and available structures. Computational methods, such as homology modeling provides reliable approaches to bridge this gap and could be a new precise tool to annotate protein functions.

RESULTS

Here, we present Active Sites Modeling and Clustering (ASMC) method, a novel unsupervised method to classify sequences using structural information of protein pockets. ASMC combines homology modeling of family members, structural alignment of modeled active sites and a subsequent hierarchical conceptual classification. Comparison of profiles obtained from computed clusters allows the identification of residues correlated to subfamily function divergence, called specificity determining positions. ASMC method has been validated on a benchmark of 42 Pfam families for which previous resolved holo-structures were available. ASMC was also applied to several families containing known protein structures and comprehensive functional annotations. We will discuss how ASMC improves annotation and understanding of protein families functions by giving some specific illustrative examples on nucleotidyl cyclases, protein kinases and serine proteases.

AVAILABILITY

http://www.genoscope.fr/ASMC/.

Dimorphic motifs in D0 and D1+D2 domains of killer cell Ig-like receptor 3DL1 combine to form receptors with high, moderate, and no avidity for the complex of a peptide derived from HIV and HLA-A*2402

Comparison of mutant killer cell Ig-like receptor (KIR) 3DL1*015 substituted at natural positions of variation showed that tryptophan/leucine dimorphism at position 283 uniquely changes receptor conformation and can strongly influence binding of the A24nef tetramer. Dimorphic motifs at positions 2, 47, and 54 in D0 and 182 and 283 in D1+D2 distinguish the two 3DL1 lineages, typified by 3DL1*005 and 3DL1*015. The interlineage recombinant, KIR3DL1*001, combines D0 of 3DL1*005 with D1+D2 of 3DL1*015 and binds A24nef more strongly than either parent. In contrast, the reciprocal recombinant with D0 from 3DL1*015 and D1+D2 from 3DL1*005 cannot bind A24nef. Thus, D0 polymorphism directly affects the avidity of the KIR3DL1 ligand binding site. From these observations, multiple sequence alignment, and homology modeling, we constructed structural models for KIR3DL1 and its complex with A24nef. In these models, D0, D1, and D2 come together to form a binding surface for A24nef, which is contacted by all three Ig-like domains. A central pocket binds arginine 83, the only Bw4 motif residue essential for KIR3DL1 interaction, similar to the binding of lysine 80 in HLA-C by KIR2DL1. Central to this interaction is a salt bridge between arginine 83 of Bw4 and glutamate 282 of 3DL1, which juxtaposes the functionally influential dimorphism at position 283. Further 3DL1 mutants were tested and shown to have A24nef-binding properties consistent with the models. A24nef was not bound by KIR3DS1, the activating counterpart of KIR3DL1. Moreover, introducing any one of three residues specific to KIR3DS1, serine 163, arginine 166, or leucine 199, into 3DL1*015, abrogated A24nef binding.

Accounting for loop flexibility during protein-protein docking

Although reliable docking can now be achieved for systems that do not undergo important induced conformational change upon association, the presence of flexible surface loops, which must adapt to the steric and electrostatic properties of a partner, generally presents a major obstacle. We report here the first docking method that allows large loop movements during a systematic exploration of the possible arrangements of the two partners in terms of position and rotation. Our strategy consists in taking into account an ensemble of possible loop conformations by a multi-copy representation within a reduced protein model. The docking process starts from regularly distributed positions and orientations of the ligand around the whole receptor. Each starting configuration is submitted to energy minimization during which the best-fitting loop conformation is selected based on the mean-field theory. Trials were carried out on proteins with significant differences in the main-chain conformation of the binding loop between isolated form and complexed form, which were docked to their partner considered in their bound form. The method is able to predict complexes very close to the crystal complex both in terms of relative position of the two partners and of the geometry of the flexible loop. We also show that introducing loop flexibility on the isolated protein form during systematic docking largely improves the predictions of relative position of the partners in comparison with rigid-body docking.

Docking macromolecules with flexible segments

We address a major obstacle to macromolecular docking algorithms by presenting a new method that takes into account the induced conformational adjustment of flexible loops situated at a protein/macromolecule interface. The method, MC2, is based on a multiple copy representation of the loops, coupled with a Monte Carlo conformational search of the relative position of the macromolecules and their side chain conformations. The selection of optimal loop conformations takes place during Monte Carlo cycling by the iterative adjustment of the weight of each copy. We describe here the parameterization of the method and trials on a protein-DNA complex of known 3-D structure, involving the Drosophila prd paired domain protein and its target oligonucleotide Wenqing, X. et al., Cell 1995, 80, 639. We demonstrate that our algorithm can correctly configure and position this protein, despite its relatively complex interactions with both grooves of DNA.