Anomalous Dynamics of Water Confined in Protein-Protein and Protein-DNA Interfaces

Dynamics of Protein-RNA Interfaces Using All-Atom Molecular Dynamics Simulations

Facing the current challenges posed by human health diseases requires the understanding of cell machinery at a molecular level. The interplay between proteins and RNA is key for any physiological phenomenon, as well protein-RNA interactions. To understand these interactions, many experimental techniques have been developed, spanning a very wide range of spatial and temporal resolutions. In particular, the knowledge of tridimensional structures of protein-RNA complexes provides structural, mechanical, and dynamical pieces of information essential to understand their functions. To get insights into the dynamics of protein-RNA complexes, we carried out all-atom molecular dynamics simulations in explicit solvent on nine different protein-RNA complexes with different functions and interface size by taking into account the bound and unbound forms. First, we characterized structural changes upon binding and, for the RNA part, the change in the puckering. Second, we extensively analyzed the interfaces, their dynamics and structural properties, and the structural waters involved in the binding, as well as the contacts mediated by them. Based on our analysis, the interfaces rearranged during the simulation time showing alternative and stable residue-residue contacts with respect to the experimental structure.

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.

Revealing the hidden dynamics of confined water in acrylate polymers: Insights from hydrogen-bond lifetime analysis

Polymers contain functional groups that participate in hydrogen bond (H-bond) with water molecules, establishing a robust H-bond network that influences bulk properties. This study utilized molecular dynamics (MD) simulations to examine the H-bonding dynamics of water molecules confined within three poly(meth)acrylates: poly(2-methoxyethyl acrylate) (PMEA), poly(2-hydroxyethyl methacrylate) (PHEMA), and poly(1-methoxymethyl acrylate) (PMC1A). Results showed that H-bonding dynamics significantly slowed as the water content decreased. Additionally, the diffusion of water molecules and its correlation with H-bond breakage were analyzed. Our findings suggest that when the H-bonds between water molecules and the methoxy oxygen of PMEA are disrupted, those water molecules persist in close proximity and do not diffuse on a picosecond time scale. In contrast, the water molecules H-bonded with the hydroxy oxygen of PHEMA and the methoxy oxygen of PMC1A diffuse concomitantly with the breakage of H-bonds. These results provide an in-depth understanding of the impact of polymer functional groups on H-bonding dynamics.

Flexibility of flanking DNA is a key determinant of transcription factor affinity for the core motif

Selective gene regulation is mediated by recognition of specific DNA sequences by transcription factors (TFs). The extremely challenging task of searching out specific cognate DNA binding sites among several million putative sites within the eukaryotic genome is achieved by complex molecular recognition mechanisms. Elements of this recognition code include the core binding sequence, the flanking sequence context, and the shape and conformational flexibility of the composite binding site. To unravel the extent to which DNA flexibility modulates TF binding, in this study, we employed experimentally guided molecular dynamics simulations of ternary complex of closely related Hox heterodimers Exd-Ubx and Exd-Scr with DNA. Results demonstrate that flexibility signatures embedded in the flanking sequences impact TF binding at the cognate binding site. A DNA sequence has intrinsic shape and flexibility features. While shape features are localized, our analyses reveal that flexibility features of the flanking sequences percolate several basepairs and allosterically modulate TF binding at the core. We also show that lack of flexibility in the motif context can render the cognate site resistant to protein-induced shape changes and subsequently lower TF binding affinity. Overall, this study suggests that flexibility-guided DNA shape, and not merely the static shape, is a key unexplored component of the complex DNA-TF recognition code.

A dynamical view of protein-protein complexes: Studies by molecular dynamics simulations

Protein-protein interactions are at the basis of many protein functions, and the knowledge of 3D structures of protein-protein complexes provides structural, mechanical and dynamical pieces of information essential to understand these functions. Protein-protein interfaces can be seen as stable, organized regions where residues from different partners form non-covalent interactions that are responsible for interaction specificity and strength. They are commonly described as a peripheral region, whose role is to protect the core region that concentrates the most contributing interactions, from the solvent. To get insights into the dynamics of protein-protein complexes, we carried out all-atom molecular dynamics simulations in explicit solvent on eight different protein-protein complexes of different functional class and interface size by taking into account the bound and unbound forms. On the one hand, we characterized structural changes upon binding of the proteins, and on the other hand we extensively analyzed the interfaces and the structural waters involved in the binding. Based on our analysis, in 6 cases out of 8, the interfaces rearranged during the simulation time, in stable and long-lived substates with alternative residue-residue contacts. These rearrangements are not restricted to side-chain fluctuations in the periphery but also affect the core interface. Finally, the analysis of the waters at the interface and involved in the binding pointed out the importance to take into account their role in the estimation of the interaction strength.

Temperature-Dependent Dynamics at Protein-Solvent Interfaces

We perform differential scanning calorimetry, broadband dielectric spectroscopy (BDS), and nuclear magnetic resonance (NMR) studies to ascertain the molecular dynamics in mixtures of ethylene glycol with elastin or lysozyme over broad temperature ranges. To focus on the protein-solvent interface, we use mixtures with about equal numbers of amino acids and solvent molecules. The elastin and lysozyme mixtures show similar glass transition steps, which extend over a broad temperature range of 157-185K. The BDS and NMR studies yield fully consistent results for the fastest process P1, which is caused by the structural relaxation of ethylene glycol between the protein molecules and follows an Arrhenius law with an activation energy of Ea=0.63eV. It involves quasi-isotropic reorientation and is very similar in the elastin and lysozyme matrices but different from the alpha and beta relaxations of bulk ethylene glycol. Two slower BDS processes P2 and P3 have protein-dependent time scales, but exhibit a similar Arrhenius-like temperature dependence with an activation energy of Ea~0.81eV. However, P2 and P3 do not have a clear NMR signature. In particular, the NMR results for the lysozyme mixture reveal that the protein backbone does not show isotropic alpha-like motion on the P2 and P3 time scales but only restricted beta-like reorientation. The different activation energies of the P1 and P2/P3 processes do not support an intimate coupling of protein and ethylene glycol dynamics. The present results are compared with previous findings for mixtures of proteins with water or glycerol, implying qualitatively different dynamical couplings at various protein-solvent interfaces.

Interfacial Water in the SARS Spike Protein: Investigating the Interaction with Human ACE2 Receptor and In Vitro Uptake in A549 Cells

The severity of global pandemic due to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has engaged the researchers and clinicians to find the key features triggering the viral infection to lung cells. By utilizing such crucial information, researchers and scientists try to combat the spread of the virus. Here, in this work, we performed in silico analysis of the protein-protein interactions between the receptor-binding domain (RBD) of the viral spike protein and the human angiotensin-converting enzyme 2 (hACE2) receptor to highlight the key alteration that happened from SARS-CoV to SARS-CoV-2. We analyzed and compared the molecular differences between spike proteins of the two viruses using various computational approaches such as binding affinity calculations, computational alanine, and molecular dynamics simulations. The binding affinity calculations showed that SARS-CoV-2 binds a little more firmly to the hACE2 receptor than SARS-CoV. The major finding obtained from molecular dynamics simulations was that the RBD-ACE2 interface is populated with water molecules and interacts strongly with both RBD and ACE2 interfacial residues during the simulation periods. The water-mediated hydrogen bond by the bridge water molecules is crucial for stabilizing the RBD and ACE2 domains. Near-ambient pressure X-ray photoelectron spectroscopy (NAP-XPS) confirmed the presence of vapor and molecular water phases in the protein-protein interfacial domain, further validating the computationally predicted interfacial water molecules. In addition, we examined the role of interfacial water molecules in virus uptake by lung cell A549 by binding and maintaining the RBD/hACE2 complex at varying temperatures using nanourchins coated with spike proteins as pseudoviruses and fluorescence-activated cell sorting (FACS) as a quantitative approach. The structural and dynamical features presented here may serve as a guide for developing new drug molecules, vaccines, or antibodies to combat the COVID-19 pandemic.

Surprising Rigidity of Functionally Important Water Molecules Buried in the Lipid Headgroup Region

Understanding water dynamics and structure is an important topic in biological systems. It is generally held in the literature that the interfacial water of hydrated phospholipids is highly mobile, in fast exchange with the bulk water ranging from the nano- to femtosecond timescale. Although nuclear magnetic resonance (NMR) is a powerful tool for structural and dynamic studies, direct probing of interfacial water in hydrated phospholipids is formidably challenging due to the extreme population difference between bulk and interfacial water. We developed a novel ¹⁷O solid-state NMR technique in combination with an ultra-high-field magnet (35.2 T) to directly probe the functionally important interfacial water. By selectively suppressing the dominant bulk water signal, we observed two distinct water species in the headgroup region of hydrated dimyristoylphosphatidylcholine (DMPC) lipid bilayers for the first time. One water species denoted as "confined water" is chemically and dynamically different from the bulk water (∼0.17 ppm downfield and a slightly shorter spin-lattice relaxation time). Another water species denoted as "bound water" has severely restricted motion and a distinct chemical shift (∼12 ppm upfield). Additionally, the bulk water is not as "free" as pure water, resulting from the fast exchange with the water molecules that weakly and transiently interact with the lipid choline groups. These new discoveries clearly indicate the existence of the interfacial water molecules that are relatively stable over the NMR timescale (on the order of milliseconds), providing an opportunity to characterize water dynamics on the millisecond or slower timescale in biomacromolecules.

Critical In-Plane Density of Polyelectrolyte Brush for the Ordered Hydrogen-Bonded Structure of Incorporated Water

A polymer electrolyte brush is a reasonable platform to confine water molecules within a nanoscopic area to study their role in the function of interacting media because of their adjustable nanospace and charge by changing the in-plane density and side chains of the brush. Here, we demonstrate how the in-plane spacing of cationic polymer brush chains, poly[2-(methacryloyloxy)ethyltrimethylammonium chloride] (PMTAC), affects the hydrogen bond configuration of incorporated water using soft X-ray emission spectroscopy. At the critical in-plane density σ = 0.30 chains/nm² of PMTAC, tetrahedrally coordinated water molecules started to melt into distorted or broken hydrogen-bonded configurations. Considering the charge on the quaternary ammonium cations, the electric field required to form a tetrahedrally coordinated hydrogen-bonded configuration was estimated as ∼500 kV cm^-1 and is effective up to ∼1 nm from the surface of the polymer chain. These findings are useful for designing specific interface properties and the resultant surface function of polyelectrolyte-based materials.

Influence of molecular rebinding on the reaction rate of complex formation

We simulated Brownian diffusion and reaction-diffusion processes to study the influence of molecular rebinding on the reaction rates of bimolecular reactions. We found that the number of rebinding events, N_reb, is proportional to the target's size and inversely proportional to the diffusion coefficient D and simulation time-step Δt. We found the proportionality constant close to π^-1/2. We confirmed that N_reb is defined as a ratio of the activation-limited rate constant k_a to the diffusion-limited rate constant, k_D. We provide the formula describing the reactivity coefficient κ, modelling the transient-native complex transition for the activation-controlled reaction rates. We show that κ is proportional to (D/Δt)^1/2. Finally, we apply our rebinding-including reaction rate model to the real reactions of photoacid dissociation and protein association. Based on literature data for both types of reactions, we found the Δt time-scale. We show that for the photodissociation of a proton, the Δt is equal to 171 ± 18 fs and the average number of rebinding events is approximately equal to 40. For proteins, Δt is of the order of 100 ps with around 20 rebinding events. In both cases the timescale is similar to the timescale of fluctuation of the solvent molecules surrounding the reactants; vibrations and bending in the case of photoacid dissociation and diffusional motion for proteins.

Effect of trehalose on protein cryoprotection: Insights into the mechanism of slowing down of hydration water

We study, with molecular dynamics simulations, a lysozyme protein immersed in a water-trehalose solution upon cooling. The aim is to understand the cryoprotectant role played by this disaccharide through the modifications that it induces on the slow dynamics of protein hydration water with its presence. The α-relaxation shows a fragile to strong crossover about 20° higher than that in the bulk water phase and 15° higher than that in lysozyme hydration water without trehalose. The protein hydration water without trehalose was found to show a second slower relaxation exhibiting a strong to strong crossover coupled with the protein dynamical transition. This slower relaxation time importantly appears enormously slowed down in our cryoprotectant solution. On the other hand, this long-relaxation in the presence of trehalose is also connected with a stronger damping of the protein structural fluctuations than that found when the protein is in contact with the pure hydration water. Therefore, this appears to be the mechanism through which trehalose manifests its cryoprotecting function.

Competition between chiral solvents and chiral monomers in the helical bias of supramolecular polymers

Solute-solvent interactions are key for the assembly and proper functioning of biomacromolecules and play important roles in many fields of organic and polymer chemistry. Despite numerous reports describing the effects of (chiral) solvents on helical conformations of (supramolecular) polymers, the combination of chiral solvents and chiral monomers is unexplored. Here we report diastereomeric differences in the supramolecular polymerization of enantiomers of chiral triphenylene-2,6,10-tricarboxamides in chiral chlorinated solvents. Competition between the preferences induced by the stereocentres of the assembled monomers and those present in the solvent molecules results in unforeseen temperature-dependent solvation effects. By combining experiments and mathematical modelling, we show that the observed differences between enantiomers originate from the combined additive entropic effects of stereocentres present in the monomer and in the solvent. Remarkably, copolymerizations show that the chiral solvent can bias the copolymer helicity and thereby overrule the helical preference of the monomers. Our results highlight the importance of cumulative solvation effects in supramolecular polymerizations.

Revisiting the Rate-Limiting Step of the ANS-Protein Binding at the Protein Surface and Inside the Hydrophobic Cavity

8-Anilino-1-naphthalenesulfonic acid (ANS) is used as a hydrophobic fluorescence probe due to its high intensity in hydrophobic environments, and also as a microenvironment probe because of its unique ability to exhibit peak shift and intensity change depending on the surrounding solvent environment. The difference in fluorescence can not only be caused by the microenvironment but can also be affected by the binding affinity, which is represented by the binding constant (K). However, the overall binding process considering the binding constant is not fully understood, which requires the ANS fluorescence binding mechanism to be examined. In this study, to reveal the rate-limiting step of the ANS-protein binding process, protein concentration-dependent measurements of the ANS fluorescence of lysozyme and bovine serum albumin were performed, and the binding constants were analyzed. The results suggest that the main factor of the binding process is the microenvironment at the binding site, which restricts the attached ANS molecule, rather than the attractive diffusion-limited association. The molecular mechanism of ANS-protein binding will help us to interpret the molecular motions of ANS molecules at the binding site in detail, especially with respect to an equilibrium perspective.

Multiresolution continuous wavelet transform for studying coupled solute-solvent vibrations via ab initio molecular dynamics

Vibrational analysis in solution and the theoretical determination of infrared and Raman spectra are of key importance in many fields of chemical interest. Vibrational band dynamics of molecules and their sensitivity to the environment can also be captured by these spectroscopies in their time dependent version. However, it is often difficult to provide an interpretation of the experimental data at the molecular scale, such as molecular mechanisms or the processes hidden behind them. In this work, we present a theoretical-computational protocol based on ab initio molecular dynamics simulations and a combination of normal-like (generalized) mode analysis of solute-solvent clusters with a wavelet transform, for the first time. The case study is the vibrational dynamics of N-methyl-acetamide (NMA) in water solution, a well-known model of hydration of peptides and proteins. Amide modes are typical bands of peptide and protein backbone, and their couplings with the environment are very challenging in terms of the accurate prediction of solvent induced intensity and frequency shifts. The contribution of water molecules surrounding NMA to the composition of generalized and time resolved modes is introduced in our vibrational analysis, showing unequivocally its influence on the amide mode spectra. It is also shown that such mode compositions need the inclusion of the first shell solvent molecules to be accurately described. The wavelet analysis is proven to be strongly recommended to follow the time evolution of the spectra, and to capture vibrational band couplings and frequency shifts over time, preserving at the same time a well-balanced time-frequency resolution. This peculiar feature also allows one to perform a combined structural-vibrational analysis, where the different strengths of hydrogen bond interactions can quantitatively affect the amide bands over time at finite temperature. The proposed method allows for the direct connection between vibrational modes and local structural changes, providing a link from the spectroscopic observable to the structure, in this case the peptide backbone, and its hydration layouts.

Heterogeneous Dynamical Environment at the Interface of a Protein-DNA Complex

Binding between protein and DNA is an essential process to regulate different biological activities. Two puzzling questions in protein-DNA recognition are (i) how the protein's binding domain identifies the DNA sequence in an aqueous solution and (ii) how the formation of the complex alters the dynamical environment around it. In this work, we present results obtained from molecular dynamics simulations of the N-terminal α-helical domain of the λ-repressor protein (in dimeric form) bound to the corresponding operator DNA. Effects of formation of the complex in modifying the microscopic dynamics of water as well as the kinetics of hydrogen bonds at the interface have been explored. Locally heterogeneous restricted water motions at the complex interface have been observed, the extent of restriction being more significant around the directly bound residues of the protein and the DNA. In particular, the calculation revealed the existence of significantly constrained motionally restricted water layer that can form either bridges around the directly bound residues of the protein and DNA or are engaged in forming water-mediated contacts between a fraction of the unbound residues. More importantly, it is observed that the restricted water motion around the complex is correlated with the hydrogen bond relaxation time scale at the interface. It is further demonstrated that the kinetics of water-water hydrogen bonds involving the bridged water are influenced more due to complex formation.

Dynamic signature of ligand binding over a protein surface

We study the motion of Zn^{2+} in the presence of ubiquitin by all-atom molecular-dynamics simulations. We observe that unlike normal diffusive liquid, metal ions show an exponential tail in the self-van Hove function (self-vHf). Moreover, the metal ions are trapped strongly by acidic residues which form a binding pocket over the protein surface. The exponential tail disappears by mutation of trapping residues, suggesting that the tail appears due to trapped motion of the ions. The mean-squared displacements, however, in all the cases show linear dependence on time. Our model establishes that ligand binding generically results in an exponential tail of self-vHf. The self-vHf may give an approach to find binding pockets over a protein surface.

Spectroscopic Analysis of Protein-Crowded Environments Using the Charge-Transfer Fluorescence Probe 8-Anilino-1-Naphthalenesulfonic Acid

The molecular behaviors of proteins under crowding conditions are crucial for understanding the protein actions in intracellular environments. Under a crowded environment, the distance between protein molecules is almost the same size as the molecular level, thus, both the excluded volume effect and short ranged soft chemical interaction on protein surface could induce the complicated influence on the protein behavior cooperatively. Recently, various kinds of analytical approaches from macroscopic to microscopic aspects have been made to evaluate the crowding effect. The method, however, has not been established to evaluate the surface specific interactions on protein surface. In this study, the analytical method to evaluate the crowding effect has been suggested by using a charge-transfer fluorescence probe, ANS. By employing the unique property of ANS attaching to charged residues on the surface of lysozyme, the crowding effect was focused, while the case was compared as a reference, in which ANS is confined in hydrophobic pockets of BSA. Consequently, the surface specific changes of fluorescence spectra were readily observed under the crowded environment, whereas the fluorescence spectra of ANS in protein inside did not change. This result suggests the fluorescence spectra of ANS binding to protein surface have the capability to estimate the crowding effect of proteins.

Nanoscopic insights of saline water in carbon nanotube appended filters using molecular dynamics simulations

Nanotube appended membranes are shown to be very promising due to their ultrafast water transport and very high salt rejection ability. Using classical molecular dynamics, the present study reports the nanoscopic assessment of various molecular events for nanotube-based desalination, which might be useful for nanoscale devices during process operation at the macroscopic scale. The characteristics of water and ion flow are discussed with varied strength of pressure gradient and salt concentration for different scales of confinement. The results revealed that the membranes comprising nanotubes of 1.0-1.1 nm diameter can be optimized for efficient water desalination with more than >95% salt rejection. Furthermore, the anomalies in water flux through nanotubes are linked with the hydration characteristics of ions inside CNTs. The results show the maximum hydration of confined ions inside the nanotubes, which indicated the minimum permeability of water due to freezing effects. Furthermore, the MD results revealed that akin to bulk phases, the mass transport through nanotubes can be linked with the component diffusivity in the medium. It has been demonstrated that not only the diffusivities of water and ions, but even the gradient of water to ion diffusivity might be utilized to predict and explore the experimental observations, which might be helpful in optimizing the operational regime in nanotube-based filtrations. Moreover, the thermodynamic characteristics of the flow are discussed in terms of the entropy of water and ions using the robust two-phase thermodynamic (2PT) method. The results reflect that the entropy of water is linked to the distortion of the hydrogen bond network inside the nanotube confinement, at the nanotube-water interface and at the bulk solution, whereas the entropy of ions seems to be majorly dominated by their oscillation. Also, the interconnection of hydration structure, mass flux and the diffusivity of water and ions along with their thermodynamic origin are discussed.

Ultrafast energy relaxation dynamics of amide I vibrations coupled with protein-bound water molecules

The influence of hydration water on the vibrational energy relaxation in a protein holds the key to understand ultrafast protein dynamics, but its detection is a major challenge. Here, we report measurements on the ultrafast vibrational dynamics of amide I vibrations of proteins at the lipid membrane/H₂O interface using femtosecond time-resolved sum frequency generation vibrational spectroscopy. We find that the relaxation time of the amide I mode shows a very strong dependence on the H₂O exposure, but not on the D₂O exposure. This observation indicates that the exposure of amide I bond to H₂O opens up a resonant relaxation channel and facilitates direct resonant vibrational energy transfer from the amide I mode to the H₂O bending mode. The protein backbone motions can thus be energetically coupled with protein-bound water molecules. Our findings highlight the influence of H₂O on the ultrafast structure dynamics of proteins.

Overhauser Dynamic Nuclear Polarization for the Study of Hydration Dynamics, Explained

We outline the physical properties of hydration water that are captured by Overhauser Dynamic Nuclear Polarization (ODNP) relaxometry and explore the insights that ODNP yields about the water and the surface that this water is coupled to. As ODNP relies on the pairwise cross-relaxation between the electron spin of a spin probe and a proton nuclear spin of water, it captures the dynamics of single-particle diffusion of an ensemble of water molecules moving near the spin probe. ODNP principally utilizes the same physics as other nuclear magnetic resonance (NMR) relaxometry (i.e., relaxation measurement) techniques. However, in ODNP, electron paramagnetic resonance (EPR) excites the electron spins probes and their high net polarization acts as a signal amplifier. Furthermore, it renders ODNP parameters highly sensitive to water moving at rates commensurate with the EPR frequency of the spin probe (typically 10GHz). Also, ODNP selectively enhances the NMR signal contributions of water moving within close proximity to the spin label. As a result, ODNP can capture ps-ns movements of hydration waters with high sensitivity and locality, even in samples with protein concentrations as dilute as 10 µM. To date, the utility of the ODNP technique has been demonstrated for two major applications: the characterization of the spatial variation in the properties of the hydration layer of proteins or other surfaces displaying topological diversity, and the identification of structural properties emerging from highly disordered proteins and protein domains. The former has been shown to correlate well with the properties of hydration water predicted by MD simulations and has been shown capable of evaluating the hydrophilicity or hydrophobicity of a surface. The latter has been demonstrated for studies of an interhelical loop of proteorhodopsin, the partial structure of α-synuclein embedded at the lipid membrane surface, incipient structures adopted by tau proteins en route to fibrils, and the structure and hydration profile of a transmembrane peptide. This chapter focuses on offering a mechanistic understanding of the ODNP measurement and the molecular dynamics encoded in the ODNP parameters. In particular, it clarifies how the electron-nuclear dipolar coupling encodes information about the molecular dynamics in the nuclear spin self-relaxation and, more importantly, the electron-nuclear spin cross-relaxation rates. The clarification of the molecular dynamics underlying ODNP should assist in establishing a connection to theory and computer simulation that will offer far richer interpretations of ODNP results in future studies.

Facilitated and Non-Gaussian Diffusion of Cholesterol in Liquid Ordered Phase Bilayers Depends on the Flip-Flop and Spatial Arrangement of Cholesterol

The diffusion of cholesterol in biological membranes is critical to cellular processes such as the formation of cholesterol-enriched domains. The cholesterol diffusion may be complicated especially when cholesterol flip-flops and/or stays at the membrane center. Understanding the diffusion mechanism of cholesterol at a molecular level should be, therefore, a topic of interest. We perform molecular dynamics simulations up to 100 μs for lipid bilayers with various concentrations of cholesterol. We find that cholesterol diffusion in the liquid ordered phase depends on whether it is within leaflets or at the bilayer center, is non-Gaussian for several microseconds, and is enhanced significantly compared to that of lipids. Cholesterol at the bilayer center diffuses fast, while cholesterol in the hydrocarbon region with upright orientation diffuses relatively slowly. Such position-dependent dynamics of cholesterol leads to facilitated and non-Gaussian diffusion.

Twist and turn: a revised structural view on the unpaired bubble of class II CPD photolyase in complex with damaged DNA

Cyclobutane pyrimidine dimer (CPD) photolyases harness the energy of blue light to repair UV-induced DNA CPDs. Upon binding, CPD photolyases cause the photodamage to flip out of the duplex DNA and into the catalytic site of the enzyme. This process, called base-flipping, induces a kink in the DNA, as well as an unpaired bubble, which are stabilized by a network of protein-nucleic acid interactions. Previously, several co-crystal structures have been reported in which the binding mode of CPD photolyases has been studied in detail. However, in all cases the internucleoside linkage of the photodamage site was a chemically synthesized formacetal analogue and not the natural phosphodiester. Here, the first crystal structure and conformational analysis via molecular-dynamics simulations of a class II CPD photolyase in complex with photodamaged DNA that contains a natural cyclobutane pyrimidine dimer with an intra-lesion phosphodiester linkage are presented. It is concluded that a highly conserved bubble-intruding region (BIR) mediates stabilization of the open form of CPD DNA when complexed with class II CPD photolyases.

Correlating interfacial water dynamics with protein-protein interaction in complex of GDF-5 and BMPRI receptors

GDF-5 mediated signal transduction regulating chondrogenesis and skeletogenesis involves three different type-I receptors viz. Act-RI, BMPRIA and BMPRIB. BMPRIA and BMPRIB generally shows temporal and spatial co-expression but some spatially different expression pattern has also been observed. BMPRIA receptor is the key receptor implicated in BMP signalling during osteogenesis and is expressed in osteoblasts during the course of bone formation. However, BMPRIB appears to be primarily expressed in mesenchymal pre-cartilage condensations and also found in differentiated osteoblast and chondrocytes. The extracellular pH affects bone cell function and it is experimentally known that mineralization of bone is affected by shift of pH in cultured osteoblast. Here we report the effect of pH on dynamics of water present at the interface of GDF-5:BMPRIA and GDF-5:BMPRIB and binding interaction energy of these complexes. Water dynamics at different pH was analysed using residence time and hydrogen bond relaxation kinetics. pH influences the interaction energy between GDF-5 and BMPRIA and BMPRIB receptors indicating the electrostatic environment modulating the activity of two receptors. This pH dependence of interaction energy is further supported by similar behaviour of hydrogen bond existence of buried water molecules at the interface. In contrast to this the slow and fast exchanging water molecules do not show similar pH dependence of hydrogen bonding relaxation kinetics. Hence; we conclude that only buried water molecule at the interface influences the protein-protein interaction and the electrostatic environment of the extracellular fluid might decide the specificity of the two receptors.

Water Dynamics in the Hydration Shells of Biomolecules

The structure and function of biomolecules are strongly influenced by their hydration shells. Structural fluctuations and molecular excitations of hydrating water molecules cover a broad range in space and time, from individual water molecules to larger pools and from femtosecond to microsecond time scales. Recent progress in theory and molecular dynamics simulations as well as in ultrafast vibrational spectroscopy has led to new and detailed insight into fluctuations of water structure, elementary water motions, electric fields at hydrated biointerfaces, and processes of vibrational relaxation and energy dissipation. Here, we review recent advances in both theory and experiment, focusing on hydrated DNA, proteins, and phospholipids, and compare dynamics in the hydration shells to bulk water.