Affinity and specificity of protein U1A-RNA complex formation based on an additive component free energy model

An MM-GBSA computational protocol was used to investigate wild-type U1A-RNA and F56 U1A mutant experimental binding free energies. The trend in mutant binding free energies compared to wild-type is well-reproduced. Following application of a linear-response-like equation to scale the various energy components, the binding free energies agree quantitatively with observed experimental values. Conformational adaptation contributes to the binding free energy for both the protein and the RNA in these systems. Small differences in DeltaGs are the result of different and sometimes quite large relative contributions from various energetic components. Residual free energy decomposition indicates differences not only at the site of mutation, but throughout the entire protein. MM-GBSA and ab initio calculations performed on model systems suggest that stacking interactions may nearly, but not completely, account for observed differences in mutant binding affinities. This study indicates that there may be different underlying causes of ostensibly similar experimentally observed binding affinities of different mutants, and thus recommends caution in the interpretation of binding affinities and specificities purely by inspection.     

Estimation of the binding affinities of FKBP12 inhibitors using a linear response method.

A series of non-immunosuppressive inhibitors of FK506 binding protein (FKBP12) are investigated using Monte Carlo statistical mechanics simulations. These small molecules may serve as scaffolds for chemical inducers of protein dimerization, and have recently been found to have FKBP12-dependent neurotrophic activity. A linear response model was developed for estimation of absolute binding free energies based on changes in electrostatic and van der Waals energies and solvent-accessible surface areas, which are accumulated during simulations of bound and unbound ligands. With average errors of 0.5 kcal/mol, this method provides a relatively rapid way to screen the binding of ligands while retaining the structural information content of more rigorous free energy calculations.     

Calculations of antibody-antigen interactions: microscopic and semi-microscopic evaluation of the free energies of binding of phosphorylcholine analogs to McPC603.

The study of antibody-antigen interactions should greatly benefit from the development of quantitative models for the evaluation of binding free energies in proteins. The present work addresses this challenge by considering the test case of the binding free energies of phosphorylcholine analogs to the murine myeloma protein McPC603. This includes the evaluation of the differential binding energy as well as the absolute binding energies and their corresponding electrostatic contributions. Four different approaches are examined: the Protein Dipoles Langevin Dipoles (PDLD) method, the semi-microscopic PDLD (PDLD/S) method, a free energy perturbation (FEP) method based on an adiabatic charging procedure and a linear response approximation that accelerates the FEP calculation. The PDLD electrostatic calculations are augmented by estimates of the relevant hydrophobic and steric contributions. The determination of the hydrophobic energy involves an approach which considers the modification of the effective surface area of the solute by local field effects. The steric contributions are analyzed in terms of the corresponding reorganization energies. This treatment, which considers the protein as a harmonic system, views the steric forces as the restoring forces for the electrostatic interactions. The FEP method is found to give unreliable results with regular cut-off radii and starts to give quantitative results only in very expensive treatment with very large cut-off radii. The PDLD and PDLD/S methods are much faster than the FEP approach and give reasonable results for both the relative and absolute binding energies. The speed and simplicity of the PDLD/S method make it an effective strategy for interactive docking studies and indeed such an option is incorporated in the program MOLARIS. A component analysis of the different energy contributions of the FEP treatment and a similar PDLD analysis indicate that electrostatic effects provide the largest contribution to the differential binding energy, while the hydrophobic and steric contributions are much smaller. This finding lends further support to the idea that electrostatic interactions play a major role in determining the antigen specificity of McPC603.     

Evaluation of catalytic free energies in genetically modified proteins

A combination of the empirical valence bond method and a free energy perturbation approach is used to simulate the activity of genetically modified enzymes. The simulations reproduce in a semiquantitative way the observed effects of mutations on the activity and binding free energies of trypsin and subtilisin. This suggests that we are approaching a stage of quantitative structure-function correlation of enzymes. The analysis of the calculations points towards the electrostatic energy of the reacting system as the key factor in enzyme catalysis. The changes in the charges of the reacting system and the corresponding changes in "solvation" free energy (generalized here as the interaction between the charges and the given microenvironment) are emphasized. It is argued that a reliable evaluation of these changes might be sufficient for correlating structure and catalysis. The use of free energy perturbation methods and thermodynamic cycles for evaluation of solvation energies and reactivity is discussed, pointing out our early contributions. The apparent elaborated nature of our treatment is clarified, explaining that such a treatment is essential for consistent calculations of chemical reactions in polar environments. The problems associated with seemingly more rigorous quantum mechanical methods are discussed, emphasizing the inconsistency associated with using gas phase charge distributions. The importance of dynamic aspects is examined by evaluating the autocorrelation of the protein "reaction field" on the reacting substrate. It is found that, at least in the present case, dynamic effects are not important. The nature of the catalytic free energy is considered, arguing that the protein provides preoriented dipoles (polarized to stabilize the transition state charge distribution) and small reorganization energy, thus reducing the activation free energy. The corresponding catalytic free energy is related to the folding free energy, which is being invested in aligning the active site dipoles.     

Mechanism and thermodynamics of binding of the polypyrimidine tract binding protein to RNA

The polypyrimidine tract binding protein (PTB) is involved in many physiological processes, including alternative splicing, internal ribosomal entry side (IRES)-mediated initiation of translation, and polyadenylation, as well as in ensuring mRNA stability. However, the role of PTB in these processes is not fully understood, and this has motivated us to undertake a computational study of the protein. PTB RNA binding domains (RBDs) 3 and 4 and their complexes with oligopyrimidine RNAs were simulated using the GROMOS simulation software using the GROMOS 45A4 force field. First, the stability and fluctuations of the tertiary fold and of the secondary structural elements in individual domains, the combined RBD34 domain, and their complexes with RNA were studied. Second, the simulation results were validated against the experimental NMR NOE data. The analysis of hydrogen bonding patterns, salt bridge networks, and stacking interactions of the RNA to the binding pockets of the protein domains showed that binding is not sequence-specific and that many RNA fragments can bind to them successfully. Further calculations of the relative free energy of binding for different polypyrimidine sequences were carried out using the thermodynamic integration (TI) and single-step perturbation (SSP) methods. It is was not possible to calculate the relative free energies with high accuracy, but the obtained results do give qualitative insights into PTB's affinity for different RNA sequences. Furthermore, the low-energy conformations of the complexes that were found provided additional information about the mechanism of binding.     

Binding diversity of the two binding sites of ricin B lectin

Ricin B is a galactose-binding protein, which contains two binding sites. We have compared the binding properties of the two binding sites of ricin B chain toward different mono- and disaccharide ligands. The free energies of binding are calculated using the free energy perturbation simulation (thermodynamic integration method) and linear interaction energy approach using CHARMM force field. The second binding site of the protein was found to be weaker compared to the first. The details of the hydrogen-bonding scheme suggested the origin of the epimeric specificity of the protein. The reason for the weaker binding capacity of the second binding site has been addressed.     

Computational Prediction of Alanine Scanning and Ligand Binding Energetics in G-Protein Coupled Receptors

Site-directed mutagenesis combined with binding affinity measurements is widely used to probe the nature of ligand interactions with GPCRs. Such experiments, as well as structure-activity relationships for series of ligands, are usually interpreted with computationally derived models of ligand binding modes. However, systematic approaches for accurate calculations of the corresponding binding free energies are still lacking. Here, we report a computational strategy to quantitatively predict the effects of alanine scanning and ligand modifications based on molecular dynamics free energy simulations. A smooth stepwise scheme for free energy perturbation calculations is derived and applied to a series of thirteen alanine mutations of the human neuropeptide Y1 receptor and series of eight analogous antagonists. The robustness and accuracy of the method enables univocal interpretation of existing mutagenesis and binding data. We show how these calculations can be used to validate structural models and demonstrate their ability to discriminate against suboptimal ones.     

Evaluation of protein-protein association energies by free energy perturbation calculations

The association energy upon binding of different amino acids in the specificity pocket of trypsin was evaluated by free energy perturbation calculations on complexes between bovine trypsin (BT) and bovine pancreatic trypsin inhibitor (BPTI). Three simulations of mutations of the primary binding residue (P(1)) were performed (P(1)-Ala to Gly, P(1)-Met to Gly and P(1)-Met to Ala) and the resulting differences in association energy (DeltaDeltaG(a)) are 2. 28, 5.08 and 2.93 kcal/mol for P(1)-Ala to Gly, P(1)-Met to Gly and to Ala with experimental values of 1.71, 4.62 and 2.91 kcal/mol, respectively. The calculated binding free energy differences are hence in excellent agreement with the experimental binding free energies. The binding free energies, however, were shown to be highly dependent on water molecules at the protein-protein interface and could only be quantitatively estimated if the correct number of such water molecules was included. Furthermore, the cavities that were formed when a large amino acid side-chain is perturbed to a smaller one seem to create instabilities in the systems and had to be refilled with water molecules in order to obtain reliable results. In addition, if the protein atoms that were perturbed away were not replaced by water molecules, the simulations dramatically overestimated the initial state of the free energy perturbations.     

Site specific point mutation changes specificity: A molecular modeling study by free energy simulations and enzyme kinetics of the thermodynamics in ribonuclease T1 substrate interactions

We have theoretically and experimentally studied the binding of two different ligands to wild-type ribonuclease T₁ (RNT1) and to a mutant of RNT1 with Glu-46 replaced by Gln. The binding of the natural substrate 3′-GMP has been compared with the binding of a fluorescent probe, 2-aminopurine 3′-monophosphate (2AP), and relative free energies of binding of these ligands to the mutant and the wild-type (wt) enzyme have been calculated by free energy perturbation methods. The free energy perturbations predict that the mutant RNT1-Gln-46 binds 2AP better than 3′GMP, in agreement with experiments on dinucleotides. Four free energy perturbations, forming a closed loop, have been performed to allow the detection of systematic errors in the simulation procedure. Because of the larger number of atoms involved, it was necessary to use a much longer simulation time for the change in the protein, i.e., the perturbation from Glu to Gln, than in the perturbation from 3′-GMP to 2AP. Finally the structure of the binding site is analyzed for understanding differences in catalytic speed and binding strength. © 1993 Wiley-Liss, Inc.     

A computational study of ion binding and protonation states in the KcsA potassium channel

We report results from microscopic molecular dynamics and free energy perturbation simulations of the KcsA potassium channel based on its experimental atomic structure. Conformational properties of selected amino acid residues as well as equilibrium positions of K(+) ions inside the selectivity filter and the internal water cavity are examined. Positions three and four (counting from the extracellular site) in the experimental structure correspond to distinctly separate binding sites for K(+) ions inside the selectivity filter. The protonation states of Glu71 and Asp80, which are close to each other and to the selectivity filter, as well as K(+) binding energies are determined using free energy perturbation calculations. The Glu71 residue which is buried inside a protein cavity is found to be most stable in the neutral form while the solvent exposed Asp80 is ionized. The channel altogether exothermically binds up to three ions, where two of them are located inside the selectivity filter and one in the internal water cavity. Ion permeation mechanisms are discussed in relation to these results.     

Molecular dynamics and free energy analysis of neuraminidase-ligand interactions

We report molecular dynamics calculations of neuraminidase in complex with an inhibitor, 4-amino-2-deoxy-2,3-didehydro-N-acetylneuraminic acid (N-DANA), with subsequent free energy analysis of binding by using a combined molecular mechanics/continuum solvent model approach. A dynamical model of the complex containing an ionized Glu119 amino acid residue is found to be consistent with experimental data. Computational analysis indicates a major van der Waals component to the inhibitor-neuraminidase binding free energy. Based on the N-DANA/neuraminidase molecular dynamics trajectory, a perturbation methodology was used to predict the binding affinity of related neuraminidase inhibitors by using a force field/Poisson-Boltzmann potential. This approach, incorporating conformational search/local minimization schemes with distance-dependent dielectric or generalized Born solvent models, correctly identifies the most potent neuraminidase inhibitor. Mutation of the key ligand four-substituent to a hydrogen atom indicates no favorable binding free energy contribution of a hydroxyl group; conversely, cationic substituents form favorable electrostatic interactions with neuraminidase. Prospects for further development of the method as an analysis and rational design tool are discussed.     

An evaluation of Poisson-Boltzmann electrostatic free energy calculations through comparison with experimental mutagenesis data

For systems involving highly and oppositely charged proteins, electrostatic forces dominate association and contribute to biomolecular complex stability. Using experimental or theoretical alanine-scanning mutagenesis, it is possible to elucidate the contribution of individual ionizable amino acids to protein association. We evaluated our electrostatic free energy calculations by comparing calculated and experimental data for alanine mutants of five protein complexes. We calculated Poisson-Boltzmann electrostatic free energies based on a thermodynamic cycle, which incorporates association in a reference (Coulombic) and solvated (solution) state, as well as solvation effects. We observe that Coulombic and solvation free energy values correlate with experimental data in highly and oppositely charged systems, but not in systems comprised of similarly charged proteins. We also observe that correlation between solution and experimental free energies is dependent on dielectric coefficient selection for the protein interior. Free energy correlations improve as protein dielectric coefficient increases, suggesting that the protein interior experiences moderate dielectric screening, despite being shielded from solvent. We propose that higher dielectric coefficients may be necessary to more accurately predict protein-protein association. Additionally, our data suggest that Coulombic potential calculations alone may be sufficient to predict relative binding of protein mutants.     

Plant derived anti-cancerous secondary metabolites as multipronged inhibitor of COX,Topo, and aromatase: molecular modeling and dynamics simulation analyses

In the present study, 300 plant derived secondary metabolites (100 each of alkaloid, flavonoid, and terpenoid), have been screened for their anti-cancerous activity through inhibition of selected key enzymatic targets, namely cyclooxygenases (COXs), topoisomerases (Topos), and aromatase by molecular docking approach. Furthermore, the stability of the complexes of top hits, from each class of secondary metabolites, with their respective enzymatic targets was analyzed using molecular dynamics (MD) simulation analyses and binding free energy calculations. Analysis of the results of the docking in light of the pharmacokinetically screened 18 alkaloids, 26 flavonoids, and 9 terpenoids, revealed that the flavonoid, curcumin, was the most potent inhibitor for all the selected enzymatic targets. The stability of the complexes of COX-1, COX-2, Topo I, Topo IIβ and aromatase with the most potent inhibitor curcumin and those of the respective drugs, namely ibuprofen, aspirin, topotecan, etoposide, and exemestane were also analyzed through MD simulation analyses which revealed better stability of curcumin complexes than those of respective drugs. Binding energy calculations of the complexes of the curcumin with all the targets, except those of Topos, exhibited lower binding energies for the curcumin complexes than those of respective drugs which corroborated with the results of molecular docking analyses. Thus, the present study affirms the versatile and multipronged nature of curcumin, the traditionally used herbal medicine, as anti-cancer molecule directed against these enzymatic targets.     

Surface salt bridges,double-mutant cycles,and protein stability: an experimental and computational analysis of the interaction of the Asp 23 side chain with the N-terminus of the N-terminal domain of the ribosomal protein l9

Experimental and theoretical double-mutant cycles have been used to investigate a salt bridge in the N-terminal domain of the protein L9. Aspartic acid 23 is the only acidic residue involved in a well-defined pairwise interaction, namely, a partially solvent-exposed salt bridge with the protonated N-terminus of the protein. Mutations were studied in which Asp 23 was substituted by alanine, asparagine, and nitrile alanine. Interactions with the N-terminus were probed by comparisons between proteins with a protonated and acetylated N-terminus. The mutants were all folded, and the structures were unchanged from wild type as judged by CD and 2-D NMR. The coupling free energy between the N-terminus and the side chain of Asp 23 measured through double-mutant cycle analysis was favorable and ranged from -0.7 to -1.7 kcal mol(-)(1), depending upon the set of mutants used. This relatively large coupling free energy for a surface salt bridge likely arises from geometric factors that reduce the entropy loss associated with salt-bridge formation and from structural relaxation in the mutants. Coupling free energies computed with continuum electrostatic calculations agreed well with the experimental values when full account was taken of all potential interactions, particularly those involving Asp 23 and the acetylated N-terminus as well as interactions with solvent. The measured and calculated coupling free energy decreased only slightly when the salt concentration was increased from 100 to 750 mM NaCl. The calculations suggest that the coupling free energy between D23 and the N-terminus measured through the experimental double-mutant cycle analysis is significantly smaller than the actual interaction free energy between the groups in the wild-type structure because of the inapplicability of assumptions frequently used to interpret double-mutant cycles.     

Probing the effect of point mutations at protein-protein interfaces with free energy calculations

We have studied the effect of point mutations of the primary binding residue (P1) at the protein-protein interface in complexes of chymotrypsin and elastase with the third domain of the turkey ovomucoid inhibitor and in trypsin with the bovine pancreatic trypsin inhibitor, using molecular dynamics simulations combined with the linear interaction energy (LIE) approach. A total of 56 mutants have been constructed and docked into their host proteins. The free energy of binding could be reliably calculated for 52 of these mutants that could unambiguously be fitted into the binding sites. We find that the predicted binding free energies are in very good agreement with experimental data with mean unsigned errors between 0.50 and 1.03 kcal/mol. It is also evident that the standard LIE model used to study small drug-like ligand binding to proteins is not suitable for protein-protein interactions. Three different LIE models were therefore tested for each of the series of protein-protein complexes included, and the best models for each system turn out to be very similar. The difference in parameterization between small drug-like compounds and protein point mutations is attributed to the preorganization of the binding surface. Our results clearly demonstrate the potential of free energy calculations for probing the effect of point mutations at protein-protein interfaces and for exploring the principles of specificity of hot spots at the interface.     

Coordination topology and stability for the native and binding conformers of chymotrypsin inhibitor 2

We demonstrate that the stabilization of the binding region is accomplished at the expense of a loss in the stability of the rest of the protein. A novel molecular mechanics (MM) approach is introduced to distinguish residue stabilities of proteins in a given conformation. As an example, the relative stabilities of folded chymotrypsin inhibitor 2 (CI2) in unbound form, and CI2 in complex with subtilisin novo is investigated. The conformation of the molecule in the two states is almost identical, with an approximately 0.6-A root-mean-square deviation (RMSD) of the Calpha atoms. On binding, the packing density changes only at the binding loop. However, residue fluctuations in the rest of the protein are greatly altered solely due to those contacts, indicating the effective propagation of perturbation and the presence of remotely controlling residues. To quantify the interplay between packing density, packing order, residue fluctuations, and residue stability, we adopt an MM approach whereby small displacements are inserted at selected residues, followed by energy minimization; the displacement of each residue in response to such perturbations are organized in a perturbation-response matrix L. We define residue stability lambda(i) = summation operator((j)L(ij))/ summation operator((j) L(ji)) as the ratio of the amount of change to which the residue is amenable, to the ability of a given residue to induce change. We then define the free energy associated with residue stability, DeltaG(lambda) = -RT ln lambda. DeltaG(lambda) intrinsically selects the residues that are in the folding core. Upon complexation, the binding loop becomes more resistant to perturbation, in contrast to the alpha-helix that favors change. Although the two forms of CI2 are structurally similar, residue fluctuations differ vastly, and the stability of many residues is altered upon binding. The decrease in entropy introduced by binding is thus compensated by these changes.     

Computational analysis of binding of P1 variants to trypsin

The binding of P1 variants of bovine pancreatic trypsin inhibitor (BPTI) to trypsin has been investigated by means of molecular dynamics simulations. The specific interaction formed between the amino acid at the primary binding (P1) position of the binding loop of BPTI and the specificity pocket of trypsin was estimated by use of the linear interaction energy (LIE) method. Calculations for 13 of the naturally occurring amino acids at the P1 position were carried out, and the results obtained were found to correlate well with the experimental binding free energies. The LIE calculations rank the majority of the 13 variants correctly according to the experimental association energies and the mean error between calculated and experimental binding free energies is only 0.38 kcal/mole, excluding the Glu and Asp variants, which are associated with some uncertainties regarding protonation and the possible presence of counter-ions. The three-dimensional structures of the complex with three of the P1 variants (Asn, Tyr, and Ser) included in this study have not at present been solved by any experimental techniques and, therefore, were modeled on the basis of experimental data from P1 variants of similar size. Average structures were calculated from the MD simulations, from which specific interactions explaining the broad variation in association energies were identified. The present study also shows that explicit treatment of the complex water-mediated hydrogen bonding network at the protein-protein interface is of crucial importance for obtaining reliable binding free energies. The successful reproduction of relative binding energies shows that this type of methodology can be very useful as an aid in rational design and redesign of biologically active macromolecules.     

Molecular dynamics,free energy,and SPR analyses of the interactions between the SH2 domain of Grb2 and ErbB phosphotyrosyl peptides

We studied the interactions between the SH2 domain of growth factor receptor binding protein 2 (Grb2) and ErbB receptor-derived phosphotyrosyl peptides using molecular dynamics, free energy calculations, and surface plasmon resonance (SPR) analysis. Binding free energies for nine phosphotyrosyl peptides were calculated using the MM-PBSA continuum solvent method, and excellent qualitative agreement with the SPR experimental data, with a correlation coefficient of 0.92, was obtained. Consistent with previous experimental findings, phosphotyrosyl peptides with the consensus sequence pYXNX showed favorable binding affinity for the Grb2. Unexpectedly, phosphotyrosyl peptides with the consensus sequence pYQQD, which had not shown any specific binding affinity for the Grb2 in earlier studies, also showed favorable binding affinity for the Grb2 in our experimental and computational analyses. Component analysis of the calculated binding free energies revealed that van der Waals interaction between the Grb2 and the phosphotyrosyl peptide was the dominant factor for specificity and binding affinity. These results indicate that current methods of estimating binding free energies are efficient for obtaining important information about protein-protein interactions, which are essential for the transmission of signals in cellular signaling pathways.