Optimal use of data in parallel tempering simulations for the construction of discrete-state Markov models of biomolecular dynamics

Parallel tempering (PT) molecular dynamics simulations have been extensively investigated as a means of efficient sampling of the configurations of biomolecular systems. Recent work has demonstrated how the short physical trajectories generated in PT simulations of biomolecules can be used to construct the Markov models describing biomolecular dynamics at each simulated temperature. While this approach describes the temperature-dependent kinetics, it does not make optimal use of all available PT data, instead estimating the rates at a given temperature using only data from that temperature. This can be problematic, as some relevant transitions or states may not be sufficiently sampled at the temperature of interest, but might be readily sampled at nearby temperatures. Further, the comparison of temperature-dependent properties can suffer from the false assumption that data collected from different temperatures are uncorrelated. We propose here a strategy in which, by a simple modification of the PT protocol, the harvested trajectories can be reweighted, permitting data from all temperatures to contribute to the estimated kinetic model. The method reduces the statistical uncertainty in the kinetic model relative to the single temperature approach and provides estimates of transition probabilities even for transitions not observed at the temperature of interest. Further, the method allows the kinetics to be estimated at temperatures other than those at which simulations were run. We illustrate this method by applying it to the generation of a Markov model of the conformational dynamics of the solvated terminally blocked alanine peptide.     

Toward biodegradable nanogel star polymers via organocatalytic ROP

Organocatalytic ring opening polymerization (OROP) is used to effect the rapid, scalable, room temperature formation of size-controlled, highly uniform, polyvalent, nanogel star polymer nanoparticles of biodegradable composition.     

Synthesis of BaZrS3 and BaHfS3 Chalcogenide Perovskite Films Using Single-Phase Molecular Precursors at Moderate Temperatures

Chalcogenide perovskites have garnered interest for applications in semiconductor devices due to their excellent predicted optoelectronic properties and stability. However, high synthesis temperatures have historically made these materials incompatible with the creation of photovoltaic devices. Here, we demonstrate the solution processed synthesis of luminescent BaZrS₃ and BaHfS₃ chalcogenide perovskite films using single-phase molecular precursors at sulfurization temperatures of 575 °C and sulfurization times as short as one hour. These molecular precursor inks were synthesized using known carbon disulfide insertion chemistry to create Group 4 metal dithiocarbamates, and this chemistry was extended to create species, such as barium dithiocarboxylates, that have never been reported before. These findings, with added future research, have the potential to yield fully solution processed thin films of chalcogenide perovskites for various optoelectronic applications.     

Electrochemically modulated liquid chromatography using a boron-doped diamond particle stationary phase

This paper reports on preliminary tests of the performance of boron-doped diamond powder (BDDP) as a stationary phase in electrochemically modulated liquid chromatography (EMLC). EMLC manipulates retention through changes in the potential applied (E(appl)) to a conductive packing. Porous graphitic carbon (PGC) has routinely been utilized as a material in EMLC separations. Herein the utility of BDDP as a stationary phase in EMLC was investigated and its stability, both compositionally and microstructurally, relative to PGC was compared. The results show that BDDP is stable over a wide range of E(appl) values (i.e., -1.2 to +1.2V vs. Ag/AgCl, sat'd NaCl). The data also reveal that electrostatics play a key role in the adsorption of the aromatic sulfonates on the BDDP stationary phase, and that these analytes are more weakly retained in comparison to the PGC support. The potential for this methodology to provide a means to advance the understanding of molecular adsorption and retention mechanisms on carbonaceous materials is briefly discussed.     

Alchemical prediction of hydration free energies for SAMPL

Hydration free energy calculations have become important tests of force fields. Alchemical free energy calculations based on molecular dynamics simulations provide a rigorous way to calculate these free energies for a particular force field, given sufficient sampling. Here, we report results of alchemical hydration free energy calculations for the set of small molecules comprising the 2011 Statistical Assessment of Modeling of Proteins and Ligands challenge. Our calculations are largely based on the Generalized Amber Force Field with several different charge models, and we achieved RMS errors in the 1.4-2.2 kcal/mol range depending on charge model, marginally higher than what we typically observed in previous studies (Mobley et al. in J Phys Chem B 111(9):2242-2254, 2007, J Chem Theory Comput 5(2):350-358, 2009, J Phys Chem B 115:1329-1332, 2011; Nicholls et al. in J Med Chem 51:769-779, 2008; Klimovich and Mobley in J Comput Aided Mol Design 24(4):307-316, 2010). The test set consists of ethane, biphenyl, and a dibenzyl dioxin, as well as a series of chlorinated derivatives of each. We found that, for this set, using high-quality partial charges from MP2/cc-PVTZ SCRF RESP fits provided marginally improved agreement with experiment over using AM1-BCC partial charges as we have more typically done, in keeping with our recent findings (Mobley et al. in J Phys Chem B 115:1329-1332, 2011). Switching to OPLS Lennard-Jones parameters with AM1-BCC charges also improves agreement with experiment. We also find a number of chemical trends within each molecular series which we can explain, but there are also some surprises, including some that are captured by the calculations and some that are not.     

Alternative expressions for energies and forces due to angle bending and torsional energy

We have derived alternative expressions for computing the energies and forces associated with angle bending and torsional energy terms commonly used in molecular mechanics and molecular dynamics computer programs. Our expressions address the problems of singularities that are intrinsic in popular angle energy functions and that occur from other chain rule derivations of force expressions. Most chain rule derivations of expressions for Cartesian forces due to angle energies make use of relations such as where ? is a bond or torsion angle, E(?) is energy, and ?/?x represents a derivative with respect to some Cartesian coordinate. This expression leads to singularities from the middle term, ?1/sin ?, when ? is 0 or π. This is a problem that prevents the use of torsional energy expressions that have phase angles, ?°, other than 0 or π, such as in E(?) = κ[1 + cos(n? ? phsi;°)]. Our derivations make use of a different, but equivalent, form of the chain rule: This form still possesses singularities for the bond angle forces since the last factor is undefined when ? is 0 or π. However, the alternate form may be used to great advantage for the torsional angle forces where no such problem arises. The new expressions are necessary if one desires the use of torsional energy expressions with general phase angles. Even for energy expressions in common use, i.e., with phase angles of 0 or π, our force expressions are as computationally efficient as the standard ones. The new expressions are applicable to all molecular simulations that employ restrained, or phase-shifted, torsional angle energy expressions.     

The melanocortin 1 receptor and the UV response of human melanocytes--a shift in paradigm

Cutaneous pigmentation is the major photoprotective mechanism against the carcinogenic and aging effects of UV. Epidermal melanocytes synthesize the pigment melanin, in the form of eumelanin or pheomelanin. Synthesis of the photoprotective eumelanin by human melanocytes is regulated mainly by the melanocortins alpha-melanocortin (alpha-MSH) and adrenocorticotropic hormone (ACTH), which bind the melanocortin 1 receptor (MC1R) and activate the cAMP pathway that is required for UV-induced tanning. Melanocortins stimulate proliferation and melanogenesis and inhibit UV-induced apoptosis of human melanocytes. Importantly, melanocortins reduce the generation of hydrogen peroxide and enhance repair of DNA photoproducts, independently of pigmentation. MC1R is a major contributor to the diversity of human pigmentation and a melanoma susceptibility gene. Certain allelic variants of this gene, namely R151C, R160W and D294H, are strongly associated with red hair phenotype and increased melanoma susceptibility. Natural expression of two of these variants sensitizes melanocytes to the cytotoxic effect of UV, and increases the burden of DNA damage and oxidative stress. We are designing potent melanocortin analogs that mimic the effects of alpha-MSH as a strategy to prevent skin cancer, particularly in individuals who express MC1R genotypes that reduce but do not abolish MC1R function, or mutations in other melanoma susceptibility genes, such as p16.