Water Dynamics Around Proteins: T- and R-States of Hemoglobin and Melittin

Cefmetazole sodium as an allosteric effector that regulates the oxygen supply efficiency of adult hemoglobin

Allosteric effectors play an important role in regulating the oxygen supply efficiency of hemoglobin for blood storage and disease treatment. However, allosteric effectors that are approved by the US FDA are limited. In this study, cefmetazole sodium (CS) was found to bind adult hemoglobin (HbA) from FDA library (1338 compounds) using surface plasmon resonance imaging high-throughput screening. Using surface plasmon resonance (SPR), the interaction between CS and HbA was verified. The oxygen dissociation curve of HbA after CS interaction showed a significant increase in P50 and theoretical oxygen-release capacity. Acid-base sensitivity (SI) exhibited a decreasing trend, although not significantly different. An oxygen dissociation assay indicated that CS accelerated HbA deoxygenation. Microfluidic modulated spectroscopy showed that CS changed the ratio of the alpha-helix to the beta-sheet of HbA. Molecular docking suggested CS bound to HbA's β-chains via hydrogen bonds, with key amino acids being N282, K225, H545, K625, K675, and V544.The results of molecular dynamics simulations (MD) revealed a stable orientation of the HbA-CS complex. CS did not significantly affect the P50 of bovine hemoglobin, possibly due to the lack of Valβ1 and Hisβ2, indicating that these were the crucial amino acids involved in HbA's oxygen affinity. Competition between the 2,3-Diphosphoglycerate (2,3-DPG) and CS in the HbA interaction was also determined by SPR, molecular docking and MD. In summary, CS could interact with HbA and regulate the oxygen supply efficiency via forming stable hydrogen bonds with the β-chains of HbA, and showed competition with 2,3-DPG.Communicated by Ramaswamy H. Sarma.

SCN as a local probe of protein structural dynamics

The dynamics of lysozyme is probed by attaching -SCN to all alanine residues. The one-dimensional infrared spectra exhibit frequency shifts in the position of the maximum absorption of 4 cm-1, which is consistent with experiments in different solvents and indicates moderately strong interactions of the vibrational probe with its environment. Isotopic substitution 12C → 13C leads to a redshift by -47 cm-1, which agrees quantitatively with experiments for CN-substituted copper complexes in solution. The low-frequency, far-infrared part of the protein spectra contains label-specific information in the difference spectra when compared with the wild type protein. Depending on the position of the labels, local structural changes are observed. For example, introducing the -SCN label at Ala129 leads to breaking of the α-helical structure with concomitant change in the far-infrared spectrum. Finally, changes in the local hydration of SCN-labeled alanine residues as a function of time can be related to the reorientation of the label. It is concluded that -SCN is potentially useful for probing protein dynamics, both in the high-frequency part (CN-stretch) and in the far-infrared part of the spectrum.

The current landscape of the antimicrobial peptide melittin and its therapeutic potential

Melittin, a main component of bee venom, is a cationic amphiphilic peptide with a linear α-helix structure. It has been reported that melittin can exert pharmacological effects, such as antitumor, antiviral and anti-inflammatory effects in vitro and in vivo. In particular, melittin may be beneficial for the treatment of diseases for which no specific clinical therapeutic agents exist. Melittin can effectively enhance the therapeutic properties of some first-line drugs. Elucidating the mechanism underlying melittin-mediated biological function can provide valuable insights for the application of melittin in disease intervention. However, in melittin, the positively charged amino acids enables it to directly punching holes in cell membranes. The hemolysis in red cells and the cytotoxicity triggered by melittin limit its applications. Melittin-based nanomodification, immuno-conjugation, structural regulation and gene technology strategies have been demonstrated to enhance the specificity, reduce the cytotoxicity and limit the off-target cytolysis of melittin, which suggests the potential of melittin to be used clinically. This article summarizes research progress on antiviral, antitumor and anti-inflammatory properties of melittin, and discusses the strategies of melittin-modification for its future potential clinical applications in preventing drug resistance, enhancing the selectivity to target cells and alleviating cytotoxic effects to normal cells.

Quantifying the Molecular Polarization Response of Liquid Water Interfaces at Heterogeneously Charged Surfaces

The hydration shells of proteins mediate interactions, such as small molecule binding, that are vital to their biological function or in some cases their dysfunction. However, even when the structure of a protein is known, the properties of its hydration environment cannot be easily predicted due to the complex interplay between protein surface heterogeneity and the collective structure of water's hydrogen bonding network. This manuscript presents a theoretical study of the influence of surface charge heterogeneity on the polarization response of the liquid water interface. We focus our attention on classical point charge models of water, where the polarization response is limited to molecular reorientation. We introduce a new computational method for analyzing simulation data that is capable of quantifying water's collective polarization response and determining the effective surface charge distribution of hydrated surfaces over atomistic length scales. To illustrate the utility of this method, we present the results of molecular dynamics simulations of liquid water in contact with a heterogeneous model surface and the CheY protein.

Local Hydration Control and Functional Implications Through S-Nitrosylation of Proteins: Kirsten Rat Sarcoma Virus (K-RAS) and Hemoglobin (Hb)

S-nitrosylation, the covalent addition of NO to the thiol side chain of cysteine, is an important post-transitional modification (PTM) that can affect the function of proteins. As such, PTMs extend and diversify protein function and thus characterizing consequences of PTM at a molecular level is of great interest. Although PTMs can be detected through various direct/indirect methods, they lack the capability to investigate the modifications with molecular detail. In the present work local and global structural dynamics, their correlation, the hydration structure, and the infrared spectroscopy for WT and S-nitrosylated Kirsten rat sarcoma virus (K-RAS) and hemoglobin (Hb) are characterized from molecular dynamics simulations. It is found that attaching NO to Cys118 in K-RAS rigidifies the protein in the Switch-I region which has functional implications, whereas for Hb, nitrosylation at Cys93 at the β₁ chain increases the flexibility of secondary structural motives for Hb in its T₀ and R₄ conformational substates. Solvent water access decreased by 40% after nitrosylation in K-RAS, similar to Hb for which, however, local hydration of the R₄SNO state is yet lower than for T₀SNO. Finally, S-nitrosylation leads to detectable peaks for the NO stretch frequency, but the congested IR spectral region will make experimental detection of these bands difficult. Overall, S-nitrosylation in these two proteins is found to influence hydration, protein flexibility, and conformational dynamics which are all eventually involved in protein regulation and function at a molecular level.

Water dynamics around T0 vs R4 of hemoglobin from local hydrophobicity analysis

The local hydration around tetrameric hemoglobin (Hb) in its T₀ and R₄ conformational substates is analyzed based on molecular dynamics simulations. Analysis of the local hydrophobicity (LH) for all residues at the α₁β₂ and α₂β₁ interfaces, responsible for the quaternary T → R transition, which is encoded in the Monod-Wyman-Changeux model, as well as comparison with earlier computations of the solvent accessible surface area, makes clear that the two quantities measure different aspects of hydration. Local hydrophobicity quantifies the presence and structure of water molecules at the interface, whereas "buried surface" reports on the available space for solvent. For simulations with Hb frozen in its T₀ and R₄ states, the correlation coefficient between LH and buried surface is 0.36 and 0.44, respectively, but it increases considerably if the 95% confidence interval is used. The LH with Hb frozen and flexible changes little for most residues at the interfaces but is significantly altered for a few select ones: Thr41α, Tyr42α, Tyr140α, Trp37β, Glu101β (for T₀) and Thr38α, Tyr42α, Tyr140α (for R₄). The number of water molecules at the interface is found to increase by ∼25% for T₀ → R₄, which is consistent with earlier measurements. Since hydration is found to be essential to protein function, it is clear that hydration also plays an essential role in allostery.

Cross-Correlated Motions in Azidolysozyme

The changes in the local and global dynamics of azide-labelled lysozyme compared with that of the wild type protein are quantitatively assessed for all alanine residues along the polypeptide chain. Although attaching -N3 to alanine residues has been considered to be a minimally invasive change in the protein it is found that depending on the location of the alanine residue, the local and global changes in the dynamics differ. For Ala92, the change in the cross-correlated motions are minimal, whereas attaching -N3 to Ala90 leads to pronounced differences in the local and global correlations as quantified by the cross-correlation coefficients of the Cα atoms. We also demonstrate that the spectral region of the asymmetric azide stretch distinguishes between alanine attachment sites, whereas changes in the low frequency, far-infrared region are less characteristic.

Dynamics of camel and human hemoglobin revealed by molecular simulations

Hemoglobin is one of the most widely studied proteins genetically, biochemically, and structurally. It is an oxygen carrying tetrameric protein that imparts the characteristic red color to blood. Each chain of hemoglobin harbors a heme group embedded in a hydrophobic pocket. Several studies have investigated structural variations present in mammalian hemoglobin and their functional implications. However, camel hemoglobin has not been thoroughly explored, especially from a structural perspective. Importantly, very little is known about how the heme group interacts with hemoglobin under varying conditions of osmolarity and temperature. Several experimental studies have indicated that the tense (T) state is more stable than the relaxed (R) state of hemoglobin under normal physiological conditions. Despite the fact that R state is less stable than the T state, no extensive structural dynamics studies have been performed to investigate global quaternary transitions of R state hemoglobin under normal physiological conditions. To evaluate this, several 500 ns all-atom molecular dynamics simulations were performed to get a deeper understanding of how camel hemoglobin behaves under stress, which it is normally exposed to, when compared to human hemoglobin. Notably, camel hemoglobin was more stable under physiological stress when compared to human hemoglobin. Additionally, when compared to camel hemoglobin, cofactor-binding regions of hemoglobin also exhibited more fluctuations in human hemoglobin under the conditions studied. Several differences were observed between the residues of camel and human hemoglobin that interacted with heme. Importantly, distal residues His58 of α hemoglobin and His63 of β hemoglobin formed more sustained interactions, especially at higher temperatures, in camel hemoglobin. These residues are important for oxygen binding to hemoglobin. Thus, this work provides insights into how camel and human hemoglobin differ in their interactions under stress.

The functional role of the hemoglobin-water interface

The interface between hemoglobin (Hb) and its environment, in particular water, is of great physiological relevance. Here, results from in vitro, in vivo, and computational experiments (molecular dynamics simulations) are summarized and put into perspective. One of the main findings from the computations is that the stability of the deoxy, ligand-free T-state (T₀) can be stabilized relative to the deoxy R-state (R₀) only in sufficiently large simulation boxes for the hydrophobic effect to manifest itself. This effect directly influences protein stability and is operative also under physiological conditions. Furthermore, molecular simulations provide a dynamical interpretation of the Perutz model for Hb function. Results from experiments using higher protein concentrations and realistic cellular environments are also discussed. One of the next great challenges for computational studies, which as we show is likely to be taken up in the near future, is to provide a molecular-level understanding of the dynamics of proteins in such crowded environments.

Molecular Dynamics of Hemoglobin Reveals Structural Alterations and Explains the Interactions Driving Sickle Cell Fibrillation

In sickle cell anemia, deoxyhemoglobin deforms RBCs by forming fibrils inside that disintegrate on oxygenation. We studied 100 ns long all-atom molecular dynamics (MD) for sickle and normal hemoglobin fibril models to understand this process, complemented by multiple 1 μs MD for a single tetramer of sickle and normal hemoglobin in deoxy and oxy states. We find that the presence of hydrophobic residues without a bulky side chain at β-6 in hemoglobin is the reason for the stability of the fibrils. Moreover, the free energy landscapes from MD of hemoglobin starting in the tensed (T) state capture the putative transition from T to relaxed (R) state, associated with oxygen binding. The three conformational wells in the landscapes are characterized by the quaternary changes where one αβ dimer rotates with respect to the other. The conformational changes from the oxygenation of sickle hemoglobin hinder the intermolecular contacts necessary for fibril formation.

Local Mutations Can Serve as a Game Changer for Global Protein Solvent Interaction

Although it is well-known that limited local mutations of enzymes, such as matrix metalloproteinases (MMPs), may change enzyme activity by orders of magnitude as well as its stability, the completely rational design of proteins is still challenging. These local changes alter the electrostatic potential and thus local electrostatic fields, which impacts the dynamics of water molecules close the protein surface. Here we show by a combined computational design, experimental, and molecular dynamics (MD) study that local mutations have not only a local but also a global effect on the solvent: In the specific case of the matrix metalloprotease MMP14, we found that the nature of local mutations, coupled with surface morphology, have the ability to influence large patches of the water hydrogen-bonding network at the protein surface, which is correlated with stability. The solvent contribution can be experimentally probed via terahertz (THz) spectroscopy, thus opening the door to the exciting perspective of rational protein design in which a systematic tuning of hydration water properties allows manipulation of protein stability and enzymatic activity.

Roles of Fe-Histidine bonds in stability of hemoglobin: Recognition of protein flexibility by Q Sepharose

Using various mutants, we investigated to date the roles of the Fe-histidine (F8) bonds in cooperative O₂ binding of human hemoglobin (Hb) and differences in roles between α- and β-subunits in the α₂β₂ tetramer. An Hb variant with a mutation in the heme cavity exhibited an unexpected feature. When the β mutant rHb (βH92G), in which the proximal histidine (His F8) of the β-subunit is replaced by glycine (Gly), was subjected to ion-exchange chromatography (Q Sepharose column) and eluted with an NaCl concentration gradient in the presence of imidazole, yielded two large peaks, whereas the corresponding α-mutant, rHb (αH87G), gave a single peak similar to Hb A. The β-mutant rHb proteins under each peak had identical isoelectric points according to isoelectric focusing electrophoresis. Proteins under each peak were further characterized by Sephadex G-75 gel filtration, far-UV CD, ¹H NMR, and resonance Raman spectroscopy. We found that rHb (βH92G) exists as a mixture of αβ-dimers and α₂β₂ tetramers, and that hemes are released from β-subunits in a fraction of the dimers. An approximate amount of released hemes were estimated to be as large as 30% with Raman relative intensities. It is stressed that Q Sepharose columns can distinguish differences in structural flexibility of proteins having identical isoelectric points by altering the exit rates from the porous beads. Thus, the role of Fe-His (F8) bonds in stabilizing the Hb tetramer first described by Barrick et al. was confirmed in this study. In addition, it was found in this study that a specific Fe-His bond in the β-subunit minimizes globin structural flexibility.

Uncertainty quantification in classical molecular dynamics

Molecular dynamics simulation is now a widespread approach for understanding complex systems on the atomistic scale. It finds applications from physics and chemistry to engineering, life and medical science. In the last decade, the approach has begun to advance from being a computer-based means of rationalizing experimental observations to producing apparently credible predictions for a number of real-world applications within industrial sectors such as advanced materials and drug discovery. However, key aspects concerning the reproducibility of the method have not kept pace with the speed of its uptake in the scientific community. Here, we present a discussion of uncertainty quantification for molecular dynamics simulation designed to endow the method with better error estimates that will enable it to be used to report actionable results. The approach adopted is a standard one in the field of uncertainty quantification, namely using ensemble methods, in which a sufficiently large number of replicas are run concurrently, from which reliable statistics can be extracted. Indeed, because molecular dynamics is intrinsically chaotic, the need to use ensemble methods is fundamental and holds regardless of the duration of the simulations performed. We discuss the approach and illustrate it in a range of applications from materials science to ligand-protein binding free energy estimation. This article is part of the theme issue 'Reliability and reproducibility in computational science: implementing verification, validation and uncertainty quantification in silico'.

Site-selective dynamics of azidolysozyme

The spectroscopic response of and structural dynamics around all azido-modified alanine residues (AlaN₃) in lysozyme are characterized. It is found that AlaN₃ is a positionally sensitive probe for the local dynamics, covering a frequency range of ∼15 cm^-1 for the center frequency of the line shape. This is consistent with findings from selective replacements of amino acids in PDZ2, which reported a frequency span of ∼10 cm^-1 for replacements of Val, Ala, or Glu by azidohomoalanine. For the frequency fluctuation correlation functions, the long-time decay constants τ₂ range from ∼1 to ∼10 ps, which compares with experimentally measured correlation times of 3 ps. Attaching azide to alanine residues can yield dynamics that decays to zero on the few ps time scale (i.e., static component Δ₀ ∼ 0 ps^-1) or to a remaining, static contribution of ∼0.5 ps^-1 (corresponding to 2.5 cm^-1), depending on the local environment on the 10 ps time scale. The magnitude of the static component correlates qualitatively with the degree of hydration of the spectroscopic probe. Although attaching azide to alanine residues is found to be structurally minimally invasive with respect to the overall protein structure, analysis of the local hydrophobicity indicates that the hydration around the modification site differs for modified and unmodified alanine residues, respectively.

Modulation of Nanoparticle Diffusion by Surface Ligand Length and Charge: Analysis with Molecular Dynamics Simulations

To help better interpret experimental measurement of nanoparticle size, it is important to understand how their diffusion depends on the physical and chemical features of surface ligands. In this study, explicit solvent molecular dynamics simulations are used to probe the effect of ligand charge and flexibility on the diffusion of small gold nanoparticles. The results suggest that despite a high bare charge (+18 e), cationic nanoparticles studied here have reduced diffusion constants compared to a hydrophobic gold nanoparticle by merely a modest amount. Increasing the ligand length by 10 CH₂ units also has a limited impact on the diffusion constant. For the three particles studied here, the difference between estimated hydrodynamic radius and radius of gyration is on the order of one solvent layer (3-5 Å), confirming that the significant discrepancies found in the size of similar nanoparticles by recent transmission electron microscopy and dynamic light scattering measurements were due to aggregation under solution conditions. The limited impact of electrostatic friction on the diffusion of highly charged nanoparticles is found to be due to the strong anticorrelation between electrostatic and van der Waals forces between nanoparticle and environment, supporting the generality of recent observation for proteins by Matyushov and co-workers. Including the first shell of solvent molecules as part of the diffusing particle has a minor impact on the total force autocorrelation function but reduces the disparity in relaxation time between the total force and its electrostatic and van der Waals components.

On the importance of statistics in molecular simulations for thermodynamics, kinetics and simulation box size

Computational simulations, akin to wetlab experimentation, are subject to statistical fluctuations. Assessing the magnitude of these fluctuations, that is, assigning uncertainties to the computed results, is of critical importance to drawing statistically reliable conclusions. Here, we use a simulation box size as an independent variable, to demonstrate how crucial it is to gather sufficient amounts of data before drawing any conclusions about the potential thermodynamic and kinetic effects. In various systems, ranging from solvation free energies to protein conformational transition rates, we showcase how the proposed simulation box size effect disappears with increased sampling. This indicates that, if at all, the simulation box size only minimally affects both the thermodynamics and kinetics of the type of biomolecular systems presented in this work.