Marcos A. Villarreal

On the Ewald Artifacts in Computer Simulations. The Test-Case of the Octaalanine Peptide With Charged Termini

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.

Binding and interactions of L-BABP to lipid membranes studied by molecular dynamic simulations.

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.

Identification of a Site in Sar1 Involved in the Interaction with the Cytoplasmic Tail of Glycolipid Glycosyltransferases

Glycolipid glycosyltransferases (GGT) are transported from the endoplasmic reticulum (ER) to the Golgi, their site of residence, via COPII vesicles. An interaction of a (R/K)X(R/K) motif at their cytoplasmic tail (CT) with Sar1 is critical for the selective concentration in the transport vesicles. In this work using computational docking, we identify three putative binding pockets in Sar1 (sites A, B, and C) involved in the interaction with the (R/K)X(R/K) motif. Sar1 mutants with alanine replacement of amino acids in site A were tested in vitro and in cells. In vitro, mutant versions showed a reduced ability to bind immobilized peptides with the CT sequence of GalT2. In cells, Sar1 mutants (Sar1D198A) specifically affect the exiting of GGT from the ER, resulting in an ER/Golgi concentration ratio favoring the ER. Neither the typical Golgi localization of GM130 nor the exiting and transport of the G protein of the vesicular stomatitis virus were affected. The protein kinase inhibitor H89 produced accumulation of Sec23, Sar1, and GalT2 at the ER exit sites; Sar1D189A also accumulated at these sites, but in this case GalT2 remained disperse along ER membranes. The results indicate that amino acids in site A of Sar1 are involved in the interaction with the CT of GGT for concentration at ER exiting sites.

Thermodynamic and structural analysis of homodimeric proteins: Model of β-lactoglobulin

The energetics of protein homo-oligomerization was analyzed in detail with the application of a general thermodynamic model. We have studied the thermodynamic aspects of protein-protein interaction employing β-lactoglobulin A from bovine milk at pH=6.7 where the protein is mainly in its dimeric form. We performed differential calorimetric scans at different total protein concentration and the resulting thermograms were analyzed with the thermodynamic model for oligomeric proteins previously developed. The thermodynamic model employed, allowed the prediction of the sign of the enthalpy of dimerization, the analysis of complex calorimetric profiles without transitions baselines subtraction and the obtainment of the thermodynamic parameters from the unfolding and the association processes and the compared with association parameters obtained with Isothermal Titration Calorimetry performed at different temperatures. The dissociation and unfolding reactions were also monitored by Fourier-transform infrared spectroscopy and the results indicated that the dimer of β-lactoglobulin (N(2)) reversibly dissociates into monomeric units (N) which are structurally distinguishable by changes in their infrared absorbance spectra upon heating. Hence, it is proposed that β-lactoglobulin follows the conformational path induced by temperature:N(2)⇌2N⇌2D. The general model was validated with these results indicating that it can be employed in the study of the thermodynamics of other homo-oligomeric protein systems.

Inter-domain interactions in charged lipid monolayers.

Phase coexistence is common in model biomembranes with the presence of domains formed by lipids in a dense phase state modulating lateral diffusion of species through hydrodynamic and electrostatic interactions. In this study, interdomain interactions in monolayers of charged surfactants were analyzed and compared with neutral systems. Interactions were investigated at different interdomain distances and by varying the ionic strength (I) of the subphase. At low percentages of condensed area (%Ac), i.e., high interdomain distances, domains were approximated as point charges or dipoles, and a comparison between the simulated and experimental results was made. At high %Ac, domains were arranged in a distorted hexagonal lattice, and the energy of a domain around its equilibrium position in the lattice was modeled using a harmonic potential and the spring constant determined. On subphases of high I, charged domains interacted in a manner similar to neutral domains with domain motion being precluded at high percentages of condensed area. At low I, a higher interdomain repulsion was observed along with a lower domain motion and, therefore, a higher apparent viscosity at comparable %Ac. Interestingly, this effect was observed at conditions where the Debye-Hückel length was still 2 orders lower than the interdomain distances.

Vinardo: A Scoring Function Based on Autodock Vina Improves Scoring, Docking, and Virtual Screening.

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.

Kinetics of lipid-membrane binding and conformational change of L-BABP

We designed an experimental approach to differentiate the kinetics of protein binding to a lipid membrane from the kinetics of the associated conformational change in the protein. We measured the fluorescence intensity of the single Trp6 in chicken liver bile acid-binding protein (L-BABP) as a function of time after mixing the protein with lipid membranes. We mixed the protein with pure lipid membranes, with lipid membranes in the presence of a soluble quencher, and with lipid membranes containing a fluorescence quencher attached to the lipid polar head group. We fitted simultaneously the experimental curves to a three-state kinetic model. We conclude that in a first step, the binding of L-BABP to the interfacial region of the anionic lipid polar head groups occurred simultaneously with a conformational change to the partly unfolded state. In a second slower step, Trp6 buried within the polar head group region, releasing contacts with the aqueous phase.

Energetic and entropic contributions to the interactions between like-charged groups in cationic peptides: A molecular dynamics simulation study

The interaction between like‐charged amino acid residues has been proposed to stabilize the folded state of peptides and proteins, and to modulate the substrate binding and the action mechanism of enzymes. We have used an alanine‐ and lysine‐based peptide as a model system to study the interaction between like charges, and we have performed a 16‐nsec molecular dynamics simulation in solution. The calculated potential of mean force for the approach of the lysines Nζ atoms showed a minimum at a distance of 0.7 nm, in agreement with the separation probabilities obtained from analysis of protein crystal structures. The analysis of the individual energy components showed that the solvent polarization pays for the approach of the like charges and that the van der Waals energies do not contribute significantly. The entropic contributions have been divided in conformational and desolvation terms. Both terms favor the formation of the charge pair. A 10‐fold increase in counterion concentration was observed—with respect to its bulk concentration—next to the peptide charges, which helps to stabilize the peptide charges at a close distance.

Interactions of the fatty acid-binding protein ReP1-NCXSQ with lipid membranes. Influence of the membrane electric field on binding and orientation

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.

A Straightforward Approach for the Determination of the Maximum Time Step for the Simulation of Nanometric Metallic Systems.

In the present work, we report on a systematic analysis to determine the maximum time step allowed in molecular dynamics simulations applied to study metal systems of current interest in nanoscience. Using the velocity Verlet integration scheme, we have found that it is possible to use a 20 fs time step for the simulation of gold nanosystems. This is roughly an order of magnitude greater than the usually employed integration step (2 to 5 fs). We also propose a general criterion to select this maximum time step for other metallic nanosystems, even in the case of bimetallic nanosystems.