Stability of Dimeric Interface in Banana Lectin: Insight from Molecular Dynamics Simulations

Banana lectin (Banlec) is a homodimeric non-glycosylated protein. It exhibits the b-prism I structure. High-temperature molecular dynamics simulations have been utilized to monitor and understand early stages of thermally induced unfolding of Banlec. The present study elucidates the behavior of the dimeric protein at four different temperatures and compares the structural and conformational changes to that of the minimized crystal structure. The process of unfolding was monitored by following the radius of gyration, the rms deviation of each residue, change in relative solvent accessibility and the pattern of inter- and intra-subunit interactions. The overall study demonstrates that the Banlec dimer is a highly stable structure, and the stability is mostly contributed by interfacial interactions. It maintains its overall conformation during high-temperature (400–500 K) simulations, with only the unstructured loop regions acquiring greater momentum under such condition. Nevertheless, at still higher temperatures (600 K) the tertiary structure is gradually lost which later extends to loss of secondary structural elements. The pattern of hydrogen bonding within the subunit and at the interface across different stages has been analyzed and has provided rationale for its intrinsic high stability.

Conformational Dynamics in Alkyl Chains of an Anchored Bilayer: A Molecular Dynamics Study

Molecular dynamics (MD) simulations are reported for an anchored bilayer formed by the intercalation of cetyl trimethyl ammonium (CTA) and CH3(CH2)15N+(CH3) ions in a layered solid, CdPS3. The intercalated CTA ions are organized with the cationic headgroups tethered to the inorganic sheet and the hydrocarbon tails arranged as bilayers. Simulations were performed at three temperatures, 65, 180, and 298 K, using an isothermal−isobaric ensemble that was subsequently switched once macroscopic parameters had converged to a canonical isothermal−isochoric ensemble. The simulations are able to reproduce the experimental features of this system, including the formation of the bilayer and layer-to-layer separation distance. An analysis of the conformation of the chains showed that at all three temperatures a fraction of the alkyl chains retained a planar all-trans conformation, and that gauche bonds occurred as part of a “kink” (gauche+−trans−gauche−) sequence and not as isolated gauche bonds. Trans−gauche isomerization rates for the alkyl chains in the anchored bilayer are slower than those in lipid bilayers at the same temperature and show a progressive increase as the torsion numbers approach the tail. A two-dimensional periodic Voronoi tessellation analysis was performed to obtain the single-molecular area of an alkyl chain in the bilayer. The single-molecular area relaxation times are an order of magnitude longer than the trans−gauche isomerization times. The results indicate that the trans−gauche isomerization is associated with the creation and annihilation of a kink defect sequence. The results of the present MD simulation explain the apparent conflicting estimates of the gauche disorder in this system as obtained from infrared and 13C nuclear magnetic resonance measurements.

Proton magnetic resonance study of molecular dynamics in (NH4)2CdI4

Proton magnetic resonance and spin-lattice relaxation studies have been carried out on (NH4)2CdI4 as a function of temperature (77–400 K) and Larmor frequency (10, 20 and 30 MHz). The T1 data indicate isotropic tumbling of ammonium ions at equivalent sites till 160 K. There is an indication of a phase transition at 265 K, the activation energy for molecular reorientation increases from 2.8 kcal/mole to 4.6 kcal/mole. The relaxation results and the linewidth data support the presence of two inequivalent sites at low temperatures, one having an environment corresponding to near-rigid-lattice limit and the other undergoing fast reorientations. The behaviour of the free induction decay with temperature below 120 K suggests a coherent motion for the faster species.

Molecular dynamics simulations of sliding in an Fe-Cu tribopair system

Experimental evidence suggests that high strain rates, stresses, strains and temperatures are experienced near sliding interfaces. The associated microstructural changes are due to several dynamic an interacting phenomena. 3D non-equilibrium molecular dynamics (MD) simulations of sliding were conducted with the aim of understanding the dynamic processes taking place in crystalline tribopairs, with a focus on plastic deformation and microstructural evolution. Embedded atom potentials were employed for simulating sliding of an Fe-Cu tribopair. Sliding velocity, crystal orientation and presence of lattice defects were some of the variables in these simulations. Extensive plastic deformation involving dislocation and twin activity, dynamic recrystallization, amorphization and/or nanocrystallization, mechanical mixing and material transfer were observed. Mechanical mixing in the vicinity of the sliding interface was observed even in the Fe-Cu system, which would cluster under equilibrium conditions, hinting at the ballistic nature of the process. Flow localization was observed at high velocities implying the possible role of adiabatic heating. The presence of preexisting defects (such as dislocations and interfaces) played a pivotal role in determining friction and microstructural evolution. The study also shed light on the relationship between adhesion and plastic deformation, and friction. Comparisons with experiments suggest that such simulations can indeed provide valuable insights that are difficult to obtain from experiments.

Neutron Scattering and Molecular Dynamics Evidence for Levitation Effect in Nanopores

Neutron Scattering and Molecular Dynamics Evidence for Levitation Effect in Nanopores ... Neutron scattering measurements and molecular dynamics simulations have been carried out on the three isomers of pentane (neopentane (neo), isopentane (iso), and n-pentane (n-)) adsorbed in zeolite NaY. ... In order to understand this surprising dependence, the dimensionless levitation parameter, γ, for atomic systems may be modified to suit molecular systems.

Can current force fields reproduce ring puckering in 2-O-sulfo-α-l-iduronic acid? A molecular dynamics simulation study

The monosaccharide 2-O-sulfo-α-l-iduronic acid (IdoA2S) is one of the major components of glycosaminoglycans. The ability of molecular mechanics force fields to reproduce ring-puckering conformational equilibrium is important for the successful prediction of the free energies of interaction of these carbohydrates with proteins. Here we report unconstrained molecular dynamics simulations of IdoA2S monosaccharide that were carried out to investigate the ability of commonly used force fields to reproduce its ring conformational flexibility in aqueous solution. In particular, the distribution of ring conformer populations of IdoA2S was determined. The GROMOS96 force field with the SPC/E water potential can predict successfully the dominant skew-boat to chair conformational transition of the IdoA2S monosaccharide in aqueous solution. On the other hand, the GLYCAM06 force field with the TIP3P water potential sampled transitional conformations between the boat and chair forms. Simulations using the GROMOS96 force field showed no pseudorotational equilibrium fluctuations and hence no inter-conversion between the boat and twist boat ring conformers. Calculations of theoretical proton NMR coupling constants showed that the GROMOS96 force field can predict the skew-boat to chair conformational ratio in good agreement with the experiment, whereas GLYCAM06 shows worse agreement. The omega rotamer distribution about the C5–C6 bond was predicted by both force fields to have torsions around 10°, 190°, and 360°.

Molecular dynamics simulation of the phosphorylation-induced conformational changes of a tau peptide fragment

Aggregation of the microtubule associated protein tau (MAPT) within neurons of the brain is the leading cause of tauopathies such as Alzheimer's disease. MAPT is a phospho-protein that is selectively phosphorylated by a number of kinases in vivo to perform its biological function. However, it may become pathogenically hyperphosphorylated, causing aggregation into paired helical filaments and neurofibrillary tangles. The phosphorylation induced conformational change on a peptide of MAPT (htau225−250) was investigated by performing molecular dynamics simulations with different phosphorylation patterns of the peptide (pThr231 and/or pSer235) in different simulation conditions to determine the effect of ionic strength and phosphate charge. All phosphorylation patterns were found to disrupt a nascent terminal β-sheet pattern (226VAVVR230 and 244QTAPVP249), replacing it with a range of structures. The double pThr231/pSer235 phosphorylation pattern at experimental ionic strength resulted in the best agreement with NMR structural characterization, with the observation of a transient α-helix (239AKSRLQT245). PPII helical conformations were only found sporadically throughout the simulations. Proteins 2014; 82:1907–1923. © 2014 Wiley Periodicals, Inc.

Ab initio molecular dynamics simulations of excited hydrogen halides and methyl halides

Using excited-state ab initio molecular dynamics simulations employing the complete-active-space self-consistent-field approach, we study the mechanism of photodissociation in terms of time evolution of structure, kinetic energy, charges and potential energy for the first excited state of hydrogen halides and methyl halides. Although the hydrogen halides and methyl halides are similar in the photodissociation mechanism, their dynamics are slightly different. The presence of the methyl group causes delay in photodissociation as compared to hydrogen halides.

A comparative molecular dynamics study of diffusion of n-decane and 3-methyl pentane in Y zeolite

Molecular dynamics simulations are reported on the structure and dynamics of n-decane and 3-methylpentane in zeolite NaY. We have calculated several properties such as the center of mass-center of mass rdf, the end-end distance distribution, bond angle distribution and dihedral angle distribution. We have also analysed trajectory to obtain diffusivity and velocity autocorrelation function (VACF). Surprisingly, the diffusivity of 3-methylpentane which is having larger cross-section perpendicular to the long molecular axis is higher than n-decane at 300 K. Activation energies have been obtained from simulations performed at 200 K, 300 K, 350 K, 400 K and 450 K in the NVE ensemble. These results can be understood in terms of the previously known levitation effect. Arrhenious plot has higher value of slope for n-decane (5 center dot 9 kJ/mol) than 3-methylpentane (3 center dot 7 kJ/mol) in agreement with the prediction of levitation effect.

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.

Conformational study of valinomycin: a molecular dynamics approach

Valinomycin is a highly flexible cyclic dodecadepsipeptide that transports ions across membranes. Such a flexibility in the conformation is required for its biological function since it has to encounter a variety of environments and liganding state. Exploration of conformational space of this molecule is therefore important and is one of the objectives of the present study that has been carried out by means of high temperature Molecular Dynamics. Further, the stability of the known bracelet-like structure of the uncomplexed valinomycin and the inherent flexibility around this structure has been investigated. The uncomplexed form of valinomycin has been simulated at 75–100 K for 1 ns in order to elucidate the average conformational properties. An alanine-analog of valinomycin has been simulated under identical conditions in order to evaluate the effect of sidechain on the conformational properties, The studies confirm the effect of sidechain on conformational equilibrium.

Insight into the early stages of thermal unfolding of peanut agglutinin by molecular dynamics simulations

Peanut agglutinin is a homotetrameric nonglycosylated protein. The protein has a unique open quaternary structure. Molecular dynamics simulations have been employed follow the atomistic details of its unfolding at different temperatures. The early events of the deoligomerization of the protein have been elucidated in the present study. Simulation trajectories of the monomer as well as those of the tetramer have been compared and the tetramer is found to be substantially more stable than its monomeric counterpart. The tetramer shows retention of most of its.. secondary structure but considerable loss of the tertiary structure at high temperature. e generation of a This observation impies the molten globule-like intermediate in the later stages of deoligomerization. The quaternary structure of the protein has weakened to a large extent, but none of the subunits are separated. In addition, the importance of the metal-binding to the stability of the protein structure has also been investigated. Binding of the metal ions not only enhances the local stability of the metal-ion binding loop, but also imparts a global stability to the overall structure. The dynamics of different interfaces vary significantly as probed through interface clusters. The differences are substantially enhanced at higher temperatures. The dynamics and the stability of the interfaces have been captured mainly by cluster analysis, which has provided detailed information on the thermal deoligomerization of the protein.

Melting of an Anchored Bilayer: Molecular Dynamics Simulations of the Structural Transition in (CnH2n+1NH3)(2)PbI4 (n=12, 14, 16, 18)

The thermally driven Structural phase transition in the organic-inorganic hybrid perovskite (CnH2n+1NH3)(2)PbI4 has been investigated using molecular dynamics (MD) simulations. This system consists of positively charged alkyl-amine chains anchored to a rigid negatively charged PbI4 sheet with the chains organized as bilayers with a herringbone arrangement. Atomistic simulations were performed using ail isothermal-isobaric ensemble over a wide temperature range from 65 to 665 K for different alkyl chain lengths, n = 12, 14, 16, and 18. The simulations are able to reproduce the essential Features of the experimental observations of this system, including the existence of a transition, the linear variation of the transition temperature with alkyl chain length, and the expansion of the bilayer thickness at the transition. By use of the distance fluctuation Criteria, it is Shown that the transition is associated With a Melting of the alkyl chains of the anchored bilayer. Ail analysis of the conformation of the alkyl chains shows increased disorder in the form of gauche defects above due melting transition. Simulations also show that the melting transition is characterized by the complete disappearance of all-trans alkyl chains in the anchored bilayer, in agreement with experimental observations. A conformationally disordered chain has a larger effective cross-sectional area, and above due transition a uniformly tilted arrangement of the anchored chains call no longer be Sustained. At the melt the angular distribution of the orientation of the chains are 110 longer uniform; the chains are splayed allowing for increased space for individual chains of the anchored bilayer. This is reflected in a sharp rise in the ratio of the mean head-to-head to tail-to-tail distance of the chains of the bilayer at the transition resulting in in expansion of the bilayer thickness. The present MD simulations provide a simple explanation as to how changes in conformation of individual alkyl-chains gives rise to the observed increase in the interlayer lattice spacing of (CnH2n+1NH3)(2)PbI4 at the melting transition.

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.