A consistent treatment of electrostatic energies is arguably the most important requirement for the realistic modeling of biological systems. An important part of electrostatic modeling is the ability to account for the polarizability of the simulated system. This can be done both macroscopically and microscopically, but the use of macroscopic models may lead to conceptual traps, which do not exist in the microscopic treatments. The present work describes the development of microscopic polarizable force fields starting with the introduction of these powerful tools and following some of the subsequent developments in the field. Special effort has been made to review a wide range of applications and emphasize cases when the use of polarizable force fields is important. Finally, a brief perspective is given on the future of this rapidly growing field.

The idea that enzymes catalyze reactions by dynamical coupling between the conformational motions and the chemical coordinates has recently attracted major experimental and theoretical interest. However, experimental studies have not directly established that the conformational motions transfer energy to the chemical coordinate, and simulating enzyme catalysis on the relevant timescales has been impractical. Here, we introduce a renormalization approach that transforms the energetics and dynamics of the enzyme to an equivalent low-dimensional system, and allows us to simulate the dynamical coupling on a ms timescale. The simulations establish, by means of several independent approaches, that the conformational dynamics is not remembered during the chemical step and does not contribute significantly to catalysis. Nevertheless, the precise nature of this coupling is a question of great importance.adenylate kinase ͉ coarse-grained model ͉ renormalization ͉ simplified model ͉ enzyme catalysis

The origin of the barrier for proton transport through the aquaporin channel is a problem of general interest. It is becoming increasingly clear that this barrier is not attributable to the orientation of the water molecules across the channel but rather to the electrostatic penalty for moving the proton charge to the center of the channel. However, the reason for the high electrostatic barrier is still rather controversial. It has been argued by some workers that the barrier is due to the so-called NPA motif and/or to the helix macrodipole or to other specific elements. However, our works indicated that the main reason for the high barrier is the loss of the generalized solvation upon moving the proton charge from the bulk to the center of the channel and that this does not reflect a specific repulsive electrostatic interaction but the absence of sufficient electrostatic stabilization. At this stage it seems that the elucidation and clarification of the origin of the electrostatic barrier can serve as an instructive test case for electrostatic models. Thus, we reexamine the free-energy surface for proton transport in aquaporins using the microscopic free-energy perturbation/umbrella sampling (FEP/US) and the empirical valence bond/umbrella sampling (EVB/US) methods as well as the semimacroscopic protein dipole Langevin dipole model in its linear response approximation version (the PDLD/S-LRA). These extensive studies help to clarify the nature of the barrier and to establish the "reduced solvation effect" as the primary source of this barrier. That is, it is found that the barrier is associated with the loss of the generalized solvation energy (which includes of course all electrostatic effects) upon moving the proton charge from the bulk solvent to the center of the channel. It is also demonstrated that the residues in the NPA region and the helix dipole cannot be considered as the main reasons for the electrostatic barrier. Furthermore, our microscopic and semimacroscopic studies clarify the problems with incomplete alternative calculations, illustrating that the effects of various electrostatic elements are drastically overestimated by macroscopic calculations that use a low dielectric constant and do not consider the protein reorganization. Similarly, it is pointed out that microscopic potential of mean force calculations that do not evaluate the electrostatic barrier relative to the bulk water cannot be used to establish the origin of the electrostatic barrier. The relationship between the present study and calculations of pK(a)s in protein interiors is clarified, pointing out that approaches that are applied to study the aquaporin barrier should be validated by pK(a)s calculations. Such calculations also help to clarify the crucial role of solvation energies in establishing the barrier in aquaporins.

The structure of quinol-dependent nitric oxide reductase (qNOR) from G. stearothermophilus, which catalyzes the reduction of NO to produce the major ozone-depleting gas N(2)O, has been characterized at 2.5 Å resolution. The overall fold of qNOR is similar to that of cytochrome c-dependent NOR (cNOR), and some structural features that are characteristic of cNOR, such as the calcium binding site and hydrophilic cytochrome c domain, are observed in qNOR, even though it harbors no heme c. In contrast to cNOR, structure-based mutagenesis and molecular dynamics simulation studies of qNOR suggest that a water channel from the cytoplasm can serve as a proton transfer pathway for the catalytic reaction. Further structural comparison of qNOR with cNOR and aerobic and microaerobic respiratory oxidases elucidates their evolutionary relationship and possible functional conversions.

Summary Paragraph Cytochrome c oxidase is a member of the heme copper oxidase superfamily (HCO)1. HCOs function as the terminal enzymes in the respiratory chain of mitochondria and aerobic prokaryotes, coupling molecular oxygen reduction to transmembrane proton pumping. Integral to the enzyme’s function is the transfer of electrons from cytochrome c to the oxidase via a transient association of the two proteins. Electron entry and exit are proposed to occur from the same site on cytochrome c2–4. Here we report the crystal structure of the caa3-type cytochrome oxidase from Thermus thermophilus, which has a covalently tethered cytochrome c domain. Crystals were grown in a bicontinuous mesophase using a synthetic short-chain monoacylglycerol as the hosting lipid. From the electron density map, at 2.36 Å resolution, a novel integral membrane subunit and a native glycoglycerophospholipid embedded in the complex were identified. Contrary to previous electron transfer mechanisms observed for soluble cytochrome c, the structure reveals the architecture of the electron transfer complex for the fused cupredoxin/cytochrome c domain which implicates different sites on cytochrome c for electron entry and exit. Support for an alternative to the classical proton gate characteristic of this HCO class is presented.

Using the nonperturbative approach to the calculation of nonlinear optical spectra developed in a foregoing paper [Mancal et al., J. Chem. Phys. 124, 234504 (2006), preceding paper], calculations of two-dimensional electronic spectra of an excitonically coupled dimer model system are presented. The dissipative exciton transfer dynamics is treated within the Redfield theory and energetic disorder within the molecular ensemble is taken into account. The manner in which the two-dimensional spectra reveal electronic couplings in the aggregate system and the evolution of the spectra in time is studied in detail. Changes in the intensity and shape of the peaks in the two-dimensional relaxation spectra are related to the coherent and dissipative dynamics of the system. It is shown that coherent electronic motion, an electronic analog of a vibrational wave packet, can manifest itself in two-dimensional optical spectra of molecular aggregate systems as a periodic modulation of both the diagonal and off-diagonal peaks.

Gaining detailed understanding of the energetics of the protonpumping process in cytochrome c oxidase (CcO) is one of the challenges of modern biophysics. Despite promising mechanistic proposals, most works have not related the activation barriers of the different assumed steps to the protein structure, and there has not been a physically consistent model that reproduced the barriers needed to create a working pump. This work reevaluates the activation barriers for the primary proton transfer (PT) steps by calculations that reflect all relevant free energy contributions, including the electrostatic energies of the generated charges, the energies of water insertion, and large structural rearrangements of the donor and acceptor. The calculations have reproduced barriers that account for the directionality and sequence of events in the primary PT in CcO. It has also been found that the PT from Glu-286 (E) to the propionate of heme a3 (Prd) provides a gate for an initial back leakage from the high pH side of the membrane. Interestingly, the rotation of E that brings it closer to Prd appears to provide a way for blocking competing pathways in the primary PT. Our study elucidates and quantifies the nature of the control of the directionality in the primary PT in CcO and provides instructive insight into the role of the water molecules in biological PT, showing that ''bridges'' of several water molecules in hydrophobic regions present a problem (rather than a solution) that is minimized in the primary PT.coupled electron transfer proton transfer ͉ proton pumps ͉ dielectric effects ͉ free energy calculations ͉ pK a in nonpolar regions

Malaria and cryptosporidiosis, caused by apicomplexan parasites, remain major drivers of global child mortality. New drugs for the treatment of malaria and cryptosporidiosis, in particular, are of high priority; however, there are few chemically validated targets. The natural product cladosporin is active against blood- and liver-stagePlasmodium falciparumandCryptosporidium parvumin cell-culture studies. Target deconvolution inP. falciparumhas shown that cladosporin inhibits lysyl-tRNA synthetase (PfKRS1). Here, we report the identification of a series of selective inhibitors of apicomplexan KRSs. Following a biochemical screen, a small-molecule hit was identified and then optimized by using a structure-based approach, supported by structures of bothPfKRS1 andC. parvumKRS (CpKRS). In vivo proof of concept was established in an SCID mouse model of malaria, after oral administration (ED90= 1.5 mg/kg, once a day for 4 d). Furthermore, we successfully identified an opportunity for pathogen hopping based on the structural homology betweenPfKRS1 andCpKRS. This series of compounds inhibitCpKRS andC. parvumandCryptosporidium hominisin culture, and our lead compound shows oral efficacy in two cryptosporidiosis mouse models. X-ray crystallography and molecular dynamics simulations have provided a model to rationalize the selectivity of our compounds forPfKRS1 andCpKRS vs. (human)HsKRS. Our work validates apicomplexan KRSs as promising targets for the development of drugs for malaria and cryptosporidiosis.