Simultaneous occupancy of the two binding sites in gramicidin A by monovalent cations is a well known property of this channel, but the energetic feasibility of this process in molecular dynamics simulations has not been established so far. Here the authors study the energetics of single and double ion occupancy in gramicidin A by constructing the potential of mean force for single and pair of ...

Recent studies reveal that the core sequences of many proteins were nearly optimized for stability by natural evolution. Surface residues, by contrast, are not so optimized, presumably because protein function is mediated through surface interactions with other molecules. Here, we sought to determine the extent to which the sequences of protein ligand-binding and enzyme active sites could be pr...

An efficient molecular simulation methodology has been developed for the evaluation of the druggability (ligandability) of a protein. Previously proposed techniques were designed to assess the druggability of crystallographic structures and cannot be tightly coupled to molecular dynamics (MD) simulations. By contrast, the present approach, JEDI (Just Exploring Druggability at protein Interfaces...

By using the lattice statistical argument, we have shown that for a protein whose subunits have the same number of neighbors, the three parameters (K(AB), K(BB), and K(S)K(t)) in the sequential theory formulated by Koshland, Nemethy, and Filmer [Biochemistry (1966) 5, 365] can be reduced to two parameters. One of the parameters, Z, measures the strength of the subunit interactions and is relate...

in this study, the binding of amlodipine (a ca ++ channel blocker) and arsenic (metalloid) to bovine serum albumin (bsa) was studied by equilibrium dialysis(ed) method in order to have an insight into their binding chemistry to bsa. free amlodipine concentration was increased due to addition of arsenic which reduced the binding of the compounds to bsa. however, the free fraction was not increas...

Recent advances in computational hardware, software and algorithms enable simulations of protein-ligand complexes to achieve timescales during which complete ligand binding and unbinding pathways can be observed. While observation of such events can promote understanding of binding and unbinding pathways, it does not alone provide information about the molecular drivers for protein-ligand assoc...

Implicit solvent hydration free energy models are an important component of most modern computational methods aimed at protein structure prediction, binding affinity prediction, and modeling of conformational equilibria. The nonpolar component of the hydration free energy, consisting of a repulsive cavity term and an attractive van der Waals solute-solvent interaction term, is often modeled usi...

The activation of BK type Ca(2+)-activated K(+) channels depends on both voltage and Ca(2+). We studied three point mutations in the putative voltage sensor S4 or S4-S5 linker regions in the mslo1 BK channels to explore the relationship between voltage and Ca(2+) in activating the channel. These mutations reduced the steepness of the open probability - voltage (P(o) - V) relation and increased ...

We recently developed a stochastic-based program that allows individual molecules in a cell signalling pathway to be simulated. This program has now been used to model the Tar complex, a multimeric signalling complex employed by coliform bacteria. This complex acts as a solid-state computational cassette, integrating and disseminating information on the presence of attractants and repellents in...

Genome-wide binding preferences of the key components of eukaryotic preinitiation complex (PIC) have been recently measured at high resolution in Saccharomyces cerevisiae by Rhee and Pugh. However, the rules determining the PIC binding specificity remain poorly understood. In this study, we show that nonconsensus protein-DNA binding significantly influences PIC binding preferences. We estimate ...