Probing the global and local dynamics of aminoacyl-tRNA synthetases using all-atom and coarse-grained simulations

Insights into the Mechanism of Tryptophan Fluorescence Quenching due to Synthetic Crowding Agents: A Combined Experimental and Computational Study

Intrinsic tryptophan fluorescence spectroscopy is an important tool for examining the effects of molecular crowding and confinement on the structure, dynamics, and function of proteins. Synthetic crowders such as dextran, ficoll, polyethylene glycols, polyvinylpyrrolidone, and their respective monomers are used to mimic crowded intracellular environments. Interactions of these synthetic crowders with tryptophan and the subsequent impact on its fluorescence properties are therefore critically important for understanding the possible interference created by these crowders. In the present study, the effects of polymer and monomer crowders on tryptophan fluorescence were assessed by using experimental and computational approaches. The results of this study demonstrated that both polymer and monomer crowders have an impact on the tryptophan fluorescence intensity; however, the molecular mechanisms of quenching were different. Using Stern-Volmer plots and a temperature variation study, a physical basis for the quenching mechanism of commonly used synthetic crowders was established. The quenching of free tryptophan was found to involve static, dynamic, and sphere-of-action mechanisms. In parallel, computational studies employing Kohn-Sham density functional theory provided a deeper insight into the effects of intermolecular interactions and solvation, resulting in differing quenching modes for these crowders. Taken together, the study offers new physical insights into the quenching mechanisms of some commonly used monomer and polymer synthetic crowders.

Exploring the intrinsic dynamics of SARS-CoV-2, SARS-CoV and MERS-CoV spike glycoprotein through normal mode analysis using anisotropic network model

COVID-19 caused by SARS-CoV-2 have become a global pandemic with serious rate of fatalities. SARS-CoV and MERS-CoV have also caused serious outbreak previously but the intensity was much lower than the ongoing SARS-CoV-2. The main infectivity factor of all the three viruses is the spike glycoprotein. In this study we have examined the intrinsic dynamics of the prefusion, lying state of trimeric S protein of these viruses through Normal Mode Analysis using Anisotropic Network Model. The dynamic modes of the S proteins of the aforementioned viruses were compared by root mean square inner product (RMSIP), spectral overlap and cosine correlation matrix. S proteins show homogenous correlated or anticorrelated motions among their domains but direction of Cα atom among the spike proteins show less similarity. SARS-CoV-2 spike shows high vertically upward motion of the receptor binding motif implying its propensity for binding with the receptor even in the lying state. MERS-CoV spike shows unique dynamical motion compared to the other two S protein indicated by low RMSIP, spectral overlap and cosine correlation value. This study will guide in developing common potential inhibitor molecules against closed state of spike protein of these viruses to prevent conformational switching from lying to standing state.

Editing Domain Motions Preorganize the Synthetic Active Site of Prolyl-tRNA Synthetase

Prolyl-tRNA synthetases (ProRSs) catalyze the covalent attachment of proline onto cognate tRNAs, an indispensable step for protein synthesis in all living organisms. ProRSs are modular enzymes and the "prokaryotic-like" ProRSs are distinguished from "eukaryotic-like" ProRSs by the presence of an editing domain (INS) inserted between motifs 2 and 3 of the main catalytic domain. Earlier studies suggested the presence of coupled-domain dynamics could contribute to catalysis; however, the role that the distal, highly mobile INS domain plays in catalysis at the synthetic active site is not completely understood. In the present study, a combination of theoretical and experimental approaches has been used to elucidate the precise role of INS domain dynamics. Quantum mechanical/molecular mechanical simulations were carried out to model catalytic Pro-AMP formation by Enterococcus faecalis ProRS. The energetics of the adenylate formation by the wild-type enzyme was computed and contrasted with variants containing active site mutations, as well as a deletion mutant lacking the INS domain. The combined results revealed that two distinct types of dynamics contribute to the enzyme's catalytic power. One set of motions is intrinsic to the INS domain and leads to conformational preorganization that is essential for catalysis. A second type of motion, stemming from the electrostatic reorganization of active site residues, impacts the height and width of the energy profile and has a critical role in fine tuning the substrate orientation to facilitate reactive collisions. Thus, motions in a distal domain can preorganize the active site of an enzyme to optimize catalysis.

Effects of Distal Mutations on Prolyl-Adenylate Formation of Escherichia coli Prolyl-tRNA Synthetase

Enzymes play important roles in many biological processes. Amino acid residues in the active site pocket of an enzyme, which are in direct contact with the substrate(s), are generally believed to be critical for substrate recognition and catalysis. Identifying and understanding how these "catalytic" residues help enzymes achieve enormous rate enhancement has been the focus of many structural and biochemical studies over the past several decades. Recent studies have shown that enzymes are intrinsically dynamic and dynamic coupling between distant structural elements is essential for effective catalysis in modular enzymes. Therefore, distal residues are expected to have impact on enzyme function. However, few studies have investigated the role of distal residues on enzymatic catalysis. In the present study, the effects of distal residue mutations on the catalytic function of an aminoacyl-tRNA synthetase, namely, prolyl-tRNA synthase, were investigated. The present study demonstrates that distal residues significantly contribute to catalysis of the modular Escherichia coli prolyl-tRNA synthetase by maintaining intrinsic protein flexibility.

Role of disordered regions in transferring tyrosine to its cognate tRNA

Aminoacyl tRNA synthetase (AARS) plays an important role in transferring each amino acid to its cognate tRNA. Specifically, tyrosyl tRNA synthetase (TyrRS) is involved in various functions including protection from DNA damage due to oxidative stress, protein synthesis and cell signaling and can be an attractive target for controlling the pathogens by early inhibition of translation. TyrRS has two disordered regions, which lack a stable 3D structure in solution, and are involved in tRNA synthetase catalysis and stability. One of the disordered regions undergoes disorder-to-order transition (DOT) upon complex formation with tRNA whereas the other remains disordered (DR). In this work, we have explored the importance of these disordered regions using molecular dynamics simulations of both free and RNA-complexed states. We observed that the DOT and DR regions of the first subunit acts as a flap and interact with the acceptor arm of the tRNA. The DOT-DR flap closes when tyrosine (TyrRS^Tyr) is present at the active site of the complex and opens in the presence of tyrosine monophosphate (TyrRS^YMP). The DOT and DR regions of the second subunit interact with the anticodon stem as well as D-loop of the tRNA, which might be involved in stabilizing the complex. The anticodon loop of the tRNA binds to the structured region present in the C-terminal of the protein, which is observed to be flexible during simulations. Detailed energy calculations also show that TyrRS^Tyr complex has stronger binding energy between tRNA and protein compared to TyrRS^YMP; on the contrary, the anticodon is strongly bound in TyrRS^YMP. The results obtained in the present study provide additional insights for understanding catalysis and the involvement of disordered regions in Tyr transfer to cognate tRNA.

Computational modeling and molecular dynamics simulations of mammalian cytoplasmic tyrosyl-tRNA synthetase and its complexes with substrates

Cytoplasmic tyrosyl-tRNA synthetase (TyrRS) is one of the key enzymes of protein biosynthesis. TyrRSs of pathogenic organisms have gained attention as potential targets for drug development. Identifying structural differences between various TyrRSs will facilitate the development of specific inhibitors for the TyrRSs of pathogenic organisms. However, there is a deficiency in structural data for mammalian cytoplasmic TyrRS in complexes with substrates. In this work, we constructed spatial structure of full-length Bos taurus TyrRS (BtTyrRS) and its complexes with substrates using the set of computational modeling techniques. Special attention was paid to BtTyrRS complexes with substrates [L-tyrosine, K⁺ and ATP:Mg²⁺] and intermediate products [tyrosyl-adenylate (Tyr-AMP), K⁺ and PP_i:Mg²⁺] with the different catalytic loop conformations. In order to analyze their dynamical properties, we performed 100 ns of molecular dynamics (MD) simulations. MD simulations revealed new structural data concerning the tyrosine activation reaction in mammalian TyrRS. Formation of strong interaction between Lys154 and γ-phosphate suggests the additional role of CP1 insertion as an important factor for ATP binding. The presence of a potassium-binding pocket within the active site of mammalian TyrRS compensates the absence of the second lysine in the KMSKS motif. Our data provide new details concerning a role of K⁺ ions at different stages of the first step of the tyrosylation reaction, including the coordination of substrates and involvement in the PP_i releasing. The results of this work suggest that differences between ATP-binding sites of mammalian and bacterial TyrRSs are meaningful and could be exploited in the drug design.

Aminoacyl-tRNA Synthetases in the Bacterial World

Aminoacyl-tRNA synthetases (aaRSs) are modular enzymes globally conserved in the three kingdoms of life. All catalyze the same two-step reaction, i.e., the attachment of a proteinogenic amino acid on their cognate tRNAs, thereby mediating the correct expression of the genetic code. In addition, some aaRSs acquired other functions beyond this key role in translation. Genomics and X-ray crystallography have revealed great structural diversity in aaRSs (e.g., in oligomery and modularity, in ranking into two distinct groups each subdivided in 3 subgroups, by additional domains appended on the catalytic modules). AaRSs show huge structural plasticity related to function and limited idiosyncrasies that are kingdom or even species specific (e.g., the presence in many Bacteria of non discriminating aaRSs compensating for the absence of one or two specific aaRSs, notably AsnRS and/or GlnRS). Diversity, as well, occurs in the mechanisms of aaRS gene regulation that are not conserved in evolution, notably between distant groups such as Gram-positive and Gram-negative Bacteria. The review focuses on bacterial aaRSs (and their paralogs) and covers their structure, function, regulation, and evolution. Structure/function relationships are emphasized, notably the enzymology of tRNA aminoacylation and the editing mechanisms for correction of activation and charging errors. The huge amount of genomic and structural data that accumulated in last two decades is reviewed, showing how the field moved from essentially reductionist biology towards more global and integrated approaches. Likewise, the alternative functions of aaRSs and those of aaRS paralogs (e.g., during cell wall biogenesis and other metabolic processes in or outside protein synthesis) are reviewed. Since aaRS phylogenies present promiscuous bacterial, archaeal, and eukaryal features, similarities and differences in the properties of aaRSs from the three kingdoms of life are pinpointed throughout the review and distinctive characteristics of bacterium-like synthetases from organelles are outlined.

Incorporating modeling and simulations in undergraduate biophysical chemistry course to promote understanding of structure-dynamics-function relationships in proteins

A project-based biophysical chemistry laboratory course, which is offered to the biochemistry and molecular biology majors in their senior year, is described. In this course, the classroom study of the structure-function of biomolecules is integrated with the discovery-guided laboratory study of these molecules using computer modeling and simulations. In particular, modern computational tools are employed to elucidate the relationship between structure, dynamics, and function in proteins. Computer-based laboratory protocols that we introduced in three modules allow students to visualize the secondary, super-secondary, and tertiary structures of proteins, analyze non-covalent interactions in protein-ligand complexes, develop three-dimensional structural models (homology model) for new protein sequences and evaluate their structural qualities, and study proteins' intrinsic dynamics to understand their functions. In the fourth module, students are assigned to an authentic research problem, where they apply their laboratory skills (acquired in modules 1-3) to answer conceptual biophysical questions. Through this process, students gain in-depth understanding of protein dynamics-the missing link between structure and function. Additionally, the requirement of term papers sharpens students' writing and communication skills. Finally, these projects result in new findings that are communicated in peer-reviewed journals.

Comparison of intrinsic dynamics of cytochrome p450 proteins using normal mode analysis

Cytochrome P450 enzymes are hemeproteins that catalyze the monooxygenation of a wide-range of structurally diverse substrates of endogenous and exogenous origin. These heme monooxygenases receive electrons from NADH/NADPH via electron transfer proteins. The cytochrome P450 enzymes, which constitute a diverse superfamily of more than 8,700 proteins, share a common tertiary fold but < 25% sequence identity. Based on their electron transfer protein partner, cytochrome P450 proteins are classified into six broad classes. Traditional methods of pro are based on the canonical paradigm that attributes proteins' function to their three-dimensional structure, which is determined by their primary structure that is the amino acid sequence. It is increasingly recognized that protein dynamics play an important role in molecular recognition and catalytic activity. As the mobility of a protein is an intrinsic property that is encrypted in its primary structure, we examined if different classes of cytochrome P450 enzymes display any unique patterns of intrinsic mobility. Normal mode analysis was performed to characterize the intrinsic dynamics of five classes of cytochrome P450 proteins. The present study revealed that cytochrome P450 enzymes share a strong dynamic similarity (root mean squared inner product > 55% and Bhattacharyya coefficient > 80%), despite the low sequence identity (< 25%) and sequence similarity (< 50%) across the cytochrome P450 superfamily. Noticeable differences in Cα atom fluctuations of structural elements responsible for substrate binding were noticed. These differences in residue fluctuations might be crucial for substrate selectivity in these enzymes.

MD Simulations of tRNA and Aminoacyl-tRNA Synthetases: Dynamics, Folding, Binding, and Allostery

While tRNA and aminoacyl-tRNA synthetases are classes of biomolecules that have been extensively studied for decades, the finer details of how they carry out their fundamental biological functions in protein synthesis remain a challenge. Recent molecular dynamics (MD) simulations are verifying experimental observations and providing new insight that cannot be addressed from experiments alone. Throughout the review, we briefly discuss important historical events to provide a context for how far the field has progressed over the past few decades. We then review the background of tRNA molecules, aminoacyl-tRNA synthetases, and current state of the art MD simulation techniques for those who may be unfamiliar with any of those fields. Recent MD simulations of tRNA dynamics and folding and of aminoacyl-tRNA synthetase dynamics and mechanistic characterizations are discussed. We highlight the recent successes and discuss how important questions can be addressed using current MD simulations techniques. We also outline several natural next steps for computational studies of AARS:tRNA complexes.