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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Insight into the kinetics and thermodynamics of the hydride transfer reactions between quinones and lumiflavin: a density functional theory study.

The kinetics and equilibrium of the hydride transfer reaction between lumiflavin and a number of substituted quinones was studied using density functional theory. The impact of electron withdrawing/donating substituents on the redox potentials of quinones was studied. In addition, the role of these substituents on the kinetics of the hydride transfer reaction with lumiflavin was investigated in detail under the transition state (TS) theory assumption. The hydride transfer reactions were found to be more favorable for an electron-withdrawing substituent. The activation barrier exhibited a quadratic relationship with the driving force of these reactions as derived under the formalism of modified Marcus theory. The present study found a significant extent of electron delocalization in the TS that is stabilized by enhanced electrostatic, polarization, and exchange interactions. Analysis of geometry, bond-orders, and energetics revealed a predominant parallel (Leffler-Hammond) effect on the TS. Closer scrutiny reveals that electron-withdrawing substituents, although located on the acceptor ring, reduce the N-H bond order of the donor fragment in the precursor complex. Carried out in the gas-phase, this is the first ever report of a theoretical study of flavin's hydride transfer reactions with quinones, providing an unfiltered view of the electronic effect on the nuclear reorganization of donor-acceptor complexes.

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.

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.

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

Coarse-grained simulations have emerged as invaluable tools for studying conformational changes in biomolecules. To evaluate the effectiveness of computationally inexpensive coarse-grained models in studying global and local dynamics of large protein systems like aminoacyl-tRNA synthetases, we have performed coarse-grained normal mode analysis, as well as principle component analysis on trajectories of all-atom and coarse-grained molecular dynamics simulations for three aminoacyl-tRNA synthetases--Escherichia coli methionyl-tRNA synthetase, Thermus thermophilus leucyl-tRNA synthetase, and Enterococcus faecium prolyl-tRNA synthetase. In the present study, comparison of predicted dynamics based on B-factor and overlap calculations revealed that coarse-grained methods are comparable to the all-atom simulations in depicting the intrinsic global dynamics of the three enzymes. However, the principal component analyses of the motions obtained from the all-atom molecular dynamics simulations provide a superior description of the local fluctuations of these enzymes. In particular, the all-atom model was able to capture the functionally relevant substrate-induced dynamical changes in prolyl-tRNA synthetase. The alteration in the coupled dynamics between the catalytically important proline-binding loop and its neighboring structural elements due to substrate binding has been characterized and reported for the first time. Taken together, the study portrays comparable and contrasting situations in studying the functional dynamics of large multi-domain aminoacyl-tRNA synthetases using coarse-grained and all-atom simulation methods.

Comparison of the intrinsic dynamics of aminoacyl-tRNA synthetases.

Aminoacyl-tRNA synthetases (AARSs) are an important family of enzymes that catalyze tRNA aminoacylation reaction (Ibba and Soll in Annu Rev Biochem 2000, 69:617-650) [1]. AARSs are grouped into two broad classes (class I and II) based on sequence/structural homology and mode of their interactions with the tRNA molecule (Ibba and Soll in Annu Rev Biochem 2000, 69:617-650) [1]. As protein dynamics play an important role in enzyme function, we explored the intrinsic dynamics of these enzymes using normal mode analysis and investigated if the two classes and six subclasses (Ia-c and IIa-c) of AARSs exhibit any distinct patterns of motion. The present study found that the intrinsic dynamics-based classification of these enzymes is similar to that obtained based on sequence/structural homology for most enzymes. However, the classification of seryl-tRNA synthetase was not straightforward; the internal mobility patterns of this enzyme are comparable to both IIa and IIb AARSs. This study revealed only a few general mobility patterns in these enzymes--(1) the insertion domain is generally engaged in anticorrelated motion with respect to the catalytic domain for both classes of AARSs and (2) anticodon binding domain dynamics are partly correlated and partly anticorrelated with respect to other domains for class I enzymes. In most of the class II AARSs, the anticodon binding domain is predominately engaged in anticorrelated motion with respect to the catalytic domain and correlated to the insertion domain. This study supports the notion that dynamic-based classification could be useful for functional classification of proteins.

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.