Inhibition of staphylococcal nuclease by benzimidazole-based Ligand: Implications in DNA-Mediated entrapment and uptake of MRSA by Macrophage-like cells.

The staphylococcal nuclease also referred as micrococcal nuclease (MNase) is a key drug target as the enzyme degrades the neutrophil extracellular trap (NET) and empowers the pathogen to subvert the host innate immune system. To this end, the current study presents a critical evaluation of MNase inhibition rendered by benzimidazole-based ligands (C1 and C2) and probes its therapeutic implications. A nuclease assay indicated that MNase inhibition rendered by C1 and C2 was âˆ¼ 55 % and âˆ¼ 72 %, respectively, at the highest tested concentration of 10 µM. Studies on enzyme kinetics revealed that C2 rendered non-competitive inhibition and significantly reduced MNase turnover number (Kcat) and catalytic efficiency (Kcat/Km) with an IC50 value of âˆ¼ 1122 nM. In CD spectroscopy, a notable perturbation in the ß-sheet content of MNase was observed in presence of C2. Fluorescence-microscope analysis indicated that MNase inhibition by C2 could restore entrapment of methicillin-resistant Staphylococcus aureus (MRSA) in calf-thymus DNA (CT-DNA). Flow cytometry and confocal microscope analysis revealed that uptake of DNA-entrapped MRSA by activated THP-1 cells was reinstated by MNase inhibition rendered by C2. Inhibition of nuclease by the non-toxic ligand C2 holds therapeutic prospect as it has the potential to bolster the DNA-mediated entrapment machinery and mitigate MRSA infections.

Extending the Stochastic Titration CpHMD to CHARMM36m.

The impact of pH on proteins is significant but often neglected in molecular dynamics simulations. Constant-pH Molecular Dynamics (CpHMD) is the state-of-the-art methodology to deal with these effects. However, it still lacks widespread adoption by the scientific community. The stochastic titration CpHMD is one of such methods that, until now, only supported the GROMOS force field family. Here, we extend this method's implementation to include the CHARMM36m force field available in the GROMACS software package. We test this new implementation with a diverse group of proteins, namely, lysozyme, Staphylococcal nuclease, and human and E. coli thioredoxins. All proteins were conformationally stable in the simulations, even at extreme pH values. The RMSE values (pKa prediction vs experimental) obtained were very encouraging, in particular for lysozyme and human thioredoxin. We have also identified a few residues that challenged the CpHMD simulations, highlighting scenarios where the method still needs improvement independently of the force field. The CHARMM36m all-atom implementation was more computationally efficient when compared with the GROMOS 54A7, taking advantage of a shorter nonbonded interaction cutoff and a less frequent neighboring list update. The new extension will allow the study of pH effects in many systems for which this force field is particularly suited, i.e., proteins, membrane proteins, lipid bilayers, and nucleic acids.

Electronic Polarization Is Essential for the Stabilization and Dynamics of Buried Ion Pairs in Staphylococcal Nuclease Mutants.

Buried charged residues play important roles in the modulation of protein stabilities and conformational dynamics and make crucial contributions to protein functions. Considering the generally nonpolar nature of protein interior, a key question concerns the contribution of electronic polarization to the stabilization and properties of buried charges. We answer this question by conducting free energy simulations using the latest polarizable CHARMM force field based on Drude oscillators for a series of Staphylococcal nuclease mutants that involve a buried Glu-Lys pair in different titration states and orientations. While a nonpolarizable model suggests that the ionized form of the buried Glu-Lys pair is more than 40 kcal/mol less stable than the charge-neutral form, the two titration states are comparable in stability when electronic polarization is included explicitly, a result better reconcilable with available experimental data. Analysis of free energy components suggests that additional stabilization of the ionized Glu-Lys pair has contributions from both the enhanced salt-bridge strength and stronger interaction between the ion-pair and surrounding protein residues and penetrated water. Despite the stronger direct interaction between Glu and Lys, the ion-pair exhibits considerably larger and faster structural fluctuations when polarization is included, due to compensation of interactions in the cavity. Collectively, observations from this work provide compelling evidence that electronic polarization is essential to the stability, hydration, dynamics, and therefore function of buried charges in proteins. Therefore, our study advocates for the explicit consideration of electronic polarization for mechanistic and engineering studies that implicate buried charged residues, such as enzymes and ion transporters.

Morphological, Gene, and Hormonal Changes in Gonads and In-Creased Micrococcal Nuclease Accessibility of Sperm Chromatin Induced by Mercury.

Mercury is one of the most dangerous environmental pollutants. In this work, we analysed the effects of exposure of Mytilus galloprovincialis to 1, 10 and 100 pM HgCl2 for 24 h on the gonadal morphology and on the expression level of three stress genes: mt10, hsp70 and πgst. In this tissue we also evaluated the level of steroidogenic enzymes 3ß-HSD and 17ß-HSD and the expression of PL protein genes. Finally, we determined difference in sperm chromatin accessibility to micrococcal nuclease. We found alterations in gonadal morphology especially after exposure to 10 and 100 pM HgCl2 and hypo-expression of the three stress genes, particularly for hsp70. Furthermore, decreased labelling with both 3ß-HSD and 17ß-HSD antibodies was observed following exposure to 1 and 10 pM HgCl2 and complete absence at 100 pM HgCl2 exposure. Gonads of mussels exposed to all HgCl2 doses showed decreased expression of PL protein genes especially for PLIII. Finally, micrococcal nuclease digestions showed that all doses of HgCl2 exposure resulted in increased sperm chromatin accessibility to this enzyme, indicative of improper sperm chromatin structure. All of these changes provide preliminary data of the potential toxicity of mercury on the reproductive health of this mussel.

Structural Evaluation of Protein/Metal Complexes via Native Electrospray Ultraviolet Photodissociation Mass Spectrometry.

Ultraviolet photodissociation (UVPD) has emerged as a promising tool to characterize proteins with regard to not only their primary sequences and post-translational modifications, but also their tertiary structures. In this study, three metal-binding proteins, Staphylococcal nuclease, azurin, and calmodulin, are used to demonstrate the use of UVPD to elucidate metal-binding regions via comparisons between the fragmentation patterns of apo (metal-free) and holo (metal-bound) proteins. The binding of staphylococcal nuclease to calcium was evaluated, in addition to a series of lanthanide(III) ions which are expected to bind in a similar manner as calcium. On the basis of comparative analysis of the UVPD spectra, the binding region for calcium and the lanthanide ions was determined to extend from residues 40-50, aligning with the known crystal structure. Similar analysis was performed for both azurin (interrogating copper and silver binding) and calmodulin (four calcium binding sites). This work demonstrates the utility of UVPD methods for determining and analyzing the metal binding sites of a variety of classes of proteins.

Molecular dynamics simulation of proteins under high pressure: Structure, function and thermodynamics.

BACKGROUND: Molecular dynamics (MD) simulation is well-recognized as a powerful tool to investigate protein structure, function, and thermodynamics. MD simulation is also used to investigate high pressure effects on proteins. For conducting better MD simulation under high pressure, the main issues to be addressed are: (i) protein force fields and water models were originally developed to reproduce experimental properties obtained at ambient pressure; and (ii) the timescale to observe the pressure effect is often much longer than that of conventional MD simulations. SCOPE OF REVIEW: First, we describe recent developments in MD simulation methodologies for studying the high-pressure structure and dynamics of protein molecules. These developments include force fields for proteins and water molecules, and enhanced simulation techniques. Then, we summarize recent studies of MD simulations of proteins in water under high pressure. MAJOR CONCLUSIONS: Recent MD simulations of proteins in solution under pressure have reproduced various phenomena identified by experiments using high pressure, such as hydration, water penetration, conformational change, helix stabilization, and molecular stiffening. GENERAL SIGNIFICANCE: MD simulations demonstrate differences in the properties of proteins and water molecules between ambient and high-pressure conditions. Comparing the results obtained by MD calculations with those obtained experimentally could reveal the mechanism by which biological molecular machines work well in collaboration with water molecules.

Toward the computational design of protein crystals with improved resolution.

Substantial advances have been made in the computational design of protein interfaces over the last 20 years. However, the interfaces targeted by design have typically been stable and high-affinity. Here, we report the development of a generic computational design method to stabilize the weak interactions at crystallographic interfaces. Initially, we analyzed structures reported in the Protein Data Bank to determine whether crystals with more stable interfaces result in higher resolution structures. We found that for 22 variants of a single protein crystallized by a single individual, the Rosetta-calculated `crystal score' correlates with the reported diffraction resolution. We next developed and tested a computational design protocol, seeking to identify point mutations that would improve resolution in a highly stable variant of staphylococcal nuclease (SNase). Using a protocol based on fixed protein backbones, only one of the 11 initial designs crystallized, indicating modeling inaccuracies and forcing us to re-evaluate our strategy. To compensate for slight changes in the local backbone and side-chain environment, we subsequently designed on an ensemble of minimally perturbed protein backbones. Using this strategy, four of the seven designed proteins crystallized. By collecting diffraction data from multiple crystals per design and solving crystal structures, we found that the designed crystals improved the resolution modestly and in unpredictable ways, including altering the crystal space group. Post hoc, in silico analysis of the three observed space groups for SNase showed that the native space group was the lowest scoring for four of six variants (including the wild type), but that resolution did not correlate with crystal score, as it did in the preliminary results. Collectively, our results show that calculated crystal scores can correlate with reported resolution, but that the correlation is absent when the problem is inverted. This outcome suggests that more comprehensive modeling of the crystallographic state is necessary to design high-resolution protein crystals from poorly diffracting crystals.

Accurate Calculation of Barnase and SNase Folding Energetics Using Short Molecular Dynamics Simulations and an Atomistic Model of the Unfolded Ensemble: Evaluation of Force Fields and Water Models.

As proteins perform most cellular functions, quantitative understanding of protein energetics is required to gain control of biological phenomena. Accurate models of native proteins can be obtained experimentally, but the lack of equally fine models of unfolded ensembles impedes the calculation of protein folding energetics from first principles. Here, we show that an atomistic unfolded ensemble model, consisting of a few dozen conformations built from a protein sequence, can be used in conjunction with an X-ray structure of its native state to calculate accurately by difference the changes in enthalpy and heat capacity of the polypeptide upon folding. The calculation is done using molecular dynamics simulations, popular force fields, and water models, and for the two model proteins studied (barnase and SNase), the results agree within error or are very close to their experimentally determined properties. The enthalpy sampling of the unfolded ensemble is done through short 2 ns simulations that do not significantly modify the representative distribution of Rg of the starting conformations. The impressive accuracy obtained opens the possibility to investigate quantitatively systems or phenomena not amenable to experiment and paves the way for addressing the calculation of protein conformational stability (i.e., the change in Gibbs energy upon folding), a central goal of structural biology. So far, these calculated enthalpy and heat capacity changes, combined with the experimentally determined melting temperatures of the corresponding protein, allow us to reproduce the stability curves of both barnase and SNase.

pH-Dependent Conformational Changes Due to Ionizable Residues in a Hydrophobic Protein Interior: The Study of L25K and L125K Variants of SNase.

Ionizable residues in the hydrophobic interior of certain proteins are known to play important roles in life processes like energy transduction and enzyme catalysis. These internal ionizable residues show experimental apparent pKa values having large shifts as compared to their values in solution. In the present work, we study the pH-dependent conformational changes undergone by two variants of staphylococcal nuclease (SNase), L25K and L125K, using pH replica exchange molecular dynamics (pH-REMD) in explicit solvent. Our results show that the observed pKa of Lys25 and Lys125 are significantly different than their pKa in solution. We observed that the internal lysine residues prefer to be water-exposed when protonated at low pH, but they remain buried within the hydrophobic pocket when deprotonated at high pH. Using thermodynamic laws, we estimate the microscopic conformation-specific pKa of the water-exposed and buried conformations of the internal lysine residues and explain their relation to the macroscopic observed pKa values. We present the differences in the microscopic mechanisms that lead to similar experimentally observed apparent pKa of Lys25 and Lys125, and explain the need of thermodynamic models of different complexities to account for our calculations. We see that L25K displays pH-dependent fluctuations throughout the entire ß barrel and the α1 helix. In contrast, pH-independent fluctuations are observed in L125K, primarily limited to the α3 helix. The present computational study offers a detailed atomistic understanding of the determinants of the observed anomalous pKa of internal ionizable residues, bolstering the experimental findings.

Universality and Structural Implications of the Boson Peak in Proteins.

The softness and rigidity of proteins are reflected in the structural dynamics, which are in turn affected by the environment. The characteristic low-frequency vibrational spectrum of a protein, known as boson peak, is an indication of the structural rigidity of the protein at a cryogenic temperature or dehydrated conditions. In this article, the effect of hydration, temperature, and pressure on the boson peak and volumetric properties of a globular protein are evaluated by using inelastic neutron scattering and molecular dynamics simulation. Hydration, pressurization, and cooling shift the boson peak position to higher energy and depress the peak intensity and decreases the protein and cavity volumes. We found the correlation between the boson peak and cavity volume in a protein. A decrease of cavity volume means the increase of rigidity, which is the origin of the boson peak shift. Boson peak is the universal property of a protein, which is rationalized by the correlation.

Single-cell chromatin immunocleavage sequencing (scChIC-seq) to profile histone modification.

Our method for analyzing histone modifications, scChIC-seq (single-cell chromatin immunocleavage sequencing), involves targeting of the micrococcal nuclease (MNase) to a histone mark of choice by tethering to a specific antibody. Cleaved target sites are then selectively PCR amplified. We show that scChIC-seq reliably detects H3K4me3 and H3K27me3 target sites in single human white blood cells. The resulting data are used for clustering of blood cell types.

Detection of micrococcal nuclease for identifying Staphylococcus aureus based on DNA templated fluorescent copper nanoclusters.

Micrococcal nuclease (MNase) is a naturally-secreted nucleic acid degrading enzyme with important role in the spread of the bacteria in an infected host. The content of MNase can be used to estimate the pathogenicity of Staphylococcus aureus (S. aureus). A fluorometric method is described here for determination of the activity of MNase and for identification of S. aureus using DNA templated fluorescent copper nanoclusters (CuNC). A double-stranded DNA (dsDNA) with AT-rich regions and protruding 3'-termini was identified as a high-selectivity substrate for MNase and as a template for CuNC. In the absence of MNase, the long AT-rich dsDNA templates the formation of CuNC that display bright yellow fluorescence, with excitation/emission peaks at 340/570 nm. However, the substrates are enzymatically digested to mononucleotides or short-oligonucleotide fragments, which fail to synthesize fluorescent CuNC. The method works in the 1.0 × 10-3 - 5.0 × 10-2 U mL-1 MNase activity range, has a 1.0 mU mL-1 detection limit, and is highly selective over other exonucleases. The assay was successfully applied to the detection of MNase secreted by S. aureus and to the identification of S. aureus. Graphical abstract A smart dsDNA, with AT-rich regions and 3'-protruding termini, is screened as the high-selectivity substrate for MNase and template for the formation of copper nanoclusters (CuNC). A method is described for determination of the activity of MNase and for identification of S. aureus via smart DNA templated formation of fluorescent CuNCs.

Theoretical Insights on Nonlinear Response Theory of Fluorescence Spectroscopy in Liquids.

Correlations between the nonequilibrium solvation dynamics upon the photon excitation of the chromophore and a system's equilibrium fluctuations are deeply studied. As the linear response of the solvent has been linked with Gaussian statistics of the energy fluctuations in the literature, we specifically explore the cases beyond the regime of the linear response theory due to deviation from Gaussian fluctuations. As a continuation of our previous work, an analytical formalism is presented to project the energy shift with various order moments, where the non-Gaussian statistics arise from the overlap of the energy basins on the perturbed potential energy surface. It is shown that the nonequilibrium dynamics still correlate with the spontaneous regressions at equilibrium and are controlled by the decay rates of those higher order components with the prevailing contributions to the energy shift. Molecular dynamics simulations were performed in the protein Staphylococcus nuclease, in which even the dynamics of the high order moments are available. The results further verify the above relationship. Our scheme is used to evaluate Stokes shift using the information on non-Gaussian statistics at equilibrium, thus presenting a broad picture on the correlation between the nonequilibrium process and equilibrium properties in liquids.

A Novel Method to Generate MNase Ladders Reveal Rapid Chromatin Remodeling upon Gametogenesis and Mating in Chlamydomonas.

To circumvent nuclei isolation for nucleosomal mapping of wild-type (cell walled) algal cells, we developed a quick and versatile methodology, by abrasion of whole cells (Chlamydomonas, Scenedesmus and yeast), allowing Micrococcal Nuclease (MNase) direct access to nuclear chromatin, in situ. Varying parameters such as bead abrasion, vortex and incubation conditions, we optimized capture of an 'early digest' which may probe chromatin differentially, based on nucleosome accessibility. A comparison of such ladders across vegetative cells, gametes and zygotes revealed an increase in the average nucleosomal repeat length (+17-34nt) upon gametogenesis, indicating a trend of chromatin compaction. Using PCR, we compared promoter enrichment in increasing orders of fractionated nucleosomal repeats (mono-, di-, up to penta-), each differing in cleavability based on chromatin accessibility. Concordant with higher gene expression (mating locus), promoters revealed an enrichment in mono-nucleosomal fractions. Interestingly, the zygote specific gene, MT0828 displayed rapid remodelling from penta-nucleosomal enrichment when completely repressed (vegetative), to intermediate states during gametogenesis (24hrs), which finally shifted to being largely mono-nucleosomal, when induced (1h zygotes). Summarizing three candidate genes from the mating locus, we conclude that the MNase based 'Chromatin Accessibility Assay' can track a range of large-scale rapid chromatin remodelling transitions within the binaries of gene expression.

Characterization of a nontypeable Haemophilus influenzae thermonuclease.

Nontypeable Haemophilus influenzae (NTHi) has been shown to form biofilms, comprised of extracellular DNA (eDNA), in the middle ear and bronchus during clinical infections. Studies in our laboratory have shown that NTHi possesses a homolog of Staphylococcus aureus thermonuclease (staphylococcal thermonuclease), NTHi nuclease (NTHi Nuc, HI_1296). This enzyme had similar size, heat stability, and divalent cation requirements to those of the staphylococcal homolog as determined by light scattering and circular dichroism spectroscopy. Small angle X-ray scattering (SAXS) analysis suggested an overall shape and substrate-binding site comparable to those of staphylococcal nuclease. However, NTHi Nuc was approximately 25-fold more active in fluorescence resonance energy transfer (FRET) activity assay than staphylococcal thermonuclease. Homology modeling implicates shorter NTHi Nuc loops near the active site for this enhanced activity.

Hydronium Ions Accompanying Buried Acidic Residues Lead to High Apparent Dielectric Constants in the Interior of Proteins.

Internal ionizable groups are known to play important roles in protein functions. A mystery that has attracted decades of extensive experimental and theoretical studies is the apparent dielectric constants experienced by buried ionizable groups, which are much higher than values expected for protein interiors. Many interpretations have been proposed, such as water penetration, conformational relaxation, local unfolding, protein intrinsic backbone fluctuations, etc. However, these interpretations conflict with many experimental observations. The virtual mixture of multiple states (VMMS) simulation method developed in our lab provides a direct approach for studying the equilibrium of multiple chemical states and can monitor p Ka values along simulation trajectories. Through VMMS simulations of staphylococcal nuclease (SNase) variants with internal Asp or Glu residues, we discovered that cations were attracted to buried deprotonated acidic groups and the presence of the nearby cations were essential to reproduce experimentally measured p Ka values. This finding, combined with structural analysis and validation simulations, suggests that the proton released from a deprotonation process stays near the deprotonated group inside proteins, possibly in the form of a hydronium ion. The existence of a proton near a buried charge has many implications in our understanding of protein functions.

Potential of Pyridine Amphiphiles as Staphylococcal Nuclease Inhibitor.

The present study explores the potential of pyridine-based synthetic amphiphiles C1 and C2 having 4-carbon and 12-carbon hydrophobic tails, respectively, as staphylococcal nuclease inhibitors. UV-visible titration with calf-thymus DNA (CT-DNA) revealed a hypochromic shift in the absorbance bands of C1 and C2, whereas fluorescence titration indicated a reduction in the emission intensity of the monomer bands of the amphiphiles. Interaction of deoxyribonuclease I (DNase 1) and micrococcal nuclease (MNase) with C1 or C2 led to a decrease in the emission intensity of tryptophan at λ=345 nm along with an increase in the monomer emission intensity of C1 and C2 at λ=375 nm for DNase I and excimer emission intensity at λ=470 nm for both DNase I and MNase. Scatchard's analysis indicated superior interaction of C2 with DNase I. Circular dichroism spectroscopy revealed major changes in the secondary structures of both DNase I and MNase upon interaction with the amphiphiles. A solution-based nuclease assay in conjunction with gel electrophoresis indicated amphiphile-mediated protection against nuclease-directed DNA cleavage. Interestingly, C2 could render inhibition of nuclease present in the culture supernatant of Staphylococcus aureus MTCC 96, which highlights the therapeutic prospect of the amphiphile against S. aureus.

Investigating Nucleosome Accessibility for MNase, FeII †Peplomycin, and Duocarmycin†B2 by Using Capillary Electrophoresis.

Capillary electrophoresis, coupled with DNA 5' Texas Red labeling, was used to investigate the ability of MNase, FeII peplomycin, and duocarmycin B2 to access the nucleosome. Distinct accessibility patterns of these species to the nucleosome were observed. MNase was completely prevented from approaching the nucleosome core and exhibited a higher site specificity for targeting DNA sites located close to the core region. Intercalation of peplomycin in the nucleosomal core region was highly suppressed, but reaction sites located at the ends of the nucleosomal core remained accessible, which implied flexibility of the core DNA end. Duocarmycin B2 was able to enter and react in the core region, although its alkylating efficiency decreased significantly.

A Coupled Ionization-Conformational Equilibrium Is Required To Understand the Properties of Ionizable Residues in the Hydrophobic Interior of Staphylococcal Nuclease.

Ionizable residues in the interior of proteins play essential roles, especially in biological energy transduction, but are relatively rare and seem incompatible with the complex and polar environment. We perform a comprehensive study of the internal ionizable residues on 21 variants of staphylococcal nuclease with internal Lys, Glu, or Asp residues. Using pH replica exchange molecular dynamics simulations, we find that, in most cases, the pKa values of these internal ionizable residues are shifted significantly from their values in solution. Our calculated results are in excellent agreement with the experimental observations of the Garcia-Moreno group. We show that the interpretation of the experimental pKa values requires the study of not only protonation changes but also conformational changes. The coupling between the protonation and conformational equilibria suggests a mechanism for efficient pH-sensing and regulation in proteins. This study provides new physical insights into how internal ionizable residues behave in the hydrophobic interior of proteins.

Anomalous Properties of Lys Residues Buried in the Hydrophobic Interior of a Protein Revealed with 15N-Detect NMR Spectroscopy.

Ionizable residues buried in hydrophobic environments in proteins are essential for many fundamental biochemical processes. These residues titrate with anomalous pKa values that are challenging to reproduce with structure-based calculations owing to the conformational reorganization coupled to their ionization. Detailed characterization of this conformational reorganization is of interest; unfortunately, the properties of buried Lys residues are difficult to study experimentally. Here we demonstrate the utility of 15N NMR spectroscopy to gain insight into the protonation state, state of hydration and conformational dynamics of the Nζ amino group of buried Lys residues. The experiments were applied to five variants of staphylococcal nuclease, with internal Lys residues that titrate with pKa values ranging from 6.2 to 8.1. Direct detection of buried Lys residues with these NMR spectroscopy methods will enable correlation between thermodynamic and structural data as well as unprecedented examination of how conformational transitions coupled to their ionization affect their pKa values.