Lemniscates and inequalities for the logarithmic capacities of continua

It is shown that if P(z) = z ⁿ + ? is a polynomial with connected lemniscate E(P) = {z: ¦P(z)¦ ≤ 1} and m critical points, then, for any n? m+1 points on the lemniscate E(P), there exists a continuum γ ? E(P) of logarithmic capacity cap γ ≤ 2^?1/n which contains these points and all zeros and critical points of the polynomial. As corollaries, estimates for continua of minimum capacity containing given points are obtained.     

Norms and perfect graphs

The weak Berge hypothesis states that a graph is perfect if and only if its complement is perfect. Previous proofs of this hypothesis have used combinatorial or polyhedral methods.In this paper, the concept of norms related to graphs is used to provide an alternative proof for the weak Berge hypothesis.This is a written account of an invited lecture delivered by the second author on occasion of the 12. Symposium on Operations Research, Passau, 9.–11. 9. 1987.     

Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.     

The Pharmacophore Hypothesis of Novel Inhibitors for Aurora A Kinase

A three-dimensional pharmacophore model was developed from a series of inhibitors of Aurora A kinase to discover new potent anti-cancer agents using the HypoGen module in the Catalyst software. The pharmacophore model was developed based on the structure of 20 currently available inhibitors, which were carefully selected from the literature. The best hypothesis （Hypo 1） was defined by four features： one hydrogen-bond donor and three hy- drophobic points, with the best correlation coefficient of 0.909, the lowest rms deviation of 1.563, and the highest cost difference of 99.075. The Hypo 1 was then validated by a test set consisting of 24 compounds and by a cross-validation of 95% confidence level through randomizing the data using the CatScramble program, which suggested that a predictive pharmacophore model had been successfully obtained.     

COSMO-RS: A novel view to physiological solvation and partition questions

Both, dielectric continuum solvation models as well as surface or group based methods using polarity and lipophilicity parameters have been proven to be useful tools for the analysis of solvation and partition questions. For the first time, COSMO-RS provides an integrated theory, which combines the aspects of continuum solvation and surface interactions, and which ends up with chemical potentials of molecules in almost arbitrary solvents and mixtures. Due to its sound theoretical basis, COSMO-RS does not only provide a new quantitative access to solvation and partition properties in well defined solvents, but it also opens a novel view and gives a better understanding of the general problem of solvation. Finally, this allows for a generalisation of COSMO-RS to sophisticatedphysiological partition problems involving as complex phases as blood, brain, or cell membranes. The use of COSMO-RS for drug discovery and design is demonstrated by applications to blood-brain partition coefficients, and water solubility.     

Calculation of affinities of peptides for proteins

Several methodologies were employed to calculate the Gibbs standard free energy of binding for a collection of protein-ligand complexes, where the ligand is a peptide and the protein is representative for various protein families. Almost 40 protein-ligand complexes were employed for a continuum approach, which considers the protein and the peptide at the atomic level, but includes solvent as a polarizable continuum. Five protein-ligand complexes were employed for an all-atom approach that relies on a combination of the double decoupling method with thermodynamic integration and molecular dynamics. These affinities were also computed by means of the linear interaction energy method. Although it generally proved rather difficult to predict the absolute free energies correctly, for some protein families the experimental ranking order was correctly reproduced by the continuum and all-atom approach. Considerable attention has also been given to correctly analyze the affinities of charged peptides, where it is required to judge the effect of one or more ions that are being decoupled in an all-atom approach to preserve electroneutrality. The various methods are further judged upon their merits.     

Some Comments on the Matter Wave-light Wave Hypothesis

From the de Broglie matter wave hypothesis and Planck’s energy quantization law, and assuming conservation of energy in the absorption of a photon and its consequent conversion to kinetic energy of motion by a material particle initially at rest, one can deduce a simple mathematical relationship between the wavelength λ (or frequency ν), of the photon absorbed by the particle at rest, and the resulting de Broglie matter wave length, λ_D, of the particle with kinetic energy of motion of mv²/2. The relationship so deduced, λ_D∝√λ, suggests that visible wavelengths of light, from about 4000 ?, in the violet, to beyond about 7000 ?, in the red, on absorption by an electron at rest, lead to material electron wavelengths, λ_D, of the order of the size of the electron transfer proteins seen in the photosynthetic reaction centers of photosynthesizing organisms, at about a size of 50–100 ?. In addition to understanding the mechanism of photosynthesis as a material wave mediated phenomenon, further areas of importance of the relations pointed out in this paper are in the design of experiments to gain a deeper understanding of the basic tenets of wave mechanics, and in the use of tunable lasers to probe various properties of material waves, and to precisely control their properties for applications including lithography.     

Solvent effects on 17O nuclear magnetic shielding: N-methylformamide in polar and apolar solutions

The solvent effect on ¹⁷O isotropic shielding was computed by different methods: the polarizable continuum model and a mixed approach, including a few real solvent molecules treated as the solute. The experimental data show that the behavior of protic and aprotic solvents is markedly different: we found that the continuum approach describes well the observed shielding at various dielectric constants for aprotic solvents, while the mixed procedure is needed when hydrogen bonds to the magnetic centre are present.Contribution to the Jacopo Tomasi Honorary Issue     

Effective atom volumes for implicit solvent models: comparison between Voronoi volumes and minimum fluctuation volumes

An essential element of implicit solvent models, such as the generalized Born method, is a knowledge of the volume associated with the individual atoms of the solute. Two approaches for determining atomic volumes for the generalized Born model are described; one is based on Voronoi polyhedra and the other, on minimizing the fluctuations in the overall volume of the solute. Volumes to be used with various parameter sets for protein and nucleic acids in the CHARMM force field are determined from a large set of known structures. The volumes resulting from the two different approaches are compared with respect to various parameters, including the size and solvent accessibility of the structures from which they are determined. The question of whether to include hydrogens in the atomic representation of the solute volume is examined. Copyright 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1857-1879, 2001     

Improving the efficiency and convergence of geometry optimization with the polarizable continuum model: new energy gradients and molecular surface tessellation

New equations are derived and implemented for efficient and accurate computation of solvation energy derivatives for the conductor-like polarizable continuum model (C-PCM) and the isotropic integral equation formalism polarizable continuum model (IEF-PCM). Two new molecular surface tessellation procedures GEPOL-RT and GEPOL-AS that generate near continuous potential energy surfaces are proposed for PCM geometry optimization. The combined use of these new techniques leads to efficient and convergent geometry optimizations with the PCMs.