A New Multi-Objective Approach for Molecular Docking Based on RMSD and Binding Energy

Ligand-protein docking is an optimization problem based on predicting the position of a ligand with the lowest binding energy in the active site of the receptor. Molecular docking problems are traditionally tackled with single-objective, as well as with multi-objective approaches, to minimize the binding energy. In this paper, we propose a novel multi-objective formulation that considers: the Root Mean Square Deviation (RMSD) difference in the coordinates of ligands and the binding (intermolecular) energy, as two objectives to evaluate the quality of the ligand-protein interactions. To determine the kind of Pareto front approximations that can be obtained, we have selected a set of representative multi-objective algorithms such as NSGA-II, SMPSO, GDE3, and MOEA/D. Their performances have been assessed by applying two main quality indicators intended to measure convergence and diversity of the fronts. In addition, a comparison with LGA, a reference single-objective evolutionary algorithm for molecular docking (AutoDock) is carried out. In general, SMPSO shows the best overall results in terms of energy and RMSD (value lower than 2A for successful docking results). This new multi-objective approach shows an improvement over the ligand-protein docking predictions that could be promising in in silico docking studies to select new anticancer compounds for therapeutic targets that are multidrug resistant.

Energy bands and acceptor binding energies of GaN

The energy bands of zinc-blende and wurtzite GaN are calculated with the empirical pseudopotential method, and the pseudopotential parameters for Ga and N atoms are-given. The calculated energy bands are in agreement with those obtained by the ab initio method. The effective-mass theory for the semiconductors of wurtzite structure is established, and the effective-mass parameters of GaN for both structures are given The binding energies of acceptor states are calculated by solving strictly the effective-mass equations. The binding energies of donor and acceptor are 24 and 142 meV for the zinc-blende structure, 20 and 131, and 97 meV for the wurtzite structure, respectively, which are consistent with recent experimental results. It is proposed that there are two kinds of acceptor in wurtzite GaN. One kind is the general acceptor such as C, which substitutes N, which satisfies the effective-mass theory. The other kind of acceptor includes Mg, Zn, Cd, etc., the binding energy of these accepters is deviated from that given by the effective mass theory. In this report, wurtzite GaN is grown by the molecular-beam epitaxy method, and the photoluminescence spectra were measured. Three main peaks are assigned to the donor-acceptor transitions from two kinds of accepters. Some of the transitions were identified as coming from the cubic phase of GaN, which appears randomly within the predominantly hexagonal material. [S0163-1829(99)15915-0].

Energy landscape theory, funnels, specificity, and optimal criterion of biomolecular binding

We study the nature of biomolecular binding. We found that in general there exists several thermodynamic phases: a native binding phase, a non-native phase, and a glass or local trapping phase. The quantitative optimal criterion for the binding specificity is found to be the maximization of the ratio of the binding transition temperature versus the trapping transition temperature, or equivalently the ratio of the energy gap of binding between the native state and the average non-native states versus the dispersion or variance of the non-native states. This leads to a funneled binding energy landscape.

Binding properties of peptidic affinity ligands for plasmid DNA capture and detection

Peptides constructed from α-helical subunits of the Lac repressor protein (LacI) were designed then tailored to achieve particular binding kinetics and dissociation constants for plasmid DNA purification and detection. Surface plasmon resonance was employed for quantification and characterization of the binding of double stranded Escherichia coli plasmid DNA (pUC19) via the lac operon (lacO) to "biomimics" of the DNA binding domain of LacI. Equilibrium dissociation constants (K D), association (k a), and dissociation rates (k d) for the interaction between a suite of peptide sequences and pUC19 were determined. K D values measured for the binding of pUC19 to the 47mer, 27mer, 16mer, and 14mer peptides were 8.8 ± 1.3 × 10 -10 M, 7.2 ± 0.6 × 10 -10 M, 4.5 ± 0.5 × 10 -8 M, and 6.2 ± 0.9 × 10 -6 M, respectively. These findings show that affinity peptides, composed of subunits from a naturally occurring operon-repressor interaction, can be designed to achieve binding characteristics suitable for affinity chromatography and biosensor devices.

Binding of nucleobases with single-walled carbon nanotubes: Theory and experiment

We report the binding energy of various nucleobases (guanine (G), adenine (A), thymine (T) and cytosine (C)) with (5,5) single-walled carbon nanotube (SWNT) calculated using first-principle Hartre–Fock method (HF) together with classical force field. The binding energy without including the solvation effects of water decreases in the order G>A>T>C. The inclusion of solvation energy changes the order of binding preference to be G>T>A>C. Using isothermal titration (micro) calorimetry experiments, we also show the relative binding affinity to be T>A>C, in agreement with our calculations.

Clustering and solute-vacancy binding energies in Al-4.4% Zn alloy with 0.01% ternary additions

Isochronal and isothermal ageing experiments have been carried out to determine the influence of 0.01 at. % addition of a second solute on the clustering rate in the quenched Al-4,4 a/o Zn alloy. The influence of quenching and ageing temperatures has been interpreted to obtain the apparent vacancy formation and vacancy migration energies in the various ternary alloys. Using a vacancy-aided clustering model the following values of binding free energy have been evaluated: Ce-0.18; Dy-0.24; Fe-0.18; Li-0.25; Mn-0.27; Nb-0.18; Pt-0.23; Sb-0.21; Si-0.30; Y-0.25; and Yb-0.23 (± 0.02 eV). These binding energy values refer to that between a solute atom and a single vacancy. The values of vacancy migration energy (c. 0.4 eV) and the experimental activation energy for solute diffusion (c. 1.1 eV) are unaffected by the presence of the ternary atoms in the Al-Zn alloy.

Modes of binding of guanosine monophosphates to ribonuclease T1 - A computer modelling study

Computer-modelling studies on the modes of binding of the three guanosine monophosphate inhibitors 2'-GMP, 3'-GMP, and 5'-GMP to ribonuclease (RNase) T1 have been carried out by energy minimization in Cartesian-coordinate space. The inhibitory power was found to decrease in the order 2'-GMP > 3'-GMP > 5'-GMP in agreement with the experimental observations. The ribose moiety was found to form hydrogen bonds with the protein in all the enzyme-inhibitor complexes, indicating that it contributes to the binding energy and does not merely act as a spacer between the base and the phosphate moieties as suggested earlier. 2'-GMP and 5'-GMP bind to RNase T1 in either of the two ribose puckered forms (with C3'-endo more favoured over the C2'-endo) and 3'-GMP binds to RNase T1 predominantly in C3'-endo form. The catalytically important residue His-92 was found to form hydrogen bond with the phosphate moiety in all the enzyme-inhibitor complexes, indicating that this residue may serve as a general acid group during catalysis. Such an interaction was not found in either X-ray or two-dimensional NMR studies.

Crystal Structure Engineering by Fine-Tuning the Surface Energy: The Case of CdE (E = S/Se) Nanocrystals

We prove that CdS nanocrystals can be thermodynamically stabilized in both wurtzite and zinc-blende crystallographic phases at will, just by the proper choice of the capping ligand. As a striking demonstration of this, the largest CdS nanocrystals (similar to 15 nm diameter) ever formed with the zinc-blende structure have been synthesized at a high reaction temperature of 310 degrees C, in contrast to previous reports suggesting the formation of zinc-blende CdS only in the small size limit (< 4.5 nm) or at a lower reaction temperature (<= 240 degrees C). Theoretical analysis establishes that the binding energy of trioctylphosphine molecules on the (001) surface of zinc-blende CdS is significantly larger than that for any of the wurtzite planes. Consequently, trioctylphosphine as a capping agent stabilizes the zinc-blende phase via influencing the surface energy that plays an important role in the overall energetics of a nanocrystal. Besides achieving giant zinc-blende CdS nanocrystals, this new understanding allows us to prepare CdSe and CdSe/CdS core/shell nanocrystals in the zinc-blende structure.

Excitonic energy of vertically stacked self-assmbled InAs quantum dots

We have studied the exciton states in vertically stacked self-assembled quantum disks within the effective mass approximation. The energy spectrum of the electron and hole is calculated using the transfer matrix formalism in the adiabatic approximation. The Coulomb interaction between the electron and the hole is treated accurately by the direct diagonalization of the Hamiltonian matrix. The effect of the vertical alignment of the disks on the ground energy of heavy- and light-hole exciton is presented and discussed. The binding energy is discussed in terms of the probability of the ground wave function. The ground energy of heavy- and light-hole excitons as a function of the magnetic field is presented and the effect of the disk size (the radius of disks) on the exciton energy is discussed.

GdxSi grown with mass-analyzed low energy dual ion beam epitaxy technique

Semiconducting gadolinium silicide GdxSi samples were prepared by mass-analyzed low-energy dual ion beam epitaxy technique. Auger electron spectroscopy depth profiles indicate that the gadolinium ions are implanted into the single-crystal silicon substrate and formed 20 nm thick GdxSi film. X-ray double-crystal diffraction measurement shows that there is no new phase formed. The XPS spectra show that one type of silicon peaks whose binding energy is between that of silicide and silicon dioxide, and the gadolinium peak of binding energy is between that of metal Gd and Gd2O3. All of these results indicate that an amorphous semiconductor is formed. (C) 2002 Elsevier Science B.V. All rights reserved.

Potential energy surface and moleculardynamics of MbNO: existence of an unsuspected FeON minimum

Ligands such as CO, O2, or NO are involved in the biological function of myoglobin. Here we investigate the energetics and dynamics of NO interacting with the Fe(II) heme group in native myoglobin using ab initio and molecular dynamics simulations. At the global minimum of the ab initio potential energy surface (PES), the binding energy of 23.4 kcal/mol and the Fe-NO structure compare well with the experimental results. Interestingly, the PES is found to exhibit two minima: There exists a metastable, linear Fe-O-N minimum in addition to the known, bent Fe-N-O global minimum conformation. Moreover, the T-shaped configuration is found to be a saddle point, in contrast to the corresponding minimum for NO interacting with Fe(III). To use the ab initio results for finite temperature molecular dynamics simulations, an analytical function was fitted to represent the Fe-NO interaction. The simulations show that the secondary minimum is dynamically stable up to 250 K and has a lifetime of several hundred picoseconds at 300 K. The difference in the topology of the heme-NO PES from that assumed previously (one deep, single Fe-NO minimum) suggests that it is important to use the full PES for a quantitative understanding of this system. Why the metastable state has not been observed in the many spectroscopic studies of myoglobin interacting with NO is discussed, and possible approaches to finding it are outlined.

Dopant-vacancy binding effects in Li-doped magnesium hydride

We use a combination of ab initio calculations and statistical mechanics to investigate the substitution of Li+ for Mg2+ in magnesium hydride (MgH2) accompanied by the formation of hydrogen vacancies with positive charge (with respect to the original ion at the site). We show that the binding energy between dopants and vacancy defects leads to a significant fraction of trapped vacancies and therefore a dramatic reduction in the number of free vacancies available for diffusion. The concentration of free vacancies initially increases with dopant concentration but reaches a maximum at around 1 mol % Li doping and slowly decreases with further doping. At the optimal level of doping, the corresponding concentration of free vacancies is much higher than the equilibrium concentrations of charged and neutral vacancies in pure MgH2 at typical hydrogen storage conditions. We also show that Li-doped MgH2 is thermodynamically metastable with respect to phase separation into pure magnesium and lithium hydrides at any significant Li concentration, even after considering the stabilization provided by dopant-vacancy interactions and configurational entropic effects. Our results suggest that lithium doping may enhance hydrogen diffusion hydride but only to a limited extent determined by an optimal dopant concentration and conditioned to the stability of the doped phase.

Many-body energy decomposition of hydrogen-bonded glycine clusters in gas-phase

A detailed analysis of the many-body contribution to the interaction energies of the gas-phase hydrogen-bonded glycine clusters, (Gly)(N), N = 1-4 is presented. The energetics of the hydrogen-bonded dimer, trimer and tetramer complexes have been analyzed using density-functional theory. The magnitude of the two-through four-body energy terms have been calculated and compared. The relaxation energy and the two-body energy terms are the principal contributors to the total binding energy. Four-body contribution is negligible. However, the three-body contribution is found to be sizable and the formation of the cyclic glycine trimer presents geometric strains that make it less favorable. (C) 2010 Elsevier B.V. All rights reserved.

Design of multiligand inhibitors for the swine flu H1N1 neuraminidase binding site

Viral neuraminidase inhibitors such as oseltamivir and zanamivir prevent early virus multiplication by blocking sialic acid cleavage on host cells. These drugs are effective for the treatment of a variety of influenza subtypes, including swine flu (H1N1). The binding site for these drugs is well established and they were designed based on computational docking studies. We show here that some common natural products have moderate inhibitory activity for H1N1 neuraminidase under docking studies. Significantly, docking studies using AutoDock for biligand and triligand forms of these compounds (camphor, menthol, and methyl salicylate linked via methylene bridges) indicate that they may bind in combination with high affinity to the H1N1 neuraminidase active site. These results also indicate that chemically linked biligands and triligands of these natural products could provide a new class of drug leads for the prevention and treatment of influenza. This study also highlights the need for a multiligand docking algorithm to understand better the mode of action of natural products, wherein multiple active ingredients are present.

Dimensional Compactification and Two-Particle Binding

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)