Autodock Vina is a very popular, and highly cited, open source docking program. Here we present a scoring function which we call Vinardo (Vina RaDii Optimized). Vinardo is based on Vina, and was trained through a novel approach, on state of the art datasets. We show that the traditional approach to train empirical scoring functions, using linear regression to optimize the correlation of predicted and experimental binding affinities, does not result in a function with optimal docking capabilities. On the other hand, a combination of scoring, minimization, and re-docking on carefully curated training datasets allowed us to develop a simplified scoring function with optimum docking performance. This article provides an overview of the development of the Vinardo scoring function, highlights its differences with Vina, and compares the performance of the two scoring functions in scoring, docking and virtual screening applications. Vinardo outperforms Vina in all tests performed, for all datasets analyzed. The Vinardo scoring function is available as an option within Smina, a fork of Vina, which is freely available under the GNU Public License v2.0 from http://smina.sf.net. Precompiled binaries, source code, documentation and a tutorial for using Smina to run the Vinardo scoring function are available at the same address.

The interactions of the cryoprotective agent trehalose with a lipid membrane made of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine at 323 K were studied by means of molecular dynamics simulations. It was observed that trehalose binds to the phospholipid headgroups with its main axis parallel to the membrane normal. Trehalose establishes hydrogen bonds with the carbonyl and phosphate groups and replaces water molecules from the lipid headgroup. Notably, the number of hydrogen bonds (HBs) that the membrane made with its environment was conserved after trehalose binding. The HBs between lipid and trehalose have a longer lifetime than those established between lipid and water. The binding of the sugar does not produce changes either in the lipid area or in the lipid order parameter. The effect of trehalose on the dipole potential is in agreement with experimental results. The contribution of the different components to the membrane dipole potential was analyzed. It was observed that the binding of trehalose produces changes in the different components and the sugar itself contributes to the surface potential due to the polarization of its hydroxyl in the interface.

MutS, a member of the ABC ATPases superfamily, is a mismatch DNA-binding protein constituent of the DNA post-replicative mismatch repair system (MMRS). In this work, it is shown that the ATPase activity of Pseudomonas aeruginosa and Escherichia coli MutS is inhibited by ortho- and decavanadate. Structural comparison of the region involved in the ATP binding of E.coli MutS with the corresponding region of other ABC ATPases inhibited by vanadate, including the myosin- orthovanadate-Mg complex, showed that they are highly similar. From these results it is proposed that the orthovanadate inhibition of MutS ATPase can take place by a similar mechanism to that described for other ATPases. Docking of decavanadate on the ATP-binding region of MutS showed that the energetically more favorable interaction of this compound would take place with the complex MutS- ADP-Mg, suggesting that the inhibitory effect could be produced by a steric impediment of the protein ATP/ADP exchange. Besides the effect observed on the ATPase activity, vanadate also affects the DNA-binding capability of the protein, and partially inhibits the oligomerization of MutS and the temperature-induced inactivation of the protein. From the results obtained, and considering that vanadate is an intracellular trace component, this compound could be considered as a new modulator of the MMRS.

Chicken liver bile acid-binding protein (L-BABP) is a member of the fatty acid-binding proteins super family. The common fold is a beta-barrel of ten strands capped with a short helix-loop-helix motif called portal region, which is involved in the uptake and release of non-polar ligands. Using multiple-run molecular dynamics simulations we studied the interactions of L-BABP with lipid membranes of anionic and zwitterionic phospholipids. The simulations were in agreement with our experimental observations regarding the electrostatic nature of the binding and the conformational changes of the protein in the membrane. We observed that L-BABP migrated from the initial position in the aqueous bulk phase to the interface of anionic lipid membranes and established contacts with the head groups of phospholipids through the side of the barrel that is opposite to the portal region. The conformational changes in the protein occurred simultaneously with the binding to the membrane. Remarkably, these conformational changes were observed in the portal region which is opposite to the zone where the protein binds directly to the lipids. The protein was oriented with its macrodipole aligned in the configuration of lowest energy within the electric field of the anionic membrane, which indicates the importance of the electrostatic interactions to determine the preferred orientation of the protein. We also identified this electric field as the driving force for the conformational change. For all the members of the fatty acid-binding protein family, the interactions with lipid membranes is a relevant process closely related to the uptake, release and transfer of the ligand. The observations presented here suggest that the ligand transfer might not necessarily occur through the domain that directly interacts with the lipid membrane. The interactions with the membrane electric field that determine orientation and conformational changes described here can also be relevant for other peripheral proteins.

No abstract

The treatment of electrostatic interactions in molecular simulations is of fundamental importance. Ewald and related methods are being increasingly used to the detriment of cutoff schemes, which are known to produce several artifacts. A potential drawback of the Ewald method is the spatial periodicity that is imposed to the system, which could produce artifacts when applied in the simulation of liquids. In this work we analyze the octaalanine peptide with charged termini in explicit solvent, for which severe effects due to the use of Ewald sums were predicted using continuum electrostatics. Molecular Dynamics simulations for a total of 158 nanoseconds were performed in cells of different sizes. From the comparison of the results of different system sizes, no significant periodicity-induced artifacts were observed. It is argued that in current biomolecular simulations, the incomplete sampling is likely to affect the results to a larger extent than the artifacts induced by the use of Ewald sums.

Molecular dynamics simulations of proteins are usually performed on a single molecule, and coarse-grained protein models are calibrated using single-molecule simulations, therefore ignoring intermolecular interactions. We present here a new coarse-grained force field for the study of many protein systems. The force field, which is implemented in the context of the discrete molecular dynamics algorithm, is able to reproduce the properties of folded and unfolded proteins, in both isolation, complexed forming well-defined quaternary structures, or aggregated, thanks to its proper evaluation of protein-protein interactions. The accuracy and computational efficiency of the method makes it a universal tool for the study of the structure, dynamics, and association/dissociation of proteins.

The regulatory protein of the squid nerve sodium calcium exchanger (ReP1-NCXSQ) is a 15kDa soluble, intracellular protein that regulates the activity of the Na(+)/Ca(2+) exchanger in the squid axon. It is a member of the cellular retinoic acid-binding proteins family and the fatty acid-binding proteins superfamily. It is composed of ten beta strands defining an inner cavity and a domain of two short alpha helix segments. In this work, we studied the binding and orientation of ReP1-NCXSQ in anionic and zwitterionic lipid membranes using molecular dynamics (MD) simulations. Binding to lipid membranes was also measured by filtration binding assay. ReP1-NCXSQ acquired an orientation in the anionic membranes with the positive end of the macrodipole pointing to the lipid membrane. Potential of mean force calculations, in agreement with experimental measurements, showed that the binding to the anionic interfaces in low ionic strength was stronger than the binding to anionic interfaces in high ionic strength or to zwitterionic membranes. The results of MD showed that the electrostatic binding can be mediated not only by defined patches or domains of basic residues but also by a global asymmetric distribution of charges. A combination of dipole-electric field interaction and local interactions determined the orientation of ReP1-NCXSQ in the interface.