Comparison of dielectric response models for simulating electrostatic effects in proteins

Two recent approaches for calculating pK shifts in proteins are compared. The first of these uses Coulomb's law with a distance-dependent dielectric permittivity, epsilon (r), to model the screening effects of the environment, and the second uses a finite difference approach to solve Poisson's equation. It is shown that an explicit form of epsilon (r) which has been fitted to experimentally determined values of the dielectric permittivity in a range from 1 to 21 A can be approximated by a linear form in the functionally significant range of charge separations of approximately 3-10 A, but for distances greater than 10 A the effective permittivity is strongly nonlinear. A statistical analysis of the errors in calculated pK shifts due to electrostatic interactions between charges with separations greater than 10 A shows that there are only marginal differences in reliability between using Coulomb's law with an appropriate form of epsilon (r) or the finite difference approach for solving Poisson's equation. Thus it is concluded that pK shifts can be calculated just as well, and with considerably less effort, using Coulomb's law.     

The molecular electrostatic potential and steric accessibility of A-DNA.

The molecular electrostatic potential and steric accessibility of A-DNA are computed for base sequences (dG.dC)n and (dA.dT)n. An interpretation of the results in terms of the structure of A-DNA is provided and differences with respect to other forms of DNA, namely B-DNA and Z-DNA, are discussed.     

The electrostatic potential and steric accessibility of reactive sites within Z-DNA.

A theoretical study of indices potentially useful for investigation of the reactivity of the recently discovered Z-DNA double helix is presented. The electrostatic potential minima and the steric accessibility of reactive sites are calculated. The effect of screening the phosphate groups by metal cations is investigated. The results are compared with those for the B-DNA double helix.     

Theoretical studies of enzymic reactions: Dielectric,electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme

A general method for detailed study of enzymic reactions is presented. The method considers the complete enzyme-substrate complex together with the surrounding solvent and evaluates all the different quantum mechanical and classical energy factors that can affect the reaction pathway. These factors include the quantum mechanical energies associated with bond cleavage and charge redistribution of the substrate and the classical energies of steric and electrostatic interactions between the substrate and the enzyme. The electrostatic polarization of the enzyme atoms and the orientation of the dipoles of the surrounding water molecules is simulated by a microscopic dielectric model. The solvation energy resulting from this polarization is considerable and must be included in any realistic calculation of chemical reactions involving anything more than an isolated molecule in vacuo. Without it, acidic groups can never become ionized and the charge distribution on the substrate will not be reasonable. The same dielectric model can also be used to study the reaction of the substrate in solution. In this way the reaction in solution can be compared with the enzymic reaction.In this paper we study the stability of the carbonium ion intermediate formed in the cleavage of a glycosidic bond by lysozyme. It is found that electrostatic stabilization is an important factor in increasing the rate of the reaction step that leads to the formation of the carbonium ion intermediate. Steric factors, such as the strain of the substrate on binding to lysozyme, do not seem to contribute significantly.     

Implications of electrostatic potentials on ribosomal proteins.

Potentiometric studies of ribosomal particles 30S, 50S, and 70S, were designed to investigate possible implications of the electrostatic potentials developed by the 16S and 23S rRNA fractions. Release of protons and proton titrations of these ribosomal fractions were examined as a function of Mg2+ and K+ concentrations. The effects of these cations fit the polyelectrolyte theory remarkably well and are discussed accordingly.     

Theoretical models of eumelanin protomolecules and their optical properties

The molecular structure of melanin, one of the most ubiquitous natural pigments in living organisms, is not known and its multifaceted biological role is still debated. We examine structural models for eumelanin protomolecules, based on tetramers consisting of four monomer units (hydroquinone, indolequinone, and its two tautomers), in arrangements that contain an interior porphyrin ring. These models reproduce convincingly many aspects of eumelanin's experimentally observed behavior. In particular, we present a plausible synthetic pathway of the tetramers and their further complexation through interlayer stacking, or through formation of helical superstructures, into eumelanin macromolecules. The unsaturated nature of C-C bonds in indolequinone units and the finite size of protomolecules introduce covalent bond formation between stacked layers. We employ time-dependent density functional theory to calculate the optical absorption spectrum of each molecule along the eumelanin synthesis pathway, which gradually develops into the characteristic broad-band adsorption of melanin pigment due to electron delocalization. These optical spectra may serve as signatures for identifying intermediate structures.     

Ab initio quantum mechanical models of peptide helices and their vibrational spectra

Structural parameters for standard peptide helices (alpha, 3(10), 3(1) left-handed) were fully ab initio optimized for Ac-(L-Ala)(9)-NHMe and for Ac-(L-Pro)(9)-NHMe (poly-L-proline-PLP I and PLP II-forms), in order to better understand the relative stability and minimum energy geometries of these conformers and the dependence of the ir absorption and vibrational CD (VCD) spectra on detailed variation in these conformations. Only the 3(10)-helical Ala-based conformation was stable in vacuum for this decaamide structure, but both Pro-based conformers minimized successfully. Inclusion of solvent effects, by use of the conductor-like screening solvent model (COSMO), enabled ab initio optimizations [at the DFT/B3LYP/SV(P) level] without any constraints for the alpha- and 3(10)-helical Ala-based peptides as well as the two Pro-based peptides. The geometries obtained compare well with peptide chain torsion angles and hydrogen-bond distances found for these secondary structure types in x-ray structures of peptides and proteins. For the simulation of VCD spectra, force field and intensity response tensors were obtained ab initio for the complete Ala-based peptides in vacuum, but constrained to the COSMO optimized torsional angles, due to limitations of the solvent model. Resultant spectral patterns reproduce well many aspects of the experimental spectra and capture the differences observed for these various helical types.     

The electrostatic molecular potential of tRNAPhe. IV. The potentials and steric accessibilities of sites associated with the bases.

The sites of the 76 nucleic acid bases of tRNAPhe potentially reactive towards electrophiles are studied by calculations on the associated molecular electrostatic potentials and the static steric accessibilities. Each of these sites is treated in its environment within the macromolecule. The influence of various schemes of screening by countercations of the backbone phosphates on the electrostatic potentials is investigated. The possible significance of the potentials and accessibilities in connection with observed chemical reactivities is discussed.     

Electrostatic effects in proteins: comparison of dielectric and charge models.

Two approaches for calculating electrostatic effects in proteins are compared and ana analysis is presented of the dependence of calculated properties on the model used to define the charge distribution. Changes in electrostatic free energy have been calculated using a screened Coulomb potential (SCP) with a distance-dependent effective dielectric permittivity to model bulk solvent effects and a finite difference approach to solve the Poisson-Boltzmann (FDPB) equation. The properties calculated include shifts in dissociation constants of ionizable groups, the effect of annihilating surface charges on the binding of metals, and shifts in redox potentials due to changes in the charge of ionizable groups. In the proteins considered the charged sites are separated by 3.5-12 A. It is shown that for the systems studied in this distance range the SCP yields calculated values which are at least as accurate as those obtained from solution of the FDPB equation. In addition, in the distance range 3-5 A the SCP gives substantially better results than the FDPB equation. Possible sources of this difference between the two methods are discussed. Shifts in binding constants and redox potentials were calculated with several standard charge sets, and the resulting values show a variation of 20-40% between the 'best' and 'worst' cases. From this study it is concluded that in most applications, changes in electrostatic free energies can be calculated economically and reliably using an SCP approach with a single functional form of the screening function.     

Multi-stage transitional models with random effects and their application to the Einstein aging study

Longitudinal studies of aging often gather repeated observations of cognitive status to describe the development of dementia and to assess the influence of risk factors. Clinical progression to dementia is often conceptualized by a multi-stage model of several transitions that synthesizes time-varying effects. In this study, we assess the influence of risk factors on the transitions among three cognitive status: cognitive stability (normal cognition for age), memory impairment, and clinical dementia. We have developed a shared random effects model that not only links the propensity of transitions and to the probability of informative missingness due to death, but also incorporates heterogeneous transition between subjects. We evaluate four approaches using generalized logit and four using proportional odds models to the first-order Markov transition probabilities as a function of covariates. Random effects were incorporated into these models to account for within-subject correlations. Data from the Einstein Aging Study are used to evaluate the goodness-of-fit of these models using the Akaike information criterion. The best fitting model for each type (generalized logit and proportional odds) is recommended and their results are discussed in more details.     

Somnogenic cytokines and models concerning their effects on sleep

All the sleep-promoting substances currently identified also have other biological activities. Despite years of effort, a single specific central nervous system sleep center has not been described. These observations led us to propose a biochemical model of a sleep activational system in which the effects of several sleep factors are integrated into a regulatory scheme. These sleep factors interact by altering the metabolism, production, or activity of each other and thereby result in multiple feedback loops. This web of interactions leads to sleep stability in that minor challenges to the system will not greatly alter sleep. The system, however, is responsive to strong perturbations, such as sleep deprivation and infectious disease. The sleep-promoting effects of cytokines and their interactions with prostaglandins and the neuroendocrine system are used to illustrate the functioning of a part of the sleep activational system under normal conditions and during infectious disease. Although the actions of individuals sleep factors are not specific to sleep, their interactions at various levels of the neuraxis can mediate a specific sleep response. Such a system would also be responsive to the autonomic and environmental parameters that alter sleep.     

Myristoylation of neutrophil proteins and their biological characteristics

We have studied protein acylation in neutrophils of guinea pigs using [3H]myristate. A large number of neutrophil proteins were acylated with exogenously added myristic acid. The myristoylation was detected on 110, 77, 56, 54, 52, 42, and 37 kDa proteins. These myristoylations were stronger in peripheral blood than in peritoneal cells. Myristic acid was found to be covalently linked by an amid bond to these proteins since the proteins were resistant to boiling, chloroform/methanol and hydroxylamine treatment. Most myristoylated proteins appeared to be associated with the membrane fraction, while some of the proteins such as 77 kDa one was distributed also in the cytoplasm and translocated from the cytoplasm to the plasma membrane by stimulation. Lysozyme was myristoylated in vitro by the N-hydroxysuccinimide ester of myristic acid. The myristoylated lysozyme had an ability to be associated with phospholipid liposomes, and the membrane-associated lysozyme became a substrate of the rat brain Ca2+- and phospholipid dependent protein kinase (protein kinase C). These results indicate that myristoylation in neutrophil proteins may have an important role in metabolic regulation through their membrane association.     

DNA binding properties of novel cytotoxic gold(III) complexes of terpyridine ligands: the impact of steric and electrostatic effects

Four gold(III) complexes of terpyridine derivatives 1–4 have been synthesized and characterized by spectroscopic methods. In vitro data demonstrated that all of them showed higher cytotoxicity than cisplatin against the human non-small-cell lung cancer cell line (A-549), the human stomach carcinoma cell line (SGC-7901), the human cervix carcinoma cell line (HELA), the human colon carcinoma cell line (HCT-116), the human liver carcinoma cell line (BEL-7402), the murine leukemia cell line (P-388) and the human acute promyelocytic leukemia cell line (HL-60). Complex 3 exhibits the highest activity, with growth inhibition rates of over 80% at 10⁻⁸ mol L⁻¹ against the A-549, HCT-116 and HELA tumor cell lines. Interestingly, ligands L1–L4 are also very cytotoxic against the cell lines tested. Complexes 1–4 are stable in aqueous solution for 2 days in the presence of the biological reducing agent glutathione. The inductively coupled plasma mass spectrometry data showed that DNA isolated from cells treated with complexes 1 and 3 contained gold with gold-to-nucleotide ratios of approximately 1:6,400 and 1:4,900, respectively. Fluorescence titration, UV and circular dichroism analyses proved that the steric and electrostatic effects of the ligand remarkably influence the interactions of their gold(III) complexes with DNA. The DNA binding ability of the complexes has been correlated with their cytotoxicity, which could potentially provide a new rationale for the future design of terpyridine-based metal complexes with antitumor potential.Electronic Supplementary Material Supplementary material is available to authorized users in the online version of this article at .     

Some considerations on vibrational mathematical models of the ulna

An approximate mechanical model of the ulna is proposed. It has realistic complicating factors: variable cross-sectional area and moment of inertia, and non-uniform mechanical properties.     

Erythrocyte membrane cholesterol levels and their effects on membrane proteins

The susceptibility of the band 3 protein of the erythrocyte membrane to proteolytic digestion at either surface of the membrane was not altered when the membrane cholesterol level was increased by 65–103%. Cross-linking of the major membrane proteins by o-phenanthroline · Cu, glutaraldehyde, dimethylsuberimidate and dimethyladipimidate was also unaffected.     

Vibrational Stark effects calibrate the sensitivity of vibrational probes for electric fields in proteins

Infrared spectroscopy is widely used to probe local environments and dynamics in proteins. The introduction of a unique vibration at a specific site of a protein or more complex assembly offers many advantages over observing the spectra of an unmodified protein. We have previously shown that infrared frequency shifts in proteins can arise from differences in the local electric field at the probe vibration. Thus, vibrational frequencies can be used to map electric fields in proteins at many sites or to measure the change in electric field due to a perturbation. The Stark tuning rate gives the sensitivity of a vibrational frequency to an electric field, and for it to be useful, the Stark tuning rate should be as large as possible. Vibrational Stark effect spectroscopy provides a direct measurement of the Stark tuning rate and allows a quantitative interpretation of frequency shifts. We present vibrational Stark spectra of several bond types, extending our work on nitriles and carbonyls and characterizing four additional bond types (carbon-fluorine, carbon-deuterium, azide, and nitro bonds) that are potential probes for electric fields in proteins. The measured Stark tuning rates, peak positions, and extinction coefficients provide the primary information needed to design amino acid analogues or labels to act as probes of local environments in proteins.     

A new theoretical index of biochemical reactivity combining steric and electrostatic factors. An application to yeast tRNAPhe

A new theoretical index of the chemical reactivity of sites within macromolecules is developed, which combines both steric and electrostatic factors. It is applied to the study of yeast tRNAPhe and the results obtained are compared with known experimental reactivities. A comparison indicates the superiority of the new index over the sole use of the surface accessibility.