Insights derived from molecular dynamics simulation into the molecular motions of serine protease proteinase K

Serine protease proteinase K, a member of the subtilisin family of enzymes, is of significant industrial, agricultural and biotechnological importance. Despite the wealth of structural information about proteinase K provided by static X-ray structures, a full understanding of the enzymatic mechanism requires further insight into the dynamic properties of this enzyme. Molecular dynamics simulations and essential dynamics (ED) analysis were performed to investigate the molecular motions in proteinase K. The results indicate that the internal core of proteinase K is relatively rigid, whereas the surface-exposed loops, most notably the substrate-binding regions, exhibit considerable conformational fluctuations. Further ED analysis reveals that the large concerted motions in the substrate-binding regions cause opening/closing of the substrate-binding pockets, thus supporting the proposed induced-fit mechanism of substrate binding. The distinct electrostatic/hydrogen-bonding interactions between Asp39 and His69 and between His69 and Ser224 within the catalytic triad lead to different thermal motions and orientations of these three catalytic residues, which can be related to their different functional roles in the catalytic process. Statistical analyses of the geometrical/functional properties as well as evolutionary conservation of the glycines in proteinase K-like proteins reveal that glycines may play an important role in determining the folding architecture and structural flexibility of this class of enzymes. Our simulation study complements the biochemical and structural studies and provides new insights into the dynamic structural basis of the functional properties of this class of enzymes.     

Protein interactions and fluctuations in a proteomic network using an elastic network model

A set of protein conformations are analyzed by normal mode analysis. An elastic network model is used to obtain fluctuation and cooperativity of residues with low amplitude fluctuations across different species. Slow modes that are associated with the function of proteins have common features among different protein structures. We show that the degree of flexibility of the protein is important for proteins to interact with other proteins and as the species gets more complex its proteins become more flexible. In the complex organism, higher cooperativity arises due to protein structure and connectivity.     

Large domain fluctuations on 50-ns timescale enable catalytic activity in phosphoglycerate kinase

Large-scale domain motions of enzymes are often essential for their biological function. Phosphoglycerate kinase has a wide open domain structure with a hinge near the active center between the two domains. Applying neutron spin echo spectroscopy and small-angle neutron scattering we have investigated the internal domain dynamics. Structural analysis reveals that the holoprotein in solution seems to be more compact compared to the crystal structure but would not allow the functionally important phosphoryl transfer between the substrates if the protein were static. Brownian large-scale domain fluctuation dynamics on a timescale of 50 ns was revealed by neutron spin echo spectroscopy. The dynamics observed was compared to the displacement patterns of low-frequency normal modes. The displacements along the normal-mode coordinates describe our experimental results reasonably well. In particular, the domain movements facilitate a close encounter of the key residues in the active center to build the active configuration. The observed dynamics shows that the protein has the flexibility to allow fluctuations and displacements that seem to enable the function of the protein. Moreover, the presence of the substrates increases the rigidity, which is deduced from a faster dynamics with smaller amplitude.     

Protein Interactions and Fluctuations in a Proteomic Network using an Elastic Network Model

Abstract

A set of protein conformations are analyzed by normal mode analysis. An elastic network model is used to obtain fluctuation and cooperativity of residues with low amplitude fluctuations across different species. Slow modes that are associated with the function of proteins have common features among different protein structures. We show that the degree of flexibility of the protein is important for proteins to interact with other proteins and as the species gets more complex its proteins become more flexible. In the complex organism, higher cooperativity arises due to protein structure and connectivity.     

Conservation of structural fluctuations in homologous protein kinases and its implications on functional sites

Our aim is to explore the similarities in structural fluctuations of homologous kinases. Gaussian Network Model based Normal Mode Analysis was performed on 73 active conformation structures in Ser/Thr/Tyr kinase superfamily. Categories of kinases with progressive evolutionary divergence, viz. (i) Same kinase with many crystal structures, (ii) Within‐Subfamily, (iii) Within‐Family, (iv) Within‐Group, and (v) Across‐Group, were analyzed. We identified a flexibility signature conserved in all kinases involving residues in and around the catalytic loop with consistent low‐magnitude fluctuations. However, the overall structural fluctuation profiles are conserved better in closely related kinases (Within‐Subfamily and Within‐family) than in distant ones (Within‐Group and Across‐Group). A substantial 65.4% of variation in flexibility was not accounted by variation in sequences or structures. Interestingly, we identified substructural residue‐wise fluctuation patterns characteristic of kinases of different categories. Specifically, we recognized statistically significant fluctuations unique to families of protein kinase A, cyclin‐dependent kinases, and nonreceptor tyrosine kinases. These fluctuation signatures localized to sites known to participate in protein‐protein interactions typical of these kinase families. We report for the first time that residues characterized by fluctuations unique to the group/family are involved in interactions specific to the group/family. As highlighted for Src family, local regions with differential fluctuations are proposed as attractive targets for drug design. Overall, our study underscores the importance of consideration of fluctuations, over and above sequence and structural features, in understanding the roles of sites characteristic of kinases. Proteins 2016; 84:957–978. © 2016 Wiley Periodicals, Inc.     

Comparison of full-atomic and coarse-grained models to examine the molecular fluctuations of c-AMP dependent protein kinase

Molecular fluctuations of the native conformation of c-AMP dependent protein kinase (cAPK) have been investigated with three different approaches. The first approach is the full atomic normal mode analysis (NMA) with empirical force fields. The second and third approaches are based on a coarse-grained model with a single single-parameter- harmonic potential between close residues in the crystal structure of the molecule without any residue specificity. The second method calculates only the magnitude of fluctuations whereas the third method is developed to find the directionality of the fluctuations which are essential to understand the functional importance of biological molecules. The aim, in this study, is to determine whether using such coarse-grained models are appropriate for elucidating the global dynamic characteristics of large proteins which reduces the size of the system at least by a factor of ten. The mean-square fluctuations of C(alpha) atoms and the residue cross-correlations are obtained by three approaches. These results are then compared to test the results of coarse grained models on the overall collective motions. All three of the approaches show that highly flexible regions correspond to the activation and solvent exposed loops, whereas the conserved residues (especially in substrate binding regions) exhibit almost no flexibility, adding stability to the structure. The anti-correlated motions of the two lobes of the catalytic core provide flexibility to the molecule. High similarities among the results of these methods indicate that the slowest modes governing the most global motions are preserved in the coarse grained models for proteins. This finding may suggest that the general shapes of the structures are representative of their dynamic characteristics and the dominant motions of protein structures are robust at coarse-grained levels.     

Relationship of protein flexibility to thermostability

Thermostability of proteins arises from the simultaneous effect of several forces, which in fact lead to decreased flexibility of the polypeptide chain. This is verified by flexibility indices, which are derived from normalized B-values of individual amino acids in several refined three-dimensional structures. Flexibility indices show that overall flexibility is reduced when thermostability is increased. Protein molecules require both flexibility and rigidity to function, but the higher the temperature optimum and stability the more rigid is the structure needed to compensate for increased thermal fluctuations. Flexibilities of proteins performing the same catalytic activity seem to be about the same at their temperature optima, but the more rigid thermostable proteins reach the flexibility of thermolabile proteins at higher temperatures. In several proteins such as allosteric enzymes, some local sites of flexibility are highly conserved. The relevance of reduced flexibility to overall stability of proteins is also discussed. Flexibility indices and profiles can be used in the design of more stable proteins by site-directed mutagenesis.     

Comparison of Full-atomic and Coarse-grained Models to Examine the Molecular Fluctuations of c-AMP Dependent Protein Kinase

Abstract

Molecular fluctuations of the native conformation of c-AMP dependent protein kinase (cAPK) have been investigated with three different approaches. The first approach is the full atomic normal mode analysis (NMA) with empirical force fields. The second and third approaches are based on a coarse-grained model with a single single-parameter- harmonic potential between close residues in the crystal structure of the molecule without any residue specificity. The second method calculates only the magnitude of fluctuations whereas the third method is developed to find the directionality of the fluctuations which are essential to understand the functional importance of biological molecules. The aim, in this study, is to determine whether using such coarse-grained models are appropriate for elucidating the global dynamic characteristics of large proteins which reduces the size of the system at least by a factor of ten. The mean-square fluctuations of C^α atoms and the residue cross-correlations are obtained by three approaches. These results are then compared to test the results of coarse grained models on the overall collective motions. AH three of the approaches show that highly flexible regions correspond to the activation and solvent exposed loops, whereas the conserved residues (especially in substrate binding regions) exhibit almost no flexibility, adding stability to the structure. The anti-correlated motions of the two lobes of the catalytic core provide flexibility to the molecule. High similarities among the results of these methods indicate that the slowest modes governing the most global motions are preserved in the coarse grained models for proteins. This finding may suggest that the general shapes of the structures are representative of their dynamic characteristics and the dominant motions of protein structures are robust at coarse-grained levels.     

Investigating the local flexibility of functional residues in hemoproteins

It is now widely accepted that protein function depends not only on structure, but also on flexibility. However, the way mechanical properties contribute to catalytic mechanisms remains unclear. Here, we propose a method for investigating local flexibility within protein structures that combines a reduced protein representation with Brownian dynamics simulations. An analysis of residue fluctuations during the dynamics simulation yields a rigidity profile for the protein made up of force constants describing the ease of displacing each residue with respect to the rest of the structure. This approach has been applied to the analysis of a set of hemoproteins, one of the functionally most diverse protein families. Six proteins containing one or two heme groups have been studied, paying particular attention to the mechanical properties of the active-site residues. The calculated rigidity profiles show that active site residues are generally associated with high force constants and thus rigidly held in place. This observation also holds for diheme proteins if their mechanical properties are analyzed domain by domain. We note, however, that residues other than those in the active site can also have high force constants, as in the case of residues belonging to the folding nucleus of c-type hemoproteins.     

Flexibility analysis of enzyme active sites by crystallographic temperature factors

Protein flexibility is inherent to protein structural behavior. Experimental evidence for protein flexibility is extensive both in solution and in the solid state. A major question is whether the flexibility observed in enzymes is simply an inherent property of proteins that must always be borne in mind or is essential for catalysis or substrate binding. The temperature factors or B-values, as determined crystallographically, are linearly related to the mean square displacement of an atom and give an indication of atomic flexibility in the crystalline state. In this paper, we describe the frequency distributions of the normalized B-factor (B'-factor) for the active site and non-active site residues in the selected 69 apo-enzymes. This analysis was performed over the entire sequences and for different structural subsets defined by the three-dimensional structure of proteins, as alpha-helices, beta-structures and coil conformation and buried and non-buried residues. The results show that in all cases, the active site residues predominantly occur in region of low B'-factor and the non-active site residues have a tendency to exist in the high B'-factor region. This observation suggests that the active site residues, in general, are less flexible than the non-active site residues and therefore the vibrational and the fast collective motions of the C(alpha) atoms of proteins appear not to have clear biological significance.     

Structural analysis of a glycoside hydrolase family 43 arabinoxylan arabinofuranohydrolase in complex with xylotetraose reveals a different binding mechanism compared with other members of the same family

AXHs (arabinoxylan arabinofuranohydrolases) are alpha-L-arabinofuranosidases that specifically hydrolyse the glycosidic bond between arabinofuranosyl substituents and xylopyranosyl backbone residues of arabinoxylan. Bacillus subtilis was recently shown to produce an AXH that cleaves arabinose units from O-2- or O-3-mono-substituted xylose residues: BsAXH-m2,3 (B. subtilis AXH-m2,3). Crystallographic analysis reveals a two-domain structure for this enzyme: a catalytic domain displaying a five-bladed beta-propeller fold characteristic of GH (glycoside hydrolase) family 43 and a CBM (carbohydrate-binding module) with a beta-sandwich fold belonging to CBM family 6. Binding of substrate to BsAXH-m2,3 is largely based on hydrophobic stacking interactions, which probably allow the positional flexibility needed to hydrolyse both arabinose substituents at the O-2 or O-3 position of the xylose unit. Superposition of the BsAXH-m2,3 structure with known structures of the GH family 43 exo-acting enzymes, beta-xylosidase and alpha-L-arabinanase, each in complex with their substrate, reveals a different orientation of the sugar backbone.     

Analysis of catalytic residues in enzyme active sites

We present an analysis of the residues directly involved in catalysis in 178 enzyme active sites. Specific criteria were derived to define a catalytic residue, and used to create a catalytic residue dataset, which was then analysed in terms of properties including secondary structure, solvent accessibility, flexibility, conservation, quaternary structure and function. The results indicate the dominance of a small set of amino acid residues in catalysis and give a picture of a general active site environment. It is hoped that this information will provide a better understanding of the molecular mechanisms involved in catalysis and a heuristic basis for predicting catalytic residues in enzymes of unknown function.     

Flexibility of myosin-subfragment-1 in its complex with actin as revealed by fluorescence resonance energy transfer.

The flexibility of the acto-myosin complex in rigor conditions was characterized by measuring the temperature profile of normalized fluorescence resonance energy transfer efficiency, f' [Somogyi, B., Matkó, J., Papp, S., Hevessy, J., Welch, G.R. & Damjanovich, S. (1984) Biochemistry 23, 3403-3411]. Fluorescence acceptors were introduced to the Cys374 residues of actin and the donors were covalently attached either to Cys707 in the catalytic domain or to Cys177 in the essential light-chain of myosin S1. Fluorescence resonance energy transfer measurements revealed that the protein matrix between Cys374 of actin and Cys707 of S1 is rigid. In contrast, the link between the catalytic and light-chain-binding domains in myosin S1 is flexible. We have recently shown that the positional distribution of Cys707 was narrow relative to the actin filament, while that of the Cys177 was broad. Accordingly, the broad positional distribution of Cys177 is likely to be due to the large flexibility of the link between the catalytic and light-chain-binding domains. This flexibility is probably essential for the interdomain reorganization of the myosin head during the force generation process and for accommodating the symmetry difference between actin and myosin filaments to allow the formation of cross-bridges.     

Insertion of new sequences into the catalytic domain of an enzyme

Activities of enzymes can be modified by the replacement of active-site amino acids with residues that strengthen specific interactions with substrates or that alter the specificity. The scope for engineered enzymes would be broadened if additional, new sequences could be inserted into a catalytic domain. Properly designed, these sequences could encode new ligand binding sites, be intermediates in the construction of chimeric enzymes, or alter the internal flexibility and "breathing" modes of the active-site region. As a first step toward this objective, we inserted oligopeptides of up to 14 amino acids into various locations within an 82 amino acid region of the adenylate synthesis domain of Escherichia coli methionyl-tRNA synthetase. These sites include ones that are flanked by sequences that are conserved between the proteins from E. coli and the yeast Saccharomyces cerevisiae and those that are essential for activity and stability. We found that all of the insertional mutants are stable and some have catalytic parameters for adenylate synthesis that are comparable to those of the wild-type enzyme. Thus, such an approach may provide for a variety of novel applications.     

New insights into the active site structure and catalytic mechanism of tyrosinase and its related proteins

Tyrosinases are widely distributed in nature. They are copper‐containing oxidases belonging to the type 3 copper protein family, together with catechol oxidases and haemocyanins. Tyrosinases are essential enzymes in melanin biosynthesis and therefore responsible for pigmentation of skin and hair in mammals, where two more enzymes, the tyrosinase‐related proteins (Tyrps), participate in the pathway. The structure and catalytic mechanism of mammalian tyrosinases have been extensively studied but they are not completely understood because of the lack of information on the tertiary structure. The availability of crystallographic data of one plant catechol oxidase and one bacterial tyrosinase has improved the model of the three‐dimensional structure of the active site of the enzyme. Furthermore, sequence comparison of tyrosinase and the Tyrps reveals that the three orthologue proteins share many key structural features, because of their common origin from an ancestral gene, although the specific residues responsible for their different catalytic capabilities have not been identified yet. This review summarizes our current knowledge of tyrosinase and Tyrps structure and function and describes the catalytic mechanism of tyrosinase and Dct/Tyrp2, which are better characterized.     

Mechanism of hyaluronan degradation by Streptococcus pneumoniae hyaluronate lyase. Structures of complexes with the substrate

Hyaluronate lyase enzymes degrade hyaluronan, the main polysaccharide component of the host connective tissues, predominantly into unsaturated disaccharide units, thereby destroying the normal connective tissue structure and exposing the tissue cells to various endo- and exogenous factors, including bacterial toxins. The crystal structures of Streptococcus pneumoniae hyaluronate lyase with tetra- and hexasaccharide hyaluronan substrates bound in the active site were determined at 1.52- and 2.0-A resolution, respectively. Hexasaccharide is the longest substrate segment that binds entirely within the active site of these enzymes. The enzyme residues responsible for substrate binding, positioning, catalysis, and product release were thereby identified and their specific roles characterized. The involvement of three residues in catalysis, Asn(349), His(399), and Tyr(408), is confirmed, and the details of proton acceptance and donation within the catalytic machinery are described. The mechanism of processivity of the enzyme is analyzed. The flexibility (allosteric) behavior of the enzyme may be understood in terms of the results of flexibility analysis of this protein, which identified two modes of motion that are also proposed to be involved in the hyaluronan degradation process. The first motion describes an opening and closing of the catalytic cleft located between the alpha- and beta-domains. The second motion demonstrates the mobility of a binding cleft, which may facilitate the binding of the negatively charged hyaluronan to the enzyme.     

The roles of highly conserved,non‐catalytic residues in class A β‐lactamases

Evolution minimizes the number of highly conserved amino acid residues in proteins to ensure evolutionary robustness and adaptability. The roles of all highly conserved, non‐catalytic residues, 11% of all residues, in class A β‐lactamase were analyzed by studying the effect of 146 mutations on in cell and in vitro activity, folding, structure, and stability. Residues around the catalytic residues (second shell) contribute to fine‐tuning of the active site structure. Mutations affect the structure over the entire active site and can result in stable but inactive protein. Conserved residues farther away (third shell) ensure a favorable balance of folding versus aggregation or stabilize the folded form over the unfolded state. Once folded, the mutant enzymes are stable and active and show only localized structural effects. These residues are found in clusters, stapling secondary structure elements. The results give an integral picture of the different roles of essential residues in enzymes.     

Analysis and prediction of the location of catalytic residues in enzymes

The catalytic residues of an enzyme are defined as the amino acids directly involved in chemical catalysis. They mainly act as a general acid--base, electrophilic or nucleophilic catalyst or they polarize and stabilize the transition state. An analysis of the structural features of 36 catalytic residues in 17 enzymes of known structure and with defined mechanism is reported. Residues that bind metal ions (Zn2+ and Cu2+) are considered separately. The features examined are: residue type, location in secondary structure, separation between the residues, accessibility to solvent, intra-protein electrostatic interactions, mobility as evaluated from crystallographic temperature factors, polarity of the environment and the sequence conservation between homologous enzymes of residues that were sequentially or spatially close to the catalytic residue. In general the environment of catalytic residues is similar to that of polar side chains that have low accessibility to solvent. Two algorithms have been developed to identify probable catalytic residues. Scanning an alignment of homologous enzyme sequences for peaks of sequence conservation identifies 13 out of the 16 catalytic residues with 50 residues overpredicted. When the conservation of the spatially close residues is used instead, a different set of 13 residues are identified with 47 residues overpredicted. A combination of the two algorithms identifies 11 residues with 36 residues overpredicted.     

Bioinformatics and molecular modeling in chemical enzymology. Active sites of hydrolases

Comparison and multiple alignments of amino acid sequences of a representative number of related enzymes demonstrate the existence of certain positions of amino acid residues which are permanently reproducible in all members of the whole family. The use of the bioinformatic approach revealed conservative residues in each of the related enzymes and ranked amino acid conservatism for the overall enzymatic catalysis. Glycine and aspartic acid residues were shown to be the most essential for structure and catalytic activity of enzymes. Amino acid residues forming catalytic subsite of the active site of enzymes are always highly conservative. Analysis revealed that aspartic acid carboxyl group is the most frequently employed nucleophilic (in deprotonated form) and electrophilic (in protonated form) agent involved in activation of molecules by the mechanism of general base and acidic catalyses in the catalytic sites of enzymes. Glycine is a unique amino acid possessing the highest possibilities for rotation along C–C and C–N bonds of the polypeptide chain. The conservative fixation of the glycine residue in polypeptide chains of related enzymes provides a possibility for directed assembly of amino acid residues into the catalytic subsite structure. It is possible that the conservative glycines provide known conformational mobility of the protein and the active site. Methods of molecular modeling were used for analysis of structural substitutions of conservative and non-conservative glycines and their effects on geometry of catalytic site of typical hydrolases. The substitution of glycine(s) for alanine significantly altered the catalytic site structures.     

Coupling between catalytic site and collective dynamics: a requirement for mechanochemical activity of enzymes

Growing evidence supports the view that enzymatic activity results from a subtle interplay between chemical kinetics and molecular motions. A systematic analysis is performed here to delineate the type and level of coupling between catalysis and conformational mechanics. The dynamics of a set of 98 enzymes representative of different EC classes are analyzed with the Gaussian network model (GNM) and compared with experimental data. In more than 70% of the examined enzymes, the global hinge centers predicted by the GNM are found to be colocalized with the catalytic sites experimentally identified. Low translational mobility (< 7%) is observed for the catalytic residues, consistent with the fine-tuned design of enzymes to achieve precise mechanochemical activities. Ligand binding sites, while closely neighboring catalytic sites, enjoy a moderate flexibility to accommodate the ligand binding. These findings could serve as additional criteria for assessing drug binding residues and could lessen the computational burden of substrate docking searches.