Dynamics of water clusters confined in proteins: a molecular dynamics simulation study of interfacial waters in a dimeric hemoglobin

Dynamics of Hydrogen Bonds between Water and Intrinsically Disordered and Structured Regions of Proteins

Recent studies indicate more restricted dynamics of water around intrinsically disordered proteins (IDPs) than structured proteins. We examine here the dynamics of hydrogen bonds between water molecules and two proteins, small ubiquitin-related modifier-1 (SUMO-1) and ubiquitin-conjugating enzyme E2I (UBC9), which we compare around intrinsically disordered regions (IDRs) and structured regions of these proteins. It has been recognized since some time that excluded volume effects, which influence access of water molecules to hydrogen-bonding sites, and the strength of hydrogen bonds between water and protein affect hydrogen bond lifetimes. While we find those two properties to mediate lifetimes of hydrogen bonds between water and protein residues in this study, we also find that the lifetimes are affected by the concentration of charged groups on other nearby residues. These factors are more important in determining the hydrogen bond lifetimes than whether a residue hydrogen bonding with water belongs to an IDR or to a structured region.

Energy Transport in Class B GPCRs: Role of Protein-Water Dynamics and Activation

We compute energy exchange networks (EENs) through glucagon-like peptide-1 receptor (GLP-1R), a class B G-protein-coupled receptor (GPCR), in inactive and two active states, one activated by a peptide ligand and the other by a small molecule agonist, from results of molecular dynamics simulations. The reorganized network upon activation contains contributions from structural as well as from dynamic changes and corresponding entropic contributions to the free energy of activation, which are estimated in terms of the change in rates of energy transfer across non-covalent contacts. The role of water in the EENs and in activation of GLP-1R is also investigated. The dynamics of water in contact with the central polar network of the transmembrane region is found to be significantly slower for both activated states compared to the inactive state. This result is consistent with the contribution of water molecules to activation of GLP-1R previously suggested and resembles water dynamics in parts of the transmembrane region found in earlier studies of rhodopsin-like GPCRs.

Ionic Dynamics and Vibrational Spectral Diffusion of a Protic Alkylammonium Ionic Salt through Intrinsic Cationic N-H Vibrational Probe from FPMD Simulations

We employed density functional theory (DFT)-based molecular dynamics simulations to explore the structure, dynamics, and spectral properties of the protic ionic entity trimethylammonium chloride (TMACl). Structural investigations include calculating the site-site radial distribution functions (RDFs), the distribution of constituent cations and anions in three-dimensional space, and combined distribution functions of the hydrogen-bonded pair RDF versus angle, revealing the structural characteristics of the ionic solvation and the intermolecular interactions within ions. Further, we determined the instantaneous vibrational stretching frequencies of the intrinsic N-H stretch probe modes by applying the time-series wavelet method. The associated ionic dynamics within the protic ionic compound were investigated by calculating the time-evolution of the fluctuating frequencies and the frequency-time correlation functions (FFCFs). The time scale related to the local structural relaxation process and the average hydrogen bond lifetime, ion cage dynamics, and mean squared displacement were investigated. The faster decay component of the FFCFs, depicting the intermolecular motion of intact hydrogen bonds in TMACl, is 0.07 ps for the Perdew-Burke-Ernzerhof (PBE)-based simulation and 0.06 ps for the PBE-D2 representation. The slower time scale of the longer picosecond decay time component of PBE and PBE-D2 representations are 3.13 and 2.87 ps, respectively. These picosecond time scales represent more significant fluctuations of the hydrogen-bonding partners in the ionic entity and hydrogen-bond jump events accompanied by large angular jumps. The longest picosecond time scales represent structural relaxation, including large angular jumps and ion-pair dynamics. Also, ion cage lifetimes correlate with the slowest time scale of the associated dynamics of vibrational spectral diffusion despite the type of DFT functional. This study benchmarks DFT treatments of the exchange-correlation functional with and without the van der Waals (vdW) dispersion correction scheme. The inclusion of vdW interactions to the PBE functional represents a less structured state of the ionic entity and faster dynamics of the molecular motions relative to the one predicted by the PBE system. All the results illustrate the necessity of accurately describing the Coulomb interactions, vdW dispersive interactive forces, and localized hydrogen bonds required to sustain the energetic balance in this ionic salt.

Amphiphilicity of Intricate Layered Graphene/g-C3N4 Nanosheets

The hybrid heterostructure of the tri-s-triazine form of graphitic carbon nitride (g-C₃N₄), a stable two-dimensional material, results from intricate layer formation with graphene. In this material, g-C₃N₄, an amphiphilic material, stabilizes Pickering emulsions as an emulsifier and can effectively disperse graphene. Due to the various technological applications of the hybrid nanosheets in an aqueous environment, it is essential to study the interaction of water molecules with graphene and g-C₃N₄ (Gr/g-C₃N₄)-combined heterostructure. Although few studies have been performed signifying the water orientation in the interfacial layer, we find that there is a lack of detailed studies using various dynamical and structural properties of the interfacial water molecules. The interface of the Gr/g-C₃N₄ hybrid structure, one of the rarely found amphiphilic interfaces (on the g-C₃N₃ side), is appropriate for exploring the water affinity due to the availability of heterogeneous interfacial aqueous interactions. We adopted classical molecular dynamics simulations using two models for water molecules to study the structure and dynamics of an aqueous interface. We have correlated the structural properties to dynamics and spectral properties to understand the overall behavior of the amphiphilic interface. Our results branch into two significant hydrogen bond (HB) properties in HB count and HB strength among the water molecules in the different layers. The HB counts in the different layers of water are correlated using the average distance distribution (P_rO4), tetrahedral order parameters, HB donor/acceptor count, and total HBs per water molecule. A conspicuous difference is found in the HB count and related dynamics of the system. The HB lifetime and diffusion coefficient hint at the equivalent strength of HBs in the different layers. All the findings conclude that the amphiphilicity of the Gr/g-C₃N₄ interface can help in understanding various interfacial physical and chemical processes.

Activation-Induced Reorganization of Energy Transport Networks in the β2 Adrenergic Receptor

We compute energy exchange networks (EENs) through the β₂ adrenergic receptor (β₂AR), a G-protein coupled receptor (GPCR), in inactive and active states, based on the results of molecular dynamics simulations of this membrane bound protein. We introduce a new definition for the reorganization of EENs upon activation that depends on the relative change in rates of energy transfer across noncovalent contacts throughout the protein. On the basis of the reorganized network that we obtain for β₂AR upon activation, we identify a branched pathway between the agonist binding site and the cytoplasmic region, where a G-protein binds to the receptor when activated. The pathway includes all of the motifs containing molecular switches previously identified as contributing to the allosteric transition of β₂AR upon agonist binding. EENs and their reorganization upon activation are compared with structure-based contact networks computed for the inactive and active states of β₂AR.

Dynamics of Ionic Liquid through Intrinsic Vibrational Probes Using the Dispersion-Corrected DFT Functionals

First principles molecular dynamics simulations have been utilized to study the spectral properties of the protic ionic liquid, methylammonium formate (MAF). All simulations were performed using density functional theory (DFT) and various van der Waals-corrected exchange-correlation functionals. We calculated the vibrational stretch frequency distributions, determined the time-frequency correlations of the intrinsic vibrational probes, the N-H and C-O modes in MAF, and the frequency-structure correlations. We also estimated the average hydrogen-bond lifetimes and orientation dynamics to capture the ultrafast spectral response. The spectroscopic signature of the N-H stretching vibrations using the Becke-Lee-Yang-Parr (BLYP) and Perdew-Burke-Ernzerhof (PBE) functionals displays a spectral shift in the lower frequency side, suggesting stronger hydrogen-bonding interactions represented by the gradient approximation functionals than the van der Waals (vdW)-corrected simulations. The carboxylate frequency profiles with the dispersion-corrected representations are almost similar without a significant difference in the normalized distributions. Besides, the COO stretching frequencies at the peak maxima positions of the PBE functionals exhibit a lesser deviation from the experimental data. Spectral diffusion dynamics of the intrinsic vibrational probes on the cationic and anionic sites of the ionic liquid proceed through a short time relaxation of the intact hydrogen bonds followed by an intermediate time constant and a longer time decay indicating the switchover of hydrogen bonds. Dispersion-corrected atom-centered one-electron potential (DCACP) correction added to the BLYP system slows down the picosecond time scales of frequency correlation and the time constants of rotational motion, lengthening the overall system dynamics. The observed trends in the time-dependent decays of frequency fluctuations and the orientation autocorrelation functions correlate with the structural interactions in liquid MAF and hydrogen-bond dynamics. In this study, we examine the predictions made by different density functional treatments comparing the results of the uncorrected BLYP and PBE representations with the semiempirical vdW methods of Grimme and matching our calculated data with the experimental observations.

Ultrafast coherent motion and helix rearrangement of homodimeric hemoglobin visualized with femtosecond X-ray solution scattering

Ultrafast motion of molecules, particularly the coherent motion, has been intensively investigated as a key factor guiding the reaction pathways. Recently, X-ray free-electron lasers (XFELs) have been utilized to elucidate the ultrafast motion of molecules. However, the studies on proteins using XFELs have been typically limited to the crystalline phase, and proteins in solution have rarely been investigated. Here we applied femtosecond time-resolved X-ray solution scattering (fs-TRXSS) and a structure refinement method to visualize the ultrafast motion of a protein. We succeeded in revealing detailed ultrafast structural changes of homodimeric hemoglobin involving the coherent motion. In addition to the motion of the protein itself, the time-dependent change of electron density of the hydration shell was tracked. Besides, the analysis on the fs-TRXSS data of myoglobin allows for observing the effect of the oligomeric state on the ultrafast coherent motion.

Femtosecond time-resolved X-ray solution scattering (fs-TRXSS) measurements provide information on the structural dynamics of proteins in solution. Here, the authors present a structure refinement method for the analysis of fs-TRXSS data and use it to characterise the ultrafast structural changes of homodimeric haemoglobin.

Energy Transfer across Nonpolar and Polar Contacts in Proteins: Role of Contact Fluctuations

Molecular dynamics simulations of the villin headpiece subdomain HP36 have been carried out to examine relations between rates of vibrational energy transfer across non-covalently bonded contacts and equilibrium structural fluctuations, with focus on van der Waals contacts. Rates of energy transfer across van der Waals contacts vary inversely with the variance of the contact length, with the same constant of proportionality for all nonpolar contacts of HP36. A similar relation is observed for hydrogen bonds, but the proportionality depends on contact pairs, with hydrogen bonds stabilizing the α-helices all exhibiting the same constant of proportionality, one that is distinct from those computed for other polar contacts. Rates of energy transfer across van der Waals contacts are found to be up to 2 orders of magnitude smaller than rates of energy transfer across polar contacts.

Distinctive behavior and two-dimensional vibrational dynamics of water molecules inside glycine solvation shell

We present a first principles molecular dynamics study of a deuterated aqueous solution of a single glycine moiety to explore the structure, dynamics, and two-dimensional infrared spectra of water molecules found in the solvation shell of glycine.

Protein Structural Dynamics of Wild-Type and Mutant Homodimeric Hemoglobin Studied by Time-Resolved X-Ray Solution Scattering

The quaternary transition between the relaxed (R) and tense (T) states of heme-binding proteins is a textbook example for the allosteric structural transition. Homodimeric hemoglobin (HbI) from Scapharca inaequivalvis is a useful model system for investigating the allosteric behavior because of the relatively simple quaternary structure. To understand the cooperative transition of HbI, wild-type and mutants of HbI have been studied by using time-resolved X-ray solution scattering (TRXSS), which is sensitive to the conformational changes. Herein, we review the structural dynamics of HbI investigated by TRXSS and compare the results of TRXSS with those of other techniques.

Molecules and the Eigenstate Thermalization Hypothesis

We review a theory that predicts the onset of thermalization in a quantum mechanical coupled non-linear oscillator system, which models the vibrational degrees of freedom of a molecule. A system of N non-linear oscillators perturbed by cubic anharmonic interactions exhibits a many-body localization (MBL) transition in the vibrational state space (VSS) of the molecule. This transition can occur at rather high energy in a sizable molecule because the density of states coupled by cubic anharmonic terms scales as N³, in marked contrast to the total density of states, which scales as exp(aN), where a is a constant. The emergence of a MBL transition in the VSS is seen by analysis of a random matrix ensemble that captures the locality of coupling in the VSS, referred to as local random matrix theory (LRMT). Upon introducing higher order anharmonicity, the location of the MBL transition of even a sizable molecule, such as an organic molecule with tens of atoms, still lies at an energy that may exceed the energy to surmount a barrier to reaction, such as a barrier to conformational change. Illustrative calculations are provided, and some recent work on the influence of thermalization on thermal conduction in molecular junctions is also discussed.

Dynamical crossover of water confined within the amphiphilic nanocores of aggregated amyloid β peptides

It is believed that the self-assembly of amyloid beta (Aβ) peptides in the brain is the cause of Alzheimer's disease. Atomistic molecular dynamics simulations of aqueous solutions of Aβ protofilaments of different sizes at room temperature have been carried out to explore the dynamic properties of water confined within the core and at the exterior surface of the protofilaments. Attempts have been made to understand how the non-uniform distortion of the protofilaments associated with their structural crossover influences the diffusivity and the hydrogen bonding environment of the confined water molecules. In contrast to the homogeneous solvent dynamical environment at the exterior surface, the calculations revealed heterogeneously restricted motions of water confined within the distorted cores of the protofilaments. Importantly, it is demonstrated that the structural crossover of the aggregates observed for the decamer is associated with a dynamical transition of water confined within its core. A direct one-to-one correlation between the heterogeneously restricted core water motions and the kinetics of the breaking and formation of hydrogen bonds quantitatively demonstrated that a modified hydrogen bond arrangement within the cores of higher order Aβ protofilaments is the origin behind the crossover in core water mobility. A fraction of the water molecules forming short-lived water-water hydrogen bonds within the core of the crossover protofilament decamer are believed to diffuse away easily from the core and thus play a crucial role in further growth of the protofilament by facilitating the binding of new peptide monomers.

Computational approaches for deciphering the equilibrium and kinetic properties of iron transport proteins

With the advances in three-dimensional structure determination techniques, high quality structures of the iron transport proteins transferrin and the bacterial ferric binding protein (FbpA) have been deposited in the past decade. These are proteins of relatively large size, and developments in hardware and software have only recently made it possible to study their dynamics using standard computational resources. We review computational techniques towards understanding the equilibrium and kinetic properties of iron transport proteins under different environmental conditions. At the level of detail that requires quantum chemical treatments, the octahedral geometry around iron has been scrutinized and it has been established that the iron coordinating tyrosines are in an unusual deprotonated state. At the atomistic level, both the N-lobe and the full bilobal structure of transferrin have been studied under varying conditions of pH, ionic strength and binding of other metal ions by molecular dynamics (MD) simulations. These studies have allowed questions to be answered, among others, on the function of second shell residues in iron release, the role of synergistic anions in preparing the active site for iron binding, and the differences between the kinetics of the N- and the C-lobe. MD simulations on FbpA have led to the detailed observation of the binding kinetics of phosphate to the apo form, and to the conformational preferences of the holo form under conditions mimicking the environmental niches provided by the periplasmic space. To study the dynamics of these proteins with their receptors, one must resort to coarse-grained methodologies, since these systems are prohibitively large for atomistic simulations. A study of the complex of human transferrin (hTf) with its pathogenic receptor by such methods has revealed a potential mechanistic explanation for the defense mechanism that arises in evolutionary warfare. Meanwhile, the motions in the transferrin receptor bound hTf have been shown to disfavor apo hTf dissociation, explaining why the two proteins remain in complex during the recycling process from the endosome to the cell surface. Open problems and possible technological applications related to metal ion binding-release in iron transport proteins that may be handled by hybrid use of quantum mechanical, MD and coarse-grained approaches are discussed.

Influence of Cholesterol on the Dynamics of Hydration in Phospholipid Bilayers

We investigate the dynamics of interfacial waters in dipalmitoylphosphatidylcholine (DPPC) bilayers upon the addition of cholesterol, by molecular dynamics simulations. Our data reveal that the inclusion of cholesterol modifies the membrane aqueous interfacial dynamics: waters diffuse faster, their rotational decay time is shorter, and the DPPC/water hydrogen bond dynamics relaxes faster than in the pure DPPC membrane. The observed acceleration of the translational water dynamics agrees with recent experimental results, in which, by means of NMR techniques, an increment of the surface water diffusivity is measured upon the addition of cholesterol. A microscopic analysis of the lipid/water hydrogen bond network at the interfacial region suggests that the mechanism underlying the observed water mobility enhancement is given by the rupture of a fraction of interlipid water bridge hydrogen bonds connecting two different DPPC molecules, concomitant to the formation of new lipid/solvent bonds, whose dynamics is faster than that of the former. The consideration of a simple two-state model for the decay of the hydrogen bond correlation function yielded excellent results, obtaining two well-separated characteristic time scales: a slow one (∼250 ps) associated with bonds linking two DPPC molecules, and a fast one (∼15 ps), related to DPPC/solvent bonds.

Computational study to understand the energy transfer pathways within amicyanin

Vibrational energy diffusivities between the residues present in Amicyanin copper protein are calculated and presented in form of communication map. From those results energy flow pathways from the copper metal ion to the inter protein residue Glu31 are identified. Our finding suggests many different pathways are possible and copper metal ion in oxidized and reduced state switches the pathways. Our finding also suggests the cooperative nature of surrounding residues and water molecules towards selecting the pathways. The major transport channels in the oxidised state are, Cu²⁺---> MET28---> LYS29---> TYR30---> GLU31 and Cu²⁺---> MET98---> TYR30--- > GLU31. And in the reduced state Cu⁺---> CYS9---> TYR30---> GLU31 and Cu⁺---> MET28---> LYS2---> TYR30---> GLU31. We studied further the interaction energies between the copper ion and neighbouring residues using B3LYP/QZVP method. Both the methods complement each other in predicting the energy flow pathways and the cooperative nature of residues.

A review of the structure and dynamics of nanoconfined water and ionic liquids via molecular dynamics simulation

During the past decade, the research on fluids in nanoconfined geometries has received considerable attention as a consequence of their wide applications in different fields. Several nanoconfined systems such as water and ionic liquids, together with an equally impressive array of nanoconfining media such as carbon nanotube, graphene and graphene oxide have received increasingly growing interest in the past years. Water is the first system that has been reviewed in this article, due to its important role in transport phenomena in environmental sciences. Water is often considered as a highly nanoconfined system, due to its reduction to a few layers of water molecules between the extended surface of large macromolecules. The second system discussed here is ionic liquids, which have been widely studied in the modern green chemistry movement. Considering the great importance of ionic liquids in industry, and also their oil/water counterpart, nanoconfined ionic liquid system has become an important area of research with many fascinating applications. Furthermore, the method of molecular dynamics simulation is one of the major tools in the theoretical study of water and ionic liquids in nanoconfinement, which increasingly has been joined with experimental procedures. In this way, the choice of water and ionic liquids in nanoconfinement is justified by applying molecular dynamics simulation approaches in this review article.

Interfacial water effect on cooperativity and signal communication in Scapharca dimeric hemoglobin

Cooperativity is important in controlling the biological functions of allosteric proteins.

Probing the communication of deoxythymidine triphosphate in HIV-1 reverse transcriptase by communication maps and interaction energy studies

We calculate communication maps for HIV-1 Reverse Transcriptase (RT) to elucidate energy transfer pathways between deoxythymidine triphosphate (dTTP) and other parts of the protein.

Water-Mediated Energy Dynamics in a Homodimeric Hemoglobin

We examine energy dynamics in the unliganded and liganded states of the homodimeric hemoglobin from Scapharca inaequivalvis (HbI), which exhibits cooperativity mediated by the cluster of water molecules at the interface upon ligand binding and dissociation. We construct and analyze a dynamic network in which nodes representing the residues, hemes, and water cluster are connected by edges that represent energy transport times, as well as a nonbonded network (NBN) indicating regions that respond rapidly to local strain within the protein via nonbonded interactions. One of the two largest NBNs includes the Lys30-Asp89 salt bridge critical for stabilizing the dimer. The other includes the hemes and surrounding residues, as well as, in the unliganded state, the cluster of water molecules between the globules. Energy transport in the protein appears to be controlled by the Lys30-Asp89 salt bridge critical for stabilizing the dimer, as well as the interface water cluster in the unliganded state. Possible connections between energy transport dynamics in response to local strain identified here and allosteric transitions in HbI are discussed.

Effects of protein-DNA complex formation on the intermolecular vibrational density of states of interfacial water

Single-stranded DNAs (ss-DNAs) are formed as intermediates during DNA metabolic processes. ss-DNA binding (SSB) proteins specifically bind to the single-stranded segments of the DNA and protect it from being degraded. We have performed room temperature molecular dynamics simulations of the aqueous solution of two DNA binding K homology (KH) domains (KH3 and KH4) of the far upstream element (FUSE) binding protein (FBP) complexed with two ss-DNA oligomers. Efforts have been made to explore the influence of complex formation on low-frequency vibrational density of states of the surface water molecules. It is revealed that increased back scattering of water confined around the complexed structures leads to significant blue shifts of the band corresponding to the O···O···O bending or restricted transverse motions of water, the effect being more for the bridged water molecules. Importantly, it is demonstrated that the formation of such complexed structures of a similar type may often influence the transverse and longitudinal degrees of freedom of the surrounding water molecules in a nonuniform manner.

Cooperative protein structural dynamics of homodimeric hemoglobin linked to water cluster at subunit interface revealed by time-resolved X-ray solution scattering

Homodimeric hemoglobin (HbI) consisting of two subunits is a good model system for investigating the allosteric structural transition as it exhibits cooperativity in ligand binding. In this work, as an effort to extend our previous study on wild-type and F97Y mutant HbI, we investigate structural dynamics of a mutant HbI in solution to examine the role of well-organized interfacial water cluster, which has been known to mediate intersubunit communication in HbI. In the T72V mutant of HbI, the interfacial water cluster in the T state is perturbed due to the lack of Thr72, resulting in two less interfacial water molecules than in wild-type HbI. By performing picosecond time-resolved X-ray solution scattering experiment and kinetic analysis on the T72V mutant, we identify three structurally distinct intermediates (I1, I2, and I3) and show that the kinetics of the T72V mutant are well described by the same kinetic model used for wild-type and F97Y HbI, which involves biphasic kinetics, geminate recombination, and bimolecular CO recombination. The optimized kinetic model shows that the R-T transition and bimolecular CO recombination are faster in the T72V mutant than in the wild type. From structural analysis using species-associated difference scattering curves for the intermediates, we find that the T-like deoxy I3 intermediate in solution has a different structure from deoxy HbI in crystal. In addition, we extract detailed structural parameters of the intermediates such as E-F distance, intersubunit rotation angle, and heme-heme distance. By comparing the structures of protein intermediates in wild-type HbI and the T72V mutant, we reveal how the perturbation in the interfacial water cluster affects the kinetics and structures of reaction intermediates of HbI.

Macro and nano scale modelling of water–water interactions at ambient and low temperature: relaxation and residence times

A method to predict magnitudes in quantitative agreement with experimental data has been devised and applied to model water–water interactions at low temperature.

Features of S-nitrosylation based on statistical analysis and molecular dynamics simulation: cysteine acidity, surrounding basicity, steric hindrance and local flexibility

S-Nitrosylation is involved in protein functional regulation and cellular signal transduction. Although intensive efforts have been made, the molecular mechanisms of S-nitrosylation have not yet been fully understood. In this work, we carried out a survey on 213 protein structures with S-nitrosylated cysteine sites and molecular dynamic simulations of hemoglobin as a case study. It was observed that the S-nitrosylated cysteines showed a lower pKa, a higher population of basic residues, a lower population of big-volume residues in the neighborhood, and relatively higher flexibility. The case study of hemoglobin showed that, compared to that in the T-state, Cysβ93 in the R-state hemoglobin possessed the above structural features, in agreement with the previous report that the R-state was more reactive in S-nitrosylation. Moreover, basic residues moved closer to the Cysβ93 in the dep-R-state hemoglobin, while big-volume residues approached the Cysβ93 in the dep-T-state. Using the four characteristics, i.e. cysteine acidity, surrounding basicity, steric hindrance, and local flexibility, a 3-dimensional model of S-nitrosylation was constructed to explain 61.9% of the S-nitrosylated and 58.1% of the non-S-nitrosylated cysteines. Our study suggests that cysteine deprotonation is a prerequisite for protein S-nitrosylation, and these characteristics might be useful in identifying specificity of protein S-nitrosylation.

Infrared multiple photon dissociation spectroscopy and density functional theory (DFT) studies of protonated permethylated β-cyclodextrin-water non-covalent complexes

We present infrared multiple photon dissociation (IRMPD) spectroscopy and quantum chemical calculation results for the protonated permethylated β-cyclodextrin (CD)-water non-covalent complex, the simplest β-CD non-covalent complex, in the gas-phase. The IRMPD spectrum in the region 2700-3750 cm(-1) consisted of three strong peaks at 3096, 3315, and 3490 cm(-1). These spectral features in the experimental IRMPD spectrum were compared with a large set of infrared absorption spectra predicted using density functional theory (DFT) calculations for the protonated β-CD-water complex. Complex III (see ), in which the water molecule (at the primary rim) and the proton (at the secondary rim) were separated, was found to suitably reflect the main spectral characteristics found in the experimental IRMPD spectrum. The absence of the homodromic hydrogen bond ring, due to replacement of hydroxyl groups with methoxy groups in permethylated β-CD, rendered the primary rim open compared with the unmodified β-CD 'one-gate-closed' lowest energy conformer. This study demonstrates that IRMPD studies combined with DFT theoretical calculations can be a good method for studying molecular interactions of large host-guest pairs.

Non-site-specific allosteric effect of oxygen on human hemoglobin under high oxygen partial pressure

Protein allostery is essential for vital activities. Allosteric regulation of human hemoglobin (HbA) with two quaternary states T and R has been a paradigm of allosteric structural regulation of proteins. It is widely accepted that oxygen molecules (O₂) act as a “site-specific” homotropic effector, or the successive O₂ binding to the heme brings about the quaternary regulation. However, here we show that the site-specific allosteric effect is not necessarily only a unique mechanism of O₂ allostery. Our simulation results revealed that the solution environment of high O₂ partial pressure enhances the quaternary change from T to R without binding to the heme, suggesting an additional “non-site-specific” allosteric effect of O₂. The latter effect should play a complementary role in the quaternary change by affecting the intersubunit contacts. This analysis must become a milestone in comprehensive understanding of the allosteric regulation of HbA from the molecular point of view.

Communication maps of vibrational energy transport through Photoactive Yellow Protein

We calculate communication maps for Photoactive Yellow Protein (PYP) from the purple phototropic eubacterium Halorhodospira halophile and use them to elucidate energy transfer pathways from the chromophore through the rest of the protein in the ground and excited state. The calculations reveal that in PYP excess energy from the chromophore flows mainly to regions of the surrounding residues that hydrogen bond to the chromophore. In addition, quantum mechanics/molecular mechanics and molecular dynamics (MD) simulations of the dielectric response of the protein and solvent environment due to charge rearrangement on the chromophore following photoexcitation are also presented, with both approaches yielding similar time constants for the response. Results of MD simulations indicate that the residues hydrogen bonding to the chromophore make the largest contribution to the response.

On the binding affinity of macromolecular interactions: daring to ask why proteins interact

Interactions between proteins are orchestrated in a precise and time-dependent manner, underlying cellular function. The binding affinity, defined as the strength of these interactions, is translated into physico-chemical terms in the dissociation constant (K_d), the latter being an experimental measure that determines whether an interaction will be formed in solution or not. Predicting binding affinity from structural models has been a matter of active research for more than 40 years because of its fundamental role in drug development. However, all available approaches are incapable of predicting the binding affinity of protein–protein complexes from coordinates alone. Here, we examine both theoretical and experimental limitations that complicate the derivation of structure–affinity relationships. Most work so far has concentrated on binary interactions. Systems of increased complexity are far from being understood. The main physico-chemical measure that relates to binding affinity is the buried surface area, but it does not hold for flexible complexes. For the latter, there must be a significant entropic contribution that will have to be approximated in the future. We foresee that any theoretical modelling of these interactions will have to follow an integrative approach considering the biology, chemistry and physics that underlie protein–protein recognition.