PubMed Health⌕ Search

Biomedical subjects

Eugene Demchuk

Publications and source records attributed to Eugene Demchuk.

10 recordsLinked to original sources

Nearest-neighbor nonparametric method for estimating the configurational entropy of complex molecules.

A method for estimating the configurational (i.e., non-kinetic) part of the entropy of internal motion in complex molecules is introduced that does not assume any particular parametric form for the underlying probability density function. It is based on the nearest-neighbor (NN) distances of the points of a sample of internal molecular coordinates obtained by a computer simulation of a given molecule. As the method does not make any assumptions about the underlying potential energy function, it accounts fully for any anharmonicity of internal molecular motion. It provides an asymptotically unbiased and consistent estimate of the configurational part of the entropy of the internal degrees of freedom of the molecule. The NN method is illustrated by estimating the configurational entropy of internal rotation of capsaicin and two stereoisomers of tartaric acid, and by providing a much closer upper bound on the configurational entropy of internal rotation of a pentapeptide molecule than that obtained by the standard quasi-harmonic method. As a measure of dependence between any two internal molecular coordinates, a general coefficient of association based on the information-theoretic quantity of mutual information is proposed. Using NN estimates of this measure, statistical clustering procedures can be employed to group the coordinates into clusters of manageable dimensions and characterized by minimal dependence between coordinates belonging to different clusters.

Computer Simulation↗

Development of an empirical force field for silica. Application to the quartz-water interface.

Interactions of pulverized crystalline silica with biological systems, including the lungs, cause cell damage, inflammation, and apoptosis. To allow computational atomistic modeling of these pathogenic processes, including interactions between silica surfaces and biological molecules, new parameters for quartz, compatible with the CHARMM empirical force field were developed. Parameters were optimized to reproduce the experimental geometry of alpha-quartz, ab initio vibrational spectra, and interactions between model compounds and water. The newly developed force field was used to study interactions of water with two singular surfaces of alpha-quartz, (011) and (100). Properties monitored and analyzed include the variation of the density of water molecules in the plane perpendicular to the surface, disruption of the water H-bond network upon adsorption, and space-time correlations of water oxygen atoms in terms of Van Hove self-correlation functions. The vibrational density of states spectra of water in confined compartments were also computed and compared with experimental neutron-scattering results. Both the attenuation and shifting to higher frequencies of the hindered translational peaks upon confinement are clearly reproduced by the model. However, an upshift of librational peaks under the conditions of model confinement still remains underrepresented at the current empirical level.

Adsorption↗

Impact of silanol surface density on the toxicity of silica aerosols measured by erythrocyte haemolysis.

Exposures to silica-containing dusts are associated with a risk of developing life-threatening lung diseases. However, the mechanism of silica toxicity is poorly understood. In this work the atomic structure of the surfaces of different silica polymorphs was determined, and a relationship with in vitro silica toxicity was examined. The density of geminal and single silanol groups was quantitatively estimated for different silica polymorphs using a novel molecular modeling method. An association was found between the reported haemolytic activity and modeled densities of surface geminal (but not single) silanol groups on several silica polymorphs. These findings suggest a new view of aerosol toxicity based on the estimation of surface site densities. The results can be used in the development of new toxicological assays for respirable particulates, including nanomaterials.

Aerosols↗

Immunogenetic factors in beryllium sensitization and chronic beryllium disease.

Exposure to beryllium in the workplace can cause beryllium sensitization and chronic beryllium disease. Sensitization to beryllium can be detected in the laboratory using the beryllium lymphocyte proliferation test. It was shown that anti-HLA antibodies could block the beryllium-specific response in the beryllium lymphocyte proliferation test, thereby implicating HLA genes in chronic beryllium disease. A supratypic genetic marker, HLA-DPB1*E69, was found to be strongly associated with immunologic sensitization to beryllium and chronic beryllium disease in beryllium workers. However, there are 40 HLA-DPB1 gene variants that have E69 but that also have other DNA sequence variations. The purpose of the study was to evaluate the evidence for potential differential susceptibility that may be associated with the physical characteristics of HLA protein molecules for which different HLA-DPB1*E69 variants code; that is, do some HLA-DPB1*E69 variants convey higher risk of beryllium sensitization and chronic beryllium disease than others. To do this, two approaches were pursued: first, detailed analysis of the findings from the published literature was performed, and second, computational chemistry was used to seek clues concerning the physical properties of the HLA protein molecules for which these alleles code. Among the 40 HLA-DPB1 gene variants that code for E69, molecular epidemiological studies have suggested a risk hierarchy, where some variants appear to convey low to moderate risk of chronic beryllium disease (e.g., HLA-DPB1*0201, approximately 3-fold increased risk), some convey an intermediate risk (e.g., HLA-DPB1*1901, approximately 5-fold) and others convey high risk (e.g., HLA-DPB1*1701, >10-fold). Molecular modeling has been used to further investigate a potential mechanistic basis for these observations. We found a strong correlation between the hierarchical order of risk of chronic beryllium disease associated with specific alleles and the predicted surface electrostatic potential and charge of the corresponding isotypes. Therefore, when alleles were grouped by the relative negative charge on the molecules for which they code, the data suggest that those alleles associated with the most negatively charged proteins carry the greatest risk of beryllium sensitization and disease.

Berylliosis↗

A comparative study of unrelaxed surfaces on quartz and kaolinite, using the periodic density functional theory.

To investigate surface properties of fractured silica particles, which are commonly connected to the etiology of silica toxicity, models of low-index unrelaxed surfaces of quartz and kaolinite were constructed and analyzed using the periodic density functional theory calculations. The models were used to investigate surface sites that emerge in the processes of heterolytic and homolytic cleavage of quartz. It is found that the quartz surface is stabilized by two types of interactions. One, due to a more even charge distribution of sites, was characterized by surface energies of up to 0.025 eV x A(-2) and the other, due to a more even oxygen distribution between complementary surfaces, was up to 0.036 eV x A(-2). The total specific surface energies of unrelaxed surfaces ranged from 0.161 to 0.200 eV x A(-2) for quartz and from 0.017 to 0.158 eV x A(-2) for kaolinite. For the conchoidal fracture of quartz an average specific surface energy of 0.187 eV x A(-2) was obtained. These results provide a foundation for further characterization of the surface properties of mechanically comminuted respirable silica particulate and for reduction of occupational health hazards due to pulverized silica.

Chemistry, Physical↗

Estimation of the absolute internal-rotation entropy of molecules with two torsional degrees of freedom from stochastic simulations.

A method of statistical estimation is applied to the problem of evaluating the absolute entropy of internal rotation in a molecule with two torsional degrees of freedom. The configurational part of the entropy is obtained as that of the joint probability density of an arbitrary form represented by a two-dimensional Fourier series, the coefficients of which are statistically estimated using a sample of the torsional angles of the molecule obtained by a stochastic simulation. The internal rotors in the molecule are assumed to be attached to a common frame, and their reduced moments of inertia are initially calculated as functions of the two torsional angles, but averaged over all the remaining internal degrees of freedom using the stochastic-simulation sample of the atomic configurations of the molecule. The torsional-angle dependence of the reduced moments of inertia can be also averaged out, and the absolute internal-rotation entropy of the molecule is obtained in a good approximation as the sum of the configurational entropy and a kinetic contribution fully determined by the averaged reduced moments of inertia. The method is illustrated using Monte Carlo simulations of isomers of stilbene and halogenated derivatives of propane. The two torsional angles in cis-stilbene are found to be much more strongly correlated than those in trans-stilbene, while the degree of the angular correlation in propane increases strongly on substitution of hydrogen atoms with chlorine.

Journal Article↗

Structure-activity models for contact sensitization.

Allergic contact dermatitis (ACD) is a widespread cause of workers' disabilities. Although some substances found in the workplace are rigorously tested, the potential of the vast majority of chemicals to cause skin sensitization remains unknown. At the same time, exhaustive testing of all chemicals in workplaces is costly and raises ethical concerns. New approaches to developing information for risk assessment based on computational (quantitative) structure-activity relationship [(Q)SAR] methods may be complementary to and reduce the need for animal testing. Virtually any number of existing, de novo, and even preconceived compounds can be screened in silico at a fraction of the cost of animal testing. This work investigates the utility of ACD (Q)SAR modeling from the occupational health perspective using two leading software products, DEREK for Windows and TOPKAT, and an original method based on logistic regression methodology. It is found that the correct classification of (Q)SAR predictions for guinea pig data achieves values of 73.3, 82.9, and 87.6% for TOPKAT, DEREK for Windows, and the logistic regression model, respectively. The correct classification using LLNA data equals 73.0 and 83.2% for DEREK for Windows and the logistic regression model, respectively.

Allergens↗

Statistical thermodynamics of internal rotation in a hindering potential of mean force obtained from computer simulations.

A method of statistical estimation is applied to the problem of one-dimensional internal rotation in a hindering potential of mean force. The hindering potential, which may have a completely general shape, is expanded in a Fourier series, the coefficients of which are estimated by fitting an appropriate statistical-mechanical distribution to the random variable of internal rotation angle. The function of reduced moment of inertia of an internal rotation is averaged over the thermodynamic ensemble of atomic configurations of the molecule obtained in stochastic simulations. When quantum effects are not important, an accurate estimate of the absolute internal rotation entropy of a molecule with a single rotatable bond is obtained. When there is more than one rotatable bond, the "marginal" statistical-mechanical properties corresponding to a given internal rotational degree of freedom are reduced. The method is illustrated using Monte Carlo simulations of two public health relevant halocarbon molecules, each having a single internal-rotation degree of freedom, and a molecular dynamics simulation of an immunologically relevant polypeptide, in which several dihedral angles are analyzed.

Algorithms↗

Topographical organization of the N-terminal segment of lung pulmonary surfactant protein B (SP-B(1-25)) in phospholipid bilayers.

The location and depth of each residue of lung pulmonary surfactant protein B (SP-B(1-25)) in a phospholipid bilayer (PB) was determined by fluorescence quenching using synthesized single-residue-substituted peptides that were reconstituted into 1,2-dipalmitoyl phosphatidylcholine (DPPC)-enriched liposomes. The single-residue substitutions in peptides were either aspartate or tryptophan. The aspartate was subsequently labeled with the N-cyclohexyl-N'-(4-(dimethylamino)naphthyl)carbodiimide (NCD-4) fluorophore, whereas tryptophan is autofluorescent. Spin-labeled compounds, 5-doxylstearic acid (5-DSA), 7-doxylstearic acid (7-DSA), 12-doxylstearic acid (12-DSA), 4-(N,N-dimethyl-N-hexadecyl)ammonium-2,2,6,6-tetramethylpiperidine-1-oxyl iodide (CAT-16), and 4-trimethylammonium-2,2,6,6-tetramethylpiperidine-1-oxy iodide (CAT-1), were used in the quenching experiments. The effective quenching order is determined by the accessibility of the quencher to a fluorescent group on the peptide. The order of quenching efficiency provides information about the relative locations of individual residues in the PB. Our data indicate that residues Phe1-Pro6 are located at the surface of PB, residues Tyr7-Trp9 are embedded in PB, and residues Leu10-Ile22 are involved in an amphipathic alpha-helix with its axis parallel to the surface of PB; residues Pro23-Gly25 reside at the surface. The effects of intermolecular disulfide bond formation in the SP-B(1-25) dimer were also investigated. The experiments suggest that the SP-B helix A has to rotate at an angle to form a disulfide bond with the neighboring cysteine, which makes the hydrophobic sides of the amphipathic helices face each other, thus forming a hydrophobic domain. The detailed topographical mapping of SP-B(1-25) and its dimer in PB provides new insights into the conformational organization of the lung pulmonary surfactant proteins in the environment that mimics the native state. The environment-specific conformational flexibility of the hydrophobic domain created by SP-B folding may explain the key functional properties of SP-B including their impact on phospholipid transport between the lipid phases and in modulating the cell inflammatory response during respiratory distress syndrome.

Amino Acid Sequence↗

Electrostatic potential on human leukocyte antigen: implications for putative mechanism of chronic beryllium disease.

The pathobiology of chronic beryllium disease (CBD) involves the major histocompatibility complex class II human leukocyte antigen (HLA). Although occupational exposure to beryllium is the cause of CBD, molecular epidemiologic studies suggest that specific (Italic)HLA-DPB1(/Italic) alleles may be genetic susceptibility factors. We have studied three-dimensional structural models of HLA-DP proteins encoded by these genes. The extracellular domains of HLA-DPA1*0103/B1*1701, *1901, *0201, and *0401, and HLA-DPA1*0201/B1*1701, *1901, *0201, and *0401 were modeled from the X-ray coordinates of an HLA-DR template. Using these models, the electrostatic potential at the molecular surface of each HLA-DP was calculated and compared. These comparisons identify specific characteristics in the vicinity of the antigen-binding pocket that distinguish the different HLA-DP allotypes. Differences in electrostatics originate from the shape, specific disposition, and variation in the negatively charged groups around the pocket. The more negative the pocket potential, the greater the odds of developing CBD estimated from reported epidemiologic studies. Adverse impact is caused by charged substitutions in positions 55, 56, 69, 84, and 85, namely, the exact same loci identified as genetic markers of CBD susceptibility as well as cobalt-lung hard metal disease. These findings suggest that certain substitutions may promote an involuntary cation-binding site within a putatively metal-free peptide-binding pocket and therefore change the innate specificity of antigen recognition.

Amino Acid Sequence↗