Langevin dynamics simulations reveal biologically relevant folds arising from the incorporation of a torsional potential

The conformational behaviour of polymer chains has been examined using Langevin dynamics simulation techniques. Polymer chains were modelled as “beads” undergoing Brownian motion in a defined potential that accounted for stretching, bending and solvation energies. As expected, the competition between chain stiffness and solvent interactions was found to yield standard swollen or collapsed configurations in good or poor solvents, respectively. However, when a torsional term was introduced into the model, additional biologically relevant conformations such as helices, sheets, turns and hairpins naturally arose.     

A Conserved Isoleucine Maintains the Inactive State of Bruton's Tyrosine Kinase

Despite high level of homology among non-receptor tyrosine kinases, different kinase families employ a diverse array of regulatory mechanisms. For example, the catalytic kinase domains of the Tec family kinases are inactive without assembly of the adjacent regulatory domains, whereas the Src kinase domains are autoinhibited by the assembly of similar adjacent regulatory domains. Using molecular dynamics simulations, biochemical assays, and biophysical approaches, we have uncovered an isoleucine residue in the kinase domain of the Tec family member Btk that, when mutated to the closely related leucine, leads to a shift in the conformational equilibrium of the kinase domain toward the active state. The single amino acid mutation results in measureable catalytic activity for the Btk kinase domain in the absence of the regulatory domains. We suggest that this isoleucine side chain in the Tec family kinases acts as a “wedge” that restricts the conformational space available to key regions in the kinase domain, preventing activation until the kinase domain associates with its regulatory subunits and overcomes the energetic barrier to activation imposed by the isoleucine side chain.     

Pressure effects on the ensemble dynamics of ubiquitin inspected with molecular dynamics simulations and isotropic reorientational eigenmode dynamics

According to NMR chemical shift data, the ensemble of ubiquitin is a mixture of “open” and “closed” conformations at rapid equilibrium. Pressure perturbations provide the means to study the transition between the two conformers by imposing an additional constraint on the system's partial molar volume. Here we use nanosecond-timescale molecular dynamics simulations to characterize the network of correlated motions accessible to the conformers at low- and high-pressure conditions. Using the isotropic reorientational eigenmode dynamics formalism to analyze our simulation trajectories, we reproduce NMR relaxation data without fitting any parameters of our model. Comparative analysis of our results suggests that the two conformations behave very differently. The dynamics of the “closed” conformation are almost unaffected by pressure and are dominated by large-amplitude correlated motions of residues 23-34 in the extended α-helix. The “open” conformation under conditions of normal pressure displays increased mobility, focused on the loop residues 17-20, 46-55, and 58-59 at the bottom of the core of the structure, as well as the C-terminal residues 69-76, that directly participate in key protein-protein interactions. For the same conformation, a pressure increase induces a loss of separability between molecular tumbling and internal dynamics, while motions between different backbone sites become uncorrelated.     

Network of long-range concerted chemical shift displacements upon ligand binding to human angiogenin

Molecular recognition models of both induced fit and conformational selection rely on coupled networks of flexible residues and/or structural rearrangements to promote protein function. While the atomic details of these motional events still remain elusive, members of the pancreatic ribonuclease superfamily were previously shown to depend on subtle conformational heterogeneity for optimal catalytic function. Human angiogenin, a structural homologue of bovine pancreatic RNase A, induces blood vessel formation and relies on a weak yet functionally mandatory ribonucleolytic activity to promote neovascularization. Here, we use the NMR chemical shift projection analysis (CHESPA) to clarify the mechanism of ligand binding in human angiogenin, further providing information on long-range intramolecular residue networks potentially involved in the function of this enzyme. We identify two main clusters of residue networks displaying correlated linear chemical shift trajectories upon binding of substrate fragments to the purine- and pyrimidine-specific subsites of the catalytic cleft. A large correlated residue network clusters in the region corresponding to the V₁ domain, a site generally associated with the angiogenic response and structural stability of the enzyme. Another correlated network (residues 40–42) negatively affects the catalytic activity but also increases the angiogenic activity. ¹⁵N-CPMG relaxation dispersion experiments could not reveal the existence of millisecond timescale conformational exchange in this enzyme, a lack of flexibility supported by the very low-binding affinities and catalytic activity of angiogenin. Altogether, the current report potentially highlights the existence of long-range dynamic reorganization of the structure upon distinct subsite binding events in human angiogenin.     

Intricacies of interactions,dynamics and solvent effects in enzyme catalysis: a computational perspective on cyclophilin A

The essential role of enzymes in biological processes has continually ignited sparks of interest in their mechanism of action. Fully understanding the mechanism of enzymes has broad implications in protein engineering and drug design. The more than five order of magnitude speed-up in the rate of peptidyl–prolyl cis–trans isomerisation by cyclophilin A (CypA) has been the target of intense research. CypA serves as a tractable model system, because it reversibly catalyses the rotation around peptidyl–prolyl bonds without any bond breakage or formation. Here, we discuss the results of recent computational approaches used to study the mechanism of CypA. We highlight the critical role of enzyme and substrate conformational dynamics in the developing interactions as the substrate approaches the transition state that results in an astonishing enhancement of isomerisation rate. The rate of isomerisation is affected by the intricate coupling between the dynamics of the substrate, enzyme and solvent. CypA binds its substrates via conformational selection, where rearrangements of key active site residues are necessary for substrate recognition. The conformational plasticity of the active site allows the enzyme to accommodate the most favourable interactions with the transition state that can be exploited for structure-based drug design.     

The crystal structure analysis of group B Streptococcus sortase C1: a model for the "lid" movement upon substrate binding

A unique feature of the class-C-type sortases, enzymes essential for Gram-positive pilus biogenesis, is the presence of a flexible “lid” anchored in the active site. However, the mechanistic details of the “lid” displacement, suggested to be a critical prelude for enzyme catalysis, are not yet known. This is partly due to the absence of enzyme-substrate and enzyme-inhibitor complex crystal structures. We have recently described the crystal structures of the Streptococcus agalactiae SAG2603 V/R sortase SrtC1 in two space groups (type II and type III) and that of its “lid” mutant and proposed a role of the “lid” as a protector of the active-site hydrophobic environment. Here, we report the crystal structures of SAG2603 V/R sortase C1 in a different space group (type I) and that of its complex with a small-molecule cysteine protease inhibitor. We observe that the catalytic Cys residue is covalently linked to the small-molecule inhibitor without lid displacement. However, the type I structure provides a view of the sortase SrtC1 lid displacement while having structural elements similar to a substrate sorting motif suitably positioned in the active site. We propose that these major conformational changes seen in the presence of a substrate mimic in the active site may represent universal features of class C sortase substrate recognition and enzyme activation.     

Guanylate kinase, induced fit, and the allosteric spring probe

Since the introduction of the induced-fit theory by D. E. Koshland Jr., it has been established that conformational motion invariably accompanies the execution of protein function. The catalytic activity of kinases, specifically, is associated with large conformational changes (∼1 nm amplitude). In the case of guanylate kinase, upon substrate binding, the LID and nucleotide-monophosphate-binding domains are brought together and toward the CORE with large concerted movements about the α3 (helix 3) axis. However, whether the change in conformation mostly affects the catalytic rate or mostly increases binding affinities for one or the other substrate is unclear. We investigate this question using a nanotechnology approach based on mechanical stress. Using an “allosteric spring probe”, we bias conformational states in favor of the “open” (substrate-free) conformation of the enzyme; the result is that the binding constant for the substrate guanosine monophosphate (GMP) is reduced by up to a factor of 10, whereas the binding constant for adenosine triphosphate (ATP) and the catalytic rate are essentially unaffected. The results show that the GMP-induced conformational change, which promotes catalysis, does not promote ATP binding, consistent with previous mutagenesis studies. Furthermore, they show that this conformational change is of the induced-fit type with respect to GMP binding (but not ATP binding). We elaborate on this point by proposing a quantitative criterion for the classification of conformational changes with respect to the induced-fit theory. More generally, these results show that the allosteric spring probe can be used to affect enzymatic activity in a continuously controlled manner, and also to affect specific steps of the reaction mechanism while leaving others unaffected. It is presumed that this will enable informative comparisons with the results of future molecular dynamics or statistical mechanics computations.     

Zipping and Unzipping of Adenylate Kinase: Atomistic Insights into the Ensemble of Open ↔ Closed Transitions

Adenylate kinase (AdK), a phosphotransferase enzyme, plays an important role in cellular energy homeostasis. It undergoes a large conformational change between an open and a closed state, even in the absence of substrate. We investigate the apo-AdK transition at the atomic level both with free-energy calculations and with our new dynamic importance sampling (DIMS) molecular dynamics method. DIMS is shown to sample biologically relevant conformations as verified by comparing an ensemble of hundreds of DIMS transitions to AdK crystal structure intermediates. The simulations reveal in atomic detail how hinge regions partially and intermittently unfold during the transition. Conserved salt bridges are seen to have important structural and dynamic roles; in particular, four ionic bonds that open in a sequential, zipper-like fashion and, thus, dominate the free-energy landscape of the transition are identified. Transitions between the closed and open conformations only have to overcome moderate free-energy barriers. Unexpectedly, the closed state and the open state encompass broad free-energy basins that contain conformations differing in domain hinge motions by up to 40°. The significance of these extended states is discussed in relation to recent experimental Förster resonance energy transfer measurements. Taken together, these results demonstrate how a small number of cooperative key interactions can shape the overall dynamics of an enzyme and suggest an “all-or-nothing” mechanism for the opening and closing of AdK. Our efficient DIMS molecular dynamics computer simulation approach can provide a detailed picture of a functionally important macromolecular transition and thus help to interpret and suggest experiments to probe the conformational landscape of dynamic proteins such as AdK.     

Effect of molecular confinement on internal enzyme dynamics: frequency domain fluorometry and molecular dynamics simulation studies

The tryptophanyl emission decay of the mesophilic beta-galactosidase from Aspergillus oryzae free in buffer and entrapped in agarose gel is investigated as a function of temperature and compared to that of the hyperthermophilic enzyme from Sulfolobus solfataricus. Both enzymes are tetrameric proteins with a large number of tryptophanyl residues, so the fluorescence emission can provide information on the conformational dynamics of the overall protein structure rather than that of the local environment. The tryptophanyl emission decays are best fitted by bimodal Lorentzian distributions. The long-lived component is ascribed to close, deeply buried tryptophanyl residues with reduced mobility; the short-lived one arises from tryptophanyl residues located in more flexible external regions of each subunit, some of which are involved in forming the catalytic site. The center of both lifetime distribution components at each temperature increases when going from the free in solution mesophilic enzyme to the gel-entrapped and hyperthermophilic enzyme, thus indicating that confinement of the mesophilic enzyme in the agarose gel limits the freedom of the polypeptide chain. A more complex dependence is observed for the distribution widths. Computer modeling techniques are used to recognize that the catalytic sites are similar for the mesophilic and hyperthermophilic beta-galactosidases. The effect due to gel entrapment is considered in dynamic simulations by imposing harmonic restraints to solvent-exposed atoms of the protein with the exclusion of those around the active site. The temperature dependence of the tryptophanyl fluorescence emission decay and the dynamic simulation confirm that more rigid structures, as in the case of the immobilized and/or hyperthermophilic enzyme, require higher temperatures to achieve the requisite conformational dynamics for an effective catalytic action and strongly suggest a link between conformational rigidity and enhanced thermal stability.     

Evolutionary variation and adaptation in a conserved protein kinase allosteric network: Implications for inhibitor design

The activation of protein kinases involves conformational changes in key functional regions of the kinase domain, a detailed understanding of which is essential for the design of selective protein kinase inhibitors. Through statistical analysis of protein kinase sequences and crystal structures from diverse organisms, we recently proposed that the activation of protein kinases involves a hidden strain switch in the catalytic loop. Specifically, we demonstrated that the backbone torsion-angles of residues in the catalytic loop switch from a “relaxed” to “strained” conformation upon kinase activation and the strained geometry results in a network of hydrogen bonds involving conserved non-catalytic residues in the ATP and substrate binding lobes. Here, we further explore this activation mechanism by analyzing families that lack the canonical hydrogen bonding interactions with the strained backbone. We find that alternative mechanisms have evolved to maintain catalytic loop strain. In PIM kinase, for example, two water molecules account for the lack of a conserved aspartate in the substrate binding by hydrogen bonds to the strained backbone. We discuss the relevance of these findings in the design of family-specific allosteric inhibitors, and in predicting the structural and functional impact of cancer mutations that alter the strain associated hydrogen bonding network. This article is part of a Special Issue entitled: Inhibitors of Protein Kinases (2012).     

Src kinase activation: A switched electrostatic network

Src tyrosine kinases are essential in numerous cell signaling pathways, and improper functioning of these enzymes has been implicated in many diseases. The activity of Src kinases is regulated by conformational activation, which involves several structural changes within the catalytic domain (CD): the orientation of two lobes of CD; rearrangement of the activation loop (A-loop); and movement of an alpha-helix (alphaC), which is located at the interface between the two lobes, into or away from the catalytic cleft. Conformational activation was investigated using biased molecular dynamics to explore the transition pathway between the active and the down-regulated conformation of CD for the Src-kinase family member Lyn kinase, and to gain insight into the interdependence of these changes. Lobe opening is observed to be a facile motion, whereas movement of the A-loop motion is more complex requiring secondary structure changes as well as communication with alphaC. A key result is that the conformational transition involves a switch in an electrostatic network of six polar residues between the active and the down-regulated conformations. The exchange between interactions links the three main motions of the CD. Kinetic experiments that would demonstrate the contribution of the switched electrostatic network to the enzyme mechanism are proposed. Possible implications for regulation conferred by interdomain interactions are also discussed.     

A new clustering of antibody CDR loop conformations

Previous analyses of the complementarity-determining regions (CDRs) of antibodies have focused on a small number of “canonical” conformations for each loop. This is primarily the result of the work of Chothia and coworkers, most recently in 1997. Because of the widespread utility of antibodies, we have revisited the clustering of conformations of the six CDR loops with the much larger amount of structural information currently available. In this work, we were careful to use a high-quality data set by eliminating low-resolution structures and CDRs with high B-factors or high conformational energies. We used a distance function based on directional statistics and an effective clustering algorithm with affinity propagation. With this data set of over 300 nonredundant antibody structures, we were able to cover 28 CDR-length combinations (e.g., L1 length 11, or “L1-11” in our CDR-length nomenclature) for L1, L2, L3, H1, and H2. The Chothia analysis covered only 20 CDR-lengths. Only four of these had more than one conformational cluster, of which two could easily be distinguished by gene source (mouse/human; κ/λ) and one could easily be distinguished purely by the presence and the positions of Pro residues (L3-9). Thus, using the Chothia analysis does not require the complicated set of “structure-determining residues” that is often assumed. Of our 28 CDR-lengths, 15 have multiple conformational clusters, including 10 for which the Chothia analysis had only one canonical class. We have a total of 72 clusters for non-H3 CDRs; approximately 85% of the non-H3 sequences can be assigned to a conformational cluster based on gene source and/or sequence. We found that earlier predictions of “bulged” versus “nonbulged” conformations based on the presence or the absence of anchor residues Arg/Lys94 and Asp101 of H3 have not held up, since all four combinations lead to a majority of conformations that are bulged. Thus, the earlier analyses have been significantly enhanced by the increased data. We believe that the new classification will lead to improved methods for antibody structure prediction and design.     

Mapping the mutual information network of enzymatic families in the protein structure to unveil functional features

Amino acids committed to a particular function correlate tightly along evolution and tend to form clusters in the 3D structure of the protein. Consequently, a protein can be seen as a network of co-evolving clusters of residues. The goal of this work is two-fold: first, we have combined mutual information and structural data to describe the amino acid networks within a protein and their interactions. Second, we have investigated how this information can be used to improve methods of prediction of functional residues by reducing the search space. As a main result, we found that clusters of co-evolving residues related to the catalytic site of an enzyme have distinguishable topological properties in the network. We also observed that these clusters usually evolve independently, which could be related to a fail-safe mechanism. Finally, we discovered a significant enrichment of functional residues (e.g. metal binding, susceptibility to detrimental mutations) in the clusters, which could be the foundation of new prediction tools.     

Contribution of buried distal amino acid residues in horse liver alcohol dehydrogenase to structure and catalysis

The dynamics of enzyme catalysis range from the slow time scale (～ms) for substrate binding and conformational changes to the fast time (～ps) scale for reorganization of substrates in the chemical step. The contribution of global dynamics to catalysis by alcohol dehydrogenase was tested by substituting five different, conserved amino acid residues that are distal from the active site and located in the hinge region for the conformational change or in hydrophobic clusters. X‐ray crystallography shows that the structures for the G173A, V197I, I220 (V, L, or F), V222I, and F322L enzymes complexed with NAD⁺ and an analogue of benzyl alcohol are almost identical, except for small perturbations at the sites of substitution. The enzymes have very similar kinetic constants for the oxidation of benzyl alcohol and reduction of benzaldehyde as compared to the wild‐type enzyme, and the rates of conformational changes are not altered. Less conservative substitutions of these amino acid residues, such as G173(V, E, K, or R), V197(G, S, or T), I220(G, S, T, or N), and V222(G, S, or T) produced unstable or poorly expressed proteins, indicating that the residues are critical for global stability. The enzyme scaffold accommodates conservative substitutions of distal residues, and there is no evidence that fast, global dynamics significantly affect the rate constants for hydride transfers. In contrast, other studies show that proximal residues significantly participate in catalysis.     

Thermodynamic analysis shows conformational coupling and dynamics confer substrate specificity in fructose-1,6-bisphosphate aldolase

Conformational flexibility is emerging as a central theme in enzyme catalysis. Thus, identifying and characterizing enzyme dynamics are critical for understanding catalytic mechanisms. Herein, coupling analysis, which uses thermodynamic analysis to assess cooperativity and coupling between distal regions on an enzyme, is used to interrogate substrate specificity among fructose-1,6-(bis)phosphate aldolase (aldolase) isozymes. Aldolase exists as three isozymes, A, B, and C, distinguished by their unique substrate preferences despite the fact that the structures of the active sites of the three isozymes are nearly identical. While conformational flexibility has been observed in aldolase A, its function in the catalytic reaction of aldolase has not been demonstrated. To explore the role of conformational dynamics in substrate specificity, those residues associated with isozyme specificity (ISRs) were swapped and the resulting chimeras were subjected to steady-state kinetics. Thermodynamic analyses suggest cooperativity between a terminal surface patch (TSP) and a distal surface patch (DSP) of ISRs that are separated by >8.9 A. Notably, the coupling energy (DeltaGI) is anticorrelated with respect to the two substrates, fructose 1,6-bisphosphate and fructose 1-phosphate. The difference in coupling energy with respect to these two substrates accounts for approximately 70% of the energy difference for the ratio of kcat/Km for the two substrates between aldolase A and aldolase B. These nonadditive mutational effects between the TSP and DSP provide functional evidence that coupling interactions arising from conformational flexibility during catalysis are a major determinant of substrate specificity.     

Functional and structural dynamics of NhaA,a prototype for Na and H antiporters,which are responsible for Na and H homeostasis in cells

The crystal structure of down-regulated NhaA crystallized at acidic pH 4 [21] has provided the first structural insights into the antiport mechanism and pH regulation of a Na⁺/H⁺ antiporter [22]. On the basis of the NhaA crystal structure [21] and experimental data (reviewed in [2,22,38] we have suggested that NhaA is organized into two functional regions: (i) a cluster of amino acids responsible for pH regulation (ii) a catalytic region at the middle of the TM IV/XI assembly, with its unique antiparallel unfolded regions that cross each other forming a delicate electrostatic balance in the middle of the membrane. This unique structure contributes to the cation binding site and allows the rapid conformational changes expected for NhaA. Extended chains interrupting helices appear now a common feature for ion binding in transporters. However the NhaA fold is unique and shared by ASBTNM [30] and NapA [29]. Computation [13], electrophysiology [69] combined with biochemistry [33,47] have provided intriguing models for the mechanism of NhaA. However, the conformational changes and the residues involved have not yet been fully identified. Another issue which is still enigma is how energy is transduced “in this ‘nano-machine.’” We expect that an integrative approach will reveal the residues that are crucial for NhaA activity and regulation, as well as elucidate the pHand ligand-induced conformational changes and their dynamics. Ultimately, integrative results will shed light on the mechanism of activity and pH regulation of NhaA, a prototype of the CPA2 family of transporters. This article is part of a Special Issue entitled: 18th European Bioenergetic Conference.     

Altered dynamics upon oligomerization corresponds to key functional sites

It is known that over half of the proteins encoded by most organisms function as oligomeric complexes. Oligomerization confers structural stability and dynamics changes in proteins. We investigate the effects of oligomerization on protein dynamics and its functional significance for a set of 145 multimeric proteins. Using coarse‐grained elastic network models, we inspect the changes in residue fluctuations upon oligomerization and then compare with residue conservation scores to identify the functional significance of these changes. Our study reveals conservation of about ½ of the fluctuations, with ¼ of the residues increasing in their mobilities and ¼ having reduced fluctuations. The residues with dampened fluctuations are evolutionarily more conserved and can serve as orthosteric binding sites, indicating their importance. We also use triosephosphate isomerase as a test case to understand why certain enzymes function only in their oligomeric forms despite the monomer including all required catalytic residues. To this end, we compare the residue communities (groups of residues which are highly correlated in their fluctuations) in the monomeric and dimeric forms of the enzyme. We observe significant changes to the dynamical community architecture of the catalytic core of this enzyme. This relates to its functional mechanism and is seen only in the oligomeric form of the protein, answering why proteins are oligomeric structures. Proteins 2017; 85:1422–1434. © 2017 Wiley Periodicals, Inc.     

Emergence of novel enzyme quasi-species depends on the substrate matrix

Current research on enzyme evolution has shown that many enzymes are promiscuous and have activities with alternative substrates. Mutagenesis tends to relax substrate selectivity, and evolving enzymes can be regarded (summed over evolutionary time) as clusters of enzyme variants, or “quasi-species,” tested against a “substrate matrix” defined by all chemical substances to which the evolvants are exposed.In this investigation, the importance of the substrate matrix for identification of evolvable clusters of enzymes was evaluated by random sampling of variants from a library of glutathione transferase (GST) mutants. The variant GSTs were created by DNA shuffling of homologous Alpha class sequences. The substrate matrix was an array of alternative substrates used under defined experimental conditions. The measured enzyme activities produced a rectangular matrix, in which the rows can be projected as enzyme vectors in substrate-activity space and, reciprocally, the columns can be projected as alternative substrate vectors in enzyme-activity space. Multivariate analysis of the catalytic activities demonstrated that the enzyme vectors formed two primary clusters or functional “molecular quasi-species.” These quasi-species serve as the raw material from which more specialized enzymes eventually could evolve. The substrate vectors similarly formed two major groups. Identification of separate quasi-species of GSTs in a mutant library was critically dependent on the nature of the substrate matrix. When substrates from just one of the two groups were used, only one cluster of enzymes could be recognized. On the other hand, expansion of the substrate matrix to include additional substrates showed the presence of a third quasi-species among the GST variants already analyzed. Thus, the portrayal of the functional quasi-species is intimately linked to the effective substrate matrix. In natural evolution, the substrates actually encountered therefore play a pivotal role in determining whether latent catalytic abilities become manifest in novel enzymes.     

Structural and stability effects of phosphorylation: Localized structural changes in phenylalanine hydroxylase

Phosphorylation of phenylalanine hydroxylase (PAH) at Ser16 by cAMP-dependent protein kinase increases the basal activity of the enzyme and its resistance to tryptic proteolysis. The modeled structures of the full-length phosphorylated and unphosphorylated enzyme were subjected to molecular dynamics simulations, and we analyzed the energy of charge-charge interactions for individual ionizable residues in the final structures. These calculations showed that the conformational changes induced by incorporation of phosphate were localized and limited mostly to the region around the phosphoserine (Arg13-Asp17) and a region around the active site in the catalytic domain that includes residues involved in the binding of the iron and the substrate L-Phe (Arg270 and His285). The absence of a generalized conformational change was confirmed by differential scanning calorimetry, thermal-dependent circular dichroism, fluorescence spectroscopy, and limited chymotryptic proteolysis of the phosphorylated and unphosphorylated PAH. Our results explain the effect of phosphorylation of PAH on both the resistance to proteolysis specifically by trypsin-like enzymes and on the increase in catalytic efficiency.