Protein docking using continuum electrostatics and geometric fit

The computer program DOT quickly finds low-energy docked structures for two proteins by performing a systematic search over six degrees of freedom. A novel feature of DOT is its energy function, which is the sum of both a Poisson-Boltzmann electrostatic energy and a van der Waals energy, each represented as a grid-based correlation function. DOT evaluates the energy of interaction for many orientations of the moving molecule and maintains separate lists scored by either the electrostatic energy, the van der Waals energy or the composite sum of both. The free energy is obtained by summing the Boltzmann factor over all rotations at each grid point. Three important findings are presented. First, for a wide variety of protein-protein interactions, the composite-energy function is shown to produce larger clusters of correct answers than found by scoring with either van der Waals energy (geometric fit) or electrostatic energy alone. Second, free-energy clusters are demonstrated to be indicators of binding sites. Third, the contributions of electrostatic and attractive van der Waals energies to the total energy term appropriately reflect the nature of the various types of protein-protein interactions studied.     

Molecular mechanical perspective on halogen bonding

The nature and strength of halogen bonding in halo molecule-Lewis base complexes were studied in terms of molecular mechanics using our recently developed positive extra-point (PEP) approach, in which the σ-hole on the halogen atom is represented by an extra point of positive charge. The contributions of the σ-hole (i.e., positively charged extra point) and the halogen atom to the strength of this noncovalent interaction were clarified using the atomic parameter contribution to the molecular interaction (APCtMI) approach. The molecular mechanical results revealed that the halogen bond is electrostatic and van der Waals in nature, and its strength depends on three types of interaction: (1) the attractive electrostatic interaction between the σ-hole and the Lewis base, (2) the repulsive electrostatic interaction between the negative halogen atom and the Lewis base, and (3) the repulsive/attractive van der Waals interactions between the halogen atom and the Lewis base. The strength of the halogen bond increases with increasing σ-hole size (i.e., magnitude of the extra-point charge) and increasing halogen atom size. The van der Waals interaction's contribution to the halogen bond strength is most favorable in chloro complexes, whereas the electrostatic interaction is dominant in iodo complexes. The idea that the chloromethane molecule can form a halogen bond with a Lewis base was revisited in terms of quantum mechanics and molecular mechanics. Although chloromethane does produce a positive region along the C-Cl axis, basis set superposition error corrected second-order M?ller-Plesset calculations showed that chloromethane-Lewis base complexes are unstable, producing halogen-Lewis base contacts longer than the sum of the van der Waals radii of the halogen and O/N atoms. Molecular mechanics using the APCtMI approach showed that electrostatic interactions between chloromethane and a Lewis base are unfavorable owing to the high negative charge on the chlorine atom, which overcomes the corresponding favorable van der Waals interactions.     

Many-atom van der Waals interactions lead to direction-sensitive interactions of covalent bonds

Strict physical theory and numerical calculations show that a specific coupling of many-atom van der Waals interactions with covalent bonding can significantly (half as much) increase the strength of attractive dispersion interactions when the direction of interaction coincides with the direction of the covalent bond, and decrease this strength when the direction of interaction is perpendicular to the direction of the covalent bond. The energy effect is comparable to that caused by the replacement of atoms (e.g. N by C or O) in conventional pairwise van der Waals interactions. Analysis of protein structures shows that they bear an imprint of this effect. This means that many-atom van der Waals interactions cannot be ignored in refinement of protein structures, in simulations of their folding, and in prediction of their binding affinities.     

Influence of cation-pi interactions in different folding types of membrane proteins

Cation-pi interactions play an important role to the stability of protein structures. In this work, we analyze the influence of cation-pi interactions in three-dimensional structures of membrane proteins. We found that transmembrane strand (TMS) proteins have more number of cation-pi interactions than transmembrane helical (TMH) proteins. In TMH proteins, both the positively charged residues Lys and Arg equally experience favorable cation-pi interactions whereas in TMS proteins, Arg is more likely than Lys to be in such interactions. There is no relationship between number of cation-pi interactions and number of residues in TMH proteins whereas a good correlation was observed in TMS proteins. The average cation-pi interaction energy for TMH proteins is -16 kcal/mol and that for TMS proteins is -27 kcal/mol. The pair-wise cation-pi interaction energy between aromatic and positively charged residues showed that Lys-Trp energy is stronger in TMS proteins than TMH proteins; Arg-Phe, Arg-Tyr and Lys-Phe have higher energy in TMH proteins than TMS proteins. The decomposition of energies into electrostatic and van der Waals revealed that the contribution from electrostatic energy is twice as that from van der Waals energy in both TMH and TMS proteins. The results obtained in the present study would be helpful to understand the contribution of cation-pi interactions to the stability of membrane proteins.     

Protein-nucleic acid recognition: statistical analysis of atomic interactions and influence of DNA structure

We analyzed structural features of 11,038 direct atomic contacts (either electrostatic, H-bonds, hydrophobic, or other van der Waals interactions) extracted from 139 protein-DNA and 49 protein-RNA nonhomologous complexes from the Protein Data Bank (PDB). Globally, H-bonds are the most frequent interactions (approximately 50%), followed by van der Waals, hydrophobic, and electrostatic interactions. From the protein viewpoint, hydrophilic amino acids are over-represented in the interaction databases: Positively charged amino acids mainly contact nucleic acid phosphate groups but can also interact with base edges. From the nucleotide point of view, DNA and RNA behave differently: Most protein-DNA interactions involve phosphate atoms, while protein-RNA interactions involve more frequently base edge and ribose atoms. The increased participation of DNA phosphate involves H-bonds rather than salt bridges. A statistical analysis was performed to find the occurrence of amino acid-nucleotide pairs most different from chance. These pairs were analyzed individually. Finally, we studied the conformation of DNA in the interaction sites. Despite the prevalence of B-DNA in the database, our results suggest that A-DNA is favored in the interaction sites.     

Probing the electronic structures and properties of neutral and anionic ScSi n (0,−1) (n = 1–6) clusters using ccCA-TM and G4 theory

We apply molecular docking, molecular dynamics (MD) simulation, and binding free energy calculation to investigate and reveal the binding mechanism between five xanthine inhibitors and DPP-4. The electrostatic and van der Waals interactions of the five inhibitors with DPP-4 are analyzed and discussed. The computed binding free energies using MM-PBSA method are in qualitatively agreement with experimental inhibitory potency of five inhibitors. The hydrogen bonds of inhibitors with Ser630 and Asp663 can stabilize the inhibitors in binding sites. The van der Waals interactions, especially the key contacts with His740, Asn710, Trp629, and Tyr666 have larger contributions to the binding free energy and play important roles in distinguishing the variant bioactivity of five inhibitors.     

Contribution of enthalpy to the energetics of complex formation of aromatic ligands with DNA

The energy contributions of electrostatic, van der Waals interactions, hydrogen bonds, and interactions of charge transfer type to the enthalpy of complex formation of the double-stand DNA with the antitumor antibiotics daunomycin, nogalamycin, and novantron, as well as the mutagens ethidium bromide and proflavine have been calculated. According to the calculations, the van der Waals component (except for nogalamycin) is energetically favorable during complex formation of the antibiotics with DNA, and the contributions of H bonds and electrostatic interactions are unfavorable, with the probability of charge transfer in the complexes being low. It has been shown that the relatively low value of the experimental enthalpy of binding is the sum of components greater in absolute value and different in the sign, which is the cause of large errors in estimating the total enthalpy of complex formation of aromatic ligands with DNA.     

Conformational transitions in RNA single uridine and adenosine bulge structures: a molecular dynamics free energy simulation study

Extra unmatched nucleotides (single base bulges) are common structural motifs in folded RNA molecules and can participate in RNA-ligand binding and RNA tertiary structure formation. Often these processes are associated with conformational transitions in the bulge region such as flipping out of the bulge base from an intrahelical stacked toward a looped out state. Knowledge of the flexibility of bulge structures and energetics of conformational transitions is an important prerequisite to better understand the function of this RNA motif. Molecular dynamics simulations were performed on single uridine and adenosine bulge nucleotides at the center of eight basepair RNA molecules and indicated larger flexibility of the bulge bases compared to basepaired regions. The umbrella sampling method was applied to study the bulge base looping out process and accompanying conformational and free energy changes. Looping out toward the major groove resulted in partial disruption of adjacent basepairs and was found to be less favorable compared to looping out toward the minor groove. For both uridine and adenosine bulges, a positive free energy change for full looping out was obtained which was approximately 1.5 kcal mol-1 higher in the case of the adenosine compared to the uridine bulge system. The simulations also indicated stable partially looped out states with the bulge bases located in the RNA minor groove and forming base triples with 5'-neighboring basepairs. In the case of the uridine bulge this state was more stable than the intrahelical stacked bulge structure. Induced looping out toward the minor groove involved crossing of an energy barrier of approximately 3.5 kcal mol-1 before reaching the base triple state. A continuum solvent analysis of intermediate bulge states indicated that electrostatic interactions stabilize looped out and base triple states, whereas van der Waals interactions and nonpolar contributions favor the stacked bulge conformation.     

Stability of DNA duplexes with Watson-Crick base pairs: a predicted model

The conformational stability (difference between the free energies of the folded and unfolded states, DeltaG degrees ) of a DNA duplex is considered as a function of component energy terms, hydrophobic, base stacking, hydrogen bonding, van der Waals, and electrostatic, and a trinucleotide-level helix stiffness parameter measured in terms of its Young's modulus. Hydrophobic and base stacking energy components were determined with the use of the crystal structure data of 30 DNA duplexes judicially selected within a resolution of 1.5 A, and hydrogen bonding, van der Waals and electrostatic terms were determined through an extensive review of experimental and theoretical studies. The stiffness indices for the trinucleotides were the ones realized by M. M. Gromiha [(2000) J. Biol. Phys. 26, 43-50] using the crystal structure data of 70 DNA duplexes. The unfolded state was treated in the classical way to determine its stability. Thermodynamically determined DeltaG degrees values for 111 DNA duplexes, with the number of base pairs ranging from 4 to 16, were selected in two sets, and the regression equation formed with one set was used to predict the stabilities of the other set, taking the energy components and the stiffness parameter to be independent variables. The computed energy terms indicate that the base stacking and hydrogen bonding forces are the dominant and the hydrophobic and electrostatic forces the weak partners in imparting stability to the duplexes. This model predicts DeltaG degrees values for DNA duplexes examined with a level of accuracy similar to that used for predictions made by the widely used nearest-neighbor models. The uniqueness of this model is that it combines the crystal and thermodynamic data for interpretation of conformational stability.     

Energetics of complex formation of the DNA hairpin structure d(GCGAAGC) with aromatic ligands

The energy contributions of various physical interactions to the total Gibbs energy of complex formation of the biologically important DNA hairpin d(GCGAAGC) with aromatic antitumor antibiotics daunomycin and novantron and the mutagens ethidium and proflavine have been calculated. It has been shown that the relatively small value of the total energy of binding of the ligands to the hairpin is the sum of components great in absolute value and different in sign. The contributions of van der Waals interactions and both intra- and intermolecular hydrogen bonds and bonds with aqueous environment have been studied. According to the calculations, the hydrophobic and van der Waals components are energetically favorable in complex formation of the ligands with the DNA pairpin d(GCGAAGC), whereas the electrostatic (with consideration of hydrogen bonds) and entropic components are unfavorable.     

Energetics of complex formation of the dna hairpin structure d(GCGAAGC) with aromatic ligands

The energy contributions of various physical interactions to the total Gibbs energy of complex formation of the biologically important DNA hairpin d(GCGAAGC) with aromatic antitumor antibiotics daunomycin and novantrone and mutagens ethidium bromide and proflavin have been calculated. It has been shown that the relatively small value of the total energy of binding of the ligands to the hairpin is the sum of components great in absolute value and different in sign. The contributions of van der Waals interactions and both intra- and intermolecular hydrogen bonds and bonds with aqueous environment have been studied. According to the calculations, the hydrophobic and van der Waals components are energetically favorable in complex formation of the ligands with the hairpin d(GCGAAGC), whereas the electrostatic (with consideration of hydrogen bonds) and entropic components are unfavorable.     

On the importance of van der Waals interaction in the groove binding of DNA with ligands: restrained molecular dynamics study

Comparison of interaction energy between an oligonucleotide and a DNA-binding ligand in the minor and major groove modes was made by use of restrained molecular dynamics. Distortion in DNA was found for the major groove mode whereas less significant changes for both ligand and DNA were detected for the minor groove binding after molecular dynamics simulation. The conformation of the ligand obtained from the major groove mode resembles that computed with the ligand soaked in water. The van der Waals contact energy was found to be as significant as electrostatic energy and more important for difference in binding energy between these two binding modes. The importance of van der Waals force in groove binding was supported by computations on the complex formed by the repressor peptide fragment from the bacteriophage 434 and its operator oligonucleotide.     

Energetic analysis of binding of progesterone and 5 beta-androstane-3,17-dione to anti-progesterone antibody DB3 using molecular dynamics and free energy calculations.

Molecular dynamics simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) free energy calculations were used to study the energetics of the binding of progesterone (PRG) and 5 beta-androstane-3,17-dione (5AD) to anti-PRG antibody DB3. Although the two steroids bind to DB3 in different orientations, their binding affinities are of the same magnitude, 1 nM for PRG and 8 nM for 5AD. The calculated relative binding free energy of the steroids, 8.8 kJ/mol, is in fair agreement with the experimental energy, 5.4 kJ/mol. In addition, computational alanine scanning was applied to study the role of selected amino acid residues of the ligand-binding site on the steroid cross-reactivity. The electrostatic and van der Waals components of the total binding free energies were found to favour more the binding of PRG, whereas solvation energies were more favourable for the binding of 5AD. The differences in the free energy components are due to the binding of the A rings of the steroids to different binding pockets: PRG is bound to a pocket in which electrostatic antibody-steroid interactions are dominating, whereas 5AD is bound to a pocket in which van der Waals and hydrophobic interactions dominate.     

Molecular docking and dynamics simulations on the interaction of cationic porphyrin–anthraquinone hybrids with DNA G-quadruplexes

A series of cationic porphyrin–anthraquinone hybrids bearing either pyridine, imidazole, or pyrazole rings at the meso-positions have been investigated for their interaction with DNA G-quadruplexes by employing molecular docking and molecular dynamics simulations. Three types of DNA G-quadruplexes were utilized, which comprise parallel, antiparallel, and mixed hybrid topologies. The porphyrin hybrids have a preference to bind with parallel and mixed hybrid structures compared to the antiparallel structure. This preference arises from the end stacking of porphyrin moiety following G-stem and loop binding of anthraquinone tail, which is not found in the antiparallel due to the presence of diagonal and lateral loops that crowd the G-quartet. The binding to the antiparallel, instead, occurred with poorer affinity through both the loop and wide groove. All sites of porphyrin binding were confirmed by 6 ns molecular dynamics simulation, as well as by the negative value of the total binding free energies that were calculated using the MMPBSA method. Free energy analysis shows that the favorable contribution came from the electrostatic term, which supposedly originated from the interaction of either cationic pyridinium, pyrazole, or imidazole groups and the anionic phosphate backbone, and also from the van der Waals energy, which primarily contributed through end stacking interaction.     

The impact of <Emphasis Type="Italic">Trichoderma reesei</Emphasis> Cel7A carbohydrate binding domain mutations on its binding to a cellulose surface: a molecular dynamics free energy study

A critical role of the Family 7 cellobiohydrolase (Cel7A) carbohydrate binding domain (CBD) is to bind to a cellulose surface and increase the enzyme concentration on the surface. Several residues of Trichoderma reesei Cel7A CBD, including Y5, N29, Y31, Y32 and Q34, contribute to cellulose binding, as revealed by early experimental studies. To investigate the interactions between these important residues and cellulose, we applied a thermodynamic integration method to calculate the cellulose–Cel7A CBD binding free energy changes caused by Y5A, N29A, Y31A, Y32A and Q34A mutations. The experimental binding trend was successfully predicted, proving the effectiveness of the complex model. For the two polar residue mutants N29A and Q34A, the changes in the electrostatics are comparable to those of van der Waals, while for three Y to A mutants, the free energy differences mainly come from van der Waals interactions. However, in both cases, the electrostatics dominates the interactions between individual residues and cellulose. The side chains of these residues are rigidified after the complex is formed. The binding free energy changes for the two mutants Y5W and Y31W were also determined, and for these the van der Waals interaction was strengthened but the electrostatics was weakened.     

Comparative binding energy analysis of the substrate specificity of haloalkane dehalogenase from Xanthobacter autotrophicus GJ10

Comparative binding energy (COMBINE) analysis was conducted for 18 substrates of the haloalkane dehalogenase from Xanthobacter autotrophicus GJ10 (DhlA): 1-chlorobutane, 1-chlorohexane, dichloromethane, 1,2-dichloroethane, 1,2-dichloropropane, 2-chloroethanol, epichlorohydrine, 2-chloroacetonitrile, 2-chloroacetamide, and their brominated analogues. The purpose of the COMBINE analysis was to identify the amino acid residues determining the substrate specificity of the haloalkane dehalogenase. This knowledge is essential for the tailoring of this enzyme for biotechnological applications. Complexes of the enzyme with these substrates were modeled and then refined by molecular mechanics energy minimization. The intermolecular enzyme-substrate energy was decomposed into residue-wise van der Waals and electrostatic contributions and complemented by surface area dependent and electrostatic desolvation terms. Partial least-squares projection to latent structures analysis was then used to establish relationships between the energy contributions and the experimental apparent dissociation constants. A model containing van der Waals and electrostatic intermolecular interaction energy contributions calculated using the AMBER force field explained 91% (73% cross-validated) of the quantitative variance in the apparent dissociation constants. A model based on van der Waals intermolecular contributions from AMBER and electrostatic interactions derived from the Poisson-Boltzmann equation explained 93% (74% cross-validated) of the quantitative variance. COMBINE models predicted correctly the change in apparent dissociation constants upon single-point mutation of DhlA for six enzyme-substrate complexes. The amino acid residues contributing most significantly to the substrate specificity of DhlA were identified; they include Asp124, Trp125, Phe164, Phe172, Trp175, Phe222, Pro223, and Leu263. These residues are suitable targets for modification by site-directed mutagenesis.     

Swelling of phospholipids by monovalent salt

Critical to biological processes such as membrane fusion and secretion, ion-lipid interactions at the membrane-water interface still raise many unanswered questions. Using reconstituted phosphatidylcholine membranes, we confirm here that multilamellar vesicles swell in salt solutions, a direct indication that salt modifies the interactions between neighboring membranes. By varying sample histories, and by comparing with data from ion carrier-containing bilayers, we eliminate the possibility that swelling is an equilibration artifact. Although both attractive and repulsive forces could be modified by salt, we show experimentally that swelling is driven primarily by weakening of the van der Waals attraction. To isolate the effect of salt on van der Waals interactions, we focus on high salt concentrations at which any possible electrostatic interactions are screened. By analysis of X-ray diffraction data, we show that salt does not alter membrane structure or bending rigidity, eliminating the possibility that repulsive fluctuation forces change with salt. By measuring changes in interbilayer separation with applied osmotic stress, we have determined, using the standard paradigm for bilayer interactions, that 1 M concentrations of KBr or KCl decrease the van der Waals strength by 50%. By weakening van der Waals attractions, salt increases energy barriers to membrane contact, possibly affecting cellular communication and biological signaling.     

Crystal packing of pyrimidine nucleosides--a study in terms of empirical van der Waals' potential energy

Molecular interactions in 32 crystal structures of pyrimidine nucleosides and nucleoside derivatives have been studied in terms of empirical van der Waals' energy. The stacking patterns have been classified into three types according to the symmetry operations involved within the pattern. Halogen and sulphur atoms were included in the calculations where necessary. Model columns of stacked molecules were set up in a computer calculation giving the energy for different amounts of base overlap. There was a tendency for the crystal positions of the molecules to agree with the minima of the energy maps, which showed the predominance of the van der Waals' forces in determining the base-stacking arrangement. The position of the energy minimum was correlated with the molecular conformation defined by the C-N glycosidic torsion angle.     

Simulation of the different biological activities of diethylstilbestrol (DES) on estrogen receptor alpha and estrogen-related receptor gamma

We describe the application of the molecular dynamics (MD) and molecular mechanics-generalized Born/surface area (MM-GB/SA) approaches to the simulation of the different biological activity of diethylstilbestrol (DES) on two highly homologous nuclear receptors-estrogen receptor alpha (ER-alpha) and estrogen-related receptor gamma (ERR-gamma). DES exerts an agonistic effect against ER-alpha and an antagonistic effect against ERR-gamma. Using the x-ray crystal structures of ER-alpha in the canonical agonist bound form (PDB code: 3ERD) and antagonist bound form (PDB code: 3ERT), ERR-gamma homology models have been constructed for the receptor in two different conformations. MM-GB/SA binding free energy calculations of DES in the ER-alpha and ERR-gamma structures suggest that DES exhibits a greater free energy of binding in the agonist bound conformation of ER-alpha, while the antagonist bound conformation is preferred for ERR-gamma. Further dissection of the free energy contributions coupled with calculation of the ligand binding pocket volume suggests that the van der Waals interactions for DES within the smaller binding pocket volume of ERR-gamma are less favorable and this is the main factor for DES antagonism in ERR-gamma. This approach has potential general applicability to the prediction of the biological activity of nuclear receptor ligands.     

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.