Polyethylene Glycol Impacts Conformation and Dynamics of Escherichia coli Prolyl-tRNA Synthetase Via Crowding and Confinement Effects.

Polyethylene glycol (PEG) is a flexible, nontoxic polymer commonly used in biological and medical research, and it is generally regarded as biologically inert. PEG molecules of variable sizes are also used as crowding agents to mimic intracellular environments. A recent study with PEG crowders revealed decreased catalytic activity of Escherichia coli prolyl-tRNA synthetase (Ec ProRS), where the smaller molecular weight PEGs had the maximum impact. The molecular mechanism of the crowding effects of PEGs is not clearly understood. PEG may impact protein conformation and dynamics, thus its function. In the present study, the effects of PEG molecules of various molecular weights and concentrations on the conformation and dynamics of Ec ProRS were investigated using a combined experimental and computational approach including intrinsic tryptophan fluorescence spectroscopy, atomic force microscopy, and atomistic molecular dynamic simulations. Results of the present study suggest that lower molecular weight PEGs in the dilute regime have modest effects on the conformational dynamics of Ec ProRS but impact the catalytic function primarily via the excluded volume effect; they form large clusters blocking the active site pocket. In contrast, the larger molecular weight PEGs in dilute to semidilute regimes have a significant impact on the protein's conformational dynamics; they wrap on the protein surface through noncovalent interactions. Thus, lower-molecular-weight PEG molecules impact protein dynamics and function via crowding effects, whereas larger PEGs induce confinement effects. These results have implications for the development of inhibitors for protein targets in a crowded cellular environment.

Crowder-Induced Conformational Ensemble Shift in Escherichia coli Prolyl-tRNA Synthetase.

The effect of molecular crowding on the structure and function of Escherichia coli prolyl-transfer RNA synthetase (Ec ProRS), a member of the aminoacyl-transfer RNA synthetase family, has been investigated using a combined experimental and theoretical method. Ec ProRS is a multidomain enzyme; coupled-domain dynamics are essential for efficient catalysis. To gain insight into the mechanistic detail of the crowding effect, kinetic studies were conducted with varying concentrations and sizes of crowders. In parallel, spectroscopic and quantum chemical studies were employed to probe the "soft interactions" between crowders and protein side chains. Finally, the dynamics of the dimeric protein was examined in the presence of crowders using a long-duration (70 ns) classical molecular dynamic simulations. The results of the simulations revealed a shift in the conformational ensemble, which is consistent with the preferential exclusion of cosolutes. The "soft interactions" model of the crowding effect also explained the alteration in kinetic parameters. In summary, the study found that the effects of molecular crowding on both conformational dynamics and catalytic function are correlated in the multidomain Ec ProRS, an enzyme that is central to protein synthesis in all living cells. This study affirmed that large and small cosolutes have considerable impacts on the structure, dynamics, and function of modular proteins and therefore must be considered for stabilizing protein-based pharmaceuticals and industrial enzymes.

Effect of stacking interactions on the thermodynamics and kinetics of lumiflavin: a study with improved density functionals and density functional tight-binding protocol.

The π-π stacking interaction between lumiflavin and a number of π-electron-rich molecules has been studied by density functional theory using several new-generation density functionals. Six known lumiflavin-aromatic adducts were used and the models were evaluated by comparing the geometry and energetics with experimental results. The study found that dispersion-corrected and hybrid functionals with larger (>50%) Hartree-Fock exchanges produced superior results in modeling thermodynamic characteristics of these complexes. The functional producing the best energetics for these model systems was used to study the stacking interactions of lumiflavin with biologically relevant aromatic groups. Additionally, the reduction of flavin-in the presence of both a hydride donor and a nondonor π-electronic system was also studied. Weak interactions were observed in the stacked lumiflavin complexes of benzene, phenol, and indole, mimicking phenyl alanine, tryptophan, and tyrosine side chains, respectively, of an enzyme. The stacked complex of naphthalene and flavin showed little change in flavin's redox potential indicating insignificant effect on the thermodynamics of the hydride transfer reaction. In contrast, the hydride transfer reaction with the hydride donor N-methyl nicotinamide tells a different story, as the transition state was found to be strongly impacted by the stacking interactions. A comparison of performance between the density functional theory (DFT) and the computationally less expensive dispersion-corrected self-consistent density functional tight-binding (SCC-DFTB-D) theory revealed that the latter produces consistent energetics for this hydride transfer reaction and additional DFT-computed perturbative corrections could significantly improve these results.

Strictly conserved lysine of prolyl-tRNA Synthetase editing domain facilitates binding and positioning of misacylated tRNA(Pro.).

To ensure high fidelity in translation, many aminoacyl-tRNA synthetases, enzymes responsible for attaching specific amino acids to cognate tRNAs, require proof-reading mechanisms. Most bacterial prolyl-tRNA synthetases (ProRSs) misactivate alanine and employ a post-transfer editing mechanism to hydrolyze Ala-tRNA(Pro). This reaction occurs in a second catalytic site (INS) that is distinct from the synthetic active site. The 2'-OH of misacylated tRNA(Pro) and several conserved residues in the Escherichia coli ProRS INS domain are directly involved in Ala-tRNA(Pro) deacylation. Although mutation of the strictly conserved lysine 279 (K279) results in nearly complete loss of post-transfer editing activity, this residue does not directly participate in Ala-tRNA(Pro) hydrolysis. We hypothesized that the role of K279 is to bind the phosphate backbone of the acceptor stem of misacylated tRNA(Pro) and position it in the editing active site. To test this hypothesis, we carried out pKa, charge neutralization, and free-energy of binding calculations. Site-directed mutagenesis and kinetic studies were performed to verify the computational results. The calculations revealed a considerably higher pKa of K279 compared to an isolated lysine and showed that the protonated state of K279 is stabilized by the neighboring acidic residue. However, substitution of this acidic residue with a positively charged residue leads to a significant increase in Ala-tRNA(Pro) hydrolysis, suggesting that enhancement in positive charge density in the vicinity of K279 favors tRNA binding. A charge-swapping experiment and free energy of binding calculations support the conclusion that the positive charge at position 279 is absolutely necessary for tRNA binding in the editing active site.

Writing a literature review as a class project in an upper-level undergraduate biochemistry course.

A literature review is an important part of conducting academic research. Knowing how to conduct a literature search and write a high-quality literature review is a valuable skill. Herein, the authors describe the method of introducing a literature review writing exercise in an upper-level biochemistry course. Since 2020, authors have collaborated with numerous undergraduates writing literature reviews on topics in biochemistry that resulted in peer-reviewed publications. Authors believe that this unique idea of providing a course-based undergraduate research experience (CURE) to many undergraduates, especially those who otherwise do not receive collaborative research experience through traditional research paths, must be shared with other instructors.

Virtual Screening of Molecules via Neural Fingerprint-based Deep Learning Technique.

A machine learning-based drug screening technique has been developed and optimized using convolutional neural network-derived fingerprints. The optimization of weights in the neural network-based fingerprinting technique was compared with fixed Morgan fingerprints in regard to binary classification on drug-target binding affinity. The assessment was carried out using six different target proteins using randomly chosen small molecules from the ZINC15 database for training. This new architecture proved to be more efficient in screening molecules that less favorably bind to specific targets and retaining molecules that favorably bind to it. Scientific contribution: We have developed a new neural fingerprint-based screening model that has a significant ability to capture hits. Despite using a smaller dataset, this model is capable of mapping chemical space similar to other contemporary algorithms designed for molecular screening. The novelty of the present algorithm lies in the speed with which the models are trained and tuned before testing its predictive capabilities and hence is a significant step forward in the field of machine learning-embedded computational drug discovery.

Multiple pathways promote dynamical coupling between catalytic domains in Escherichia coli prolyl-tRNA synthetase.

Aminoacyl-tRNA synthetases are multidomain enzymes that catalyze covalent attachment of amino acids to their cognate tRNA. Cross-talk between functional domains is a prerequisite for this process. In this study, we investigate the molecular mechanism of site-to-site communication in Escherichia coli prolyl-tRNA synthetase (Ec ProRS). Earlier studies have demonstrated that evolutionarily conserved and/or co-evolved residues that are engaged in correlated motion are critical for the propagation of functional conformational changes from one site to another in modular proteins. Here, molecular simulation and bioinformatics-based analysis were performed to identify dynamically coupled and evolutionarily constrained residues that form contiguous pathways of residue-residue interactions between the aminoacylation and editing domains of Ec ProRS. The results of this study suggest that multiple pathways exist between these two domains to maintain the dynamic coupling essential for enzyme function. Moreover, residues in these interaction networks are generally highly conserved. Site-directed changes of on-pathway residues have a significant impact on enzyme function and dynamics, suggesting that any perturbation along these pathways disrupts the native residue-residue interactions that are required for effective communication between the two functional domains. Free energy analysis revealed that communication between residues within a pathway and cross-talk between pathways are important for coordinating functions of different domains of Ec ProRS for efficient catalysis.

Polyethylene Glycol 20k. Does It Fluoresce?

Polyethylene glycol (PEG) is a polyether compound commonly used in biological research and medicine because it is biologically inert. This simple polymer exists in variable chain lengths (and molecular weights). As they are devoid of any contiguous π-system, PEGs are expected to lack fluorescence properties. However, recent studies suggested the occurrence of fluorescence properties in non-traditional fluorophores like PEGs. Herein, a thorough investigation has been conducted to explore if PEG 20k fluoresces. Results of this combined experimental and computational study suggested that although PEG 20k could exhibit "through-space" delocalization of lone pairs of electrons in aggregates/clusters, formed via intermolecular and intramolecular interactions, the actual contributor of fluorescence between 300 and 400 nm is the stabilizer molecule, i.e., 3-tert-butyl-4-hydroxyanisole present in the commercially available PEG 20k. Therefore, the reported fluorescence properties of PEG should be taken with a grain of salt, warranting further investigation.

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.

Role of coupled dynamics in the catalytic activity of prokaryotic-like prolyl-tRNA synthetases.

Prolyl-tRNA synthetases (ProRSs) have been shown to activate both cognate and some noncognate amino acids and attach them to specific tRNA(Pro) substrates. For example, alanine, which is smaller than cognate proline, is misactivated by Escherichia coli ProRS. Mischarged Ala-tRNA(Pro) is hydrolyzed by an editing domain (INS) that is distinct from the activation domain. It was previously shown that deletion of the INS greatly reduced cognate proline activation efficiency. In this study, experimental and computational approaches were used to test the hypothesis that deletion of the INS alters the internal protein dynamics leading to reduced catalytic function. Kinetic studies with two ProRS variants, G217A and E218A, revealed decreased amino acid activation efficiency. Molecular dynamics studies showed motional coupling between the INS and protein segments containing the catalytically important proline-binding loop (PBL, residues 199-206). In particular, the complete deletion of INS, as well as mutation of G217 or E218 to alanine, exhibited significant effects on the motion of the PBL. The presence of coupled dynamics between neighboring protein segments was also observed through in silico mutations and essential dynamics analysis. Altogether, this study demonstrates that structural elements at the editing domain-activation domain interface participate in coupled motions that facilitate amino acid binding and catalysis by bacterial ProRSs, which may explain why truncated or defunct editing domains have been maintained in some systems, despite the lack of catalytic activity.

Pre-Existing Oxidative Stress Creates a Docking-Ready Conformation of the SARS-CoV-2 Receptor-Binding Domain.

The redox-dependent changes on the binding between the receptor-binding domain of the severe acute respiratory syndrome-coronavirus-2 spike protein and the peptidase domain of the human cell surface receptor angiotensin-converting enzyme II were investigated by performing molecular dynamics simulations. The reduced states of the protein partners were generated in silico by converting the disulfides to thiols. The role of redox transformation on the protein-protein binding affinity was assessed from the time-evolved structures after 200 ns simulations using electrostatic field calculations and implicit solvation. The present simulations revealed that the bending motion at the protein-protein interface is significantly altered when the disulfides are reduced to thiols. In the native complex, the presence of disulfide bonds preserves the structural complementarity of the protein partners and maintains the intrinsic conformational dynamics. Also, the study demonstrates that when already bound, the disulfide-to-thiol conversion of the receptor-binding domain has a limited impact on the binding of the spike protein to the receptor. However, if the reduction occurs before binding to the receptor, a spectacular conformational change of the receptor-binding domain occurs that fully impairs the binding. In other words, the formation of disulfide bonds, prevalent during oxidative stress, creates a conformation ready to bind to the receptor. Taken together, the present study demonstrates the role of pre-existing oxidative stress in elevating the binding affinity of the spike protein for the human receptor, offering future clues for alternate therapeutic possibilities.

Evolution of Stronger SARS-CoV-2 Variants as Revealed Through the Lens of Molecular Dynamics Simulations.

Using molecular dynamics simulations, the protein-protein interactions of the receptor-binding domain of the wild-type and seven variants of the severe acute respiratory syndrome coronavirus 2 spike protein and the peptidase domain of human angiotensin-converting enzyme 2 were investigated. These variants are alpha, beta, gamma, delta, eta, kappa, and omicron. Using 100 ns simulation data, the residue interaction networks at the protein-protein interface were identified. Also, the impact of mutations on essential protein dynamics, backbone flexibility, and interaction energy of the simulated protein-protein complexes were studied. The protein-protein interface for the wild-type, delta, and omicron variants contained several stronger interactions, while the alpha, beta, gamma, eta, and kappa variants exhibited an opposite scenario as evident from the analysis of the inter-residue interaction distances and pair-wise interaction energies. The study reveals that two distinct residue networks at the central and right contact regions forge stronger binding affinity between the protein partners. The study provides a molecular-level insight into how enhanced transmissibility and infectivity by delta and omicron variants are most likely tied to a handful of interacting residues at the binding interface, which could potentially be utilized for future antibody constructs and structure-based antiviral drug design.

Vitamin D and COVID-19: A review on the role of vitamin D in preventing and reducing the severity of COVID-19 infection.

Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) is a pathogenic coronavirus causing COVID-19 infection. The interaction between the SARS-CoV-2 spike protein and the human receptor angiotensin-converting enzyme 2, both of which contain several cysteine residues, is impacted by the disulfide-thiol balance in the host cell. The host cell redox status is affected by oxidative stress due to the imbalance between the reactive oxygen/nitrogen species and antioxidants. Recent studies have shown that Vitamin D supplementation could reduce oxidative stress. It has also been proposed that vitamin D at physiological concentration has preventive effects on many viral infections, including COVID-19. However, the molecular-level picture of the interplay of vitamin D deficiency, oxidative stress, and the severity of COVID-19 has remained unclear. Herein, we present a thorough review focusing on the possible molecular mechanism by which vitamin D could alter host cell redox status and block viral entry, thereby preventing COVID-19 infection or reducing the severity of the disease.

Impact of Thiol-Disulfide Balance on the Binding of Covid-19 Spike Protein with Angiotensin-Converting Enzyme 2 Receptor.

The novel coronavirus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has led to an ongoing pandemic of coronavirus disease (COVID-19), which started in 2019. This is a member of Coronaviridae family in the genus Betacoronavirus, which also includes SARS-CoV and Middle East respiratory syndrome coronavirus (MERS-CoV). The angiotensin-converting enzyme 2 (ACE2) is the functional receptor for SARS-CoV and SARS-CoV-2 to enter the host cells. In particular, the interaction of viral spike proteins with ACE2 is a critical step in the viral replication cycle. The receptor-binding domain of the viral spike proteins and ACE2 have several cysteine residues. In this study, the role of thiol-disulfide balance on the interactions between SARS-CoV/CoV-2 spike proteins and ACE2 was investigated using molecular dynamics simulations. The study revealed that the binding affinity was significantly impaired when all of the disulfide bonds of both ACE2 and SARS-CoV/CoV-2 spike proteins were reduced to thiol groups. The impact on the binding affinity was less severe when the disulfide bridges of only one of the binding partners were reduced to thiols. This computational finding possibly provides a molecular basis for the differential COVID-19 cellular recognition due to the oxidative stress.

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.

Role of Oxidative Stress on SARS-CoV (SARS) and SARS-CoV-2 (COVID-19) Infection: A Review.

Novel coronavirus disease 2019 (COVID-19) has resulted in a global pandemic and is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Several studies have suggested that a precise disulfide-thiol balance is crucial for viral entry and fusion into the host cell and that oxidative stress generated from free radicals can affect this balance. Here, we reviewed the current knowledge about the role of oxidative stress on SARS-CoV and SARS-CoV-2 infections. We focused on the impact of antioxidants, like NADPH and glutathione, and redox proteins, such as thioredoxin and protein disulfide isomerase, that maintain the disulfide-thiol balance in the cell. The possible influence of these biomolecules on the binding of viral protein with the host cell angiotensin-converting enzyme II receptor protein as well as on the severity of COVID-19 infection was discussed.

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.

Critical role of substrate conformational change in the proton transfer process catalyzed by 4-oxalocrotonate tautomerase.

4-Oxalocrotonate tautomerase enzyme (4-OT) catalyzes the isomerization of 2-oxo-4-hexenedioate to 2-oxo-3-hexenedioate. The chemical process involves two proton transfers, one from a carbon of the substrate to the nitrogen of Pro1 and another from this nitrogen atom to a different carbon of the substrate. In this paper the isomerization has been studied using the combined quantum mechanical and molecular mechanical method with a dual-level treatment of the quantum subsystem employing the MPW1BK density functional as the higher level. Exploration of the potential energy surface shows that the process is stepwise, with a stable intermediate state corresponding to the deprotonated substrate and a protonated proline. The rate constant of the overall process has been evaluated using ensemble-averaged variational transition state theory, including the quantized vibrational motion of a primary zone of active-site atoms and a transmission coefficient based on an ensemble of optimized reaction coordinates to account for recrossing trajectories and optimized multidimensional tunneling. The two proton-transfer steps have similar free energy barriers, but the transition state associated with the first proton transfer is found to be higher in energy. The calculations show that reaction progress is coupled to a conformational change of the substrate, so it is important that the simulation allows this flexibility. The coupled conformational change is promoted by changes in the electron distribution of the substrate that take place as the proton transfers occur.

Theoretical determination of the redox potentials of NRH:quinone oxidoreductase 2 using quantum mechanical/molecular mechanical simulations.

NRH:quinone oxidoreductase 2 (NQO2) is a flavoenzyme that catalyzes a one-step two-electron reduction of quinones. During this enzyme catalysis, the 7,8-dimethyl isoalloxazine (flavin) ring of the enzyme-bound cofactor, flavin adenine dinucleotide (FAD), shuttles between reduced and oxidized states as the enzyme passes through multiple cycles of binding/release of alternate substrates. These redox changes in NQO2, however, lead to unequal charge separation between the flavin ring and the active site, which must be stabilized by reorganization of the surrounding protein matrix. In this study, we have used a combined quantum mechanical/molecular mechanical method to simulate the electron and proton addition reactions of the flavin-bound NQO2. We have computed the redox potentials and pK(a)'s of the enzyme-bound flavin. The present work demonstrates that upon reduction, the NQO2 active site stabilizes the flavin anionic hydroquinone state. Simulation data has also allowed quantitative estimation of the electrostatic contributions of active site residues. Their significance in oscillatory redox transition of this flavoenzyme is discussed.

Cyclic Changes in Active Site Polarization and Dynamics Drive the 'Ping-pong' Kinetics in NRH:Quinone Oxidoreductase 2: An Insight from QM/MM Simulations.

Quinone reductases belong to the family of flavin-dependent oxidoreductases. With the redox active cofactor, flavin adenine dinucleotide, quinone reductases are known to utilize a 'ping-pong' kinetic mechanism during catalysis in which a hydride is bounced back and forth between flavin and its two substrates. However, the continuation of this catalytic cycle requires product displacement steps, where the product of one redox half-cycle is displaced by the substrate of the next half-cycle. Using improved hybrid quantum mechanical/molecular mechanical simulations, both the catalytic hydride transfer and the product displacement reactions were studied in NRH:quinone oxidoreductase 2. Initially, the self-consistent charge-density functional tight binding theory was used to describe flavin ring and the substrate atoms, while embedded in the molecular mechanically-treated solvated active site. Then, for each step of the catalytic cycle, a further improvement of energetics was made using density functional theory-based corrections. The present study showcases an integrated interplay of solvation, protonation, and protein matrix-induced polarization as the driving force behind the thermodynamic wheel of the 'ping-pong' kinetics. Reported here is the first-principles model of the 'ping-pong' kinetics that portrays how cyclic changes in the active site polarization and dynamics govern the oscillatory hydride transfer and product displacement in this enzyme.