Enzymes: An integrated view of structure, dynamics and function

Microbes utilize enzymes to perform a variety of functions. Enzymes are biocatalysts working as highly efficient machines at the molecular level. In the past, enzymes have been viewed as static entities and their function has been explained on the basis of direct structural interactions between the enzyme and the substrate. A variety of experimental and computational techniques, however, continue to reveal that proteins are dynamically active machines, with various parts exhibiting internal motions at a wide range of time-scales. Increasing evidence also indicates that these internal protein motions play a role in promoting protein function such as enzyme catalysis. Moreover, the thermodynamical fluctuations of the solvent, surrounding the protein, have an impact on internal protein motions and, therefore, on enzyme function. In this review, we describe recent biochemical and theoretical investigations of internal protein dynamics linked to enzyme catalysis. In the enzyme cyclophilin A, investigations have lead to the discovery of a network of protein vibrations promoting catalysis. Cyclophilin A catalyzes peptidyl-prolyl cis/trans isomerization in a variety of peptide and protein substrates. Recent studies of cyclophilin A are discussed in detail and other enzymes (dihydrofolate reductase and liver alcohol dehydrogenase) where similar discoveries have been reported are also briefly discussed. The detailed characterization of the discovered networks indicates that protein dynamics plays a role in rate-enhancement achieved by enzymes. An integrated view of enzyme structure, dynamics and function have wide implications in understanding allosteric and co-operative effects, as well as protein engineering of more efficient enzymes and novel drug design

Synthetic, Switchable Enzymes

The construction of switchable, radiation-controlled, aptameric enzymes – “swenzymes” – is, in principle, feasible. We propose a strategy to make such catalysts from 2 (or more) aptamers each selected to bind specifically to one of the substrates in, for example, a 2-substrate reaction. Construction of a combinatorial library of candidate swenzymes entails selecting a set of a million aptamers that bind one substrate and a second set of a million aptamers that bind the second substrate; the aptamers in these sets are then linked pairwise by a linker, thus bringing together the substrates. In the presence of the substrates, some linked aptamer pairs catalyze the reaction when exposed to external energy in the form of a specific frequency of low intensity, nonionizing electromagnetic or acoustic radiation. Such swenzymes are detected via a separate product-capturing aptamer that changes conformation on capturing the product; this altered conformation allows it (1) to bind to every potential swenzyme in its vicinity (thereby giving a higher probability of capture to the swenzymes that generate the product) and (2) to bind to a sequence on a magnetic bead (thereby permitting purification of the swenzyme plus product-capturing aptamer by precipitation). Attempts to implement the swenzyme strategy may help elucidate fundamental problems in enzyme catalysis

Nucleotide substrate binding characterization in human pancreatic-type ribonucleases

Human genome contains a group of more than a dozen similar genes with diverse biological functions including antiviral, antibacterial and angiogenesis activities. The characterized gene products of this group show significant sequence similarity and a common structural fold associated with binding and cleavage of ribonucleic acid (RNA) substrates. Therefore, these proteins have been categorized as members of human pancreatic-type ribonucleases (hRNases). hRNases differ in cell/tissue localization and display distinct substrate binding preferences and a wide range of ribonucleolytic catalytic efficiencies. Limited information is available about structural and dynamical properties that influence this diversity among these homologous RNases. Here, we use computer simulations to characterize substrate interactions, electrostatics and dynamical properties of hRNases 1–7 associated with binding to two nucleotide substrates (ACAC and AUAU). Results indicate that even with complete conservation of active-site catalytic triad associated with ribonucleolytic activity, these enzymes show significant differences in substrate interactions. Detailed characterization suggests that in addition to binding site electrostatic and van der Waals interactions, dynamics of distal regions may also play a role in binding. Another key insight is that a small difference in temperature of 300 K (used in experimental studies) and 310 K (physiological temperature) shows significant changes in enzyme-substrate interactions

Future COVID19 surges prediction based on SARS-CoV-2 mutations surveillance

COVID19 has aptly revealed that airborne viruses such as SARS-CoV-2 with the ability to rapidly mutate, combined with high rates of transmission and fatality can cause a deadly world-wide pandemic in a matter of weeks.1 Apart from vaccines and post-infection treatment options, strategies for preparedness will be vital in responding to the current and future pandemics. Therefore, there is wide interest in approaches that allow predictions of increase in infections ('surges') before they occur. We describe here real time genomic surveillance particularly based on mutation analysis, of viral proteins as a methodology for a priori determination of surge in number of infection cases. The full results are available for SARS-CoV-2 at http://pandemics.okstate.edu/covid19/, and are updated daily as new virus sequences become available. This approach is generic and will also be applicable to other pathogens

Discovering Conformational Sub-States Relevant to Protein Function

Background: Internal motions enable proteins to explore a range of conformations, even in the vicinity of native state. The role of conformational fluctuations in the designated function of a protein is widely debated. Emerging evidence suggests that sub-groups within the range of conformations (or sub-states) contain properties that may be functionally relevant. However, low populations in these sub-states and the transient nature of conformational transitions between these substates present significant challenges for their identification and characterization. Methods and Findings: To overcome these challenges we have developed a new computational technique, quasianharmonic analysis (QAA). QAA utilizes higher-order statistics of protein motions to identify sub-states in the conformational landscape. Further, the focus on anharmonicity allows identification of conformational fluctuations that enable transitions between sub-states. QAA applied to equilibrium simulations of human ubiquitin and T4 lysozyme reveals functionally relevant sub-states and protein motions involved in molecular recognition. In combination with a reaction pathway sampling method, QAA characterizes conformational sub-states associated with cis/trans peptidyl-prolyl isomerization catalyzed by the enzyme cyclophilin A. In these three proteins, QAA allows identification of conformational sub-states, with critical structural and dynamical features relevant to protein function. Conclusions: Overall, QAA provides a novel framework to intuitively understand the biophysical basis of conformational diversity and its relevance to protein function. © 2011 Ramanathan et al

Evolutionarily Conserved Linkage between Enzyme Fold, Flexibility, and Catalysis

Proteins are intrinsically flexible molecules. The role of internal motions in a protein's designated function is widely debated. The role of protein structure in enzyme catalysis is well established, and conservation of structural features provides vital clues to their role in function. Recently, it has been proposed that the protein function may involve multiple conformations: the observed deviations are not random thermodynamic fluctuations; rather, flexibility may be closely linked to protein function, including enzyme catalysis. We hypothesize that the argument of conservation of important structural features can also be extended to identification of protein flexibility in interconnection with enzyme function. Three classes of enzymes (prolyl-peptidyl isomerase, oxidoreductase, and nuclease) that catalyze diverse chemical reactions have been examined using detailed computational modeling. For each class, the identification and characterization of the internal protein motions coupled to the chemical step in enzyme mechanisms in multiple species show identical enzyme conformational fluctuations. In addition to the active-site residues, motions of protein surface loop regions (>10 Å away) are observed to be identical across species, and networks of conserved interactions/residues connect these highly flexible surface regions to the active-site residues that make direct contact with substrates. More interestingly, examination of reaction-coupled motions in non-homologous enzyme systems (with no structural or sequence similarity) that catalyze the same biochemical reaction shows motions that induce remarkably similar changes in the enzyme–substrate interactions during catalysis. The results indicate that the reaction-coupled flexibility is a conserved aspect of the enzyme molecular architecture. Protein motions in distal areas of homologous and non-homologous enzyme systems mediate similar changes in the active-site enzyme–substrate interactions, thereby impacting the mechanism of catalyzed chemistry. These results have implications for understanding the mechanism of allostery, and for protein engineering and drug design

Unfolding Simulations of Holomyoglobin from Four Mammals: Identification of Intermediates and β-Sheet Formation from Partially Unfolded States

Myoglobin (Mb) is a centrally important, widely studied mammalian protein. While much work has investigated multi-step unfolding of apoMb using acid or denaturant, holomyoglobin unfolding is poorly understood despite its biological relevance. We present here the first systematic unfolding simulations of holoMb and the first comparative study of unfolding of protein orthologs from different species (sperm whale, pig, horse, and harbor seal). We also provide new interpretations of experimental mean molecular ellipticities of myoglobin intermediates, notably correcting for random coil and number of helices in intermediates. The simulated holoproteins at 310 K displayed structures and dynamics in agreement with crystal structures (R g ~1.48-1.51 nm, helicity ~75%). At 400 K, heme was not lost, but some helix loss was observed in pig and horse, suggesting that these helices are less stable in terrestrial species. At 500 K, heme was lost within 1.0-3.7 ns. All four proteins displayed exponentially decaying helix structure within 20 ns. The C- and F-helices were lost quickly in all cases. Heme delayed helix loss, and sperm whale myoglobin exhibited highest retention of heme and D/E helices. Persistence of conformation (RMSD), secondary structure, and ellipticity between 2-11 ns was interpreted as intermediates of holoMb unfolding in all four species. The intermediates resemble those of apoMb notably in A and H helices, but differ substantially in the D-, E- and F-helices, which interact with heme. The identified mechanisms cast light on the role of metal/cofactor in poorly understood holoMb unfolding. We also observed β-sheet formation of several myoglobins at 500 K as seen experimentally, occurring after disruption of helices to a partially unfolded, globally disordered state; heme reduced this tendency and sperm-whale did not display any sheet propensity during the simulations