Hybrid Quantum Mechanical/Molecular Mechanical Methods For Studying Energy Transduction in Biomolecular Machines.

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods have become indispensable tools for the study of biomolecules. In this article, we briefly review the basic methodological details of QM/MM approaches and discuss their applications to various energy transduction problems in biomolecular machines, such as long-range proton transports, fast electron transfers, and mechanochemical coupling. We highlight the particular importance for these applications of balancing computational efficiency and accuracy. Using several recent examples, we illustrate the value and limitations of QM/MM methodologies for both ground and excited states, as well as strategies for calibrating them in specific applications. We conclude with brief comments on several areas that can benefit from further efforts to make QM/MM analyses more quantitative and applicable to increasingly complex biological problems.

DFTB+, a software package for efficient approximate density functional theory based atomistic simulations.

DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green's functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives.

Conversion of a light-driven proton pump into a light-gated ion channel.

Interest in microbial rhodopsins with ion pumping activity has been revitalized in the context of optogenetics, where light-driven ion pumps are used for cell hyperpolarization and voltage sensing. We identified an opsin-encoding gene (CsR) in the genome of the arctic alga Coccomyxa subellipsoidea C-169 that can produce large photocurrents in Xenopus oocytes. We used this property to analyze the function of individual residues in proton pumping. Modification of the highly conserved proton shuttling residue R83 or its interaction partner Y57 strongly reduced pumping power. Moreover, this mutation converted CsR at moderate electrochemical load into an operational proton channel with inward or outward rectification depending on the amino acid substitution. Together with molecular dynamics simulations, these data demonstrate that CsR-R83 and its interacting partner Y57 in conjunction with water molecules forms a proton shuttle that blocks passive proton flux during the dark-state but promotes proton movement uphill upon illumination.

Playing Tic-Tac-Toe with a Sugar-Based Molecular Computer.

Today, molecules can perform Boolean operations and circuits at a level of higher complexity. However, concatenation of logic gates and inhomogeneous inputs and outputs are still challenging tasks. Novel approaches for logic gate integration are possible when chemical programming and software programming are combined. Here it is shown that a molecular finite automaton based on the concatenated implication function (IMP) of a fluorescent two-component sugar probe via a wiring algorithm is able to play tic-tac-toe.

Creatine in mouse models of neurodegeneration and aging.

The supplementation of creatine has shown a marked neuroprotective effect in mouse models of neurodegenerative diseases (Parkinson's disease, Huntington's disease, amyotrophic lateral sclerosis). This has been assigned to the known bioenergetic, anti-apoptotic, anti-excitotoxic and anti-oxidant properties of creatine. As aging and neurodegeneration share pathophysiological pathways, we investigated the effect of oral creatine supplementation on aging in 162 aged wild-type C57Bl/6J mice. The median healthy life span of creatine-fed mice was 9% higher than in their control littermates, and they performed significantly better in neurobehavioral tests. In brains of creatine-treated mice, there was a trend toward a reduction of reactive oxygen species and significantly lower accumulation of the "aging pigment" lipofuscin. Expression profiling showed an upregulation of genes implicated in neuronal growth, neuroprotection, and learning. These data showed that creatine improves health and longevity in mice. Creatine may, therefore, be a promising food supplement to promote healthy human aging. However, the strong neuroprotective effects in animal studies of creatine have not been reproduced in human clinical trials (that have been conducted in Parkinson's disease, Huntington's disease, and amyotrophic lateral sclerosis). The reasons for this translational gap are discussed. One obvious cause seems to be that all previous human studies may have been underpowered. Large phase III trials over long time periods are currently being conducted for Parkinson's disease and Huntington's disease, and will possibly solve this issue.

Charge transport through biomolecular wires in a solvent: bridging molecular dynamics and model Hamiltonian approaches.

We present a hybrid method based on a combination of classical molecular dynamics simulations, quantum-chemical calculations, and a model Hamiltonian approach to describe charge transport through biomolecular wires with variable lengths in presence of a solvent. The core of our approach consists in a mapping of the biomolecular electronic structure, as obtained from density-functional based tight-binding calculations of molecular structures along molecular dynamics trajectories, onto a low-dimensional model Hamiltonian including the coupling to a dissipative bosonic environment. The latter encodes fluctuation effects arising from the solvent and from the molecular conformational dynamics. We apply this approach to the case of pG-pC and pA-pT DNA oligomers as paradigmatic cases and show that the DNA conformational fluctuations are essential in determining and supporting charge transport.

The mitochondrial proteome database: MitoP2.

Defining the mitochondrial proteome is a prerequisite for fully understanding the organelles function as well as mechanisms underlying mitochondrial pathology. The core functions of mitochondria include oxidative phosphorylation, amino acid metabolism, fatty acid oxidation, and ion homeostasis. In addition to these well-known functions, many crucial properties in cell signaling, cell differentiation and cell death are only now being elucidated, and with them the proteins involved. With the wealth of information arriving from single protein studies and sophisticated genome-wide approaches, MitoP2 was designed and is maintained to consolidate knowledge on mitochondrial proteins in one comprehensive database, thus making all pertinent data readily accessible (http://www.mitop2.de). Although the identification of the human mitochondrial proteome is ultimately the prime objective, integration of other species includes Saccharomyces cerevisiae, mouse, Arabidopsis thaliana, and Neurospora crassa so orthology between these species can be interrogated. Data from genome-wide studies can be individually retrieved and are also processed by a support vector machine (SVM) to generate a score that indicates the likelihood of a candidate protein having a mitochondrial location. Manually validated proteins constitute the reference set of the database that contains over 590 yeast, 920 human, and 1020 mouse entries, and that is used for benchmarking the SVM score. Multiple search options allow for the interrogation of the reference set, candidates, disease related proteins, chromosome locations as well as availability of mouse models. Taken together, MitoP2 is a valuable tool for basic scientists, geneticists, and clinicians who are investigating mitochondrial physiology and dysfunction.

Reversible cerebral vasoconstriction syndrome: a complicated clinical course treated with intra-arterial application of nimodipine.

Thunderclap headache (TCH) is a neurological emergency that warrants immediate and comprehensive diagnostic determination. When no pathology can be identified the condition is classified as primary TCH, which is considered benign and self-limiting. TCH has also been reported as the initial symptom of reversible cerebral vasoconstriction syndrome (RCVS), which subsumes a variety of conditions, inconsistently coined Call-Flemming syndrome, benign angiopathy of the central nervous system, drug-induced arteritis, or migrainous vasospasm. Serious complications such as borderline ischaemic stroke have been reported. Although no standardized treatment regime exists, one commonly described but unproven therapy is parenteral or oral application of the calcium channel blocker nimodipine. Here, we report on a case of RCVS, where a progressive course prompted intra-arterial application of nimodipine, which resolved vasoconstriction immediately. We discuss the use of intra-arterial nimodipine application as a potential emergency treatment for a complicated or treatment-refractory course of RCVS.

Mitochondrial 12S rRNA susceptibility mutations in aminoglycoside-associated and idiopathic bilateral vestibulopathy.

The mitochondrial 12S rRNA is considered a hotspot for mutations associated with nonsyndromic (NSHL) and aminoglycoside-induced hearing loss (AIHL). Although aminoglycoside ototoxicity is the most common cause of bilateral vestibular dysfunction, the conceivable role of 12S rRNA mutations has never been systematically investigated. We sequenced the 12S rRNA of 66 patients with bilateral vestibulopathy (BV) with (n=15) or without (n=51) prior exposure to aminoglycosides, as well as 155 healthy controls with intact vestibular function (sport pilots), and compared these to 2704 published sequences (Human Mitochondrial Genome Database). No mutations with a confirmed pathogenicity were found (A1555G, C1494T), but four mutations with a hitherto tentative status were detected (T669C, C960del, C960ins, T961G). Due to their predominant occurrence in patients without aminoglycoside exposure, their detection in controls and a weak evolutionary conservation, their pathogenic role in vestibulocochlear dysfunction remains provisional.

[Migraine recurrence due to intracranial metastasis of a thyroid carcinoma]. / Migränerezidiv durch intrakranielle Metastasierung eines Schilddrüsenkarzinoms.

We report a 57-year-old female with a history of migraine without aura in her early adulthood who complained about new migraine attacks after being free of them for 30 years. As a possible trigger, an intracranial metastasis of a thyroid cancer was found which also caused elevated serum prolactin. The mechanism of a para- or endocrinal effect of the tumour is discussed, showing the relevance of intracranial tumours as a human headache model. The recurrence of a primary headache syndrome after long latency should result in the exclusion of a pathological cause.

SCC-DFTB: what is the proper degree of self-consistency?

The approximate SCC-DFTB method (Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, Th.; Suhai, S.; Seifert, G. Phys. Rev. B 1998, 58, 7260) is derived from DFT by a second-order expansion of the total energy expression. In this article, basic approximations and assumptions underlying the DFTB method are discussed in detail, and further extensions to include third-order terms are proposed. Further, the SCC-DFTB and semiempirical NDDO formalisms are compared to elucidate similarities and differences.

Individual dopaminergic neurons show raised iron levels in Parkinson disease.

OBJECTIVE: Evidence suggests that abnormal iron metabolism is associated with Parkinson disease (PD), with raised iron levels found in pathologically affected areas in PD. It is unknown if this elevated iron is actually associated with neurons or reactive glia, and we therefore addressed this issue by determining if raised iron was present in single dopaminergic neurons. METHODS: We used unfixed frozen sections from postmortem tissue of PD patients and elderly normal individuals to avoid metal contamination and translocation. Levels of iron and other elements were measured using sensitive and specific wavelength dispersive electron probe x-ray microanalysis coupled with cathodoluminescence spectroscopy in individual substantia nigra dopaminergic neurons. RESULTS: We identified raised intraneuronal iron in single defined substantia nigra neurons in PD (mean neuronal iron 2,838 vs 1,611, p < 0.0001) but not in other movement disorders such as Huntington disease. These findings were unrelated to the density of remaining neurons. CONCLUSIONS: Primary changes in neuronal iron could lead to neurodegeneration in Parkinson disease.

Creatine supplementation in Parkinson disease: a placebo-controlled randomized pilot trial.

Mitochondrial dysfunction plays a major role in the pathogenesis of Parkinson disease (PD). Creatine (Cr) is an ergogenic compound that exerts neuroprotective effects in animal models of PD. We conducted a 2-year placebo-controlled randomized clinical trial on the effect of Cr in 60 patients with PD. Cr improved patient mood and led to a smaller dose increase of dopaminergic therapy but had no effect on overall Unified Parkinson's Disease Rating Scale scores or dopamine transporter SPECT.

Toward theoretical analysis of long-range proton transfer kinetics in biomolecular pumps.

Motivated by the long-term goal of theoretically analyzing long-range proton transfer (PT) kinetics in biomolecular pumps, researchers made a number of technical developments in the framework of quantum mechanics-molecular mechanics (QM/MM) simulations. A set of collective reaction coordinates is proposed for characterizing the progress of long-range proton transfers; unlike previous suggestions, the new coordinates can describe PT along highly nonlinear three-dimensional pathways. Calculations using a realistic model of carbonic anhydrase demonstrated that adiabatic mapping using these collective coordinates gives reliable energetics and critical geometrical parameters as compared to minimum energy path calculations, which suggests that the new coordinates can be effectively used as reaction coordinate in potential of mean force calculations for long-range PT in complex systems. In addition, the generalized solvent boundary potential was implemented in the QM/MM framework for rectangular geometries, which is useful for studying reactions in membrane systems. The resulting protocol was found to produce water structure in the interior of aquaporin consistent with previous studies including a much larger number of explicit solvent and lipid molecules. The effect of electrostatics for PT through a membrane protein was also illustrated with a simple model channel embedded in different dielectric continuum environments. The encouraging results observed so far suggest that robust theoretical analysis of long-range PT kinetics in biomolecular pumps can soon be realized in a QM/MM framework.

Calculating absorption shifts for retinal proteins: computational challenges.

Rhodopsins can modulate the optical properties of their chromophores over a wide range of wavelengths. The mechanism for this spectral tuning is based on the response of the retinal chromophore to external stress and the interaction with the charged, polar, and polarizable amino acids of the protein environment and is connected to its large change in dipole moment upon excitation, its large electronic polarizability, and its structural flexibility. In this work, we investigate the accuracy of computational approaches for modeling changes in absorption energies with respect to changes in geometry and applied external electric fields. We illustrate the high sensitivity of absorption energies on the ground-state structure of retinal, which varies significantly with the computational method used for geometry optimization. The response to external fields, in particular to point charges which model the protein environment in combined quantum mechanical/molecular mechanical (QM/MM) applications, is a crucial feature, which is not properly represented by previously used methods, such as time-dependent density functional theory (TDDFT), complete active space self-consistent field (CASSCF), and Hartree-Fock (HF) or semiempirical configuration interaction singles (CIS). This is discussed in detail for bacteriorhodopsin (bR), a protein which blue-shifts retinal gas-phase excitation energy by about 0.5 eV. As a result of this study, we propose a procedure which combines structure optimization or molecular dynamics simulation using DFT methods with a semiempirical or ab initio multireference configuration interaction treatment of the excitation energies. Using a conventional QM/MM point charge representation of the protein environment, we obtain an absorption energy for bR of 2.34 eV. This result is already close to the experimental value of 2.18 eV, even without considering the effects of protein polarization, differential dispersion, and conformational sampling.

Simulation of water cluster assembly on a graphite surface.

The assembly of small water clusters (H2O)n, n = 1-6, on a graphite surface is studied using a density functional tight-binding method complemented with an empirical van der Waals force correction, with confirmation using second-order Møller-Plesset perturbation theory. It is shown that the optimized geometry of the water hexamer may change its original structure to an isoenergy one when interacting with a graphite surface in some specific orientation, while the smaller water cluster will maintain its cyclic or linear configurations (for the water dimer). The binding energy of water clusters interacting with graphite is dependent on the number of water molecules that form hydrogen bonds, but is independent of the water cluster size. These physically adsorbed water clusters show little change in their IR peak position and leave an almost perfect graphite surface.

A global investigation of excited state surfaces within time-dependent density-functional response theory.

This work investigates the capability of time-dependent density functional response theory to describe excited state potential energy surfaces of conjugated organic molecules. Applications to linear polyenes, aromatic systems, and the protonated Schiff base of retinal demonstrate the scope of currently used exchange-correlation functionals as local, adiabatic approximations to time-dependent Kohn-Sham theory. The results are compared to experimental and ab initio data of various kinds to attain a critical analysis of common problems concerning charge transfer and long range (nondynamic) correlation effects. This analysis goes beyond a local investigation of electronic properties and incorporates a global view of the excited state potential energy surfaces.

Quantum mechanics simulation of protein dynamics on long timescale.

Protein structure and dynamics are the keys to a wide range of problems in biology. In principle, both can be fully understood by using quantum mechanics as the ultimate tool to unveil the molecular interactions involved. Indeed, quantum mechanics of atoms and molecules have come to play a central role in chemistry and physics. In practice, however, direct application of quantum mechanics to protein systems has been prohibited by the large molecular size of proteins. As a consequence, there is no general quantum mechanical treatment that not only exceeds the accuracy of state-of-the-art empirical models for proteins but also maintains the efficiency needed for extensive sampling in the conformational space, a requirement mandated by the complexity of protein systems. Here we show that, given recent developments in methods, a general quantum mechanical-based treatment can be constructed. We report a molecular dynamics simulation of a protein, crambin, in solution for 350 ps in which we combine a semiempirical quantum-mechanical description of the entire protein with a description of the surrounding solvent, and solvent-protein interactions based on a molecular mechanics force field. Comparison with a recent very high-resolution crystal structure of crambin (Jelsch et al., Proc Natl Acad Sci USA 2000;102:2246-2251) shows that geometrical detail is better reproduced in this simulation than when several alternate molecular mechanics force fields are used to describe the entire system of protein and solvent, even though the structure is no less flexible. Individual atomic charges deviate in both directions from "canonical" values, and some charge transfer is found between the N and C-termini. The capability of simulating protein dynamics on and beyond the few hundred ps timescale with a demonstrably accurate quantum mechanical model will bring new opportunities to extend our understanding of a range of basic processes in biology such as molecular recognition and enzyme catalysis.

Tetragonal crystalline carbon nitrides: theoretical predictions.

New tetragonal phases of crystalline carbon nitride (CN) and their atomic structures have been identified using a self-consistent-charge density-functional tight-binding method. A tetragonal rocksalt structure provides theoretical support to recent experimental evidence for a stoichiometric CN phase with tetragonal symmetry. A body-centered tetragonal CN phase with 1:1 stoichiometry is predicted to be highly stable and of interesting atomic structure with complicated C-C and N-N dimerizations along the c axis. The cubic-to-tetragonal transitions are carefully examined to understand the underlying mechanism.