Ion channels in icosahedral virus: a comparative analysis of the structures and binding sites at their fivefold axes.

An analysis of the crystallographically determined structures of the icosahedral protein coats of Tomato Bushy Stunt Virus, Southern Bean Mosaic Virus, Satellite Tobacco Necrosis Virus, Human Rhinovirus 14 and Mengovirus around their fivefold axes is presented. Accessibilities surfaces, electrostatic energy profile calculations, ion-protein interaction energy calculations, free energy perturbation methods and comparisons with structures of chelating agents are used in this study. It is concluded that the structures built around the viral fivefold axes would be adequate for ion binding and transport. Relative ion preferences are derived for the binding sites, using free energy perturbation methods, which are consistent with the experimental data when available. In the cases where crystallographic studies determined the existence of ions on the fivefold axes, our results indicate that they would correspond to ions in crystallization or purification buffers. The environment of the fivefold axes are rich in polar residues in all icosahedral viral structures whose atomic coordinates are available, including some that are not being analyzed in detail in this work. The fivefold channel-like structures have most of the basic properties expected for real ion channels including a funnel at the entrance, a polar internal environment with frequent alternation of acidic and basic residues, ion binding sites, the capability to induce ion dehydration and ion transit from the external viral surface to the binding sites.     

Benchmark of different charges for prediction of the partitioning coefficient through the hydrophilic/lipophilic index

A few different theoretical methods for assigning the partial atomic charges were benchmarked for calculation of the hydrophilic/lipophilic index (HLI). The coefficients were selected to produce the best correlation of the HLI values with the experimental octanol-water partition. Different parameters were checked in calculations of partial charges to get the best performance of the HLI values obtained. Thus, four partitioning schemes (Coulson, Mulliken, Merz-Kollman, Ford-Wang) were benchmarked for calculations of atomic charges with six semiempirical methods (AM1, PM3, RM1, PM6, PM6-D3H4, PM7). Moreover, five distinct types of partial atomic charges (Mulliken, Hirshfeld, Löwdin, CHELPG, NPA), obtained at the Hartree–Fock and DFT levels of theory with three basis sets, were tested for their ability to produce the HLI values with the best correlation to experimental logP coefficients of 50 mono-charged organic anions. In the case of the semiempirical methods, the best correlation between the HLI and logP values (the correlation coefficient r?=?0.9216) was obtained with the AM1 Ford–Wang parametric electrostatic potential charges. The Mulliken and Coulson charges calculated with the PM7 method can be used as an alternative to AM1, with the r values of 0.9107 and 0.8984, respectively. In the case of the DFT, the PBE/def2-TZVP natural population analysis charges produce the best correlation (r?=?0.9220). Nevertheless, in spite of a marginally lower performance (r?=?0.9159), the NPA charges computed at the PBE/def2-SVP level are more robust and can be regarded as the optimum choice for calculating the HLI values.

Open image in new window

Graphical abstract The hydrophilic/lipophilic index (HLI)

     

On contribution of known atomic partial charges of protein backbone in electrostatic potential density maps

Partial charges of atoms in a molecule and electrostatic potential (ESP) density for that molecule are known to bear a strong correlation. In order to generate a set of point‐field force field parameters for molecular dynamics, Kollman and coworkers have extracted atomic partial charges for each of all 20 amino acids using restrained partial charge‐fitting procedures from theoretical ESP density obtained from condensed‐state quantum mechanics. The magnitude of atomic partial charges for neutral peptide backbone they have obtained is similar to that of partial atomic charges for ionized carboxylate side chain atoms. In this study, the effect of these known atomic partial charges on ESP is examined using computer simulations and compared with the experimental ESP density recently obtained for proteins using electron microscopy. It is found that the observed ESP density maps are most consistent with the simulations that include atomic partial charges of protein backbone. Therefore, atomic partial charges are integral part of atomic properties in protein molecules and should be included in model refinement.     

Calculations of the charge distribution in dodecyltrimethylammonium: a quantum chemical investigation

In this work we report the atomic partial charges evaluated on dodecyltrimethylammonium ion. The values obtained from 17 quantum methods [CHELP, CHELPG, MK, NPA at (HF, LDA, PBE, B3LYP)//6-31G++(d,p) level and APT at B3LYP//6-31G++(d,p)] on the molecule optimised at B3LYP/6-31G++(d,p) level were compared to two semiempirical methods (Gasteiger and Qeq) and the commercial force field PCFF. All methods based on quantum calculation gave a positive charge delocalised on at least the first four alkyl groups of the tail. However, those deriving partial charges from the electrostatic potential gave an unrealistic set of alternative positive and negative alkyl group charges along the tail. In comparison, the NPA and APT methods lead to a steady decrease in the partial charges from the third alkyl group, and agreed closely with the representation of the electrostatic potential mapped onto the 0.002 au isodensity surface. The choice of the exchange correlation treatment does not drastically influence the atomic partial charges.     

A model binding site for testing scoring functions in molecular docking

Prediction of interaction energies between ligands and their receptors remains a major challenge for structure-based inhibitor discovery. Much effort has been devoted to developing scoring schemes that can successfully rank the affinities of a diverse set of possible ligands to a binding site for which the structure is known. To test these scoring functions, well-characterized experimental systems can be very useful. Here, mutation-created binding sites in T4 lysozyme were used to investigate how the quality of atomic charges and solvation energies affects molecular docking. Atomic charges and solvation energies were calculated for 172,118 molecules in the Available Chemicals Directory using a semi-empirical quantum mechanical approach by the program AMSOL. The database was first screened against the apolar cavity site created by the mutation Leu99Ala (L99A). Compared to the electronegativity-based charges that are widely used, the new charges and desolvation energies improved ranking of known apolar ligands, and better distinguished them from more polar isosteres that are not observed to bind. To investigate whether the new charges had predictive value, the non-polar residue Met102, which forms part of the binding site, was changed to the polar residue glutamine. The structure of the resulting Leu99Ala and Met102Gln double mutant of T4 lysozyme (L99A/M102Q) was determined and the docking calculation was repeated for the new site. Seven representative polar molecules that preferentially docked to the polar versus the apolar binding site were tested experimentally. All seven bind to the polar cavity (L99A/M102Q) but do not detectably bind to the apolar cavity (L99A). Five ligand-bound structures of L99A/M102Q were determined by X-ray crystallography. Docking predictions corresponded to the crystallographic results to within 0.4A RMSD. Improved treatment of partial atomic charges and desolvation energies in database docking appears feasible and leads to better distinction of true ligands. Simple model binding sites, such as L99A and its more polar variants, may find broad use in the development and testing of docking algorithms.     

The study on the interaction between seryl-histidine dipeptide and proteins by circular dichroism and molecular modeling

The selective cleavage of proteins is very important in key biological processes. Chemical (nonenzymatic) reagents such as cyanogen bromide and transition metal complexes are used extensively with great defects. In this paper, the binding of seryl-histidine dipeptide (abbreviated as SH) with bovine serum albumin (BSA) and lysozyme were investigated by the circular dichroism spectroscopy (CD) at 298K, molecular docking studies and quantum chemical calculations based on the previous results of polyacrylamide gel electrophoresis (PAGE). From the studies of CD, it showed that SH interacted strongly with BSA and lysozyme. The change percentages of the secondary structures of BSA and lysozyme were calculated. The contents of the beta-sheets decreased remarkably. It indicated that the interactions between SH and proteins could break the hydrogen bonds of beta-sheets selectively. The docking studies between SH and BSA showed that the position of the oxygen atom of the hydroxyl group of SH (O(12)) was in favor of a nucleophilic attack on carbon atom of the amide bond of a beta-sheet (C(34)) because the distance between O(12) and C(34) was 3.37A. Natural charges, natural atomic hybrid percentages and square sums of HOMO coefficients calculated by the NBO and population analysis at HF/6-31G* supported the suggested mechanism. And so SH may be an interesting agent for the therapeutic use.     

Large variation in electrostatic contours upon addition of steric parameters and the effect of charge calculation schemes in CoMFA on mutagenicity of MX analogues

Partial atomic charge is a useful concept to describe physico-chemical properties of a molecule. For this, various schemes have been devised to get reasonable values. Mutagen X is an ideal set to test the effect of partial atomic charge variation. Therefore, we collected data from previous reports and studied various charge schemes. Our systematic study covers 26 charge calculation schemes along with a broad range of levels of theory. Charge calculation schemes include charges from charge equalisation, electrostatic potential fitting, molecular orbital and atomic polar tensor. Calculation levels span from empirical, semi-empirical, Hartree–Fock, density functional and Møller–Plesset 2. We also used two validation statistics for internal prediction. To observe the electrostatic effect accurately during comparative molecular field analysis (CoMFA) modelling, we first studied isolated electrostatic parameters to avoid interaction effect with steric parameters. The results clearly show that adding steric parameters does change statistical conclusions as well as CoMFA maps. Although there was a weak trend that quantum mechanical (QM)-derived charges gave better statistical values, it is not apparent statistically (alpha = 0.05). Particularly, Mülliken population analysis (MPA) did not produce better results. Therefore, when we excluded MPA schemes from QM calculation, the QM-derived charges were found to be significant, i.e. sophisticated charge schemes other than MPA with QM methods were found to be superior to simple empirical charge schemes. In addition, we demonstrated that in order to test charge schemes properly, excluding steric parameter is more important. This work exemplifies Occam's theorem of parsimony. A simpler model is a better model.     

Partial atomic charges of amino acids in proteins

Using a semiempirical quantum mechanical procedure (FCPAC) we have calculated the partial atomic charges of amino acids from 494 high-resolution protein structures. To analyze the influence of the protein's environment, we considered each residue under two conditions: either as the center of a tripeptide with PDB structure geometry (free) or as the center of 13-16 amino acid clusters extracted from the PDB structure (buried). The partial atomic charges from residues in helices and in sheets were separated. The FCPAC partial atomic charges of the Cbeta and Calpha of most residues correlate with their helix propensity, positively for Cbeta and negatively for Calpha (r2 = 0.76 and 0.6, respectively). The main consequence of burying residues in proteins is the polarization of the backbone C=O bond, which is more pronounced in helices than in sheets. The average shift of the oxygen partial charges that results from burying is -0.120 in helix and -0.084 in sheet with the charge of the proton as unit. Linear correlations are found between the average NMR chemical shifts and the average FCPAC partial charges of Calpha (r2 = 0.8-0.85), N (r3 = 0.67-0.72), and Cbeta (r2 = 0.62) atoms. Correlations for helix and beta-sheet FCPAC partial charges show parallel regressions, suggesting that the charge variations due to burying in proteins differentiate between the dihedral angle effects and the polarization of backbone atoms.     

Spectroscopic investigation and hydrogen-bonding analysis of triazinones

NIR FT-Raman, FTIR and UV-vis spectra of the herbicide metamitron were recorded and analyzed. The aromaticities, equilibrium geometries, bonding features, electrostatic potentials, and harmonic vibrational wavenumbers of the monomers and dimers of triazinone derivatives were also investigated with the aid of BLYP/6-311 G(df, p) density functional theory. Features in the vibrational spectra were assigned with the aid of the VEDA.4 program. The calculated results were a good match to the experimental data obtained from FTIR, Raman, and electronic absorption spectra. Mulliken population analysis was performed on the atomic charges and the HOMO-LUMO energies were also calculated. NBO analysis highlighted the intra- and intermolecular N-H…O and C-H…O hydrogen bonds in the crystal structures of the triazinones. The solvent effect was calculated using time-dependent density functional theory in combination with the polarizable continuum model.     

Structural organization of DNA in chlorella viruses

Chlorella viruses have icosahedral capsids with an internal membrane enclosing their large dsDNA genomes and associated proteins. Their genomes are packaged in the particles with a predicted DNA density of ca. 0.2 bp nm(-3). Occasionally infection of an algal cell by an individual particle fails and the viral DNA is dynamically ejected from the capsid. This shows that the release of the DNA generates a force, which can aid in the transfer of the genome into the host in a successful infection. Imaging of ejected viral DNA indicates that it is intimately associated with proteins in a periodic fashion. The bulk of the protein particles detected by atomic force microscopy have a size of ～60 kDa and two proteins (A278L and A282L) of about this size are among 6 basic putative DNA binding proteins found in a proteomic analysis of DNA binding proteins packaged in the virion. A combination of fluorescence images of ejected DNA and a bioinformatics analysis of the DNA reveal periodic patterns in the viral DNA. The periodic distribution of GC rich regions in the genome provides potential binding sites for basic proteins. This DNA/protein aggregation could be responsible for the periodic concentration of fluorescently labeled DNA observed in ejected viral DNA. Collectively the data indicate that the large chlorella viruses have a DNA packaging strategy that differs from bacteriophages; it involves proteins and share similarities to that of chromatin structure in eukaryotes.     

Binding site characteristics of 17beta-estradiol imprinted polymers

The variety of applications utilizing molecularly imprinted polymers (MIPs) requires synthetic strategies yielding different MIP formats including films, irregular particles, or spheres, along with precise knowledge on the specific material characteristics, such as binding capacity and binding efficiency of these materials. In response to this demand, MIPs are prepared in different formats by variation of the polymerization methodology. It is commonly agreed that micro- and sub-microspheres are particularly advantageous MIP formats, due to their monodispersity and facile synthesis procedures in contrast to conventional imprinted polymers prepared by bulk polymerization. However, the differences in actual rebinding characteristics of different MIP formats based on molecular interactions under a variety of binding/rebinding conditions have not been studied in detail to date. Consequently, the present work details an analytical strategy generically applicable to MIP systems for rebinding studies including equilibrium binding, non-equilibrium binding, and release experiments enabling more profound understanding on the molecular interactions between the imprinted materials and the template molecules. In this study, three MIP formats were considered for the same template molecule, 17beta-estradiol: irregularly shaped particulate polymers prepared by bulk polymerization and grinding, microspheres, and sub-microspheres. The latter two formats were synthesized via precipitation polymerization using different processing strategies. The morphologies and porosities of the resulting imprinted materials were characterized by scanning electron microscopy (SEM) and Brunauer-Emmett-Teller (BET) analysis, respectively. The obtained results indicate that microspheres prepared by precipitation polymerization provide superior rebinding properties during equilibrium binding in contrast to bulk polymers and sub-microspheres, and that the rebinding properties are different during equilibrium binding versus non-equilibrium binding. The median binding affinity constant determined during non-equilibrium rebinding is higher than the values obtained from equilibrium rebinding. Furthermore, the binding site distribution appears more homogeneous thief derived from non-equilibrium rebinding, as reflected in a heterogeneity index of m=0.725. Moreover, it is hypothesized that the specific interactions between template and monomers are related to the porosity of the imprinted polymers, which implies that the amount of binding sites and the pore sized distribution of the imprinted materials are a critical factor in achieving the desired MIP performance in various analytical applications. The BET results indicate that particles prepared with lower cross-linker-to-template ratio have a reduced surface area. Furthermore, it can be expected that there are less specific binding sites available at particles with reduced surface area and pore volume given similar distribution of the binding sites, as confirmed by the equilibrium binding isotherm studies. The pore size distribution results reveal that control of the pore size in the range of 100-180 A is essential to obtain the desired retention properties and Gaussian peak shape during HPLC analysis of small molecules.     

Beyond CHELP: improved potential derived charges for sugars

The partitioning of the overall molecular charge distribution into atom centered monopole charges, while quantum mechanically ill-defined, is nevertheless a technique which finds applications in several broad classes of chemical problems. Charges derived from fits to electrostatic potentials have an intuitive appeal since, in principle, these could be derived from either theoretical or experimental data. It has been noted, however, that such potential derived charges can be conformationally dependent in ways that do not appear to reflect the changes in the molecular wavefunction. Both the algorithm used for selecting points at which the molecular electrostatic potential will be fit and the density of points used in the fit have been suggested to influence the resultant charges. Recently [Stouch TR, Williams DE (1992) J Comp Chem 13: 622–32; Stouch TR, Williams DE (1993) J Comp Chem 14: 858–66] it has been noted that numerical difficulties may make it impossible to fit all the atomic charges in a molecule. Singular value decomposition (SVD) of the linear least squares matrices used in fitting atom based monopoles to molecular electrostatic potentials provides a tool for evaluating the integrity of the calculated charges. Based on the SVD analysis for a selected group of molecules we have noted particularly that increasing the molecular size reduces the fraction of charges which can be validly assigned. Users of PD derived charges, especially those who are using those charges for tasks other than reproduction of the MEP, should be aware that there is a high probability that a significant portion of those charges are statistically unreliable. Therefore, charges in many biological molecules, such as sugars, prove to be difficult to obtain by potential derived (PD) methods such as CHELP or CHELPG. Results from the SVD can be used to both assess PD charges and to generate an improved, albeit incomplete, set. Improved PD fits are presented for a series of simple saccharides. Abbreviations: HF, Hartree-Fock; LLS, linear least squares; MEP, molecular electrostatic potential; PD, potential derived; SVD, singular value decomposition This revised version was published online in November 2006 with corrections to the Cover Date.     

Theoretical analysis of deoxycytidine-5'-phosphate protonation

The electron density distribution in deoxycytidine-5'-monophosphate (5'-dCMP) molecule and dianion has been studied by the method of CNDO/2. The comparison between the results of calculation for the neutral molecule and the data obtained by Pearlman and Kim shows that there is a linear correlation between the atomic charges calculated using quantum chemistry and those derived from X-ray results. However, partial charges for the deoxyribose fragment are correlated in a nonlinear manner. The influence of the protons added to the cytosine and phosphate residues on the atomic charges and bond orders of deoxy-cytidine-5'-monophosphate has been analyzed here. The conclusion has been drawn that the semiempirical quantum-chemical CNDO/2 technique is applicable to the mononucleotide studies.     

An improved generalized AMBER force field (GAFF) for urea

We describe an improved force field parameter set for the generalized AMBER force field (GAFF) for urea. Quantum chemical computations were used to obtain geometrical and energetic parameters of urea dimers and larger oligomers using AM1 semiempirical MO theory, density functional theory at the B3LYP/6-31G(d,p) level, MP2 and CCSD ab initio calculations with the 6-311++G(d,p), aug-cc-pVDZ, aug-cc-pVTZ, and aug-cc-pVQZ basis sets, and with the CBS-QB3 and CBS-APNO complete basis set methods. Seven different urea dimer structures were optimized at the MP2/aug-cc-pVDZ level to obtain accurate interaction energies. Atomic partial charges were calculated at the MP2/aug-cc-pVDZ level with the restrained electrostatic potential (RESP) fitting approach. The interaction energies computed with these new RESP charges in the force field are consistent with those obtained from CCSD and MP2 calculations. The linear dimer structure calculated using the force field with modified geometrical parameters and the new RESP charge set agrees well with available experimental data.     

Electronic structure and stabilities of Ni-doped germanium nanoclusters: a density functional modeling study

The present study reports the geometry, electronic structure, growth behavior and stability of neutral and ionized nickel encapsulated germanium clusters containing 1–20 germanium atoms within the framework of a linear combination of atomic orbital density functional theory (DFT) under a spin polarized generalized gradient approximation. In the growth pattern, Ni-capped Ge_n and Ni-encapsulated Ge_n clusters appear mostly as theoretical ground state at a particular size. To explain the relative stability of the ground state clusters, variation of different parameters, such as average binding energy per atom (BE), embedding energy (EE) and fragmentation energy (FE) of the clusters, were studied together with the size of the cluster. To explain the chemical stability of the clusters, different parameters, e.g., energy gap between the highest occupied and lowest unoccupied molecular orbitals (HOMO–LUMO gap), ionization energy (IP), electron affinity (EA), chemical potential (μ), chemical hardness (η), and polarizability etc. were calculated and are discussed. Finally, natural bond orbital (NBO) analysis was applied to understand the electron counting rule applied in the most stable Ge₁₀Ni cluster. The importance of the calculated results in the design of Ge-based superatoms is discussed.

Atomic charges for DNA constituents derived from single-crystal X-ray diffraction data

We have derived a complete set of atomic charges for DNA from very high resolution, low temperature, single-crystal X-ray diffraction data, collected for a variety of nucleosides and nucleotides: cytidine; deoxycytidine 5'-monophosphate; deoxythymidine; guanosine 5'-monophosphate; deoxyadenosine; adenosine. This set of charges represents the first experimentally based parameterization of an important term in the energy function used in most modeling of DNA. The resulting charges are in good agreement with chemical intuition and experimental observations. They also agree qualitatively with the theoretically derived values now commonly used, but numerous and significant quantitative differences are observed. Possible reasons for the quantitative disagreement are discussed. An averaged set of charges (derived from the experimental results), which can be used in DNA modeling calculations, is presented.     

The Influence of Charge Calculation on Molecular Dynamics Simulation of Adenine in Water

Abstract

Molecular dynamics simulations of an aqueous solution of adenine have been performed using different methods of charge calculation to evaluate the influence of the values of the atomic charges on the dynamical results and to incorporate new information about the interaction between adenine and water. Four sets of partial charges where computed using ab-initio methods. In all cases the hydration properties of adenine were similar. These results support the view that the simulations by molecular dynamics, at least for the regime of infinite dilution, are not sensitive with respect to the different sets of partial charges used. A net hydrophobic behavior of the adenine molecule, on the water was observed.     

Characterization of surface binding sites in glycoside hydrolases: A computational study

Structural properties of carbohydrate surface binding sites (SBSs) were investigated with computational methods. Eighty‐five SBSs of 44 enzymes in 119 Protein Data Bank (PDB) files were collected as a dataset. On the basis of SBSs shape, they were divided into 3 categories: flat surfaces, clefts, and cavities (types A, B, and C, respectively). Ligand varieties showed the correlation between shape of SBSs and ligands size. To reduce cut‐off differences in each SBSs with different ligand size, molecular docking were performed. Molecular docking results were used to refine SBSs classification and binding sites cut‐off. Docking results predicted putative ligands positions and displayed dependence of the ligands binding mode to the structural type of SBSs. Physicochemical properties of SBSs were calculated for all docking results with YASARA Structure. The results showed that all SBSs are hydrophilic, while their charges could vary and depended to ligand size and defined cut‐off. Surface binding sites type B had highest average values of solvent accessible surface area. Analysis of interactions showed that hydrophobic interactions occur more than hydrogen bonds, which is related to the presence of aromatic residues and carbohydrates interactions.     

Electronic structure simulations of 2,6-dimethyl-2,5-heptadien-4-one by FTIR,FT-Raman,NMR, UV–vis,NBO and density functional theory

Density functional theory (DFT) (B3LYP and B3PW91) calculations have been carried out for 2,6-dimethyl-2,5-heptadien-4-one (DMHD4O) using 6–311++ G** basis set. Complete vibrational assignment and analysis of the fundamental modes of the compound were carried out from the FTIR and FT-Raman spectral data. The theoretical electronic absorption has been calculated by using time-dependent DFT (TD-DFT) methods and compared with the experimental spectra. The theoretically computed Frontier energy gaps and TD-DFT calculations are in good agreement with the experimental UV–vis spectral absorption. The chemical hardness measured from the Frontier molecular orbital energies of DMHD4O is 0.0693 eV. Electronic stability of the compound arising from hyperconjugative interactions and charge delocalisation were also investigated based on the natural bond orbital (NBO) analysis. Effective stabilisation energy E ⁽²⁾ associated with the interactions of the π and the lone pair of electrons was determined by the NBO analysis. ¹³C and ¹H NMR chemical shifts of the compound have been calculated by means of Gauge-Invariant Atomic Orbital using B3LYP/6–311++ G** method. The partial ionic character of the carbonyl group due to resonance render a partially positive charge to the carbonyl carbon, and thus C4 chemical shift lie in the very downfield 191.6 ppm. Comparison between the experimental and the theoretical results indicates that B3LYP method is able to provide satisfactory results for predicting vibrational, electronic and NMR properties.     

A free energy based computational pathway from chemical templates to lead compounds: a case study of COX-2 inhibitors

Automation of lead compound design in silico given the structure of the protein target and a definition of its active site vies for the top of the wish list in any drug discovery programme. We present here an enumeration of steps starting from chemical templates and propose a solution at the state of the art, in the form of a system independent comprehensive computational pathway. This methodology is illustrated with cyclooxygenase-2 (COX-2) as a target. We built candidate molecules including a few Non Steroidal Anti-inflammatory Drugs (NSAIDs) from chemical templates, passed them through empirical filters to assess drug-like properties, optimized their geometries, derived partial atomic charges via quantum calculations, performed Monte Carlo docking, carried out molecular mechanics and developed free energy estimates with Molecular Mechanics Generalized Born Solvent Accessibility (MMGBSA) methodology for each of the candidate molecules. For the case of aspirin, we also conducted molecular dynamics on the enzyme, the drug and the complex with explicit solvent followed by binding free energy analysis. Collectively, the results obtained from the above studies viz. sorting of drugs from non-drugs, semi-quantitative estimates of binding free energies, amply demonstrate the viability of the strategy proposed for lead selection/design for biomolecular targets.