A many-body model to study proteins. II. Incidence of many-body polarization effects on the interaction of the calmodulin protein with four Ca2+ dications and with a target enzyme peptide

Combining a polarizable force-field and a coarse-grained polarizable solvent model: application to long dynamics simulations of bovine pancreatic trypsin inhibitor

The dynamic coupling between a polarizable protein force field and a particle-based implicit solvent model is described. The polarizable force field, TCPEp, developed recently to simulate protein systems, is characterized by a reduced number of polarizable sites, with a substantial gain in efficiency for an equal chemical accuracy. The Polarizable Pseudo-Particle (PPP) solvent model represents the macroscopic solvent polarization by induced dipoles placed on mobile Lennard-Jones pseudo-particles. The solvent-induced dipoles are sensitive to the solute electric field, but not to each other, so that the computational cost of solvent-solvent interactions is basically negligible. The solute and solvent induced dipoles are determined self-consistently and the equations of motion are solved using an efficient iterative multiple time step procedure. The solvation cost with respect to vacuum simulations is shown to decrease with solute size: the estimated multiplicative factor is 2.5 for a protein containing about 1000 atoms, and as low as 1.15 for 8000 atoms. The model is tested for six 20 ns molecular dynamics trajectories of a traditional benchmark system: the hydrated Bovine Pancreatic Trypsin Inhibitor (BPTI). Even though the TCPEp parameters have not been refined to be used with the solvent PPP model, we observe a good conservation of the BPTI structure along the trajectories. Moreover, our approach is able to provide a description of the protein solvation thermodynamic at the same accuracy as the standard Poisson-Boltzman continuum methods. It provides in addition a good description of the microscopic structural aspects concerning the solute/solvent interaction.

Anisotropic, Polarizable Molecular Mechanics Studies of Inter- and Intramolecular Interactions and Ligand-Macromolecule Complexes. A Bottom-Up Strategy

We present an overview of the SIBFA polarizable molecular mechanics procedure, which is formulated and calibrated on the basis of quantum chemistry (QC). It embodies nonclassical effects such as electrostatic penetration, exchange-polarization, and charge transfer. We address the issues of anisotropy, nonadditivity, and transferability by performing parallel QC computations on multimolecular complexes. These encompass multiply H-bonded complexes and polycoordinated complexes of divalent cations. Recent applications to the docking of inhibitors to Zn-metalloproteins are presented next, namely metallo-beta-lactamase, phosphomannoisomerase, and the nucleocapsid of the HIV-1 retrovirus. Finally, toward third-generation intermolecular potentials based on density fitting, we present the development of a novel methodology, the Gaussian electrostatic model (GEM), which relies on ab initio-derived fragment electron densities to compute the components of the total interaction energy. As GEM offers the possibility of a continuous electrostatic model going from distributed multipoles to densities, it allows an inclusion of short-range quantum effects in the molecular mechanics energies. The perspectives of an integrated SIBFA/GEM/QM procedure are discussed.

Binding of calcium and other metal ions to the EF-hand loops of calmodulin studied by quantum chemical calculations and molecular dynamics simulations

Calcium ion binding by the four EF-hand motifs of the protein calmodulin (CaM) is a central event in Ca2+-based cellular signaling. To understand molecular details of this complex process, isolated Ca2+-binding loops can be studied, by use of both experiments and calculations. In this work, we explore the metal specificity of the four Ca2+-binding loops of CaM using density functional theory (DFT) quantum chemical calculations and molecular dynamics simulations. We study CaM complexes with the physiologically important ions of calcium (Ca2+) and magnesium (Mg2+) and also with two other ions, strontium (Sr2+) and lanthanum (La3+). The former is of interest in the area of radioactive waste bioremediation, whereas the latter is often used as a probe of Ca2+-binding sites. We obtain intrinsic metal ion-loop binding energies as well as their components: vacuum, charge-transfer, solvation, entropy, and deformation terms. A detailed analysis of the results reveals that the total binding energy depends on a delicate balance among these energy components. They, in turn, are determined by the cation's charge and size as well as the amino acid composition and flexibility of the loops and the identity of the metal-chelating residues.

Structural studies of human alkaline phosphatase in complex with strontium: implication for its secondary effect in bones

Strontium is used in the treatment of osteoporosis as a ranelate compound, and in the treatment of painful scattered bone metastases as isotope. At very high doses and in certain conditions, it can lead to osteomalacia characterized by impairment of bone mineralization. The osteomalacia symptoms resemble those of hypophosphatasia, a rare inherited disorder associated with mutations in the gene encoding for tissue-nonspecific alkaline phosphatase (TNAP). Human alkaline phosphatases have four metal binding sites--two for zinc, one for magnesium, and one for calcium ion--that can be substituted by strontium. Here we present the crystal structure of strontium-substituted human placental alkaline phosphatase (PLAP), a related isozyme of TNAP, in which such replacement can have important physiological implications. The structure shows that strontium substitutes the calcium ion with concomitant modification of the metal coordination. The use of the flexible and polarizable force-field TCPEp (topological and classical polarization effects for proteins) predicts that calcium or strontium has similar interaction energies at the calcium-binding site of PLAP. Since calcium helps stabilize a large area that includes loops 210-228 and 250-297, its substitution by strontium could affect the stability of this region. Energy calculations suggest that only at high doses of strontium, comparable to those found for calcium, can strontium substitute for calcium. Since osteomalacia is observed after ingestion of high doses of strontium, alkaline phosphatase is likely to be one of the targets of strontium, and thus this enzyme might be involved in this disease.

Modulating the affinity and the selectivity of engineered calmodulin EF-Hand peptides for lanthanides

A set of engineered peptides (33 amino acids long) corresponding to the helix-turn-helix (EF-Hand) motif of the metal-binding site I of the protein calmodulin from paramecium tetraurelia have been synthesized. A disulfide bridge has been introduced in the native sequence in order to stabilize a native-like conformation. The calcium-binding carboxylate residues in positions 20, 22, 24, and 31 were mutated into other amino acids and the influence of such mutations on the binding affinity of the peptides for calcium and lanthanides have been studied. It was shown that the binding affinity for terbium ions can be modulated with dissociation constants ranging from 40 nmolar to 40 mmolar. The study of the influence of the mutations on the terbium affinity showed that the residue in position 24 played a key role on the capability of the peptides to bind lanthanides and that the affinity could be enhanced by mutations on non-coordinating positions. Such peptides with high affinity for lanthanides may facilitate the development of new highly sensitive biosensors to monitor the metal pollution in the environment.