Steered molecular dynamics simulations of ligand-receptor interaction in lipocalins

Insights from Structure-Based Simulations into the Persulfidation of Uridine Diphosphate-Glycosyltransferase71c5 Facilitating the Reversible Inactivation of Abscisic Acid

The action of abscisic acid (ABA) is closely related to its level in plant tissues. Uridine diphosphate-glycosyltransferase71c5 (UGT71C5) was characterized as a major UGT enzyme to catalyze the formation of the ABA-glucose ester (ABA-GE), a reversible inactive form of free ABA in Arabidopsis thaliana (thale cress). UGTs function in a mode where the catalytic base deprotonates an acceptor to allow a nucleophilic attack at the anomeric center of the donor, achieving the transfer of a glucose moiety. The proteomic data revealed that UGT71C5 can be persulfidated. Herein, an experimental method was employed to detect the persulfidation site of UGT71C5, and the computational methods were further used to identify the yet unknown molecular basis of ABA glycosylation as well as the regulatory role of persulfidation in this process. Our results suggest that the linker and the U-shaped loop are regulatory structural elements: the linker is associated with the binding of uridine diphosphate glucose (UPG) and the U-shaped loop is involved in binding both UPG and ABA.It was also found that it is through tuning the dynamics of the U-shaped loop that is accompanied by the movement of tyrosine (Y388) that the persulfidation of cysteine (C311) leads to the catalytic residue histidine (H16) being in place, preparing for the deprotonation of ABA, and then reorientates UPG and deprotonated ABA closer to the 'Michaelis' complex, facilitating the transfer of a glucose moiety. Ultimately, the persulfidation of UGT71C5 is in favor of ABA glycosylation. Our results provide insights into the molecular details of UGT71C5 recognizing substrates and insights concerning persulfidation as a possible mechanism for hydrogen sulfide (H2S) to modulate the content of ABA, which helps us understand how modulating ABA level strengthens plant tolerance.

The interaction analysis between human serum albumin with chlorpyrifos and its derivatives through sub-atomic docking and molecular dynamics simulation techniques

Chlorpyrifos (CPF) is an extensively used organophosphate pesticide for crop protection. However, there are concerns about it contaminating the environment and human health, with estimated three lakh deaths annually. The molecular modeling protocol was assisted in redesigning thirteen well-known CPF linkers and inserting them at five selectable CPF (R1-R5) positions of CPF to get 258 CPF derivatives. CPF and its derivatives were optimized using LigPrep and docked to a grid centralized on Trp214 using extra precision glide docking. The Binding free energy of complexes was calculated using molecular mechanics/generalized born surface area (MM-GBSA). CPF and CPFD-225 have glide scores of - 3.08 and - 6.152 kcal/mol, respectively, with human serum albumin and ΔG bind for CPF (- 33.041817 kcal/mol) (- 52.825 kcal/mol) for CPF-D225. The top ten CPF derivatives showed at least ninefold better binding free energy than the CPF proposed for polyclonal antibody production. Subsequently, molecular docking studies revealed that CPF and its derivatives could bind to human serum albumin (HSA). Furthermore, using the Desmond package, a 100-ns molecular dynamics (MD) simulation was performed on the potential binding site. The final systems of CPF-HSA and CPF-222D complexes consist of 76,014 and 76,026 atoms, respectively. The physical stability of both the systems (CPF-HSA and CPF-222D) was analyzed by considering the overall potential energy, RMSF, RMSD, Hydrophobic interactions, and water-mediated patterns, which showed total energy of - 141,610 kcal/mol and - 140,150 kcal/mol, respectively.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03344-7.

Protective effects of a small molecule inhibitor ligand against hyperphosphorylated tau-induced mitochondrial and synaptic toxicities in Alzheimer disease

The purpose of our study is to understand the protective effects of small molecule ligands for phosphorylated tau (p-tau) in Alzheimer's disease (ad) progression. Many reports show evidence that p-tau is reported to be an important contributor to the formation of paired helical filaments (PHFs) and neurofibrillary tangles (NFTs) in ad neurons. In ad, glycogen synthase kinase-3 beta (GSK3β), cyclin-dependent kinase- 5 (CDK5) and dual specificity tyrosine-phosphorylation-regulated kinase 1A (DYRK1A), are the three important kinases responsible for tau hyperphosphorylation. Currently, there are no drugs and/or small molecules that reduce the toxicity of p-tau in ad. In the present study, we rationally selected and validated small molecule ligands that binds to the phosphorylated tau at SER23 (Ser 285). We also assessed the molecular dynamics and validated molecular docking sites for the three best ligands. Based on the best docking scores -8.09, -7.9 and - 7.8 kcal/mol, we found that ligand 1 binds to key hyperphosphorylation residues of p-tau that inhibit abnormal PHF-tau, DYRK1A, and GKS3β that reduce p-tau levels in ad. Using biochemical, molecular, immunoblotting, immunofluorescence, and transmission electron microscopy analyses, we studied the ligand 1 inhibition as well as mitochondrial and synaptic protective effects in immortalized primary hippocampal neuronal (HT22) cells. We found interactions between NAT10-262501 (ligand 1) and p-tau at key phosphorylation sites and these ligand-based inhibitions decreased PHF-tau, DYRK1A and GSK3β levels. We also found increased mitochondrial biogenesis, mitochondrial fusion and synaptic activities and reduced mitochondrial fission in ligand 1-treated mutant tau HT22 cells. Based on these results, we cautiously conclude that p-tau NAT10-262501 (ligand 1) reduces hyperphosphorylation of tau based GKS3β and CDK5 kinase regulation in ad, and aids in the maintenance of neuronal structure, mitochondrial dynamics, and biogenesis with a possible therapeutic drug target for ad.

Binding characteristics of staphylococcal protein A and streptococcal protein G for fragment crystallizable portion of human immunoglobulin G

In the wide array of physiological processes, protein-protein interactions and their binding are the most basal activities for achieving adequate biological metabolism. Among the studies on binding proteins, the examination of interactions between immunoglobulin G (IgG) and natural immunoglobulin-binding ligands, such as staphylococcal protein A (spA) and streptococcal protein G (spG), is essential in the development of pharmaceutical science, biotechnology, and affinity chromatography. The widespread utilization of IgG-spA/spG binding characteristics has allowed researchers to investigate these molecular interactions. However, the detailed binding strength of each ligand and the corresponding binding mechanisms have yet to be fully investigated. In this study, the authors analyzed the binding strengths of IgG-spA and IgG-spG complexes and identified the mechanisms enabling these bindings using molecular dynamics simulation, steered molecular dynamics, and advanced Poisson-Boltzmann Solver simulations. Based on the presented data, the binding strength of the spA ligand was found to significantly exceed that of the spG ligand. To find out which non-covalent interactions or amino acid sites have a dominant role in the tight binding of these ligands, further detailed analyses of electrostatic interactions, hydrophobic bonding, and binding free energies have been performed. In investigating their binding affinity, a relatively independent and different unbinding mechanism was found in each ligand. These distinctly different mechanisms were observed to be highly correlated to the protein secondary and tertiary structures of spA and spG ligands, as explicated from the perspective of hydrogen bonding.

Using steered molecular dynamics to study the interaction between ADP and the nucleotide-binding domain of yeast Hsp70 protein Ssa1

Genetics experiments have identified six mutations located in the subdomain IA (A17V, R23H, G32D, G32S, R34K, V372I) of Ssa1 that influence propagation of the yeast [PSI⁺] prion. However, the underlining molecular mechanisms of these mutations are still unclear. The six mutation sites are present in the IA subdomain of the nucleotide-binding domain (NBD). The ATPase subdomain IA is a critical mediator of inter-domain allostery in Hsp70 molecular chaperones, so the mutation and changes in this subdomain may influence the function of the substrate-binding domain. In addition, ADP release is a rate-limiting step of the ATPase cycle and dysregulation of the ATPase cycle influences the propagation of the yeast [PSI⁺] prion. In this work, steered molecular dynamics (SMD) simulations were performed to explore the interaction between ADP and NBD. Results suggest that during the SMD simulations, hydrophobic interactions are predominant and variations in the binding state of ADP within the mutants is a potential reason for in vivo effects on yeast [PSI⁺] prion propagation. Additionally, we identify the primary residues in the ATPase domain that directly constitute the main hydrophobic interaction network and directly influence the ADP interaction state with the NBD of Ssa1. Furthermore, this in silico analysis reaffirms the importance of previously experimentally-determined residues in the Hsp70 ATPase domain involved in ADP binding and also identifies new residues potentially involved in this process.

An overview of recent molecular dynamics applications as medicinal chemistry tools for the undruggable site challenge

Molecular dynamics (MD) has become increasingly popular due to the development of hardware and software solutions and the improvement in algorithms, which allowed researchers to scale up calculations in order to speed them up. MD simulations are usually used to address protein folding issues or protein-ligand complex stability through energy profile analysis over time. In recent years, the development of new tools able to deeply explore a potential energy surface (PES) has allowed researchers to focus on the dynamic nature of the binding recognition process and binding-induced protein conformational changes. Moreover, modern approaches have been demonstrated to be effective and reliable in calculating some kinetic and thermodynamic parameters behind the host-guest recognition process. Starting from all of these considerations, several efforts have been made in order to integrate MD within the virtual screening process in drug discovery. Knowledge retrieved from MD can, in fact, be exploited as a starting point to build pharmacophores or docking constraints in the early stage of the screening campaign as well as to define key features, in order to unravel hidden binding modes and help the optimisation of the molecular structure of a lead compound. Based on these outcomes, researchers are nowadays using MD as an invaluable tool to discover and target previously considered undruggable binding sites, including protein-protein interactions and allosteric sites on a protein surface. As a matter of fact, the use of MD has been recognised as vital to the discovery of selective protein-protein interaction modulators. The use of a dynamic overview on how the host-guest recognition occurs and of the relative conformational modifications induced allows researchers to optimise small molecules and small peptides capable of tightly interacting within the cleft between two proteins. In this review, we aim to present the most recent applications of MD as an integrated tool to be used in the rational design of small molecules or small peptides able to modulate undruggable targets, such as allosteric sites and protein-protein interactions.

Molecular dynamics simulation and steered molecular dynamics simulation on irisin dimers

Irisin is found closely associated with promoting the browning of beige fat cells in white adipose tissue. The crystal structure reveals that irisin forms a continuous inter-subunit β-sheet dimer. Here, molecular dynamics (MD) simulation and steered molecular dynamics (SMD) simulation were performed to investigate the dissociation process and the intricate interactions between the two irisin monomers. In the process of MD, the interactions between the monomers were roughly analyzed through the average numbers of both hydrophobic contacts and H-bonds. Then, SMD was performed to investigate the accurate interaction energy between the monomers. By the analysis of dissociation energy, the van der Waals (vdW) force was identified as the major energy to maintain the dimer structure, which also verified the results of MD simulation. Meanwhile, 11 essential residues were discovered by the magnitude of rupture force during dissociation. Among them, residues Arg75, Glu79, Ile77, Ala88, and Trp90 were reported in a previous study using the method of mutagenesis and size exclusion chromatography, and several new important residues (Arg72, Leu74, Phe76, Gln78, Val80, and Asp91) were also identified. Interestingly, the new important residues that we discovered and the important residues that were reported are located in the opposite side of the β-sheet of the dimer.

Steered molecular dynamics simulation of the binding of the bovine auxilin J domain to the Hsc70 nucleotide-binding domain

The Hsp70 and Hsp40 chaperone machine plays critical roles in protein folding, membrane translocation, and protein degradation by binding and releasing protein substrates in a process that utilizes ATP. The activities of the Hsp70 family of chaperones are recruited and stimulated by the J domains of Hsp40 chaperones. However, structural information on the Hsp40-Hsp70 complex is lacking, and the molecular details of this interaction are yet to be elucidated. Here we used steered molecular dynamics (SMD) simulations to investigate the molecular interactions that occur during the dissociation of the auxilin J domain from the Hsc70 nucleotide-binding domain (NBD). The changes in energy observed during the SMD simulation suggest that electrostatic interactions are the dominant type of interaction. Additionally, we found that Hsp70 mainly interacts with auxilin through the surface residues Tyr866, Arg867, and Lys868 of helix II, His874, Asp876, Lys877, Thr879, and Gln881 of the HPD loop, and Phe891, Asn895, Asp896, and Asn903 of helix III. The conservative residues Tyr866, Arg867, Lys868, His874, Asp876, Lys877, and Phe891 were also found in a previous study to be indispensable to the catalytic activity of the DnaJ J domain and the binding of it with the NBD of DnaK. The in silico identification of the importance of auxilin residues Asn895, Asp896, and Asn903 agrees with previous mutagenesis and NMR data suggesting that helix III of the J domain of the T antigen interacts with Hsp70. Furthermore, our data indicate that Thr879 and Gln881 from the HPD loop are also important as they mediate the interaction between the bovine auxilin J domain and Hsc70.

A steered molecular dynamics mediated hit discovery for histone deacetylases

The inhibitors of class I histone deacetylases (HDACIs) have gained significant interest in cancer therapeutics. Virtual high throughput screening (vHTS) is one of the popular approaches used in the identification of novel scaffolds of HDACIs. However, an accurate description of ligand-protein flexibilities in the vHTS remains challenging. In this work, we implement an integrated approach, which combines the vHTS with the 'state-of-the-art' steered molecular dynamics (SMD). This approach serves as an efficient tool to identify potential hits and characterize their binding potencies against the class I HDACs in a flexible solvent environment. A hybrid pharmacophore-based and structure-based vHTS method identifies the hits with more favourable physico-chemical features against the class I HDACs. Our pharmacophore-based screening enhanced the quality of the vHTS outcomes. Further, the molecular interactions between the hits and the HDACs are investigated using the SMD-driven force profiles, which in turn resulted in filtering the hits with higher binding potencies against the HDACs. Our results, therefore, reveal that vHTS and SMD can be a complementary and effective analytical tool for accelerating the hit identification phase in structure-based drug design.

In silico insights into protein-protein interactions and folding dynamics of the saposin-like domain of Solanum tuberosum aspartic protease

The plant-specific insert is an approximately 100-residue domain found exclusively within the C-terminal lobe of some plant aspartic proteases. Structurally, this domain is a member of the saposin-like protein family, and is involved in plant pathogen defense as well as vacuolar targeting of the parent protease molecule. Similar to other members of the saposin-like protein family, most notably saposins A and C, the recently resolved crystal structure of potato (Solanum tuberosum) plant-specific insert has been shown to exist in a substrate-bound open conformation in which the plant-specific insert oligomerizes to form homodimers. In addition to the open structure, a closed conformation also exists having the classic saposin fold of the saposin-like protein family as observed in the crystal structure of barley (Hordeum vulgare L.) plant-specific insert. In the present study, the mechanisms of tertiary and quaternary conformation changes of potato plant-specific insert were investigated in silico as a function of pH. Umbrella sampling and determination of the free energy change of dissociation of the plant-specific insert homodimer revealed that increasing the pH of the system to near physiological levels reduced the free energy barrier to dissociation. Furthermore, principal component analysis was used to characterize conformational changes at both acidic and neutral pH. The results indicated that the plant-specific insert may adopt a tertiary structure similar to the characteristic saposin fold and suggest a potential new structural motif among saposin-like proteins. To our knowledge, this acidified PSI structure presents the first example of an alternative saposin-fold motif for any member of the large and diverse SAPLIP family.

Ligand binding-dependent functions of the lipocalin NLaz: an in vivo study in Drosophila

Lipocalins are small extracellular proteins mostly described as lipid carriers. The Drosophila lipocalin NLaz (neural Lazarillo) modulates the IIS pathway and regulates longevity, stress resistance, and behavior. Here, we test whether a native hydrophobic pocket structure is required for NLaz to perform its functions. We use a point mutation altering the binding pocket (NLaz(L130R)) and control mutations outside NLaz binding pocket. Tryptophan fluorescence titration reveals that NLaz(L130R) loses its ability to bind ergosterol and the pheromone 7(z)-tricosene but retains retinoic acid binding. Using site-directed transgenesis in Drosophila, we test the functionality of the ligand binding-altered lipocalin at the organism level. NLaz-dependent life span reduction, oxidative stress and starvation sensitivity, aging markers accumulation, and deficient courtship are rescued by overexpression of NLaz(WT), but not of NLaz(L130R). Transcriptional responses to aging and oxidative stress show a large set of age-responsive genes dependent on the integrity of NLaz binding pocket. Inhibition of IIS activity and modulation of oxidative stress and infection-responsive genes are binding pocket-dependent processes. Control of energy metabolites on starvation appears to be, however, insensitive to the modification of the NLaz binding pocket.

Genome-based approaches to develop epitope-driven subunit vaccines against pathogens of infective endocarditis

Infective endocarditis (IE) has emerged as a public health problem due to changes in the etiologic spectrum and due to involvement of resistant bacterial strains with increased virulence. Developing potent vaccine is an important strategy to tackle IE. Complete genome sequences of eight selected pathogens of IE paved the way to design common T-cell driven subunit vaccines. Comparative genomics and subtractive genomic analysis were applied to identify adinosine tri phosphate (ATP)-binding cassette (ABC) transporter ATP-binding protein from Streptococcus mitis (reference organism) as common vaccine target. Reverse vaccinology technique was implemented using computational tools such as ProPred, SYFPEITHI, and Immune epitope database. Twenty-one T-cell epitopes were predicted from ABC transporter ATP-binding protein. Multiple sequence alignment of ABC transporter ATP-binding protein from eight selected IE pathogens was performed to identify six conserved T-cell epitopes. The six selected T-cell epitopes were further evaluated at structure level for HLA-DRB binding through homology modeling and molecular docking analysis using Maestro v9.2. The proposed six T-cell epitopes showed better binding affinity with the selected HLA-DRB alleles. Subsequently, the docking complexes of T-cell epitope and HLA-DRBs were ranked based on XP Gscore. The T-cell epitope (208-LNYITPDVV-216)-HLA-DRB1(∗)0101 (1T5 W) complex having the best XP Gscore (-13.25 kcal/mol) was assessed for conformational stability and interaction stability through molecular dynamic simulation for 10 ns using Desmond v3.2. The simulation results revealed that the HLA-DRB-epitope complex was stable throughout the simulation time. Thus, the epitope would be ideal candidate for T-cell driven subunit vaccine design against infective endocarditis.

Para-(benzoyl)-phenylalanine as a potential inhibitor against LpxC of Leptospira spp.: homology modeling, docking, and molecular dynamics study

Leptospira interrogans, a Gram-negative bacterial pathogen is the main cause of human leptospirosis. Lipid A is a highly immunoreactive endotoxic center of lipopolysaccharide (LPS) that anchors LPS into the outer membrane of Leptospira. Discovery of compounds inhibiting lipid-A biosynthetic pathway would be promising for dissolving the structural integrity of membrane leading to cell lysis and death of Leptospira. LpxC, a unique enzyme of lipid-A biosynthetic pathway was identified as common drug target of Leptospira. Herein, homology modeling, docking, and molecular dynamics (MD) simulations were employed to discover potential inhibitors of LpxC. A reliable tertiary structure of LpxC in complex with inhibitor BB-78485 was constructed in Modeller 9v8. A data-set of BB-78485 structural analogs were docked with LpxC in Maestro v9.2 virtual screening workflow, which implements three stage Glide docking protocol. Twelve lead molecules with better XP Gscore compared to BB-78485 were proposed as potential inhibitors of LpxC. Para-(benzoyl)-phenylalanine - that showed lowest XP Gscore (-10.35 kcal/mol) - was predicted to have best binding affinity towards LpxC. MD simulations were performed for LpxC and para-(benzoyl)-phenylalanine docking complex in Desmond v3.0. Trajectory analysis showed the docking complex and inter-molecular interactions was stable throughout the entire production part of MD simulations. The results indicate para-(benzoyl)-phenylalanine as a potent drug molecule against leptospirosis. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:10.

Using steered molecular dynamics to predict and assess Hsp70 substrate-binding domain mutants that alter prion propagation

Genetic screens using Saccharomyces cerevisiae have identified an array of cytosolic Hsp70 mutants that are impaired in the ability to propagate the yeast [PSI⁺] prion. The best characterized of these mutants is the Ssa1 L483W mutant (so-called SSA1-21), which is located in the substrate-binding domain of the protein. However, biochemical analysis of some of these Hsp70 mutants has so far failed to provide major insight into the specific functional changes in Hsp70 that cause prion impairment. In order to gain a better understanding of the mechanism of Hsp70 impairment of prions we have taken an in silico approach and focused on the Escherichia coli Hsp70 ortholog DnaK. Using steered molecular dynamics simulations (SMD) we demonstrate that DnaK variant L484W (analogous to SSA1-21) is predicted to bind substrate more avidly than wild-type DnaK due to an increase in numbers of hydrogen bonds and hydrophobic interactions between chaperone and peptide. Additionally the presence of the larger tryptophan side chain is predicted to cause a conformational change in the peptide-binding domain that physically impairs substrate dissociation. The DnaK L484W variant in combination with some SSA1-21 phenotypic second-site suppressor mutations exhibits chaperone-substrate interactions that are similar to wild-type protein and this provides a rationale for the phenotypic suppression that is observed. Our computational analysis fits well with previous yeast genetics studies regarding the functionality of the Ssa1-21 protein and provides further evidence suggesting that manipulation of the Hsp70 ATPase cycle to favor the ADP/substrate-bound form impairs prion propagation. Furthermore, we demonstrate how SMD can be used as a computational tool for predicting Hsp70 peptide-binding domain mutants that impair prion propagation.

Direct non-covalent interactions between protein substrates and molecular chaperones play crucial roles in the protein folding process. [PSI⁺] is a prion of the yeast Saccharomyces cerevisiae, which is formed by mis-folding of the native Sup35 protein in a process analogous to formation of prions in mammals. While much genetic data exists showing a clear role for Hsp70 in prion propagation, biochemical data has yet to provide a clear link to how Hsp70 functions in prion propagation or how some Hsp70 mutants successfully impair in vivo propagation of prions. This paper employs a novel simulation method termed “steered molecular dynamics” to explore the different types and amounts of non-covalent interactions between wild type and mutated Hsp70, with a model substrate. Extrapolating the in silico data allowed us to decipher how a mutant Hsp70 impairs yeast prion propagation and allows us to predict other Hsp70 mutants that should behave in the same manner and to test these predictions in a yeast-based system. Our computational data shows that increasing the binding affinity of Hsp70 for substrate is one way of impairing prion propagation, a proposal that correlates very well with previous experimental genetic data.