Bifunctional abietadiene synthase: mutual structural dependence of the active sites for protonation-initiated and ionization-initiated cyclizations

Abietadiene synthase from grand fir catalyzes two sequential, mechanistically distinct cyclizations, of geranylgeranyl diphosphate and of copalyl diphosphate, in the formation of a mixture of abietadiene isomers as the committed step of diterpenoid resin acid biosynthesis. Each reaction is independently conducted at a separate active site residing in what were considered to be structurally distinct domains typical of terpene cyclases. Despite the presence of an unusual 250-residue N-terminal insertional element, a tandem pair of charged residues distal to the insertion was shown to form a functional part of the C-terminal active site. Because abietadiene synthase resembles the ancestral plant terpene cyclase, this observation suggests an early evolutionary origin of catalytically important positively charged residues at the N-terminus of enzymes of this general class. A series of N- and C-terminal truncations of this enzyme were constructed and characterized, both alone and as mixtures of adjacent polypeptide pairs, to assess the proposed domain architecture, the function of the insertional element, and the role of presumptive interdomain contacts. These studies indicated a requirement for the insertional element in functional folding and allowed definition of the minimum primary structure of N- and C-terminal active site peptides. Most importantly, the results showed that, although the two active sites of abietadiene synthase are catalytically independent, substantial contact between the two regions is essential for the functional competence of this enzyme. Thus, the two cyclization sites of abietadiene synthase cannot be dissected into catalytically distinct domains, and, therefore, abietadiene synthase is unlikely to have arisen by fusion of two previously independent genes.     

Automatic prediction of catalytic residues by modeling residue structural neighborhood

Background  

Prediction of catalytic residues is a major step in characterizing the function of enzymes. In its simpler formulation, the problem can be cast into a binary classification task at the residue level, by predicting whether the residue is directly involved in the catalytic process. The task is quite hard also when structural information is available, due to the rather wide range of roles a functional residue can play and to the large imbalance between the number of catalytic and non-catalytic residues.     

Physical basis of structural and catalytic Zn-binding sites in proteins

Zn²⁺, an element that is essential to all life forms, can play a catalytic or a solely structural role. Previous works have shown that Zn²⁺ binds preferentially to water molecules and His in catalytic sites, but to Cys⁻ instructural sites, but the molecular basis for the observed ligand preference is unclear. Here, we show that the different Zn²⁺ roles are also reflected in the different bond distances to Zn²⁺ in structural and catalytic sites. We reveal the physical basis for the observed differences between structural and catalytic Zn sites: In most catalytic sites, water is found bound to Zn²⁺ as it transfers the least charge to Zn²⁺ and is less bulky compared to the protein ligands, enabling Zn²⁺ to serve as a Lewis acid in catalysis. In most structural sites, however, ≥ 2 Cys⁻ are found bound to Zn²⁺, as Cys⁻ transfers the most charge to Zn²⁺ and reduces the Zn charge to such an extent that Zn²⁺ can no longer act as a Lewis acid; furthermore, steric repulsion among the bulky Cys(S⁻) prevents Zn²⁺ from accommodating another ligand. Based on the observed ligand preference and Zn-ligand distance differences between structural and catalytic Zn sites, we present a simple method for distinguishing the two types of sites and for verifying the catalytic role of Zn²⁺. Finally, we discuss how the physical bases revealed aid in designing potential drug molecules that target Zn proteins.     

Redesign of artificial globins: effects of residue replacements at hydrophobic sites on the structural properties

Artificial sequences of the 153 amino acids have been designed to fit the main-chain framework of the sperm whale myoglobin (Mb) structure based on a knowledge-based 3D-1D compatibility method. The previously designed artificial globin (DG1) folded into a monomeric, compact, highly helical and globular form with overall dimensions similar to those of the target structure, but it lacked structural uniqueness at the side-chain level [Isogai, Y., Ota, M., Fujisawa, T. , Izuno, H., Mukai, M., Nakamura, H., Iizuka, T., and Nishikawa, K. (1999) Biochemistry 38, 7431-7443]. In this study, we redesigned hydrophobic sites of DG1 to improve the structural specificity. Several Leu and Met residues in DG1 were replaced with beta-branched amino acids, Ile and Val, referring to the 3D profile of DG1 to produce three redesigned globins, DG2-4. These residue replacements resulted in no significant changes of their compactness and alpha-helical contents in the absence of denaturant, whereas they significantly affected the dependence of the secondary structure on the concentration of guanidine hydrochloride. The analyses of the denaturation curves revealed higher global stabilities of the designed globins than that of natural apoMb. Among DG1-4, DG3, in which 11 Leu residues of DG1 are replaced with seven Ile and four Val residues, and one Met residue is replaced with Val, displayed the lowest stability but the most cooperative folding-unfolding transition and the most dispersed NMR spectrum with the smallest line width. The present results indicate that the replacements of Leu (Met) with the beta-branched amino acids at appropriate sites reduce the freedom of side-chain conformation and improve the structural specificity at the expense of stability.     

Beyond structural genomics: computational approaches for the identification of ligand binding sites in protein structures

Structural genomics projects have revealed structures for a large number of proteins of unknown function. Understanding the interactions between these proteins and their ligands would provide an initial step in their functional characterization. Binding site identification methods are a fast and cost-effective way to facilitate the characterization of functionally important protein regions. In this review we describe our recently developed methods for binding site identification in the context of existing methods. The advantage of energy-based approaches is emphasized, since they provide flexibility in the identification and characterization of different types of binding sites.     

Functional identification of catalytic metal ion binding sites within RNA

The viability of living systems depends inextricably on enzymes that catalyze phosphoryl transfer reactions. For many enzymes in this class, including several ribozymes, divalent metal ions serve as obligate cofactors. Understanding how metal ions mediate catalysis requires elucidation of metal ion interactions with both the enzyme and the substrate(s). In the Tetrahymena group I intron, previous work using atomic mutagenesis and quantitative analysis of metal ion rescue behavior identified three metal ions (MA, MB, and MC) that make five interactions with the ribozyme substrates in the reaction's transition state. Here, we combine substrate atomic mutagenesis with site-specific phosphorothioate substitutions in the ribozyme backbone to develop a powerful, general strategy for defining the ligands of catalytic metal ions within RNA. In applying this strategy to the Tetrahymena group I intron, we have identified the pro-SP phosphoryl oxygen at nucleotide C262 as a ribozyme ligand for MC. Our findings establish a direct connection between the ribozyme core and the functionally defined model of the chemical transition state, thereby extending the known set of transition-state interactions and providing information critical for the application of the recent group I intron crystallographic structures to the understanding of catalysis.     

M42 aminopeptidase catalytic site: the structural and functional role of a strictly conserved aspartate residue

The M42 aminopeptidases are a family of dinuclear aminopeptidases widely distributed in Prokaryotes. They are potentially associated to the proteasome, achieving complete peptide destruction. Their most peculiar characteristic is their quaternary structure, a tetrahedron-shaped particle made of twelve subunits. The catalytic site of M42 aminopeptidases is defined by seven conserved residues. Five of them are involved in metal ion binding which is important to maintain both the activity and the oligomeric state. The sixth conserved residue, a glutamate, is the catalytic base deprotonating the water molecule during peptide bond hydrolysis. The seventh residue is an aspartate whose function remains poorly understood. This aspartate residue, however, must have a critical role as it is strictly conserved in all MH clan enzymes. It forms some kind of catalytic triad with the histidine residue and the metal ion of the M2 binding site. We assess its role in TmPep1050, an M42 aminopeptidase of Thermotoga maritima, through a mutational approach. Asp-62 was substituted with alanine, asparagine, or glutamate residue. The Asp-62 substitutions completely abolished TmPep1050 activity and impeded dodecamer formation. They also interfered with metal ion binding as only one cobalt ion is bound per subunit instead of two. The structure of Asp62Ala variant was solved at 1.5 Å showing how the substitution has an impact on the active site fold. We propose a structural role for Asp-62, helping to stabilize a crucial loop in the active site and to position correctly the catalytic base and a metal ion ligand of the M1 site.     

Structure of the catalytic fragment of translation initiation factor 2B and identification of a critically important catalytic residue

Eukaryotic initiation factor (eIF) 2B catalyzes the nucleotide activation of eIF2 to its active GTP-bound state. The exchange activity has been mapped to the C terminus of the eIF2Bepsilon subunit. We have determined the crystal structure of residues 544-704 from yeast eIF2Bepsilon at 2.3-A resolution, and this fragment is an all-helical protein built around the conserved aromatic acidic (AA) boxes also found in eIF4G and eIF5. The eight helices are organized in a manner similar to HEAT repeats. The molecule is highly asymmetric with respect to surface charge and conservation. One area in the N terminus is proposed to be directly involved in catalysis. In agreement with this hypothesis, mutation of glutamate 569 is shown to be lethal. An acidic belt and a second area in the C terminus containing residues from the AA boxes are important for binding to eIF2. Two mutations causing the fatal human genetic disease leukoencephalopathy with vanishing white matter are buried and appear to disrupt the structural integrity of the catalytic domain rather than interfering directly with catalysis or binding of eIF2.     

Comparison of allosteric signaling in DnaK and BiP using mutual information between simulated residue conformations

The heat shock protein 70 kDa (Hsp70) chaperone system serves as a critical component of protein quality control across a wide range of prokaryotic and eukaryotic organisms. Divergent evolution and specialization to particular organelles have produced numerous Hsp70 variants which share similarities in structure and general function, but differ substantially in regulatory aspects, including conformational dynamics and activity modulation by cochaperones. The human Hsp70 variant BiP (also known as GRP78 or HSPA5) is of therapeutic interest in the context of cancer, neurodegenerative diseases, and viral infection, including for treatment of the pandemic virus SARS-CoV-2. Due to the complex conformational rearrangements and high sequential variance within the Hsp70 protein family, it is in many cases poorly understood which amino acid mutations are responsible for biochemical differences between protein variants. In this study, we predicted residues associated with conformational regulation of human BiP and Escherichia coli DnaK. Based on protein structure networks obtained from molecular dynamics simulations, we analyzed the shared information between interaction timelines to highlight residue positions with strong conformational coupling to their environment. Our predictions, which focus on the binding processes of the chaperone's substrate and cochaperones, indicate residues filling potential signaling roles specific to either DnaK or BiP. By combining predictions of individual residues into conformationally coupled chains connecting ligand binding sites, we predict a BiP specific secondary signaling pathway associated with substrate binding. Our study sheds light on mechanistic differences in signaling and regulation between Hsp70 variants, which provide insights relevant to therapeutic applications of these proteins.     

Assessing the structural conservation of protein pockets to study functional and allosteric sites: implications for drug discovery

Background

Nonribosomal peptide synthetases (NRPSs) are multienzymatic, multidomain megasynthases involved in the biosynthesis of pharmaceutically important nonribosomal peptides. The peptaibol synthetase from Trichoderma virens (TPS) is an important member of the NRPS family that exhibits antifungal properties. The majority of the NRPSs terminate peptide synthesis with the thioesterase (TE) domain, which either hydrolyzes the thioester linkage, releasing the free peptic acid, or catalyzes the intramolecular macrocyclization to produce a macrolactone product. TPS is an important NRPS that does not encompass a TE domain, but rather a reductase domain (R domain) to release the mature peptide product reductively with the aid of a NADPH cofactor. However, the catalytic mechanism of the reductase domain has not yet been elucidated.

Results

We present here a three-dimensional (3D) model of the reductase domain based on the crystal structure of vestitone reductase (VR). VR belongs to the short-chain dehydrogenase/reductase (SDR) superfamily and is responsible for the nicotinamide dinucleotide phosphate (NADPH)-dependent reduction of the substrate into its corresponding secondary alcohol product. The binding sites of the probable linear substrates, alamethicin, trichotoxin, antiamoebin I, chrysopermin C and gramicidin, were identified within the modeled R domain using multiple docking approaches. The docking results of the ligand in the active site of the R domain showed that reductase side chains have a high affinity towards ligand binding, while the thioester oxygen of each substrate forms a hydrogen bond with the OH group of Tyr176 and the thiol group of the substrate is closer to the Glu220. The modeling and docking studies revealed the reaction mechanism of reduction of thioester into a primary alcohol.

Conclusion

Peptaibol biosynthesis incorporates a single R domain, which appears to catalyze the four-electron reduction reaction of a peptidyl carrier protein (PCP)-bound peptide to its corresponding primary alcohol. Analysis of R domains present in the non-redundant (nr) database of the NCBI showed that the R domain always resides in the last NRPS module and is involved in either a two or four-electron reduction reaction.     

The structural basis of the recognition of phenylalanine and pterin cofactors by phenylalanine hydroxylase: implications for the catalytic mechanism

Phenylalanine hydroxylase (PAH) is a tetrahydrobiopterin and non-heme iron-dependent enzyme that hydroxylates L-Phe to l-Tyr using molecular oxygen as additional substrate. A dysfunction of this enzyme leads to phenylketonuria (PKU). The conformation and distances to the catalytic iron of both L-Phe and the cofactor analogue L-erythro-7,8-dihydrobiopterin (BH2) simultaneously bound to recombinant human PAH have been estimated by (1)H NMR. The resulting bound conformers of both ligands have been fitted into the crystal structure of the catalytic domain by molecular docking. In the docked structure L-Phe binds to the enzyme through interactions with Arg270, Ser349 and Trp326. The mode of coordination of Glu330 to the iron moiety seems to determine the amino acid substrate specificity in PAH and in the homologous enzyme tyrosine hydroxylase. The pterin ring of BH2 pi-stacks with Phe254, and the N3 and the amine group at C2 hydrogen bond with the carboxylic group of Glu286. The ring also establishes specific contacts with His264 and Leu249. The distance between the O4 atom of BH2 and the iron (2.6(+/-0.3) A) is compatible with coordination, a finding that is important for the understanding of the mechanism of the enzyme. The hydroxyl groups in the side-chain at C6 hydrogen bond with the carbonyl group of Ala322 and the hydroxyl group of Ser251, an interaction that seems to have implications for the regulation of the enzyme by substrate and cofactor. Some frequent mutations causing PKU are located at residues involved in substrate and cofactor binding. The sites for hydroxylation, C4 in L-Phe and C4a in the pterin are located at a distance of 4.2 and 4.3 A from the iron moiety, respectively, and at 6.3 A from each other. These distances are adequate for the intercalation of iron-coordinated molecular oxygen, in agreement with a mechanistic role of the iron moiety both in the binding and activation of dioxygen and in the hydroxylation reaction.     

Role of alphaPhe-291 residue in the phosphate-binding subdomain of catalytic sites of Escherichia coli ATP synthase

The role of αPhe-291 residue in phosphate binding by Escherichia coli F₁F₀-ATP synthase was examined. X-ray structures of bovine mitochondrial enzyme suggest that this residue resides in close proximity to the conserved βR246 residue. Herein, we show that mutations αF291D and αF291E in E. coli reduce the ATPase activity of F₁F₀ membranes by 350-fold. Yet, significant oxidative phosphorylation activity is retained. In contrast to wild-type, ATPase activities of mutants were not inhibited by MgADP-azide, MgADP-fluoroaluminate, or MgADP-fluoroscandium. Whereas, 7-chloro-4-nitrobenzo-2-oxa-1,3-diazole (NBD-Cl) inhibited wild-type ATPase essentially completely, ATPase in mutants was inhibited maximally by ∼75%, although reaction still occurred at residue βTyr-297, proximal to αPhe-291 in the phosphate-binding pocket. Inhibition characteristics supported the conclusion that NBD-Cl reacts in βE (empty) catalytic sites, as shown previously by X-ray structure analysis. Phosphate protected against NBD-Cl inhibition in wild-type but not in mutants. In addition, our data suggest that the interaction of αPhe-291 with phosphate during ATP hydrolysis or synthesis may be distinct.     

The sites of tryptic cleavage in bovine secretory component: structural and functional implications

Tryptic fragments from bovine secretory component and sIgA have been separated by HPLC and/or SDS polyacrylamide gel electrophoresis and electroblotting. Their N-terminal amino acid sequences have been determined and their positions in the secretory component molecule deduced by homology with the amino acid sequences of human secretory component and rabbit polyimmunoglobulin receptor. Taken in conjunction with the known binding affinities of the tryptic fragments, the results imply that the three most N-terminal domains of secretory component are directly involved in binding IgM and IgA dimers. The results also favour the concept of an extended 'zig-zag' structure for the secretory component molecule.     

Networks of coevolving sites in structural and functional domains of serpin proteins

Amino acids do not occur randomly in proteins; rather, their occurrence at any given site is strongly influenced by the amino acid composition at other sites, the structural and functional aspects of the region of the protein in which they occur, and the evolutionary history of the protein. The goal of our research study is to identify networks of coevolving sites within the serpin proteins (serine protease inhibitors) and classify them as being caused by structural-functional constraints or by evolutionary history. To address this, a matrix of pairwise normalized mutual information (NMI) values was computed among amino acid sites for the serpin proteins. The NMI matrix was partitioned into orthogonal patterns of amino acid variability by factor analysis. Each common factor pattern was interpreted as having phylogenetic and/or structural-functional explanations. In addition, we used a bootstrap factor analysis technique to limit the effects of phylogenetic history on our factor patterns. Our results show an extensive network of correlations among amino acid sites in key functional regions (reactive center loop, shutter, and breach). Additionally, we have discovered long-range coevolution for packed amino acids within the serpin protein core. Lastly, we have discovered a group of serpin sites which coevolve in the hydrophobic core region (s5B and s4B) and appear to represent sites important for formation of the "native" instead of the "latent" serpin structure. This research provides a better understanding on how protein structure evolves; in particular, it elucidates the selective forces creating coevolution among protein sites.     

COLLAPSE: A representation learning framework for identification and characterization of protein structural sites

The identification and characterization of the structural sites which contribute to protein function are crucial for understanding biological mechanisms, evaluating disease risk, and developing targeted therapies. However, the quantity of known protein structures is rapidly outpacing our ability to functionally annotate them. Existing methods for function prediction either do not operate on local sites, suffer from high false positive or false negative rates, or require large site‐specific training datasets, necessitating the development of new computational methods for annotating functional sites at scale. We present COLLAPSE (Compressed Latents Learned from Aligned Protein Structural Environments), a framework for learning deep representations of protein sites. COLLAPSE operates directly on the 3D positions of atoms surrounding a site and uses evolutionary relationships between homologous proteins as a self‐supervision signal, enabling learned embeddings to implicitly capture structure–function relationships within each site. Our representations generalize across disparate tasks in a transfer learning context, achieving state‐of‐the‐art performance on standardized benchmarks (protein–protein interactions and mutation stability) and on the prediction of functional sites from the prosite database. We use COLLAPSE to search for similar sites across large protein datasets and to annotate proteins based on a database of known functional sites. These methods demonstrate that COLLAPSE is computationally efficient, tunable, and interpretable, providing a general‐purpose platform for computational protein analysis.     

Looking at enzymes from the inside out: the proximity of catalytic residues to the molecular centroid can be used for detection of active sites and enzyme-ligand interfaces

Analysis of the distances of the exposed residues in 175 enzymes from the centroids of the molecules indicates that catalytic residues are very often found among the 5% of residues closest to the enzyme centroid. This property of catalytic residues is implemented in a new prediction algorithm (named EnSite) for locating the active sites of enzymes and in a new scheme for re-ranking enzyme-ligand docking solutions. EnSite examines only 5% of the molecular surface (represented by surface dots) that is closest to the centroid, identifying continuous surface segments and ranking them by their area size. EnSite ranks the correct prediction 1-4 in 97% of the cases in a dataset of 65 monomeric enzymes (rank 1 for 89% of the cases) and in 86% of the cases in a dataset of 176 monomeric and multimeric enzymes from all six top-level enzyme classifications (rank 1 in 74% of the cases). Importantly, identification of buried or flat active sites is straightforward because EnSite "looks" at the molecular surface from the inside out. Detailed examination of the results indicates that the proximity of the catalytic residues to the centroid is a property of the functional unit, defined as the assembly of domains or chains that form the active site (in most cases the functional unit corresponds to a single whole polypeptide chain). Using the functional unit in the prediction further improves the results. The new property of active sites is also used for re-evaluating enzyme-inhibitor unbound docking results. Sorting the docking solutions by the distance of the interface to the centroid of the enzyme improves remarkably the ranks of nearly correct solutions compared to ranks based on geometric-electrostatic-hydrophobic complementarity scores.