Disruption and formation of surface salt bridges are coupled to DNA binding by the integration host factor: a computational analysis

Modeling coronavirus spike protein dynamics: implications for immunogenicity and immune escape

The ongoing COVID-19 pandemic is a global public health emergency requiring urgent development of efficacious vaccines. While concentrated research efforts have focused primarily on antibody-based vaccines that neutralize SARS-CoV-2, and several first-generation vaccines have either been approved or received emergency use authorization, it is forecasted that COVID-19 will become an endemic disease requiring updated second-generation vaccines. The SARS-CoV-2 surface spike (S) glycoprotein represents a prime target for vaccine development because antibodies that block viral attachment and entry, i.e., neutralizing antibodies, bind almost exclusively to the receptor-binding domain. Here, we develop computational models for a large subset of S proteins associated with SARS-CoV-2, implemented through coarse-grained elastic network models and normal mode analysis. We then analyze local protein domain dynamics of the S protein systems and their thermal stability to characterize structural and dynamical variability among them. These results are compared against existing experimental data and used to elucidate the impact and mechanisms of SARS-CoV-2 S protein mutations and their associated antibody binding behavior. We construct a SARS-CoV-2 antigenic map and offer predictions about the neutralization capabilities of antibody and S mutant combinations based on protein dynamic signatures. We then compare SARS-CoV-2 S protein dynamics to SARS-CoV and MERS-CoV S proteins to investigate differing antibody binding and cellular fusion mechanisms that may explain the high transmissibility of SARS-CoV-2. The outbreaks associated with SARS-CoV, MERS-CoV, and SARS-CoV-2 over the last two decades suggest that the threat presented by coronaviruses is ever-changing and long term. Our results provide insights into the dynamics-driven mechanisms of immunogenicity associated with coronavirus S proteins and present a new, to our knowledge, approach to characterize and screen potential mutant candidates for immunogen design, as well as to characterize emerging natural variants that may escape vaccine-induced antibody responses.

Lysine Side-Chain Dynamics in the Binding Site of Homeodomain/DNA Complexes As Observed by NMR Relaxation Experiments and Molecular Dynamics Simulations

An important but poorly characterized contribution to the thermodynamics of protein-DNA interactions is the loss of entropy that occurs from restricting the conformational freedom of amino acid side chains. The effect of restricting the flexibility of several side chains at a protein-DNA interface may be comparable in many cases to the other factors that determine the binding thermodynamics and may, therefore, play a key role in dictating the binding affinity and/or specificity. Because the entropic contributions, including the presence and influence of side-chain dynamics, are especially difficult to estimate based on structural information, it is important to pursue experimental and theoretical studies that can provide direct information regarding these issues. We report on studies of a model system, the homeodomain/DNA complex, focusing on the Lys50 class of homeodomains where a key lysine residue in position 50 was shown previously to be critical for binding site specificity. NMR methodology was employed for determining the dynamics of lysine side-chain amino groups via 15N relaxation measurements in the Lys50-class homeodomains from the Drosophila protein Bicoid and the human protein Pitx2. In the case of Pitx2, complexes with both a consensus and a nonconsensus DNA binding site were examined. NMR-derived order parameters indicated moderate to substantial conformational freedom for the lysine NH3+ group in the complexes studied. To complement the experimental NMR measurements, molecular dynamics simulations were performed for the consensus complexes to gain further, detailed insights regarding the dynamics of the Lys50 side chain and other important residues in the protein-DNA interface.

Wing 1 of protein HOP2 is as important as helix 3 in DNA binding by MD simulation

The repair of programmed DNA double-strand breaks through recombination is required for proper association and disjunction of the meiotic homologous chromosomes. Meiosis-specific protein HOP2 plays essential roles in recombination by promoting recombinase activities. The N-terminal domain of HOP2 interacts with DNA through helix 3 (H3) and wing 1 (W1). Mutations in wing 1 (Y65A/K67A/Q68A) slightly weakened the binding but mutations in helices 2 and 3 (Q30A/K44A/K49A) nearly abolished the binding. To better understand such differential effects at atomic level, molecular dynamics simulations were employed. Despite losing some hydrogen bonds, the W1-mutant DNA complex was rescued by stronger hydrophobic interactions. For the wild type and W1-mutant, the protein was found to slide along the DNA grooves as the DNA rolls along its double-helix axis. This motion could be functionally important to facilitate the precise positioning of the single-stranded DNA with the homologous double-stranded DNA. The sliding motion was reduced in the W1-mutant. The H-mutant nearly lost all intermolecular interactions. Moreover, an additional mutation in wing 1 (Y65A/K67A/Q68A/K69A) also caused complete complex dissociation. Therefore, both wing 1 and helix 3 make important contribution to the DNA binding, which could be important to the strand invasion function of HOP2 homodimer and HOP2-MND1 heterodimer. Similar to cocking a medieval crossbow with the archer's foot placed in the stirrup, wing 1 may push the minor groove to cause distortion while helix 3 grabs the major groove.

Structural and evolutionary analyses reveal determinants of DNA binding specificities of nucleoid-associated proteins HU and IHF

Nucleoid-associated proteins (NAPs) are chromosome-organizing factors, which affect the transcriptional landscape of a bacterial cell. HU is an NAP, which binds to DNA with a broad specificity while homologous IHF (Integration Host Factor), binds DNA with moderately higher specificity. Specificity and differential binding affinity of HU/IHF proteins towards their target binding sites play a crucial role in their regulatory dynamics. Decades of biochemical and genomic studies have been carried out for HU and IHF like proteins. Yet, questions related to their DNA binding specificity, and differential ability to bend DNA thus affecting the binding site length remained unanswered. In addition, the problem has not been investigated from an evolutionary perspective. Our phylogenetic analysis revealed three major clades belonging to HU, IHFα and IHFβ like proteins with reference to E. coli. We carried out a comparative analysis of three-dimensional structures of HU/IHF proteins to gain insight into the structural basis of clade division. The present study revealed three major features which contribute to differential DNA binding specificity of HU/IHF proteins, (I) conformational restriction of DNA binding residues due to salt-bridge formation, (II) the enrichment of alanine in the DNA binding site increasing conformational space of flexible side chains in its vicinity and (III) nature of DNA binding residue (Arg to Lys bias in different clades) which interacts differentially to DNA bases. We observed an extended electropositive surface at the DNA draping site for IHF clade proteins compared to HU, which stabilizes the DNA bend. Differences in the dimer stabilization strategies between HU and IHF were also observed. Our analysis reveals a comprehensive evolutionary picture, which rationalizes the origin of multi-specificity of HU/IHF proteins using sequence and structure-based determinants, which could also be applied to understand differences in binding specificities of other nucleic acid binding proteins.

Theoretical design of a new chimeric protein for the treatment of breast cancer

p28 and NRC peptides are two anticancer peptides with various mechanisms have shown to be effective against breast cancer. Therefore, it seems that construction of a chimeric protein containing the two peptides might cause synergistic cytotoxic effects. However, since the two peptides bear opposite charges, production of a chimeric protein in which the two moieties do not intervene each other is difficult. In this study, our goal was to find a suitable peptide linker for the new chimeric protein in a manner that none of the peptides intervene the other's function. We selected some linkers with different characteristics and lengths and created a small library of the chimeric proteins harboring these linkers. Homology modeling and molecular dynamic simulation revealed that (PA)5P and (EAAAK)3 linkers can separate the p28 and NRC peptides effectively. Thus, the chimeric protein linked with (PA)5P or (EAAAK)3 linkers might show synergistic and stronger anticancer effects than the separate peptide moieties because they could exert their cytotoxic effects freely which is not influenced by the other part.

Interplay of Electrostatics and Hydrophobic Effects in the Metamorphic Protein Human Lymphotactin

The human lymphotactin (hLtn) is a protein that features two native states both of which are physiologically relevant: it is a monomer (hLtn10) at 10 °C with 200 mM salt and a dimer (hLtn40) at 40 °C and without salt. Here we focus on the networks of electrostatic and hydrophobic interactions that display substantial changes upon the conversion from hLtn10 to hLtn40 since they are expected to modulate the relative stability of the two folds. In addition to the Arg 23-Arg 43 interaction discussed in previous work, we find several other like-charge pairs that are likely important to the stability of hLtn10. Free energy perturbation calculations are carried out to explicitly evaluate the contribution of the Arg 23-Arg 43 interaction to the hLtn10 stability. hLtn40 features a larger number of salt bridges, and a set of hydrophobic residues undergo major changes in the solvent accessible surface area between hLtn10 and hLtn40, pointing to their importance to the relative stability of the two folds. We also discuss the use of explicit and implicit solvent simulations for characterizing the conformational ensembles under different solution conditions.

β-Arm flexibility of HU fromStaphylococcus aureusdictates the DNA-binding and recognition mechanism

HU, one of the major nucleoid-associated proteins, interacts with the minor groove of DNA in a nonspecific manner to induce DNA bending or to stabilize bent DNA. In this study, crystal structures are reported for both free HU fromStaphylococcus aureusMu50 (SHU) and SHU bound to 21-mer dsDNA. The structures, in combination with electrophoretic mobility shift assays (EMSAs), isothermal titration calorimetry (ITC) measurements and molecular-dynamics (MD) simulations, elucidate the overall and residue-specific changes in SHU upon recognizing and binding to DNA. Firstly, structural comparison showed the flexible nature of the β-sheets of the DNA-binding domain and that the β-arms bend inwards upon complex formation, whereas the other portions are nearly unaltered. Secondly, it was found that the disruption and formation of salt bridges accompanies DNA binding. Thirdly, residue-specific free-energy analyses using the MM-PBSA method with MD simulation data suggested that the successive basic residues in the β-arms play a central role in recognizing and binding to DNA, which was confirmed by the EMSA and ITC analyses. Moreover, residue Arg55 resides in the hinge region of the flexible β-arms, exhibiting a remarkable role in their flexible nature. Fourthly, EMSAs with various DNAs revealed that SHU prefers deformable DNA. Taken together, these data suggest residue-specific roles in local shape and base readouts, which are primarily mediated by the flexible β-arms consisting of residues 50–80.

Salt bridge exchange binding mechanism between streptavidin and its DNA aptamer--thermodynamics and spectroscopic evidences

Protein-nucleic acids binding driven by electrostatic interactions typically are characterized by the release of counter ions, and the salt-inhibited binding association constant (K(a)) and the magnitude of exothermic binding enthalpy (ΔH). Here, we report a non-classical thermodynamics of streptavidin (SA)-aptamer binding in NaCl (140-350 mM) solutions near room temperatures (23-27 °C). By using isothermal titration calorimetry (ITC) and circular dichroism (CD)/fluorescence spectroscopy, we found that the binding was enthalpy driven with a large entropy cost (ΔH -20.58 kcal mol(-1), TΔS -10.99 kcal mol(-1), and K(a) 1.08 × 10(7)  M(-1) at 140 mM NaCl 25 °C). With the raise of salt concentrations, the ΔH became more exothermic, yet the K(a) was almost unchanged (ΔH -26.29 kcal mol(-1) and K(a) 1.50 × 10(7)  M(-1) at 350 mM NaCl 25 °C). The data suggest that no counter Na(+) was released in the binding. Spectroscopy data suggest that the binding, with a stoichiometry of 2, was accompanied with substantial conformational changes on SA, and the changes were insensitive to the variation of salt concentrations. To account for the non-classical results, we propose a salt bridge exchange model. The intramolecular binding-site salt bridge(s) of the free SA and the charged phosphate group of aptamers re-organize to form the binding complex by forming a new intermolecular salt bridge(s). The salt bridge exchange binding process requires minimum amount of counter ions releasing but dehydration of the contacting surface of SA and the aptamer. The energy required for dehydration is reduced in the case of binding solution with higher salt concentration and account for the higher binding exothermic mainly.

A simple but effective modeling strategy for structural properties of non-heme Fe(II) sites in proteins: test of force field models and application to proteins in the AlkB family

To facilitate computational study of proteins in the AlkB family and related α-ketoglutarate/Fe(II)-dependent dioxygenases, we have tested a simple modeling strategy for the non-heme Fe(II) site in which the iron is represented by a simple +2 point charge with Lennard-Jones parameters. Calculations for an AlkB active site model in the gas phase and ∼150 ns molecular dynamics (MD) simulations for two enzyme-dsDNA complexes (E. coli AlkB-dsDNA and ABH2-dsDNA) suggest that this simple modeling strategy provides a satisfactory description of structural properties of the Fe(II) site in AlkB enzymes, provided that care is exercised to control the binding mode of carboxylate (Asp) to the iron. MD simulations using the model for AlkB-dsDNA and ABH2-dsDNA systems find that although the structural features for the latter are overall in good agreement with the crystal structure, the dsDNA, and AlkB-dsDNA interface undergo substantial changes during the MD simulations from the crystal structure. Even for ABH2, new interactions form between a long loop region and dsDNA upon structural relaxation of the loop, supporting the role of this loop in DNA binding despite the lack of interactions between them in the crystal structure. Analysis of DNA backbone torsional distributions helps identify regions that adopt strained conformations. Collectively, the results highlight that crystal packing may have a significant impact on the structure of protein-DNA complexes; the simulations also provide additional insights regarding why AlkB and ABH2 prefer single-strand and double-strand DNA, respectively, as substrate.

Investigating the structural impacts of I64T and P311S mutations in APE1-DNA complex: a molecular dynamics approach

Background

Elucidating the molecular dynamic behavior of Protein-DNA complex upon mutation is crucial in current genomics. Molecular dynamics approach reveals the changes on incorporation of variants that dictate the structure and function of Protein-DNA complexes. Deleterious mutations in APE1 protein modify the physicochemical property of amino acids that affect the protein stability and dynamic behavior. Further, these mutations disrupt the binding sites and prohibit the protein to form complexes with its interacting DNA.

Principal Findings

In this study, we developed a rapid and cost-effective method to analyze variants in APE1 gene that are associated with disease susceptibility and evaluated their impacts on APE1-DNA complex dynamic behavior. Initially, two different in silico approaches were used to identify deleterious variants in APE1 gene. Deleterious scores that overlap in these approaches were taken in concern and based on it, two nsSNPs with IDs rs61730854 (I64T) and rs1803120 (P311S) were taken further for structural analysis.

Significance

Different parameters such as RMSD, RMSF, salt bridge, H-bonds and SASA applied in Molecular dynamic study reveals that predicted deleterious variants I64T and P311S alters the structure as well as affect the stability of APE1-DNA interacting functions. This study addresses such new methods for validating functional polymorphisms of human APE1 which is critically involved in causing deficit in repair capacity, which in turn leads to genetic instability and carcinogenesis.