The saposins are small, membrane-active proteins that exist in both soluble and lipid-bound states. Saposin A has roles in sphingolipid catabolism and transport and is required for the breakdown of galactosylceramide by β-galactosylceramidase. In the absence of lipid, saposin A adopts a closed monomeric apo conformation typical of this family. To study a lipid-bound state of this protein, we determined the crystal structure of saposin A in the presence of detergent to 1.9 Å resolution. The structure reveals two chains of saposin A in an open conformation encapsulating 40 internally bound detergent molecules organized in a highly ordered bilayer-like hydrophobic core. The complex provides a high-resolution view of a discoidal lipoprotein particle in which all of the internalized acyl chains are resolved. Saposin A lipoprotein discs exhibit limited selectivity with respect to the incorporated lipid, and can solubilize phospholipids, sphingolipids, and cholesterol into discrete, monodisperse particles with mass of approximately 27 kDa. These discs may be the smallest possible lipoprotein structures that are stabilized by lipid self-assembly.

Uracil (U) can be found in DNA as a mismatch paired either to adenine (A) or to guanine (G). Removal of U from DNA is performed by a class of enzymes known as uracilDNAglycosylases (UDG). Recent studies suggest that recognition of UA and UG mismatches by UDG takes place via an extra-helical mechanism. In this work, we use molecular dynamics simulations to analyze the structure, dynamics and open base pair kinetics of UA base pairs relative to their natural TA counterpart in 12 dodecamers. Our results show that the presence of U does not alter the local conformation of B-DNA. Breathing dynamics and base pair closing kinetics are only weakly dependent on the presence of U versus T, with open TA and UA pairs lifetimes in the nanosecond timescale. Additionally, we observed spontaneous base flipping in UA pairs. We analyze the structure and dynamics for this event and compare the results to available crystallographic data of open base pair conformations. Our results are in agreement with both structural and kinetic data derived from NMR imino proton exchange measurements, providing the first detailed description at the molecular level of elusive events such as spontaneous base pair opening and flipping in mismatched UA sequences in DNA. Based on these results, we propose that base pair flipping can occur spontaneously at room temperature via a 3-step mechanism with an open base pair intermediate. Implications for the molecular basis of U recognition by UDG are discussed.

We used molecular dynamics (MD) simulations in conjunction with the Kirkwood-Buff theory to compute the partial molar volumes for a number of small solutes of various chemical natures. We repeated our computations using modified pair potentials, first, in the absence of the Coulombic term and, second, in the absence of the Coulombic and the attractive Lennard-Jones terms. Comparison of our results with experimental data and the volumetric results of Monte Carlo simulation with hard sphere potentials and scaled particle theory-based computations led us to conclude that, for small solutes, the partial molar volume computed with the Lennard-Jones potential in the absence of the Coulombic term nearly coincides with the cavity volume. On the other hand, MD simulations carried out with the pair interaction potentials containing only the repulsive Lennard-Jones term produce unrealistically large partial molar volumes of solutes that are close to their excluded volumes. Our simulation results are in good agreement with the reported schemes for parsing partial molar volume data on small solutes. In particular, our determined interaction volumes() and the thickness of the thermal volume for individual compounds are in good agreement with empirical estimates. This work is the first computational study that supports and lends credence to the practical algorithms of parsing partial molar volume data that are currently in use for molecular interpretations of volumetric data.

The capacity to form -sheet structure and to self-organize into amyloid aggregates is a property shared by many proteins. Severe neurodegenerative pathologies such as Alzheimer's disease are thought to involve the interaction of amyloidogenic protein oligomers with neuronal membranes. To understand the experimentally observed catalysis of amyloid formation by lipid membranes and other water-hydrophobic interfaces, we examine the physico-chemical basis of peptide adsorption and aggregation in a model membrane using atomistic molecular simulations. Blocked octapeptides with simple, repetitive sequences, (Gly-Ala), and (Gly-Val), are used as models of -sheet-forming polypeptide chains found in the core of amyloid fibrils. In the presence of an n-octane phase mimicking the core of lipid membranes, the peptides spontaneously partition at the octane-water interface. The adsorption of nonpolar sidechains displaces the peptides' conformational equilibrium from a heterogeneous ensemble characterized by a high degree of structural disorder toward a more ordered ensemble favoring -hairpins and elongated -strands. At the interface, peptides spontaneously aggregate and rapidly evolve -sheet structure on a 10 to 100 ns time scale, while aqueous aggregates remain amorphous. Catalysis of -sheet formation results from the combination of the hydrophobic effect and of reduced conformational entropy of the polypeptide chain. While the former drives interfacial partition and displaces the conformational equilibrium of monomeric peptides, the planar interface further facilitates -sheet organization by increasing peptide concentration and reducing the dimensionality of self-assembly from three to two. These findings suggest a general mechanism for the formation of -sheets on the surface of globular proteins and for amyloid self-organization at hydrophobic interfaces.

Cytochrome c oxidase, the terminal enzyme of the respiratory chain, utilizes the reduction of dioxygen into water to pump protons across the mitochondrial inner membrane. The principal pathway of proton uptake into the enzyme, the D channel, is a 2.5 nm long channel-like cavity named after a conserved, negatively charged aspartic acid (D) residue thought to help recruiting protons to its entrance (D132 in the first subunit of the S. sphaeroides enzyme). The single-point mutation of D132 to asparagine (N), a neutral residue, abolishes enzyme activity. Conversely, replacing conserved N139, one-third into the D channel, by D, induces a decoupled phenotype, whereby oxygen reduction proceeds but not proton pumping. Intriguingly, the double mutant D132N/N139D, which conserves the charge of the D channel, restores the wild-type phenotype. We use molecular dynamics simulations and electrostatic calculations to examine the structural and physical basis for the coupling of proton pumping and oxygen chemistry in single and double N139D mutants. The potential of mean force for the conformational isomerization of N139 and N139D side chains reveals the presence of three rotamers, one of which faces the channel entrance. This out-facing conformer is metastable in the wild-type and in the N139D single mutant, but predominant in the double mutant thanks to the loss of electrostatic repulsion with the carboxylate group of D132. The effects of mutations and conformational isomerization on the pKa of E286, an essential proton-shuttling residue located at the top of the D channel, are shown to be consistent with the electrostatic control of proton pumping proposed recently (Fadda et al 2008 Biochim. Biophys. Acta 1777 277-84). Taken together, these results suggest that preserving the spatial distribution of charges at the entrance of the D channel is necessary to guarantee both the uptake and the relay of protons to the active site of the enzyme. These findings highlight the interplay of long-range electrostatic forces and local structural fluctuations in the control of proton movement and provide a physical explanation for the restoration of proton pumping activity in the double mutant.

Nanoparticles are promising drug delivery systems whose selection and optimization can be gainfully conducted by theoretical methods. This review is targeted to experimentalists who are interested in enhancing their time and cost efficiency through the incorporation of theoretical approaches. This review thus begins with a brief overview of theoretical approaches available to the development of contem- porary drug delivery systems. Approaches include solubility parameters, Flory-Huggins theory, analytical predictions of partition coeffi- cients, and molecular simulations. These methods are then compared as they relate to the optimization of drug-material pairs using important performance-related parameters including the size of the delivery particles, their surface properties, and the compatibility of the materials with the drug to be sequestered. Next, this review explores contemporary efforts to optimize a selection of existing nanoparticle platforms, including nanoemulsions, linear and star-shaped block co-polymer micelles, and dendrimers. The review concludes with an outlook on the challenges remaining in the successful application of these theoretical methods to the development of new drug formulations.

Enzymatic reactions involving bilayer lipids occur in an environment with strict physical and topological constraints. The integral membrane enzyme PagP transfers a palmitoyl group from a phospholipid to lipid A in order to assist Escherichia coli in evading host immune defenses during infection. PagP measures the palmitoyl group with an internal hydrocarbon ruler that is formed in the interior of the eight-stranded antiparallel barrel. The access and egress of the palmitoyl group is thought to take a lateral route from the bilayer phase to the barrel interior. Molecular dynamics, mutagenesis, and a 1.4 A crystal structure of PagP in an SDS 2-methyl-2,4-pentanediol (MPD) cosolvent system reveal that phospholipid access occurs at the crenel present between strands F and G of PagP. In this way, the phospholipid head group can remain exposed to the cell exterior while the lipid acyl chain remains in a predominantly hydrophobic environment as it translocates to the protein interior.

Magnesium translocation across cell membranes is essential for numerous physiological processes. Three recently reported crystal structures of the CorA magnesium transport system revealed a surprising architecture, with a bundle of giant alpha-helices forming a 60-A-long pore that extends beyond the membrane before widening into a funnel-shaped cytosolic domain. The presence of divalent cations in putative intracellular regulation sites suggests that these structures correspond to the closed conformation of CorA. To examine the nature of the conduction pathway, we performed 110-ns molecular-dynamics simulations of two of these structures in a lipid bilayer with and without regulatory ions. The results show that a 15-A-long hydrophobic constriction straddling the membrane-cytosol interface constitutes a steric bottleneck whose location coincides with an electrostatic barrier opposing cation translocation. In one of the simulations, structural relaxation after the removal of regulatory ions led to concerted changes in the tilt of the pore helices, resulting in iris-like dilation and spontaneous hydration of the hydrophobic neck. This simple and robust mechanism is consistent with the regulation of pore opening by intracellular magnesium concentration, and explains the unusual architecture of CorA.

Protein disorder is abundant in proteomes throughout all kingdoms of life and serves many biologically important roles. Disordered states of proteins are challenging to study experimentally due to their structural heterogeneity and tendency to aggregate. Computer simulations, which are not impeded by these properties, have recently emerged as a useful tool to characterize the conformational ensembles of intrinsically disordered proteins. In this review, we provide a survey of computational studies of protein disorder with an emphasis on the interdisciplinary nature of these studies. The application of simulation techniques to the study of disordered states is described in the context of experimental and bioinformatics approaches. Experimental data can be incorporated into simulations, and simulations can provide predictions for experiment. In this way, simulations have been integrated into the existing methodologies for the study of disordered state ensembles. We provide recent examples of simulations of disordered states from the literature and our own work. Throughout the review, we emphasize important predictions and biophysical understanding made possible through the use of simulations. This review is intended as both an overview and a guide for structural biologists and theoretical biophysicists seeking accurate, atomic-level descriptions of disordered state ensembles.

We use principal component analysis (PCA) to detect functionally interesting collective motions in molecular-dynamics simulations of membrane-bound gramicidin A. We examine the statistical and structural properties of all PCA eigenvectors and eigenvalues for the backbone and side-chain atoms. All eigenvalue spectra show two distinct power-law scaling regimes, quantitatively separating large from small covariance motions. Time trajectories of the largest PCs converge to Gaussian distributions at long timescales, but groups of small-covariance PCs, which are usually ignored as noise, have subdiffusive distributions. These non-Gaussian distributions imply anharmonic motions on the free-energy surface. We characterize the anharmonic components of motion by analyzing the mean-square displacement for all PCs. The subdiffusive components reveal picosecond-scale oscillations in the mean-square displacement at frequencies consistent with infrared measurements. In this regime, the slowest backbone mode exhibits tilting of the peptide planes, which allows carbonyl oxygen atoms to provide surrogate solvation for water and cation transport in the channel lumen. Higher-frequency modes are also apparent, and we describe their vibrational spectra. Our findings expand the utility of PCA for quantifying the essential features of motion on the anharmonic free-energy surface made accessible by atomistic molecular-dynamics simulations.

The inhibition of specific SH2 domain mediated protein-protein interactions as an effective chemotherapeutic approach in the treatment of diseases remains a challenge. That different conformations of peptide-ligands are preferred by different SH2 domains is an underappreciated observation from the structural analysis of phosphotyrosine peptide binding to SH2 domains that may aid in future drug design. To explore the nature of ligand binding, we use simulated annealing (SA) to sample the conformational space of phosphotyrosine-containing peptides complexed with the Src SH2 domain. While in good agreement with the crystallographic and NMR studies of high-affinity phosphopeptide-SH2 domain complexes, the results suggest that the structural basis for phopsphopeptide- Src SH2 interactions is more complex than the two-pronged plug two-hole socket model. A systematic study of peptides of type pYEEX, where pY is phosphotyrosine and X is a hydrophobic residue, indicates that these peptides can assume two conformations, one extended and one helical, representing the balance between the interaction of residue X with the hydrophobic hole on the surface of the Src SH2 domain, and its contribution to the inherent tendency of the two glutamic acids to form an -helix. In contrast, a -turn conformation, almost identical to that observed in the crystal structure of pYVNV bound to the Grb2 SH2 domain, predominates for pYXNX peptides, even in the presence of isoleucine at the third position. While peptide binding affinities, as measured by fluorescence polarization, correlate with the relative proportion of extended peptide conformation, these results suggest a model where all three residues C-terminal to the phosphotyrosine determine the conformation of the bound phosphopeptide. The information obtained in this work can be used in the design of specific SH2 domain inhibitors.

Star copolymers (SCPs) have recently attracted considerable attention due to their unique applicability in a wide range of biomedical fields. With the intention of rationally designing a stable unimolecular SCP, atomistic molecular dynamics simulations of thirteen SCPs are conducted. The SCPs each have six identical arms of methoxypoly(ethylene glycol)-b-polycaprolactone (MePEGx-b-PCLy) and systematically vary in terms of total molecular weight and ratio of hydrophobic to hydrophilic block length. For all hydrated SCPs, the simulations predict a densely packed hydrophobic PCL core that excludes water and is phase separated from a highly mobile hydrophilic PEG corona. The radii of the hydrophobic PCL core and the PEG blocks are independent of each other and can be predicted over a broad molecular weight range. A linear relationship between the hydration and the molecular weight of the PEG blocks is observed with the average number of water molecules bound per PEG repeat unit within the range of that determined experimentally. As well, a quantitative relationship relates the water accessed surface area of the hydrophobic PCL core to the molecular weights of PCL and PEG moieties. We postulate that the propensity for aggregation of SCPs into multimolecular micelles is correlated with the partial hydration of the hydrophobic core of unimers. Our results suggest that SCPs with a hydrophobic PCL core ≤2 kDa per arm are fully protected from water when the hydrophilic PEG blocks approach 14.6 kDa per arm. We therefore predict that SCPs of this composition yield unimolecular micelles that are thermodynamically stable at low concentrations.

Cytochrome c oxidase couples the reduction of dioxygen to proton pumping against an electrochemical gradient. The D-channel, a 25-A-long cavity, provides the principal pathway for the uptake of chemical and pumped protons. A water chain is thought to mediate the relay of protons via a Grotthuss mechanism through the D-channel, but it is interrupted at N139 in all available crystallographic structures. We use free-energy simulations to examine the proton uptake pathway in the wild type and in single-point mutants N139V and N139A, in which redox and pumping activities are compromised. We present a general approach for the calculation of water occupancy in protein cavities and demonstrate that combining efficient sampling algorithms with long simulation times (hundreds of nanoseconds) is required to achieve statistical convergence of equilibrium properties in the protein interior. The relative population of different conformational and hydration states of the D-channel is characterized. Results shed light on the role of N139 in the mechanism of proton uptake and clarify the physical basis for inactive phenotypes. The conformational isomerization of the N139 side chain is shown to act as a gate controlling the formation of a functional water chain or "proton wire." In the closed state of N139, the spatial distribution of water in the D-channel is consistent with available crystallographic models. However, a metastable state of N139 opens up a narrow bottleneck in which 50% occupancy by a water molecule establishes a proton pathway throughout the D-channel. Results for N139V suggest that blockage of proton uptake resulting from persistent interruption of the water pathway is the cause of this mutant's marginal oxidase activity. In contrast, results for N139A indicate that the D-channel is a continuously hydrated cavity, implying that the decoupling of oxidase activity from proton pumping measured in this mutant is not due to interruption of the proton relay chain.

Generalized-ensemble algorithms in temperature space have become popular tools to enhance conformational sampling in biomolecular simulations. A random walk in temperature leads to a corresponding random walk in potential energy, which can be used to cross over energetic barriers and overcome the problem of quasi-nonergodicity. In this paper, we introduce two novel methods: simulated tempering distributed replica sampling (STDR) and virtual replica exchange (VREX). These methods are designed to address the practical issues inherent in the replica exchange (RE), simulated tempering (ST), and serial replica exchange (SREM) algorithms. RE requires a large, dedicated, and homogeneous cluster of CPUs to function efficiently when applied to complex systems. ST and SREM both have the drawback of requiring extensive initial simulations, possibly adaptive, for the calculation of weight factors or potential energy distribution functions. STDR and VREX alleviate the need for lengthy initial simulations, and for synchronization and extensive communication between replicas. Both methods are therefore suitable for distributed or heterogeneous computing platforms. We perform an objective comparison of all five algorithms in terms of both implementation issues and sampling efficiency. We use disordered peptides in explicit water as test systems, for a total simulation time of over 42 μs. Efficiency is defined in terms of both structural convergence and temperature diffusion, and we show that these definitions of efficiency are in fact correlated. Importantly, we find that ST-based methods exhibit faster temperature diffusion and correspondingly faster convergence of structural properties compared to RE-based methods. Within the RE-based methods, VREX is superior to both SREM and RE. On the basis of our observations, we conclude that ST is ideal for simple systems, while STDR is well-suited for complex systems.

Oxygen and Ni(II) are ideal paramagnetic species for NMR studies of immersion depth since they establish prominent concentration gradients across the membrane-water interface of either bilayers or micelles. Corresponding gradients of paramagnetic shifts and relaxation rates are observed by NMR for membrane embedded amphiphiles. Specifically, upon dissolution of oxygen at a partial pressure of 20 bar or more, (13)C NMR spectra of membrane embedded amphiphiles reveal chemical shift perturbations which depend sensitively on average immersion depth in the membrane. Similarly, depth-dependent enhancements of spin-lattice relaxation rates can be detected by (1)H NMR. Generally, such paramagnetic effects depend both on steric or accessibility factors and on the local concentration of the paramagnet. The steric terms can be factored out by combining paramagnetic rates arising from O(2) and Ni, in the form of a ratio, R(1P)(O(2))/R(1P)(Ni). The natural logarithm of this ratio is shown to depend linearly on immersion depth in a micelle. The analysis is verified using molecular dynamics simulations of dodecylphosphocholine in a detergent micelle, while thorough consideration of the paramagnetic rate data also allows for the determination of the orientation of imipramine in the micelle. Thus, a complete picture of topology arises from this approach which is readily applicable to studies of drugs and amphiphiles in fast-tumbling bicelles, small unilamellar vesicles, and micelles.

to US and demonstrate how may enhance As our methodology, we utilize distributed function (DRPE) introduced by et al.

As part of the mitochondrial respiratory chain, cytochrome c oxidase utilizes the energy produced by the reduction of O2 to water to fuel vectorial proton transport. The mechanism coupling proton pumping to redox chemistry is unknown. Recent advances have provided evidence that each of the four observable transitions in the complex catalytic cycle consists of a similar sequence of events. However, the physico-chemical basis underlying this recurring sequence has not been identified. We identify this recurring pattern based on a comprehensive model of the catalytic cycle derived from the analysis of oxygen chemistry and available experimental evidence. The catalytic cycle involves the periodic repetition of a sequence of three states differing in the spatial distribution of charge in the active site: [0|1], [1|0], and [1|1], where the total charge of heme a and the binuclear center appears on the left and on the right, respectively. This sequence recurs four times per turnover despite differences in the redox chemistry. This model leads to a simple, robust, and reproducible sequence of electron and proton transfer steps and rationalizes the pumping mechanism in terms of electrostatic coupling of proton translocation to redox chemistry. Continuum electrostatic calculations support the proposed mechanism and suggest an electrostatic origin for the decoupled and inactive phenotypes of ionic mutants in the principal proton-uptake pathway.

Distributed replica sampling [T. Rodinger et al., J. Chem. Theory Comput. 2, 725 (2006)] is a simple and general scheme for Boltzmann sampling of conformational space by computer simulation in which multiple replicas of the system undergo a random walk in reaction coordinate or temperature space. Individual replicas are linked through a generalized Hamiltonian containing an extra potential energy term or bias which depends on the distribution of all replicas, thus enforcing the desired sampling distribution along the coordinate or parameter of interest regardless of free energy barriers. In contrast to replica exchange methods, efficient implementation of the algorithm does not require synchronicity of the individual simulations. The algorithm is inherently suited for large-scale simulations using shared or heterogeneous computing platforms such as a distributed network. In this work, we build on our original algorithm by introducing Boltzmann-weighted jumping, which allows moves of a larger magnitude and thus enhances sampling efficiency along the reaction coordinate. The approach is demonstrated using a realistic and biologically relevant application; we calculate the standard binding free energy of benzene to the L99A mutant of T4 lysozyme. Distributed replica sampling is used in conjunction with thermodynamic integration to compute the potential of mean force for extracting the ligand from protein and solvent along a nonphysical spatial coordinate. Dynamic treatment of the reaction coordinate leads to faster statistical convergence of the potential of mean force than a conventional static coordinate, which suffers from slow transitions on a rugged potential energy surface.

The structural basis of lipid acyl-chain selection by membrane-intrinsic enzymes is poorly understood because most integral membrane enzymes of lipid metabolism have proven refractory to structure determination; however, robust enzymes from the outer membranes of gram-negative bacteria are now providing a first glimpse at the underlying mechanisms. The methylene unit resolution of the phospholipid:lipid A palmitoyltransferase PagP is determined by the hydrocarbon ruler, a 16-carbon saturated acyl-chain-binding pocket buried within the transmembrane beta-barrel structure. Substitution of Gly88 lining the floor of the hydrocarbon ruler with Ala or Met makes the enzyme select specifically 15- or 12-carbon saturated acyl chains, respectively, indicating that hydrocarbon ruler depth determines acyl-chain selection. However, the Gly88Cys PagP resolution does not diminish linearly because it selects both 14- and 15-carbon saturated acyl chains. We discovered that an exciton, emanating from a buried Tyr26-Trp66 phenol-indole interaction, is extinguished by a local structural perturbation arising from the proximal Gly88Cys PagP sulfhydryl group. Site-specific S-methylation of the single Cys afforded Gly88Cys-S-methyl PagP, which reasserted both the exciton and methylene unit resolution by specifically selecting 13-carbon saturated acyl chains for transfer to lipid A. Unlike the other Gly88 substitutions, the Cys sulfhydryl group recedes from the hydrocarbon ruler floor and locally perturbs the subjacent Tyr26 and Trp66 aromatic rings. The resulting hydrocarbon ruler expansion thus occurs at the exciton's expense and accommodates an extra methylene unit in the selected acyl chain. The hydrocarbon ruler-exciton juxtaposition endows PagP with a molecular gauge for probing the structural basis of lipid acyl-chain selection in a membrane-intrinsic environment.

Elastin provides extensible tissues, including arteries and skin, with the propensity for elastic recoil, whereas amyloid fibrils are associated with tissue-degenerative diseases, such as Alzheimer's. Although both elastin-like and amyloid-like materials result from the self-organization of proteins into fibrils, the molecular basis of their differing physical properties is poorly understood. Using molecular simulations of monomeric and aggregated states, we demonstrate that elastin-like and amyloid-like peptides are separable on the basis of backbone hydration and peptide-peptide hydrogen bonding. The analysis of diverse sequences, including those of elastin, amyloids, spider silks, wheat gluten, and insect resilin, reveals a threshold in proline and glycine composition above which amyloid formation is impeded and elastomeric properties become apparent. The predictive capacity of this threshold is confirmed by the self-assembly of recombinant peptides into either amyloid or elastin-like fibrils. Our findings support a unified model of protein aggregation in which hydration and conformational disorder are fundamental requirements for elastomeric function.

It has been shown that a small number of FPGAs can significantly accelerate certain computing tasks by up to two or three orders of magnitude. However, particularly intensive large-scale computing applications, such as molecular dynamics simulations of biological systems, underscore the need for even greater speedups to address relevant length and time scales. In this work, we propose an architecture for a scalable computing machine built entirely using FPGA computing nodes. The machine enables designers to implement large-scale computing applications using a heterogeneous combination of hardware accelerators and embedded microprocessors spread across many FPGAs, all interconnected by a flexible communication network. Parallelism at multiple levels of granularity within an application can be exploited to obtain the maximum computational throughput. By focusing on applications that exhibit a high computation-to-communication ratio, we narrow the extent of this investigation to the development of a suitable communication infrastructure for our machine, as well as an appropriate programming model and design flow for implementing applications. By providing a simple, abstracted communication interface with the objective of being able to scale to thousands of FPGA nodes, the proposed architecture appears to the programmer as a unified, extensible FPGA fabric. A programming model based on the MPI message-passing standard is also presented as a means for partitioning an application into independent computing tasks that can be implemented on our architecture. Finally, we demonstrate the first use of our design flow by developing a simple molecular dynamics simulation application for the proposed machine, which runs on a small platform of development boards

Many different methods exist for computing free energy changes from molecular simulations. Recent advances have led to improvements in the theoretical framework underlying these calculations, as well as in the accuracy and sampling efficiency of the algorithms. Novel methods combining the advantages afforded by various existing approaches offer promising strategies and open up new perspectives to help elucidate the physical basis of important biological processes.

An optimized technique for calculating the excess chemical potential of small molecules in dense liquids and the binding affinity of molecular ligands to biomolecules is reported. In this method, a molecular species is coupled to the system of interest via a nonphysical fourth spatial dimension w through which insertion or extraction can be carried out [R. Pomes, E. Eisenmesser, C. B. Post et al., J. Chem. Phys. 111, 3387 (1999)]. Molecular simulations are used to compute the potential of mean force (PMF) acting on the solute molecule in the fourth dimension. The excess chemical potential of that molecule is obtained as the difference in the PMF between fully coupled and fully decoupled systems. The simplicity, efficiency, and generality of the method are demonstrated for the calculation of the hydration free energies of water and methanol as well as sodium, cesium, and chloride ions. A significant advantage over other methods is that the 4D-PMF approach provides a single effective and general route for decoupling all nonbonded interactions (i.e., both Lennard-Jones and Coulombic) at once for both neutral and charged solutes. Direct calculation of the mean force from thermodynamic integration is shown to be more computationally efficient than calculating the PMF from umbrella sampling. Statistical error analysis suggests a simple strategy for optimizing sampling. The detailed analysis of systematic errors arising from the truncation of Coulombic interactions in a solvent droplet of finite size leads to straightforward corrections to ionic hydration free energies.

Multiple ligand co-recognition of 3'-sulfogalactosylceramide (SGC) and sulfotyrosine initiated the comparison of SGC and sulfotyrosine and, subsequently, phosphotyrosine (pY) binding. SGC is a receptor for ligands involved in cell adhesion/microbial pathology. pY forms a Src homology domain 2 recognition motif in intracellular signaling. Using hsp70, anti-SGC, and anti-pY antibodies, ligand binding is retained following phosphate/sulfate and tyrosine/galactose substitution in SGC and sulfate/phosphate exchange in pY. Remarkable lipid-dependent binding to phosphatidylethanolamine-conjugated sulfotyrosine suggests "microenvironmental" modulation of sulfotyrosine-containing receptors, similar to glycosphingolipids. Based on an aryl substrate-bound co-crystal of arylsulfatase A, a sulfogalactose and phosphotyrosine esterase, modeling provides a solvation basis for co-recognition. c-Src/Src homology domain binding confirms our hypothesis, heralding a carbohydrate-based approach to regulation of phosphotyrosine-mediated recognition.

Detailed comparative studies of proton relay in native and chemically modified gramicidin channels provide a unique opportunity to uncover the structural basis of biological proton transport. The function of ion channels hinges on their ability to provide surrogate solvation in narrow pore filters so as to overcome the dielectric barrier presented by biological membranes. In the potassium channel KcsA and in the cation channel gramicidin, permeant selectivity and mobility are determined by the proteinaceous matrix via hydrogen bonding, charge-dipole, and dipole-dipole interactions. In particular, main-chain carbonyl groups in these pore interiors play an essential role in the solvation of alkali ions and of protons. In this study, molecular dynamics simulations reveal how the translocation of H+ is controlled by nanosecond conformational transitions exchanging distorted states of the peptidic backbone in the single-file region of a dioxolane-linked analogue of the gramicidin dimer. These results underline the functional role of channel dynamics and provide a mechanism for the modulation of proton currents by fluctuating dipoles.