CIRSE: a solvation energy estimator compatible with flexible protein docking and design applications

New biochemical insights into the dynamics of water molecules at the GMP or IMP binding site of human GMPR enzyme: A molecular dynamics study

Human GMP reductase (hGMPR) enzyme is involved in a cellular metabolic pathway, converting GMP into IMP, and also it is an important target for anti-leukemic agents. Present computational investigations explain dynamical behavior of water molecules during the conformational transition process from GMP to IMP using molecular dynamics simulations. Residues at substrate-binding site of cancerous protein (PDB Id. 2C6Q) are mostly more dynamic in nature than the normal protein (PDB Id. 2BLE). Nineteen conserved water molecules are identified at the GMP/IMP binding site and are classified as (i) conserved stable dynamic and (ii) infrequent dynamic. Water molecules W11, W14, and W16 are classified as conserved stable dynamic due to their immobile character, whereas remaining water molecules (W1, W2, W3, W4, W5, W7, W8, W9, W10, W12, W13, W15, W17, W18, and W19) are infrequent with dynamic nature. Entrance or displacement of these infrequent water molecules at GMP/IMP sites may occur due to forward and backward movement of reference residues involving ligands. Four water molecules of hGMPR-I and nine water molecules of hGMPR-II are observed in repetitive transitions from GMP to IMP pathway, which indicates discrimination between two isoforms of hGMPRs. Water molecules in cancerous protein are more dynamic and unstable compared to normal protein. These water molecules execute rare dynamical events at GMP binding site and could assist in detailed understanding of conformational transitions that influence the hGMPR's biological functionality. The present study should be of interest to the experimental community engaged in leukemia research and drug discovery for CML cancer.

Structural and Dynamical Impact of Water Molecules at Substrate- or Product-Binding Sites in Human GMPR Enzyme: A Study by Molecular Dynamics Simulations

Human guanosine monophosphate reductase (hGMPR) enzyme maintains the intracellular balance between adenine and guanine nucleotide pools, and it is an excellent target for the design of isoform-specific antileukemic agents. In the present study, we have investigated solvation properties of substrate GMP or product inosine-5'-monophosphate (IMP)-binding pocket of hGMPR by employing molecular dynamics simulations on conformations A (substrate GMP), B [substrate GMP with cofactor nicotinamide adenine dinucleotide phosphate (NDP)], C (product IMP with cofactor NDP), and D (product IMP). Nineteen water sites are identified precisely; they are responsible for the catalytic activity of this site, control structural and dynamical integrity, and electronic consequences of GMP or IMP in the binding site of hGMPR. The water sites of category-1 (W1, W4, W5, W6, W13, and W15) in normal protein and category-2 (W2, W3, W7, W8, W10, W17, and W18) in cancerous protein are unique and stabilize the guanosine or inosine group of GMP or IMP for participation in the enzymatic reaction, whereas the remaining water centers either stabilize pentose sugar ribose or the phosphate group of GMP or IMP. Furthermore, water sites of category-4 (W11, W14, and W16) appear to be conserved in all conformations during the entire simulation. The GMP-binding site in cancerous protein 2C6Q is significantly expanded, and its dynamics are very different from normal protein 2BLE. Furthermore, unique interactions of GMP(N1)···W2···Asp129/Asn158, IMP(N1)···W3···Glu289, and IMP(O6)···W10···Ser270 might be used in a water mimic drug design for hGMPR-II. In this context, water finding probability, relative interaction energy (J) associated with water site W, entropy, and topologies of these three water sites are thermodynamically acceptable for the water displacement method by the modified ligand. Hence, their positions in the catalytic pocket may also facilitate future drug discovery for chronic myelogenous leukemia by the design of appropriately oriented chemical groups that may displace these water molecules to mimic their structural, electronic, and thermodynamic properties.

Improving the solubility, activity, and stability of reteplase using in silico design of new variants

Reteplase (recombinant plasminogen activator, r-PA) is a thrombolytic agent recombined from tissue-type plasminogen activator (t-PA), which has several prominent features such as strong thrombolytic ability and E. coli expressibility. Despite these outstanding features, it demonstrates reduced fibrin binding affinity, reduced stimulation of protease activity, and lower solubility, hence higher aggregation propensity, compared to t-PA. The present study was devoted to design r-PA variants with comparable structural stability, enhanced biological activity, and high solubility. For this purpose, computational molecular modeling techniques were utilized. The supercharging technique was applied for r-PA to designing new species of the protein. Based on the results from in silico evaluation of selected mutations in comparison to the wild-type r-PA, the designed supercharged mutant (S7 variant) exhibited augmented stability, decreased solvation energy, as well as enhanced binding affinity to fibrin. The data also implied increased plasminogen cleavage activity of the new variant. These findings have implications to therapies which involve removal of intravascular blood clots, including the treatment of acute myocardial infarction.

Rapid context-dependent ligand desolvation in molecular docking

In structure-based screens for new ligands, a molecular docking algorithm must rapidly score many molecules in multiple configurations, accounting for both the ligand's interactions with receptor and its competing interactions with solvent. Here we explore a context-dependent ligand desolvation scoring term for molecular docking. We relate the Generalized-Born effective Born radii for every ligand atom to a fractional desolvation and then use this fraction to scale an atom-by-atom decomposition of the full transfer free energy. The fractional desolvation is precomputed on a scoring grid by numerically integrating over the volume of receptor proximal to a ligand atom, weighted by distance. To test this method's performance, we dock ligands versus property-matched decoys over 40 DUD targets. Context-dependent desolvation better enriches ligands compared to both the raw full transfer free energy penalty and compared to ignoring desolvation altogether, though the improvement is modest. More compellingly, the new method improves docking performance across receptor types. Thus, whereas entirely ignoring desolvation works best for charged sites and overpenalizing with full desolvation works well for neutral sites, the physically more correct context-dependent ligand desolvation is competitive across both types of targets. The method also reliably discriminates ligands from highly charged molecules, where ignoring desolvation performs poorly. Since this context-dependent ligand desolvation may be precalculated, it improves docking reliability with minimal cost to calculation time and may be readily incorporated into any physics-based docking program.

An insight to the dynamics of conserved water molecular triad in IMPDH II (human): recognition of cofactor and substrate to catalytic Arg 322

Inosine 5' monophosphate dehydrogenase (IMPDH II) is a key enzyme involved in the de novo biosynthesis pathway of purine nucleotides and is also considered to be an excellent target for cancer inhibitor design. The conserve R 322 residue (in human) is thought to play some role in the recognition of inhibitor and cofactor through the catalytic D 364 and N 303. The 15 ns simulation and the water dynamics of the three different PDB structures (1B3O, 1NF7, and 1NFB) of human IMPDH by CHARMM force field have clearly indicated the involvement of three conserved water molecules (W(L), W(M), and W(C)) in the recognition of catalytic residues (R 322, D 364, and N 303) to inhibitor and cofactor. Both the guanidine nitrogen atoms (NH1 and NH 2) of the R 322 have anchored the di- and mono-nucleotide (cofactor and inhibitor) binding domains via the conserved W(C) and W(L) water molecules. Another conserved water molecule WM seems to bridge the two domains including the R 322 and also the W(C) and W(L) through seven centers H-bonding coordination. The conserved water molecular triad (W(C)-W(M)-W(L)) in the protein complex may thought to play some important role in the recognition of inhibitor and cofactor to the protein through R 322 residue.

An improved pairwise decomposable finite-difference Poisson-Boltzmann method for computational protein design

Our goal is to develop accurate electrostatic models that can be implemented in current computational protein design protocols. To this end, we improve upon a previously reported pairwise decomposable, finite difference Poisson-Boltzmann (FDPB) model for protein design (Marshall et al., Protein Sci 2005, 14, 1293). The improvement involves placing generic sidechains at positions with unknown amino acid identity and explicitly capturing two-body perturbations to the dielectric environment. We compare the original and improved FDPB methods to standard FDPB calculations in which the dielectric environment is completely determined by protein atoms. The generic sidechain approach yields a two to threefold increase in accuracy per residue or residue pair over the original pairwise FDPB implementation, with no additional computational cost. Distance dependent dielectric and solvent-exclusion models were also compared with standard FDPB energies. The accuracy of the new pairwise FDPB method is shown to be superior to these models, even after reparameterization of the solvent-exclusion model.

Designed protein-protein association

The analysis of natural contact interfaces between protein subunits and between proteins has disclosed some general rules governing their association. We have applied these rules to produce a number of novel assemblies, demonstrating that a given protein can be engineered to form contacts at various points of its surface. Symmetry plays an important role because it defines the multiplicity of a designed contact and therefore the number of required mutations. Some of the proteins needed only a single side-chain alteration in order to associate to a higher-order complex. The mobility of the buried side chains has to be taken into account. Four assemblies have been structurally elucidated. Comparisons between the designed contacts and the results will provide useful guidelines for the development of future architectures.

Solvent reaction field potential inside an uncharged globular protein: a bridge between implicit and explicit solvent models?

The solvent reaction field potential of an uncharged protein immersed in simple point charge/extended explicit solvent was computed over a series of molecular dynamics trajectories, in total 1560 ns of simulation time. A finite, positive potential of 13-24 kbTec(-1) (where T=300 K), dependent on the geometry of the solvent-accessible surface, was observed inside the biomolecule. The primary contribution to this potential arose from a layer of positive charge density 1.0 A from the solute surface, on average 0.008 ec/A3, which we found to be the product of a highly ordered first solvation shell. Significant second solvation shell effects, including additional layers of charge density and a slight decrease in the short-range solvent-solvent interaction strength, were also observed. The impact of these findings on implicit solvent models was assessed by running similar explicit solvent simulations on the fully charged protein system. When the energy due to the solvent reaction field in the uncharged system is accounted for, correlation between per-atom electrostatic energies for the explicit solvent model and a simple implicit (Poisson) calculation is 0.97, and correlation between per-atom energies for the explicit solvent model and a previously published, optimized Poisson model is 0.99.