Conserved water-mediated H-bonding dynamics of catalytic Asn 175 in plant thiol protease

The role of invariant water molecules in the activity of plant cysteine protease is ubiquitous in nature. On analysing the 11 different Protein DataBank (PDB) structures of plant thiol proteases, the two invariant water molecules W I and W2 (W220 and W222 in the template 1PPN structure) were observed to form H-bonds with the Ob atom of Asn 175. Extensive energy minimization and molecular dynamics simulation studies up to 2 ns on all the PDB and solvated structures clearly revealed the involvement of the H-bonding association of the two water molecules in fixing the orientation of the asparagine residue of the catalytic triad. From this study, it is suggested that H-bonding of the water molecule at the W1 invariant site better stabilizes the Asn residue at the active site of the catalytic triad.

Brownian particles in stationary and moving traps: The mean and variance of the heat distribution function

A recent theoretical model developed by Imparato et al. Phys of the experimentally measured heat and work effects produced by the thermal fluctuations of single micron-sized polystyrene beads in stationary and moving optical traps has proved to be quite successful in rationalizing the observed experimental data. The model, based on the overdamped Brownian dynamics of a particle in a harmonic potential that moves at a constant speed under a time-dependent force, is used to obtain an approximate expression for the distribution of the heat dissipated by the particle at long times. In this paper, we generalize the above model to consider particle dynamics in the presence of colored noise, without passing to the overdamped limit, as a way of modeling experimental situations in which the fluctuations of the medium exhibit long-lived temporal correlations, of the kind characteristic of polymeric solutions, for instance, or of similar viscoelastic fluids. Although we have not been able to find an expression for the heat distribution itself, we do obtain exact expressions for its mean and variance, both for the static and for the moving trap cases. These moments are valid for arbitrary times and they also hold in the inertial regime, but they reduce exactly to the results of Imparato et al. in appropriate limits. DOI: 10.1103/PhysRevE.80.011118 PACS.

Interaction of substrate uridyl 3',5'-adenosine with ribonuclease A:a molecular dynamics study

A wealth of information available from x-ray crystallographic structures of enzyme-ligand complexes makes it possible to study interactions at the molecular level. However, further investigation is needed when i) the binding of the natural substrate must be characterized, because ligands in the stable enzyme-ligand complexes are generally inhibitors or the analogs of substrate and transition state, and when ii) ligand binding is in part poorly characterized. We have investigated these aspects i? the binding of substrate uridyl 3',5'-adenosine (UpA) to ribonuclease A (RNase A). Based on the systematically docked RNase A-UpA complex resulting from our previous study, we have undertaken a molecular dynamics simulation of the complex with solvent molecules. The molecular dynamics trajectories of this complex are analyzed to provide structural explanations for varied experimental observations on the ligand binding at the B2 subsite of ribonuclease A. The present study suggests that B2 subsite stabilization can be effected by different active site groups, depending on the substrate conformation. Thus when adenosine ribose pucker is O4'-endo, Gln69 and Glu111 form hydrogen-bonding contacts with adenine base, and when it is C2'-endo, Asn71 is the only amino acid residue in direct contact with this base. The latter observation is in support of previous mutagenesis and kinetics studies. Possible roles for the solvent molecules in the binding subsites are described. Furthermore, the substrate conformation is also examined along the simulation pathway to see if any conformer has the properties of a transition state. This study has also helped us to recognize that small but concerted changes in the conformation of the substrate can result in substrate geometry favorable for 2',3' cyclization. The identified geometry is suitable for intraligand proton transfer between 2'-hydroxyl and phosphate oxygen atom. The possibility of intraligand proton transfer as suggested previously and the mode of transfer before the formation of cyclic intermediate during transphosphorylation are discussed.

Structures and molecular-dynamics studies of three active-site mutants of bovine pancreatic phospholipase A(2)

Phospholipase A(2) hydrolyzes phospholipids at the sn-2 position to cleave the fatty-acid ester bond of L-glycerophospholipids. The catalytic dyad (Asp99 and His48) along with a nucleophilic water molecule is responsible for enzyme hydrolysis. Furthermore, the residue Asp49 in the calcium-binding loop is essential for controlling the binding of the calcium ion and the catalytic action of phospholipase A2. To elucidate the structural role of His48 and Asp49, the crystal structures of three active-site single mutants H48N, D49N and D49K have been determined at 1.9 angstrom resolution. Although the catalytically important calcium ion is present in the H48N mutant, the crystal structure shows that proton transfer is not possible from the catalytic water to the mutated residue. In the case of the Asp49 mutants, no calcium ion was found in the active site. However, the tertiary structures of the three active-site mutants are similar to that of the trigonal recombinant enzyme. Molecular-dynamics simulation studies provide a good explanation for the crystallographic results.

A molecular dynamics study of a PbO.SiO2 glass and melt

The structures of a PbO.SiO2 glass and melt have been studied using molecular dynamics simulation employing Born-Mayer-Huggins pair potentials. Various pair distribution functions are presented and discussed. Pb-Pb correlations persist in the melt, in agreement with experimental observations. The calculated and experimental radial distribution functions are compared.

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.

The stability of transmembrane helices: A molecular dynamics study on the isolated helices of bacteriorhodopsin

Bacteriorhodopsin (bR) continues to be a proven testing ground for the study of integral membrane proteins (IMPs). It is important to study the stability of the individual helices of bR, as they are postulated to exist as independently stable transmembmne helices (TMHs) and also for their utility as templates for modeling other IMPs with the postulated seven-helix bundle topology. Toward this purpose, the seven helices of bR have been studied by molecular dynamics simulation in this study. The suitability of using the backbone-dependent rotamer library of side-chain conformations arrived at from the data base of globular protein structures in the case TMHs has been tested by another set of ? helix simulations with the side-chain orientations taken from this library. The influence of the residue's net charge oil the helix stability was examined by simulating the helices III, IV, and VI (from both of the above sets of helices) with zero net charge on the side chains. The results of these 20 simulations demonstrate in general the stability of the isolated helices of bR in conformity with the two-stage hypothesis of IMP folding. However, the helices I, II, V, and VII are more stable than the other three helices. The helical nature of certain regions of III, IV, and VI are influenced by factors such as the net charge and orientation of several residues. It is seen that the residues Arg, Lys, Asp, and Glu (charged residues), and Ser, Thr, Gly, and Pro, play a crucial role in the stability of the helices of bR. The backbone-dependent rotamer library for the side chains is found to be suitable for the study of TMHs in IMP. (C) 1996 John Wiley & Sons, Inc.

Structure and silver ion transport in AgI-Ag2MoO4 glasses: A molecular dynamics study

The structure and dynamics of silver ion conducting AgI-Ag2MoO4 glasses have been simulated by molecular dynamics simulation over a wide range of compositions. Formation of silver iodide like aggregates have been identified only in the AgI rich glasses. Increase in silver ion conductivity with an increase in AgI content in the glass is seen as in experiments. The dynamics of ion transport suggests that Ag+ ion transport occurs largely through paths connected by silver ion sites of mixed iodide-oxide coordination. The Van Hove correlation functions indicate that Ag+ ions prefer migration along the pathways formed with connected sites of similar coordination.

Slow orientational dynamics of water molecules at a micellar surface

The dynamics of water molecules near an aqueous micellar interface is studied in an atomistic molecular dynamics simulation of cesium pentadecafluorooctanoate (CsPFO) in water. The dipolar orientational time correlation function (tcf) and the translational diffusion of the water molecules are investigated. Results show that both the reorientational and the translational motion of water molecules near the micelle are restricted. In particular, the orientational tcf exhibits a very slow component in the long time which is slower than its bulk value by 2 orders of magnitude. This slow decay seems to be related to the slow decay often observed in experiments. The origin of the slow decay is analyzed.

Effect of orientational motion of mobile chromophores on the dynamics of Forster energy transfer in polymers

Brownian dynamics (BD) simulations have been carried out to explore the effects of the orientational motion of the donor-acceptor (D-A) chromophore pair on the Forster energy transfer between the D-A pair embedded in a polymer chain in solution. It is found that the usually employed orientational averaging (that is, replacing the orientational factor, kappa, by kappa (2) = 2/3) may lead to an error in the estimation of the rate of the reaction by about 20%. In the limit of slow orientational relaxation, the preaveraging of the orientational factor leads to an overestimation of the rate, while in the opposite limit of very fast orientational relaxation, the usual scheme underestimates the rate. The latter results from an interesting interplay between reaction and diffusion. On the other hand, when one of the chromophores is fixed, the preaveraged rate is found to be fairly reliable if the rotational relaxation of the chromophore is sufficiently fast. The present study also reveals a power law dependence of the FRET rate on the chain length (rate proportional to N- alpha, with alpha approximate to 2.6).

Hydrophobic hydration driven self-assembly of curcumin in water: Similarities to nucleation and growth under large metastability, and an analysis of water dynamics at heterogeneous surfaces

As the beneficial effects of curcumin have often been reported to be limited to its small concentrations, we have undertaken a study to find the aggregation properties of curcumin in water by varying the number of monomers. Our molecular dynamics simulation results show that the equilibrated structure is always an aggregated state with remarkable structural rearrangements as we vary the number of curcumin monomers from 4 to 16 monomers. We find that the curcumin monomers form clusters in a very definite pattern where they tend to aggregate both in parallel and anti-parallel orientation of the phenyl rings, often seen in the formation of beta-sheet in proteins. A considerable enhancement in the population of parallel alignments is observed with increasing the system size from 12 to 16 curcumin monomers. Due to the prevalence of such parallel alignment for large system size, a more closely packed cluster is formed with maximum number of hydrophobic contacts. We also follow the pathway of cluster growth, in particular the transition from the initial segregated to the final aggregated state. We find the existence of a metastable structural intermediate involving a number of intermediate-sized clusters dispersed in the solution. We have constructed a free energy landscape of aggregation where the metatsable state has been identified. The course of aggregation bears similarity to nucleation and growth in highly metastable state. The final aggregated form remains stable with the total exclusion of water from its sequestered hydrophobic core. We also investigate water structure near the cluster surface along with their orientation. We find that water molecules form a distorted tetrahedral geometry in the 1st solvation layer of the cluster, interacting rather strongly with the hydrophilic groups at the surface of the curcumin. The dynamics of such quasi-bound water molecules near the surface of curcumin cluster is considerably slower than the bulk signifying a restricted motion as often found in protein hydration layer. (C) 2014 AIP Publishing LLC.

Impact of glycosylation on stability, structure and unfolding of soybean agglutinin (SBA): an insight from thermal perturbation molecular dynamics simulations

Glycosylation has been recognized as one of the most prevalent and complex post-translational modifications of proteins involving numerous enzymes and substrates. Its effect on the protein conformational transitions is not clearly understood yet. In this study, we have examined the effect of glycosylation on protein stability using molecular dynamics simulation of legume lectin soybean agglutinin (SBA). Its glycosylated moiety consists of high mannose type N-linked glycan (Man(9)GlcNAc(2)). To unveil the structural perturbations during thermal unfolding of these two forms, we have studied and compared them to the experimental results. From the perspective of dynamics, our simulations revealed that the nonglycosylated monomeric form is less stable than corresponding glycosylated form at normal and elevated temperatures. Moreover, at elevated temperature thermal destabilization is more prominent in solvent exposed loops, turns and ends of distinct beta sheets. SBA maintains it folded structure due to some important saltbridges, hydrogen bonds and hydrophobic interactions within the protein. The reducing terminal GlcNAc residues interact with the protein residues VAL161, PRO182 and SER225 via hydrophobic and via hydrogen bonding with ASN 9 and ASN 75. Our simulations also revealed that single glycosylation (ASN75) has no significant effect on corresponding cis peptide angle orientation. This atomistic description might have important implications for understanding the functionality and stability of Soybean agglutinin.

Diffusion of flexible, charged, nanoscopic molecules in solution: Size and pH dependence for PAMAM dendrimer

In order to understand self-diffusion (D) of a charged, flexible, and porous nanoscopic molecule in water, we carry out very long, fully atomistic molecular dynamics simulation of PAMAM dendrimer up to eight generations in explicit salt water under varying pH. We find that while the radius of gyration (R-g) varies as N-1/3, the self-diffusion constant (D) scales, surprisingly, as N-alpha, with alpha=0.39 at high pH and 0.5 at neutral pH, indicating a dramatic breakdown of Stokes-Einstein relation for diffusion of charged nanoscopic molecules. The variation in D as a function of radius of gyration demonstrates the importance of treating water and ions explicitly in the diffusion process of a flexible nanoscopic molecule. In agreement with recent experiments, the self-diffusion constant increases with pH, revealing the importance of dielectric friction in the diffusion process. The shape of a dendrimer is found to fluctuate on a nanosecond time scale. We argue that this flexibility (and also the porosity) of the dendrimer may play an important role in determining the mean square displacement of the dendrimer and the breakdown of the Stokes-Einstein relation between diffusion constant and the radius.

Correlation between conductivity or diffusivity and activation energy in amorphous solids

There exist many investigations of ionic transport in a variety of glasses. These studies exhibit strong correlation between ionic conductivity and activation energy: Typically, it is found that higher conductivity is associated with lower activation energies and vice versa. Although there are explanations for this at a phenomenological level, there is no consistent physical picture to explain the correlation between conductivity and activation energy. We have carried out molecular dynamics simulation as a function of the size of the impurity atom or diffusant (both neutral and charged) in a host amorphous matrix. We find that there is a maximum in self-diffusivity as a function of the size of the impurity atom suggesting that there is an appropriate size for which the diffusivity is maximum. The activation energy is found to be the lowest for this size of the impurity. A similar maximum has been previously found in other condensed phases, such as confined fluids and dense liquids, and has its origin in the levitation effect. The implications of this result for understanding ionic conductivity in glasses are discussed. Our results suggest that there is a relation between microscopic structure of the amorphous solid, diffusivity or conductivity, and activation energy. The nature of this relationship is discussed in terms of the levitation parameter showing that diffusivity is maximum when the size of the neck or doorway radius is comparable with the size of the diffusant. Our computational results here are in excellent agreement with independent experimental results of Nascimento et al. [Braz. J. Phys. 35, 626 (2005)] that structural features of the glass are important in determining the ionic conductivity.

The 5-Methyl Group in Thymine Dynamically Influences the Structure of A-Tracts in DNA at the Local and Global Level

DNA sequences containing a stretch of several A:T basepairs without a 5'-TA-3' step are known as A-tracts and have been the subject of extensive investigation because of their unique structural features such as a narrow minor groove and their crucial role in several biological processes. One of the aspects under investigation has been the influence of the 5-methyl group of thymine on the properties of A-tracts. Detailed molecular dynamics simulation studies of the sequences d(CGCAAAUUUGCG) and d(CGCAAATTTGCG) indicate that the presence of the 5-methyl group in thymine increases the frequency of a narrow minor groove conformation, which could facilitate its specific recognition by proteins, and reduce its susceptibility to cleavage by DNase I. The bias toward a wider minor groove in the absence of the thymine 5-methyl group is a static structural feature. Our results also indicate that the presence of the thymine 5-methyl group is necessary for calibrating the backbone conformation and the basepair and dinucleotide step geometry of the core A-tract as well as the flanking CA/TG and the neighboring GC/GC steps, as observed in free and protein-bound DNA. As a consequence, it also fine-tunes the curvature of the longer DNA fragment in which the A-tract is embedded.