Studies of ice melting using molecular dynamics

In this paper, the melting of ice 1h is studied using molecular dynamics (MD). Common potential functions employed in the MD simulations include SPC/E, TIP4P, TIP5P, TIP4/ice and TIP5P/E. We first conducted melting of ice bulks and then studied the melting speed of the ice/water interface during ice melting. It is found that various potentials result in similar ice-melting phenomena. The result is compared with the analytical solution for phase change problem. We also studied size effects and temperature gradient effects on ice melting.     

A Water‐Soluble Cu Complex as Molecular Catalyst for Electrocatalytic CO2 Reduction on Graphene‐Based Electrodes

A structurally simple molecular 1,10‐phenanthroline‐Cu complex on a mesostructured graphene matrix that can be active and selective toward CO₂ reduction over H₂ evolution in an aqueous solution is reported. The active sites consist of Cu(I) center in a distorted trigonal bipyramidal geometry, which enables the adsorption of CO₂ with η¹‐COO‐like configuration to commence the catalysis, with a turnover frequency of ≈45 s^?1 at ?1 V versus reversible hydrogen electrode. Using in situ infrared spectroelectrochemical investigation, it is demonstrated that the Cu complex can be reversibly heterogenized near the graphene surface via potential control. An increase of electron density in the complex is observed as a result of the interaction from the electric field, which further tunes the electron distribution in the neighboring CO₂. It is also found that the mesostructure of graphene matrix favored CO₂ reduction on the Cu center over hydrogen evolution by limiting mass transport from the bulk solution to the electrode surface.     

Multi‐spectroscopic and molecular dynamics simulations investigation of the binding mechanism of polybrominated diphenyl ethers to hen egg white lysozyme

Three PBDEs (BDE25, BDE47, and BDE154) were selected to investigate the interactions between PBDEs and hen egg white lysozyme (HEWL) by molecular modeling, fluorescence spectroscopy, and FT‐IR spectra. The docking results showed that hydrogen bonds were formed between BDE25 and residue TRP63 and between BDE47 and TRP63 with bond lengths of 2.178 Å and 2.146 Å, respectively. The molecular dynamics simulations indicated that van der Waals forces played a predominant role in the binding of three PBDEs to HEWL. The observed fluorescence quenching can be attributed to the formation of complexes between HEWL and PBDEs, and the quenching mechanism is a static quenching. According to Förster's non‐radiative energy transfer theory, the binding distances r were < 7 nm, indicating a high probability of energy transfer from HEWL to the three PBDEs. The synchronous fluorescence showed that the emission maximum wavelength of tryptophan (TRP) residues emerged a red‐shift. FT‐IR spectra indicated that BDE25, BDE47 and BDE154 induced the α‐helix percentage of HEWL decreased from 32.70% ± 1.64% to 28.27% ± 1.41%, 27.50% ± 1.38% and 29.78% ± 1.49%, respectively, whereas the percentage of random coil increased from 26.67% ± 1.33% to 27.60% ± 1.38%, 29.18% ± 1.46% and 30.59% ± 1.53%, respectively.     

Structure-function relationship studies of PTH(1-11) analogues containing sterically hindered dipeptide mimetics.

The N-terminal 1-34 fragment of parathyroid hormone (PTH) is fully active in vitro and in vivo and reproduces all biological responses characteristic of the native intact PTH. In order to develop safer and non-parenteral PTH-like bone anabolic agents, we have studied the effect of introducing conformationally constrained dipeptide mimetics into the N-terminal portion of PTH in an effort to generate miniaturized PTH-mimetics. To this end, we have synthesized and conformationally and biologically characterized PTH(1-11) analogues containing 3R-carboxy-6S-amino-7,5-bicyclic thiazolidinlactam (7,5-bTL), a rigidified dipeptide mimetic unit. The wild type sequence of PTH(1-11) is H-Ser-Val-Ser-Glu-Ile-Gln-Leu-Met-His-Asn-Leu-NH(2). The following pseudo-undecapeptides were prepared: [Ala(1), 7,5-bTL(3, 4), Nle(8), Arg(11)]hPTH(1-11)NH(2) (I); [Ala(1), 7,5-bTL(6, 7), Nle(8), Arg(11)]hPTH(1-11)NH(2) (II); [Ala(1), Nle(8), 7,5-bTL(9, 10), Arg(11)]hPTH(1-11)NH(2) (III). In aqueous solution containing 20% TFE, only analogue I exhibited the typical CD pattern of the alpha-helical conformation. NMR experiments and molecular dynamics calculations located the alpha-helical stretch in the sequence Ile(5)-His(9). The dipeptide mimetic unit 7,5-bTL induces a type III beta-turn, occupying the positions i - 1 and i of the turn. Analogue II exhibited an equilibrium between a type I beta-turn and an alpha-helix, and analogue III did not show any ordered structure. Biological tests revealed poor activity for all analogues (EC(50) > 0.1 mM). Apparently, the relative side-chain orientation of Val(2), Ile(5) and Met(8) can be critical for effective analogue-receptor interaction. Considering helicity as an essential property to obtain active PTH agonists, one must decorate the correctly positioned dipeptide mimetic azabicycloalkane scaffold with substitutions corresponding to the displaced amino acids.     

Detection and identification of ‘Candidatus Phytoplasma asteris’ in Brassica rapa subsp. pekinensis in Poland

Severe growth abnormalities, including leaf yellowing, sprout proliferation and flower virescence and phyllody, were found on Brassica rapa subsp. pekinensis plants in Poland. The presence of phytoplasma in naturally infected plants was demonstrated by polymerase chain reaction assay employing phytoplasma universal P1/P7 followed by R16F2n/R16R2 primer pairs. The detected phytoplasma was identified using restriction fragment length polymorphism analysis (RFLP) of the 16S rRNA gene fragment with AluI, HhaI, MseI and RsaI endonucleases. After enzymatic digestion, all tested samples showed restriction pattern similar to that of ‘Candidatus phytoplasma asteris’. Nested PCR‐amplified products, obtained with primers R16F2n/R16R2, were sequenced. Sequences of the 16S rDNA gene fragment of analysed phytoplasma isolates were nearly identical. They revealed high nucleotide sequence identity (>98%) with corresponding sequences of other phytoplasma isolates from subgroup 16SrI‐B, and they were classified as members of ‘Candidatus phytoplasma asteris’. This is the first report of the natural occurrence of phytoplasma‐associated disease in plants of Chinese cabbage.     

Psb30 contributes to structurally stabilise the Photosystem II complex in the thermophilic cyanobacterium Thermosynechococcus elongatus

A deletion mutant that lacks the Psb30 protein, one of the small subunits of Photosystem II, was constructed in a Thermosynechococcus elongatus strain in which the D1 protein is expressed from the psbA₃ gene (WT*). The ΔPsb30 mutant appears more susceptible to photodamage, has a cytochrome b₅₅₉ that is converted into the low potential form, and probably also lacks the PsbY subunit. In the presence of an inhibitor of protein synthesis, the ?Psb30 lost more rapidly the water oxidation function than the WT* under the high light conditions. These results suggest that Psb30 contributes to structurally and functionally stabilise the Photosystem II complex in preventing the conversion of cytochrome b₅₅₉ into the low potential form. Structural reasons for such effects are discussed.     

Solution NMR Studies of Chlorella Virus DNA Ligase-adenylate

DNA ligases are essential guardians of genome integrity by virtue of their ability to recognize and seal 3′-OH/5′-phosphate nicks in duplex DNA. The substrate binding and three chemical steps of the ligation pathway are coupled to global and local changes in ligase structure, involving both massive protein domain movements and subtle remodeling of atomic contacts in the active site. Here we applied solution NMR spectroscopy to study the conformational dynamics of the Chlorella virus DNA ligase (ChVLig), a minimized eukaryal ATP-dependent ligase consisting of nucleotidyltransferase, OB, and latch domains. Our analysis of backbone ¹⁵N spin relaxation and ¹⁵N,¹H residual dipolar couplings of the covalent ChVLig-AMP intermediate revealed conformational sampling on fast (picosecond to nanosecond) and slow timescales (microsecond to millisecond), indicative of interdomain and intradomain flexibility. We identified local and global changes in ChVLig-AMP structure and dynamics induced by phosphate. In particular, the chemical shift perturbations elicited by phosphate were clustered in the peptide motifs that comprise the active site. We hypothesize that phosphate anion mimics some of the conformational transitions that occur when ligase-adenylate interacts with the nick 5′-phosphate.     

Loss of intramolecular electrostatic interactions and limited conformational ensemble may promote self‐association of cis–tau peptide

Self‐association of proteins can be triggered by a change in the distribution of the conformational ensemble. Posttranslational modification, such as phosphorylation, can induce a shift in the ensemble of conformations. In the brain of Alzheimer's disease patients, the formation of intra‐cellular neurofibrillary tangles deposition is a result of self‐aggregation of hyper‐phosphorylated tau protein. Biochemical and NMR studies suggest that the cis peptidyl prolyl conformation of a phosphorylated threonine‐proline motif in the tau protein renders tau more prone to aggregation than the trans isomer. However, little is known about the role of peptidyl prolyl cis/trans isomerization in tau aggregation. Here, we show that intra‐molecular electrostatic interactions are better formed in the trans isomer. We explore the conformational landscape of the tau segment containing the phosphorylated‐Thr²³¹‐Pro²³² motif using accelerated molecular dynamics and show that intra‐molecular electrostatic interactions are coupled to the isomeric state of the peptidyl prolyl bond. Our results suggest that the loss of intra‐molecular interactions and the more restricted conformational ensemble of the cis isomer could favor self‐aggregation. The results are consistent with experiments, providing valuable complementary atomistic insights and a hypothetical model for isomer specific aggregation of the tau protein. Proteins 2015; 83:436–444. © 2014 Wiley Periodicals, Inc.     

Zipping and Unzipping of Adenylate Kinase: Atomistic Insights into the Ensemble of Open ↔ Closed Transitions

Adenylate kinase (AdK), a phosphotransferase enzyme, plays an important role in cellular energy homeostasis. It undergoes a large conformational change between an open and a closed state, even in the absence of substrate. We investigate the apo-AdK transition at the atomic level both with free-energy calculations and with our new dynamic importance sampling (DIMS) molecular dynamics method. DIMS is shown to sample biologically relevant conformations as verified by comparing an ensemble of hundreds of DIMS transitions to AdK crystal structure intermediates. The simulations reveal in atomic detail how hinge regions partially and intermittently unfold during the transition. Conserved salt bridges are seen to have important structural and dynamic roles; in particular, four ionic bonds that open in a sequential, zipper-like fashion and, thus, dominate the free-energy landscape of the transition are identified. Transitions between the closed and open conformations only have to overcome moderate free-energy barriers. Unexpectedly, the closed state and the open state encompass broad free-energy basins that contain conformations differing in domain hinge motions by up to 40°. The significance of these extended states is discussed in relation to recent experimental Förster resonance energy transfer measurements. Taken together, these results demonstrate how a small number of cooperative key interactions can shape the overall dynamics of an enzyme and suggest an “all-or-nothing” mechanism for the opening and closing of AdK. Our efficient DIMS molecular dynamics computer simulation approach can provide a detailed picture of a functionally important macromolecular transition and thus help to interpret and suggest experiments to probe the conformational landscape of dynamic proteins such as AdK.