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.

Scaling in few-body nuclear physics

The nuclear matter calculations with realistic nucleon-nucleon potentials present a general scaling between the nucleon-nucleus binding energy, the corresponding saturation density, and the triton binding energy. The Thomas-Efimov three-body effect implies in correlations among low-energy few-body and many-body observables. It is also well known that, by varying the short-range repulsion, keeping the two-nucleon information (deuteron and scattering) fixed, the four-nucleon and three-nucleon binding energies lie on a very narrow band known as a Tjon line. By looking for a universal scaling function connecting the proper scales of the few-body system with those of the many-body system, we suggest that the general nucleus-nucleon scaling mechanism is a manifestation of a universal few-body effect.

The few scales of nuclei and nuclear matter

The well-known correlations of low-energy three and four-nucleon observables with a typical three-nucleon scale (e.g., the Tjon line) is extended to light nuclei and nuclear matter. Evidence for the scaling between light nuclei binding energies and the triton one are pointed out. We argue that the saturation energy and density of nuclear matter are correlated to the triton binding energy. The available systematic nuclear matter calculations indicate a possible band structure representing these correlations. (c) 2006 Elsevier B.V. All rights reserved.

Three- and Four-Nucleon Bound States in Three Dimensions, Without PW Decomposition

A brief review of a three-dimensional (3D) numerical method to solve few-nucleon bound and scattering states, without the standard partial-wave (PW) decomposition, is presented. The approach is applied to three-and four-nucleon bound states, by considering the solutions of the corresponding Faddeev-Yakubovsky (FY) integral equations in momentum space. Realistic spin-isospin dependent 3D and PW formalism are presented for the alpha particle and the triton binding energies, with numerical results given in both schemes for comparison.

RELATIVISTIC 3-PARTICLE DYNAMICAL EQUATIONS .2. APPLICATION TO THE TRINUCLEON SYSTEM

We calculate the contribution of relativistic dynamics on the neutron-deutron scattering length and triton binding energy employing five sets trinucleon potential models and four types of three-dimensional relativistic three-body equations suggested in the preceding paper. The relativistic correction to binding energy may vary a lot and even change sign depending on the relativistic formulation employed. The deviations of these observables from those obtained in nonrelativistic models follow the general universal trend of deviations introduced by off- and on-shell variations of two- and three-nucleon potentials in a nonrelativistic model calculation. Consequently, it will be difficult to separate unambiguously the effect of off- and on-shell variations of two- and three-nucleon potentials on low-energy three-nucleon observables from the effect of relativistic dynamics. (C) 1994 Academic Press, Inc.

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.