There is, I think, something ethereal about i —the square root of minus one. I remember first hearing about it at school. It seemed an odd beast at that time—an intruder hovering on the edge of reality.

Usually familiarity dulls this sense of the bizarre, but in the case of i it was the reverse: over the years the sense of its surreal nature intensified. It seemed that it was impossible to write mathematics that described the real world in …

I and i

J. Appl. Cryst.の発刊に際して

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecu- lar simulation program. It has been developed over the last three decades with a primary focus on molecules of bio- logical interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estima- tors, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numer- ous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.

/pdf/charmm-the-biomolecular-simulation-program-4kexem0xkx.pdf

CHARMM: the biomolecular simulation program.

A frequent problem in computational modeling is the interconversion of chemical structures between different formats. While standard interchange formats exist (for example, Chemical Markup Language) and de facto standards have arisen (for example, SMILES format), the need to interconvert formats is a continuing problem due to the multitude of different application areas for chemistry data, differences in the data stored by different formats (0D versus 3D, for example), and competition between software along with a lack of vendor-neutral formats. We discuss, for the first time, Open Babel, an open-source chemical toolbox that speaks the many languages of chemical data. Open Babel version 2.3 interconverts over 110 formats. The need to represent such a wide variety of chemical and molecular data requires a library that implements a wide range of cheminformatics algorithms, from partial charge assignment and aromaticity detection, to bond order perception and canonicalization. We detail the implementation of Open Babel, describe key advances in the 2.3 release, and outline a variety of uses both in terms of software products and scientific research, including applications far beyond simple format interconversion. Open Babel presents a solution to the proliferation of multiple chemical file formats. In addition, it provides a variety of useful utilities from conformer searching and 2D depiction, to filtering, batch conversion, and substructure and similarity searching. For developers, it can be used as a programming library to handle chemical data in areas such as organic chemistry, drug design, materials science, and computational chemistry. It is freely available under an open-source license from http://openbabel.org
 .

/pdf/open-babel-an-open-chemical-toolbox-2c5goscrz3.pdf

Open Babel: An open chemical toolbox

BIOE 402. Medical Technology Assessment. 2 or 3 hours. Bioentrepreneur course. Assessment of medical technology in the context of commercialization. Objectives, competition, market share, funding, pricing, manufacturing, growth, and intellectual property; many issues unique to biomedical products. Course Information: 2 undergraduate hours. 3 graduate hours. Prerequisite(s): Junior standing or above and consent of the instructor.

“Bioinformatics” 특집을 내면서

The efficient and accurate characterization of solvent effects is a key element in the theoretical and computational study of biological problems. Implicit solvent models, particularly generalized ...

Balancing Solvation and Intramolecular Interactions: Toward a Consistent Generalized Born Force Field

Harvesting energy from the environment offers the promise of clean power for self-sustained systems1,2. Known technologies—such as solar cells, thermoelectric devices and mechanical generators—have specific environmental requirements that restrict where they can be deployed and limit their potential for continuous energy production3–5. The ubiquity of atmospheric moisture offers an alternative. However, existing moisture-based energy-harvesting technologies can produce only intermittent, brief (shorter than 50 seconds) bursts of power in the ambient environment, owing to the lack of a sustained conversion mechanism6–12. Here we show that thin-film devices made from nanometre-scale protein wires harvested from the microbe Geobacter sulfurreducens can generate continuous electric power in the ambient environment. The devices produce a sustained voltage of around 0.5 volts across a 7-micrometre-thick film, with a current density of around 17 microamperes per square centimetre. We find the driving force behind this energy generation to be a self-maintained moisture gradient that forms within the film when the film is exposed to the humidity that is naturally present in air. Connecting several devices linearly scales up the voltage and current to power electronics. Our results demonstrate the feasibility of a continuous energy-harvesting strategy that is less restricted by location or environmental conditions than other sustainable approaches. A new type of energy-harvesting device, based on protein nanowires from the microbe Geobacter sulforreducens, can generate a sustained power output by producing a moisture gradient across the nanowire film using natural humidity.

Power generation from ambient humidity using protein nanowires.

Recent advances in implicit solvent based methods for biomolecular simulations

Accurate description of the solvent environment is critical in computer simulations of protein structure and dynamics. An implicit treatment of solvent aims to capture the mean influence of water molecules on the solute via direct estimation of the solvation free energy. It has emerged as a powerful alternative to explicit solvent, and provides a favorable compromise between computational cost and level of detail. We review the current theory and techniques for implicit modeling of nonpolar solvation in the context of simulating protein folding and conformational transitions, and discuss the main directions for further development. It is demonstrated that the current surface area based nonpolar models have severe limitations, including insufficient description of the conformational dependence of solvation, over-estimation of the strength of pair-wise nonpolar interactions, and incorrect prediction of anti-cooperativity for three-body hydrophobic associations. We argue that, to improve beyond current level of accuracy of implicit solvent models, two important aspects of nonpolar solvation need to be incorporated, namely, the length-scale dependence of hydrophobic association and solvent screening of solute–solute dispersion interactions. We recognize that substantial challenges exist in constructing a sufficiently balanced, yet reasonably efficient, implicit solvent protein force field. Nonetheless, most of the fundamental problems are understood, and exciting progress has been made over the last few years. We believe that continual work along the frontiers outlined will greatly improve one’s ability to study protein folding and large conformational transitions at atomistic detail.

https://sakai.rutgers.edu/access/content/user/kparis/biomaps_513_references/04_D_02_PhysChemChemPhys10_471_nonpolar_solvation.pdf

Implicit modeling of nonpolar solvation for simulating protein folding and conformational transitions.

Recent advances in efficient and accurate treatment of solvent with the generalized Born approximation (GB) have made it possible to substantially refine the protein structures generated by various prediction tools through detailed molecular dynamics simulations. As demonstrated in a recent CASPR experiment, improvement can be quite reliably achieved when the initial models are sufficiently close to the native basin (e.g., 3-4 A C(alpha) RMSD). A key element to effective refinement is to incorporate reliable structural information into the simulation protocol. Without intimate knowledge of the target and prediction protocol used to generate the initial structural models, it can be assumed that the regular secondary structure elements (helices and strands) and overall fold topology are largely correct to start with, such that the protocol limits itself to the scope of refinement and focuses the sampling in vicinity of the initial structure. The secondary structures can be enforced by dihedral restraints and the topology through structural contacts, implemented as either multiple pair-wise C(alpha) distance restraints or a single sidechain distance matrix restraint. The restraints are weakly imposed with flat-bottom potentials to allow sufficient flexibility for structural rearrangement. Refinement is further facilitated by enhanced sampling of advanced techniques such as the replica exchange method (REX). In general, for single domain proteins of small to medium sizes, 3-5 nanoseconds of REX/GB refinement simulations appear to be sufficient for reasonable convergence. Clustering of the resulting structural ensembles can yield refined models over 1.0 A closer to the native structure in C(alpha) RMSD. Substantial improvement of sidechain contacts and rotamer states can also be achieved in most cases. Additional improvement is possible with longer sampling and knowledge of the robust structural features in the initial models for a given prediction protocol. Nevertheless, limitations still exist in sampling as well as force field accuracy, manifested as difficulty in refinement of long and flexible loops.

Jianhan Chen

Papers

Balancing Solvation and Intramolecular Interactions: Toward a Consistent Generalized Born Force Field

Power generation from ambient humidity using protein nanowires.

Recent advances in implicit solvent based methods for biomolecular simulations

Implicit modeling of nonpolar solvation for simulating protein folding and conformational transitions.

Can molecular dynamics simulations provide high-resolution refinement of protein structure?