Intra and Inter-Molecular Communications Through Protein Structure Network

Communication within and across proteins is crucial for the biological functioning of proteins. Experiments such as mutational studies on proteins provide important information on the amino acids, which are crucial for their function. However, the protein structures are complex and it is unlikely that the entire responsibility of the function rests on only a few amino acids. A large fraction of the protein is expected to participate in its function at some level or other. Thus, it is relevant to consider the protein structures as a completely connected network and then deduce the properties, which are related to the global network features. In this direction, our laboratory has been engaged in representing the protein structure as a network of non-covalent connections and we have investigated a variety of problems in structural biology, such as the identification of functional and folding clusters, determinants of quaternary association and characterization of the network properties of protein structures. We have also addressed a few important issues related to protein dynamics, such as the process of oligomerization in multimers, mechanism on protein folding, and ligand induced communications (allosteric effect). In this review we highlight some of the investigations which we have carried out in the recent past. A review on protein structure graphs was presented earlier, in which the focus was on the graphs and graph spectral properties and their implementation in the study of protein structure graphs/networks (PSN). In this article, we briefly summarize the relevant parts of the methodology and the focus is on the advancement brought out in the understanding of protein structure-function relationships through structure networks. The investigations of structural/biological problems are divided into two parts, in which the first part deals with the analysis of PSNs based on static structures obtained from x-ray crystallography. The second part highlights the changes in the network, associated with biological functions, which are deduced from the network analysis on the structures obtained from molecular dynamics simulations.

The Influence of Bilayer Composition on the Gel to Liquid Crystalline Transition

We report molecular dynamics simulations of bilayers using a united atom model with explicit solvent molecules. The bilayer consists of the single tail cationic surfactant behenyl trimethyl ammonium chloride (BTMAC) with stearyl alcohol (SA) as the cosurfactant. We study the gel to liquid crystalline transitions in the bilayer by varying the amount of water at fixed BTMAC to SA ratio as well as by varying the BTMAC to SA ratio at fixed water content. The bilayer is found to exist in the tilted, Lβ′ phase at low temperatures, and for the compositions investigated in this study, the Lβ′ to Lα melting transition occurred in the temperature range 330−338 K. For the highest BTMAC to SA composition (2:3 molar ratio), a diffuse headgroup−water interface is observed at lower temperatures, and an increase in the d-spacing occurs prior to the melting transition. This pretransition swelling is accompanied by a sharpening in the water density variation across the headgroup region of the bilayer. Signatures of this swelling effect which can be observed in the alkane density distributions, area per headgroup, and membrane thickness are attributed to the hydrophobic effect. At a fixed bilayer composition, the transition temperature (>338 K) from the Lβ′ to Lα transition obtained for the high water content bilayer (80 wt %) is similar to that obtained with low water content (54.3 wt %), confirming that the melting transition at these water contents is dominated by chain melting.

Reorientation of Water Inside Carbon Nanorings by Large Angular Jumps

Molecular dynamics simulations of the orientational dynamics of water molecules confined inside narrow carbon nanorings reveal that reorientational relaxation is mediated by large amplitude angular jumps. The distribution of waiting time between jumps peaks at about 60 fs, and has a slowly decaying exponential tail with a timescale of about 440 fs. These time scales are much faster than the mean waiting time between jumps of the water molecules in bulk.

Enhanced Tetrahedral Ordering of Water Molecules in Minor Grooves of DNA: Relative Role of DNA Rigidity, Nanoconfinement, and Surface Specific Interactions

Confinement and Surface specific interactions call induce Structures otherwise unstable at that temperature and pressure. Here we Study the groove specific water dynamics ill the nucleic acid sequences, poly-AT and poly-GC, in long B-DNA duplex chains by large scale atomistic molecular dynamics simulations, accompanied by thermodynamic analysis. While water dynamics in the major groove remains insensitive to the sequence differences, exactly the opposite is true for the minor groove water. Much slower water dynamics observed in the minor grooves (especially in the AT minor) call be attributed to all enhanced tetrahedral ordering (< t(h)>) of water. The largest value of < t(h)> in the AT minor groove is related to the spine of hydration found in X-ray Structure. The calculated configurational entropy (S-C) of the water molecules is found to be correlated with the self-diffusion coefficient of water in different region via Adam-Gibbs relation D = A exp(-B/TSC), and also with < t(h)>.

Vibrational phase relaxation of O–H stretch in bulk water: Role of large amplitude angular jumps and negative cross-correlations among the forces on the O–H bond

The ultrafast vibrational phase relaxation of O–H stretch in bulk water is investigated in molecular dynamics simulations. The dephasing time (T2) of the O–H stretch in bulk water calculated from the frequency fluctuation time correlation function (Cω(t)) is in the range of 70–80 femtosecond (fs), which is comparable to the characteristic timescale obtained from the vibrational echo peak shift measurements using infrared photon echo [W.P. de Boeij, M.S. Pshenichnikov, D.A. Wiersma, Ann. Rev. Phys. Chem. 49 (1998) 99]. The ultrafast decay of Cω(t) is found to be responsible for the ultrashort T2 in bulk water. Careful analysis reveals the following two interesting reasons for the ultrafast decay of Cω(t). (A) The large amplitude angular jumps of water molecules (within 30–40 fs time duration) provide a large scale contribution to the mean square vibrational frequency fluctuation and gives rise to the rapid spectral diffusion on 100 fs time scale. (B) The projected force, due to all the atoms of the solvent molecules on the oxygen (FO(t)) and hydrogen (FH(t)) atom of the O–H bond exhibit a large negative cross-correlation (NCC). We further find that this NCC is partly responsible for a weak, non-Arrhenius temperature dependence of the dephasing rate.

Single-File Diffusion of Confined Water Inside SWNTs: An NMR Study

We report a nuclear magnetic resonance (NMR) study of confined water inside similar to 1.4 nm diameter single-walled carbon nanotubes (SWNTs). We show that the confined water does not freeze even up to 223 K. A pulse field gradient (PFG) NMR method is used to determine the mean squared displacement (MSD) of the water molecules inside the nanotubes at temperatures below 273 K, where the bulk water outside the nanotubes freezes and hence does not contribute to the proton NMR signal. We show that the mean squared displacement varies as the square root of time, predicted for single-file diffusion in a one-dimensional channel. We propose a qualitative understanding of our results based on available molecular dynamics simulations.

Hydration of ions under confinement

Molecular dynamics simulations are used to examine the changes in water density and hydration characteristics of NaCl solutions confined in slit-shaped graphitic pores. Using a structural signature, we define the hydration limit as the salt concentration at which a sharp drop in the hydration number is observed. At small pores (H = 8.0-10 angstrom), confined water does not possess bulk-like features and remains in a layered arrangement between two surfaces. Despite this high degree of confinement, ions are able to form a quasi-2D hydration shell between two surfaces. Our results indicate the strong propensity of ions to form the first hydration shell, even under extremely confined aqueous environments. The hydration of ions is seen to weakly perturb the oxygen density distributions between two surfaces. The hydration number of Na+ reduces to about 4.15 at a pore width of H = 0.8 nm, when compared with the bulk hydration number of 6.25. At larger pore widths, above H = 16 angstrom, where bulk-like water densities are observed in the central regions of the pore, the hydration number is above 6.

Elucidation of the conformational free energy landscape in H.pylori LuxS and its implications to catalysis

Background: One of the major challenges in understanding enzyme catalysis is to identify the different conformations and their populations at detailed molecular level in response to ligand binding/environment. A detail description of the ligand induced conformational changes provides meaningful insights into the mechanism of action of enzymes and thus its function. Results: In this study, we have explored the ligand induced conformational changes in H. pylori LuxS and the associated mechanistic features. LuxS, a dimeric protein, produces the precursor (4,5-dihydroxy-2,3-pentanedione) for autoinducer-2 production which is a signalling molecule for bacterial quorum sensing. We have performed molecular dynamics simulations on H. pylori LuxS in its various ligand bound forms and analyzed the simulation trajectories using various techniques including the structure network analysis, free energy evaluation and water dynamics at the active site. The results bring out the mechanistic details such as co operativity and asymmetry between the two subunits, subtle changes in the conformation as a response to the binding of active and inactive forms of ligands and the population distribution of different conformations in equilibrium. These investigations have enabled us to probe the free energy landscape and identify the corresponding conformations in terms of network parameters. In addition, we have also elucidated the variations in the dynamics of water co-ordination to the Zn2+ ion in LuxS and its relation to the rigidity at the active sites. Conclusions: In this article, we provide details of a novel method for the identification of conformational changes in the different ligand bound states of the protein, evaluation of ligand-induced free energy changes and the biological relevance of our results in the context of LuxS structure-function. The methodology outlined here is highly generalized to illuminate the linkage between structure and function in any protein of known structure.

Influence of non-geometrical factors on intracrystalline diffusion -- Role of sorbate-zeolite interactions

Molecular dynamics simulations on Xe in NaY and Ar in NaCaA zeolite are reported. Rates of cage-to-cage crossovers in the two zeolites exhibit trends which are contrary to that expected from geometrical considerations. The results suggest the important role of the sorbate-zeolite interactions in determining the molecular sieve properties of zeolites for small sized sorbates. The results are explained in terms of the barrier height for cage-to-cage crossover in the two zeolites.

DNA Compaction by a Dendrimer

At physiological pH, a PAMAM dendrimer is positively charged and can effectively bind negatively charged DNA. Currently, there has been great interest in understanding this complexation reaction both for fundamental (as a model for complex biological reactions) as well as for practical (as a gene delivery material and probe for sensing DNA sequence) reasons. Here, we have studied the complexation between double-stranded DNA (dsDNA) and various generations of PAMAM dendrimers (G3-05) through atomistic molecular dynamics simulations in the presence of water and ions. We report the compaction of DNA on a nanosecond time scale. This is remarkable, given the fact that such a short DNA duplex with a length close to 13 nm is otherwise thought to be a rigid rod. Using several nanoseconds long MD simulations, we have observed various binding modes of dsDNA and dendrimers for various generations of PAMAM dendrimers at varying charge ratios, and it confirms some of the binding modes proposed earlier. The binding is driven by the electrostatic interaction, and the larger the dendrimer charge, the stronger the binding affinity. As DNA wraps/binds to the dendrimer, counterions originally condensed onto DNA (Na+) and the dendrimer (Cl-) get released. We calculate the entropy of counterions and show that there is gain in entropy due to counterion release during the complexation. MD simulations demonstrate that, when the charge ratio is greater than 1 (as in the case of the G5 dendrimer), the optimal wrapping of DNA is observed. Calculated binding energies of the complexation follow the trend G5 > 04 > 03, in accordance with the experimental data. For a lower-generation dendrimer, such as G3, and, to some extent, for G4 also, we see considerable deformation in the dendrimer structure due to their flexible nature. We have also calculated the various helicoidal parameters of DNA to study the effect of dendrimer binding on the structure of DNA. The B form of the DNA is well preserved in the complex, as is evident from various helical parameters, justifying the use of the PAMAM dendrimer as a suitable delivery vehicle.

Inherent structures of phase-separating binary mixtures: Nucleation, spinodal decomposition, and pattern formation

An energy landscape view of phase separation and nonideality in binary mixtures is developed by exploring their potential energy landscape (PEL) as functions of temperature and composition. We employ molecular dynamics simulations to study a model that promotes structure breaking in the solute-solvent parent binary liquid, at low temperatures. The PEL of the system captures the potential energy distribution of the inherent structures (IS) of the system and is obtained by removing the kinetic energy (including that of intermolecular vibrations). The broader distribution of the inherent structure energy for structure breaking liquid than that of the structure making liquid demonstrates the larger role of entropy in stabilizing the parent liquid of the structure breaking type of binary mixtures. At high temperature, although the parent structure of the structure breaking binary mixture is homogenous, the corresponding inherent structure is found to be always phase separated, with a density pattern that exhibits marked correlation with the energy of its inherent structure. Over a broad range of intermediate inherent structure energy, bicontinuous phase separation prevails with interpenetrating stripes as signatures of spinodal decomposition. At low inherent structure energy, the structure is largely phase separated with one interface where as at high inherent structure energy we find nucleation type growth. Interestingly, at low temperature, the average inherent structure energy (< EIS >) exhibits a drop with temperature which signals the onset of crystallization in one of the phases while the other remains in the liquid state. The nonideal composition dependence of viscosity is anticorrelated with average inherent structure energy.

Interaction of single-walled carbon nanotubes with poly(propyl ether imine) dendrimers

We study the complexation of nontoxic, native poly(propyl ether imine) dendrimers with single-walled carbon nanotubes (SWNTs). The interaction was monitored by measuring the quenching of inherent fluorescence of the dendrimer. The dendrimer-nanotube binding also resulted in the increased electrical resistance of the hole doped SWNT, due to charge-transfer interaction between dendrimer and nanotube. This charge-transfer interaction was further corroborated by observing a shift in frequency of the tangential Raman modes of SWNT. We also report the effect of acidic and neutral pH conditions on the binding affinities. Experimental studies were supplemented by all atom molecular dynamics simulations to provide a microscopic picture of the dendrimer-nanotube complex. The complexation was achieved through charge transfer and hydrophobic interactions, aided by multitude of oxygen, nitrogen, and n-propyl moieties of the dendrimer. (C) 2011 American Institute of Physics. doi:10.1063/1.3561308]

Collective Effects on Single Particle Orientational Relaxation in Slow Dipolar Liquids

Theoretical and computer simulation studies of orientational relaxation in dense molecular liquids are presented. The emphasis of the study is to understand the effects of collective orientational relaxation on the single-particle orientational dynamics. The theoretical analysis is based on a recently developed molecular hydrodynamic theory which allows a self-consistent description of both the collective and the single-particle orientational relaxation. The molecular hydrodynamic theory can be used to derive a relation between the memory function for the collective orientational correlation function and the frequency-dependent dielectric function. A novel feature of the present work is the demonstration that this collective memory function is significantly different from the single-particle rotational friction. However, a microscopic expression for the single-particle rotational friction can be derived from the molecular hydrodynamic theory where the collective memory function can be used to obtain the single-particle orientational friction. This procedure allows, us to calculate the single-particle orientational correlation function near the alpha-beta transition in the supercooled liquid. The calculated correlation function shows an interesting bimodal decay below the bifurcation temperature as the glass transition is approached from above. Brownian dynamics simulations have been carried out to check the validity of the above procedure of translating the memory function from the dielectric relaxation data. We have also investigated the following two issues important in understanding the orientational relaxation in slow liquids. First, we present an analysis of the ''orientational caging'' of translational motion. The value of the translational friction is found to be altered significantly by the orientational caging. Second, we address the question of the rank dependence of the dielectric friction using both simulation and the molecular hydrodynamic theory.

Orientational relaxation in dipolar systems: How much do we understand the role of correlations

Dipolar systems, both liquids and solids, constitute a class of naturally abundant systems that are important in all branches of natural science. The study of orientational relaxation provides a powerful method to understand the microscopic properties of these systems and, fortunately, there are many experimental tools to study orientational relaxation in the condensed phases. However, even after many years of intense research, our understanding of orientational relaxation in dipolar systems has remained largely imperfect. A major hurdle towards achieving a comprehensive understanding is the long range and complex nature of dipolar interactions which also made reliable theoretical study extremely difficult. These difficulties have led to the development of continuum model based theories, which although they provide simple, elegant expressions for quantities of interest, are mostly unsatisfactory as they totally neglect the molecularity of inter-molecular interactions. The situation has improved in recent years because of renewed studies, led by computer simulations. In this review, we shall address some of the recent advances, with emphasis on the work done in our laboratory at Bangalore. The reasons for the failure of the continuum model, as revealed by the recent Brownian dynamics simulations of the dipolar lattice, are discussed. The main reason is that the continuum model predicts too fast a decay of the torque-torque correlation function. On the other hand, a perturbative calculation, based on Zwanzig's projection operator technique, provides a fairly satisfactory description of the single particle orientational dynamics for not too strongly polar dipolar systems. A recently developed molecular hydrodynamic theory that properly includes the effects of intermolecular orientational pair correlations provides an even better description of the single-particle orientational dynamics. We also discuss the rank dependence of the dielectric friction. The other topics reviewed here includes dielectric relaxation and solvation dynamics, as they are intimately connected with orientational relaxation. Recent molecular dynamics simulations of the dipolar lattice are also discussed. The main theme of the present review is to understand the effects of intermolecular interactions on orientational relaxation. The presence of strong orientational pair correlation leads to a strong coupling between the single particle and the collective dynamics. This coupling can lead to rich dynamical properties, some of which are detailed here, while a major part remains yet unexplored.

Modeling of angiogenin-3 '-NMP complex

Angiogenin belongs to the Ribonuclease superfamily and has a weak enzymatic activity that is crucial for its biological function of stimulating blood vessel growth. Structural studies on ligand bound Angiogenin will go a long way in understanding the mechanism of the protein as well as help in designing drugs against it. In this study we present the first available structure of nucleotide ligand bound Angiogenin obtained by computer modeling. The importance of this study in itself notwithstanding, is a precursor to modeling a full dinucleotide substrate onto Angiogenin. Bovine Angiogenin, the structure of which has been solved at a high resolution, was earlier subjected to Molecular Dynamics simulations for a nanosecond. The MD structures offer better starting points for docking as they offer lesser obstruction than the crystal structure to ligand binding. The MD structure with the least serious short contacts was modeled to obtain a steric free Angiogenin - 3' mononucleotide complex structure. The structures were energetically minimized and subjected to a brief spell of Molecular Dynamics. The results of the simulation show that all the li,ligand-Angiogenin interactions and hydrogen bonds are retained, redeeming the structure and docking procedure. Further, following ligand - protein interactions in the case of the ligands 3'-CMP and 3'-UMP we were able to speculate on how Angiogenin, a predominantly prymidine specific ribonuclease prefers Cytosine to Uracil in the first base position.