Emergence of hidden phases of methylammonium lead iodide (CH3NH3PbI3) upon compression

We perform a thorough structural search with the minima hopping method (MHM) to explore low-energy structures of methylammonium lead iodide. By combining the MHM with a forcefield, we efficiently screen vast portions of the configurational space with large simulation cells containing up to 96 atoms. Our search reveals two structures of methylammonium iodide perovskite (MAPI) that are substantially lower in energy than the well-studied experimentally observed low-temperature Pnma orthorhombic phase according to density functional calculations. Both structures have not yet been reported in the literature for MAPI, but our results show that they could emerge as thermodynamically stable phases via compression at low temperatures. In terms of the electronic properties, the two phases exhibit larger band gaps than the standard perovskite-type structures. Hence, the pressure-induced phase selection at technologically achievable pressures (i.e., via thin-film strain) is a viable route towards the synthesis of several MAPI polymorph with variable band gaps

Accuracy of Molecular Simulation-Based Predictions of koff Values: A Metadynamics Study

The koff values of ligands unbinding to proteins are key parameters for drug discovery. Their predictions based on molecular simulation may under- or overestimate experiment in a system- and/or technique-dependent way. Here we use an established method-infrequent metadynamics, based on the AMBER force field-to compute the koff of the ligand iperoxo (in clinical use) targeting the muscarinic receptor M2. The ligand charges are calculated by either (i) the Amber standard procedure or (ii) B3LYP-DFT. The calculations using (i) turn out not to provide a reasonable estimation of the transition-state free energy. Those using (ii) differ from experiment by 2 orders of magnitude. On the basis of B3LYP DFT QM/MM simulations, we suggest that the observed discrepancy in (ii) arises, at least in part, from the lack of electronic polarization and/or charge transfer in biomolecular force fields. These issues might be present in other systems, such as DNA-protein complexes

Genetic-Algorithm-Based Optimization of a Peptidic Scaffold for Sequestration and Hydration of CO

Biomimicry is a strategy that makes practical use of evolution to find efficient and sustainable ways to produce chemical compounds or engineer products. Exploring the natural machinery of enzymes for the production of desired compounds is a highly profitable investment, but the design of efficient biomimetic systems remains a considerable challenge. An ideal biomimetic system self-assembles in solution, binds a desired range of substrates and catalyzes reactions with turnover rates similar to the native system. To this end, tailoring catalytic functionality in engineered peptides generally requires site-directed mutagenesis or the insertion of additional amino acids, which entails an intensive search across chemical and sequence space. Here we discuss a novel strategy for the computational design of biomimetic compounds and processes that consists of a) characterization of the wild-type and biomimetic systems; b) identification of key descriptors for optimization; c) an efficient search through sequence and chemical space to tailor the catalytic capabilities of the biomimetic system. Through this proof-of-principle study, we are able to decisively understand and identify whether a given scaffold is useful, appropriate and tailorable for a given, desired task

Experimental library screening demonstrates the successful application of computational protein design to large structural ensembles

The stability, activity, and solubility of a protein sequence are determined by a delicate balance of molecular interactions in a variety of conformational states. Even so, most computational protein design methods model sequences in the context of a single native conformation. Simulations that model the native state as an ensemble have been mostly neglected due to the lack of sufficiently powerful optimization algorithms for multistate design. Here, we have applied our multistate design algorithm to study the potential utility of various forms of input structural data for design. To facilitate a more thorough analysis, we developed new methods for the design and high-throughput stability determination of combinatorial mutation libraries based on protein design calculations. The application of these methods to the core design of a small model system produced many variants with improved thermodynamic stability and showed that multistate design methods can be readily applied to large structural ensembles. We found that exhaustive screening of our designed libraries helped to clarify several sources of simulation error that would have otherwise been difficult to ascertain. Interestingly, the lack of correlation between our simulated and experimentally measured stability values shows clearly that a design procedure need not reproduce experimental data exactly to achieve success. This surprising result suggests potentially fruitful directions for the improvement of computational protein design technology

Baffin Island Expedition, 1953: A Preliminary Field Report

Contains description and brief summaries of results of the second Baffin Island Expedition of the Arctic Institute of North America, May-Sept. 1953. The 13-man party carried out a program in the Penny Highland region of Cumberland Peninsula, centering on Pangnirtung Pass. Glaciological work on the Penny Icecap and Highway Glacier included studies of glacier physics, seismic work, and meteorological observations on bedrock geology and structure are given. Corrie formation and the geomorphology of Pangnirtung Pass were studied. Specimens of local fauna were collected and biological studies of some species were made. Comparisons of the vegetation of Penny Highland and East and North Greenland were made through extensive plant collections. Eight peaks were climbed and their altitudes determined. A list of 17 new names approved by the Canadian Board on Geographical Names is appended

Ab initio molecular dynamics using density based energy functionals: application to ground state geometries of some small clusters

The ground state geometries of some small clusters have been obtained via ab initio molecular dynamical simulations by employing density based energy functionals. The approximate kinetic energy functionals that have been employed are the standard Thomas-Fermi $(T_{TF})$ along with the Weizsacker correction $T_W$ and a combination $F(N_e)T_{TF} + T_W$. It is shown that the functional involving $F(N_e)$ gives superior charge densities and bondlengths over the standard functional. Apart from dimers and trimers of Na, Mg, Al, Li, Si, equilibrium geometries for $Li_nAl, n=1,8$ and $Al_{13}$ clusters have also been reported. For all the clusters investigated, the method yields the ground state geometries with the correct symmetries with bondlengths within 5\% when compared with the corresponding results obtained via full orbital based Kohn-Sham method. The method is fast and a promising one to study the ground state geometries of large clusters.Comment: 15 pages, 3 PS figure

Forecasting the Impacts of Silver and Bighead Carp on the Lake Erie Food Web

Nonindigenous bigheaded carps (Bighead Carp Hypophthalmichthys nobilis and Silver Carp H. molitrix; hereafter, “Asian carps” [AC]) threaten to invade and disrupt food webs and fisheries in the Laurentian Great Lakes through their high consumption of plankton. To quantify the potential effects of AC on the food web in Lake Erie, we developed an Ecopath with Ecosim (EwE) food web model and simulated four AC diet composition scenarios (high, low, and no detritus and low detritus with Walleye Sander vitreus and Yellow Perch Perca flavescens larvae) and two nutrient load scenarios (the 1999 baseline load and 2× the baseline [HP]). We quantified the uncertainty of the potential AC effects by coupling the EwE model with estimates of parameter uncertainty in AC production, consumption, and predator diets obtained using structured expert judgment. Our model projected mean ± SD AC equilibrium biomass ranging from 52 ± 34 to 104 ± 75 kg/ha under the different scenarios. Relative to baseline simulations without AC, AC invasion under all detrital diet scenarios decreased the biomass of most fish and zooplankton groups. The effects of AC in the HP scenario were similar to those in the detrital diet scenarios except that the biomasses of most Walleye and Yellow Perch groups were greater under HP because these fishes were buffered from competition with AC by increased productivity at lower trophic levels. Asian carp predation on Walleye and Yellow Perch larvae caused biomass declines among all Walleye and Yellow Perch groups. Large food web impacts of AC occurred in only 2% of the simulations, where AC biomass exceeded 200 kg/ha, resulting in biomass declines of zooplankton and planktivorous fish near the levels observed in the Illinois River. Our findings suggest that AC would affect Lake Erie's food web by competing with other planktivorous fishes and by providing additional prey for piscivores. Our methods provide a novel approach for including uncertainty into forecasts of invasive species' impacts on aquatic food webs. Received December 6, 2014; accepted July 15, 201

Four small puzzles that Rosetta doesn't solve

A complete macromolecule modeling package must be able to solve the simplest structure prediction problems. Despite recent successes in high resolution structure modeling and design, the Rosetta software suite fares poorly on deceptively small protein and RNA puzzles, some as small as four residues. To illustrate these problems, this manuscript presents extensive Rosetta results for four well-defined test cases: the 20-residue mini-protein Trp cage, an even smaller disulfide-stabilized conotoxin, the reactive loop of a serine protease inhibitor, and a UUCG RNA tetraloop. In contrast to previous Rosetta studies, several lines of evidence indicate that conformational sampling is not the major bottleneck in modeling these small systems. Instead, approximations and omissions in the Rosetta all-atom energy function currently preclude discriminating experimentally observed conformations from de novo models at atomic resolution. These molecular "puzzles" should serve as useful model systems for developers wishing to make foundational improvements to this powerful modeling suite.Comment: Published in PLoS One as a manuscript for the RosettaCon 2010 Special Collectio

Ultrafast pulse shaping modulates perceived visual brightness in living animals

Vision is usually assumed to be sensitive to the light intensity and spectrum but not to its spectral phase. However, experiments performed on retinal proteins in solution showed that the first step of vision consists in an ultrafast photoisomerization that can be coherently controlled by shaping the phase of femtosecond laser pulses, especially in the multiphoton interaction regime. The link between these experiments in solution and the biological process allowing vision was not demonstrated. Here, we measure the electric signals fired from the retina of living mice upon femtosecond multipulse and single-pulse light stimulation. Our results show that the electrophysiological signaling is sensitive to the manipulation of the light excitation on a femtosecond time scale. The mechanism relies on multiple interactions with the light pulses close to the conical intersection, like pump-dump (photoisomerization interruption) and pump-repump (reverse isomerization) processes. This interpretation is supported both experimentally and by dynamics simulations