Shape and size of simple cations in aqueous solutions: A theoretical reexamination of the hydrated ion via computer simulations

The simplest representation of monoatomic cations in aqueous solutions by means of a sphere with a radius chosen on the basis of a well-defined property (that of the bare ion or its hydrate) is reexamined considering classical molecular dynamics simulations. Two charged sphere–water interaction potentials were employed to mimic the bare and hydrated cation in a sample of 512 water molecules. Short-range interactions of trivalent cations were described by Lennard-Jones potentials which were fitted from ab initio calculations. Five statistically independent runs of 150 ps for each of the trivalent spheres in water were carried out in the microcanonical ensemble. A comparison of structural and dynamical properties of these simple ion models in solution with those of a system containing the Cr3+ hydrate ([Cr(H2O)6]3+) is made to get insight into the size and shape definition of simple ions in water, especially those that are highly charged. Advantages and shortcomings of using simple spherical approaches are discussed on the basis of reference calculations performed with a more rigorous hydrated ion model [J. Phys. Chem. B 102, 3272 (1998)]. The importance of nonspherical shape for the hydrate of highly charged ions is stressed and it is paradoxically shown that when spherical shape is retained, the big sphere representing the hydrate leads to results of ionic solution worse than those obtained with the small sphere. A low-cost method to generate hydrated ion–water interaction potentials taking into account the shape of the ionic aggregate is proposed.Dirección General de Investigación Científica y Técnica PB95-054

Quantum mechanical polar surface area

A correlation has been established between the absorbed fraction of training-set molecules after oral administration in humans and the Quantum Mechanical Polar Surface Area (QMPSA). This correlation holds for the QMPSA calculated with structures where carboxyl groups are deprotonated. The correlation of the absorbed fraction and the QMPSA calculated on the neutral gas phase optimized structures is much less pronounced. This suggests that the absorption process is mainly determined by polar interactions of the drug molecules in water solution. Rules are given to derive the optimal polar/apolar ranges of the electrostatic potential

High seroprevalence of human herpesviruses in HIV-infected individuals attending primary healthcare facilities in rural South Africa

Seroprevalence data of human herpesviruses (HHVs) are limited for sub-Saharan Africa. These are important to provide an indication of potential burden of HHV-related disease, in particular in human immunodeficiency virus (HIV)-infected individuals who are known to be at increased risk of these conditions in the Western world. In this cross-sectional study among 405 HIV-infected and antiretroviral therapy naïve individuals in rural South Africa the seroprevalence of HHVs was: herpes simplex virus type 1 (HSV-1) (98%), herpes simplex virus type 2 (HSV-2) (87%), varicella zoster virus (VZV) (89%), and 100% for both Epstein-Barr virus (EBV) and cytomegalovirus (CMV). Independent factors associated with VZV seropositivity were low educational status and having children. Lack of in-house access to drinking water was independently associated with positive HSV-1 serostatus, whereas Shangaan ethnicity was associated with HSV-2 seropositivity. Increasing age was associated with higher IgG titres to both EBV and CMV, whereas CD4 cell count was negatively associated with EBV and CMV IgG titres. Moreover, IgG titres of HSV-1 and 2, VZV and CMV, and CMV and EBV were positively correlated. The high HHV seroprevalence emphasises the importance of awareness of these viral infections in HIV-infected individuals in South Africa

Triazole Fungicides Can Induce Cross-Resistance to Medical Triazoles in Aspergillus fumigatus

Contains fulltext : 103858.pdf (publisher's version ) (Open Access)BACKGROUND: Azoles play an important role in the management of Aspergillus diseases. Azole resistance is an emerging global problem in Aspergillus fumigatus, and may develop through patient therapy. In addition, an environmental route of resistance development has been suggested through exposure to 14alpha-demethylase inhibitors (DMIs). The main resistance mechanism associated with this putative fungicide-driven route is a combination of alterations in the Cyp51A-gene (TR(34)/L98H). We investigated if TR(34)/L98H could have developed through exposure to DMIs. METHODS AND FINDINGS: Thirty-one compounds that have been authorized for use as fungicides, herbicides, herbicide safeners and plant growth regulators in The Netherlands between 1970 and 2005, were investigated for cross-resistance to medical triazoles. Furthermore, CYP51-protein homology modeling and molecule alignment studies were performed to identify similarity in molecule structure and docking modes. Five triazole DMIs, propiconazole, bromuconazole, tebuconazole, epoxiconazole and difenoconazole, showed very similar molecule structures to the medical triazoles and adopted similar poses while docking the protein. These DMIs also showed the greatest cross-resistance and, importantly, were authorized for use between 1990 and 1996, directly preceding the recovery of the first clinical TR(34)/L98H isolate in 1998. Through microsatellite genotyping of TR(34)/L98H isolates we were able to calculate that the first isolate would have arisen in 1997, confirming the results of the abovementioned experiments. Finally, we performed induction experiments to investigate if TR(34)/L98H could be induced under laboratory conditions. One isolate evolved from two copies of the tandem repeat to three, indicating that fungicide pressure can indeed result in these genomic changes. CONCLUSIONS: Our findings support a fungicide-driven route of TR(34)/L98H development in A. fumigatus. Similar molecule structure characteristics of five triazole DMIs and the three medical triazoles appear the underlying mechanism of cross resistance development. Our findings have major implications for the assessment of health risks associated with the use of triazole DMIs

The Dynamical Mechanism of Auto-Inhibition of AMP-Activated Protein Kinase

We use a novel normal mode analysis of an elastic network model drawn from configurations generated during microsecond all-atom molecular dynamics simulations to analyze the mechanism of auto-inhibition of AMP-activated protein kinase (AMPK). A recent X-ray and mutagenesis experiment (Chen, et al Nature 2009, 459, 1146) of the AMPK homolog S. Pombe sucrose non-fermenting 1 (SNF1) has proposed a new conformational switch model involving the movement of the kinase domain (KD) between an inactive unphosphorylated open state and an active or semi-active phosphorylated closed state, mediated by the autoinhibitory domain (AID), and a similar mutagenesis study showed that rat AMPK has the same auto-inhibition mechanism. However, there is no direct dynamical evidence to support this model and it is not clear whether other functionally important local structural components are equally inhibited. By using the same SNF1 KD-AID fragment as that used in experiment, we show that AID inhibits the catalytic function by restraining the KD into an unproductive open conformation, thereby limiting local structural rearrangements, while mutations that disrupt the interactions between the KD and AID allow for both the local structural rearrangement and global interlobe conformational transition. Our calculations further show that the AID also greatly impacts the structuring and mobility of the activation loop

ATP and its N6-substituted analogues: parameterization, molecular dynamics simulation and conformational analysis

In this work we used a combination of classical molecular dynamics and simulated annealing techniques to shed more light on the conformational flexibility of 12 adenosine triphosphate (ATP) analogues in a water environment. We present simulations in AMBER force field for ATP and 12 published analogues [Shah et al. (1997) Proc Natl Acad Sci USA 94: 3565–3570]. The calculations were carried out using the generalized Born (GB) solvation model in the presence of the cation Mg2+. The ion was placed at a close distance (2 Å) from the charged oxygen atoms of the beta and gamma phosphate groups of the −3 negatively charged ATP analogue molecules. Analysis of the results revealed the distribution of inter-proton distances H8–H1′ and H8–H2′ versus the torsion angle ψ (C4–N9-C1′–O4′) for all conformations of ATP analogues. There are two gaps in the distribution of torsion angle ψ values: the first is between −30 and 30 degrees and is described by cis-conformation; and the second is between 90 and 175 degrees, which mostly covers a region of anti conformation. Our results compare favorably with results obtained in experimental assays [Jiang and Mao (2002) Polyhedron 21:435–438]

Young elephants in a large herd maintain high levels of elephant endotheliotropic herpesvirus-specific antibodies and do not succumb to fatal haemorrhagic disease

Elephant endotheliotropic herpesviruses (EEHVs) have co-existed with elephants for millions of years, yet may cause fatal haemorrhagic disease (EEHV-HD), typically in elephants between 1 and 10 years of age. EEHV is omnipresent in (sub)adult elephants, and young elephants with low EEHV-specific antibody levels are at risk for EEHV-HD, suggesting that fatal disease may occur due to an insufficiently controlled primary infection. To further address this hypothesis, sera of three large elephant cohorts were subjected to a multiple EEHV species ELISA: (I) 96 Asian elephants between 0 and 57 years, including 13 EEHV-HD fatalities, from European zoo herds typically sized five to six elephants, (II) a herd of 64 orphaned elephants aged 0–15 years at the Elephant Transit Home in Sri Lanka and (III) 31 elephants aged 8–63 years, part of a large herd of 93 elephants at Pinnawala Elephant Orphanage, Sri Lanka. All Sri Lankan elephants showed high EEHV-specific antibody levels regardless of their age. While antibody levels of most European zoo elephants were comparable to those of Sri Lankan elephants, the average antibody level of the European juveniles (1–5 years of age) was significantly lower than those of age-matched Sri Lankan individuals. Moreover, the European juveniles showed a gradual decrease between 1 and 4 years of age, to be attributed to waning maternal antibodies. Maintenance of high levels of antibodies in spite of waning maternal antibodies in young Sri Lankan elephants is likely due to the larger herd size that increases the likelihood of contact with EEHV-shedding elephants. Together with the observation that low levels of EEHV-specific antibodies correlate with increased numbers of EEHV-HD fatalities, these results suggest that infection in presence of high maternal antibody levels may protect calves from developing EEHV-HD, while at the same time activating an immune response protective in future encounters with this virus

Regulation of Bestrophins by Ca2+: A Theoretical and Experimental Study

Bestrophins are a recently discovered family of Cl− channels, for which no structural information is available. Some family members are activated by increased intracellular Ca2+ concentration. Bestrophins feature a well conserved Asp-rich tract in their COOH terminus (Asp-rich domain), which is homologous to Ca2+-binding motifs in human thrombospondins and in human big-conductance Ca2+- and voltage-gated K+ channels (BKCa). Consequently, the Asp-rich domain is also a candidate for Ca2+ binding in bestrophins. Based on these considerations, we constructed homology models of human bestrophin-1 (Best1) Asp-rich domain using human thrombospondin-1 X-ray structure as a template. Molecular dynamics simulations were used to identify Asp and Glu residues binding Ca2+ and to predict the effects of their mutations to alanine. We then proceeded to test selected mutations in the Asp-rich domain of the highly homologous mouse bestrophin-2. The mutants expressed in HEK-293 cells were investigated by electrophysiological experiments using the whole-cell voltage-clamp technique. Based on our molecular modeling results, we predicted that Asp-rich domain has two defined binding sites and that D301A and D304A mutations may impact the binding of the metal ions. The experiments confirmed that these mutations do actually affect the function of the protein causing a large decrease in the Ca2+-activated Cl− current, fully consistent with our predictions. In addition, other studied mutations (E306A, D312A) did not decrease Ca2+-activated Cl− current in agreement with modeling results

Regulation of Bestrophins by Ca2+: A Theoretical and Experimental Study

Bestrophins are a recently discovered family of Cl− channels, for which no structural information is available. Some family members are activated by increased intracellular Ca2+ concentration. Bestrophins feature a well conserved Asp-rich tract in their COOH terminus (Asp-rich domain), which is homologous to Ca2+-binding motifs in human thrombospondins and in human big-conductance Ca2+- and voltage-gated K+ channels (BKCa). Consequently, the Asp-rich domain is also a candidate for Ca2+ binding in bestrophins. Based on these considerations, we constructed homology models of human bestrophin-1 (Best1) Asp-rich domain using human thrombospondin-1 X-ray structure as a template. Molecular dynamics simulations were used to identify Asp and Glu residues binding Ca2+ and to predict the effects of their mutations to alanine. We then proceeded to test selected mutations in the Asp-rich domain of the highly homologous mouse bestrophin-2. The mutants expressed in HEK-293 cells were investigated by electrophysiological experiments using the whole-cell voltage-clamp technique. Based on our molecular modeling results, we predicted that Asp-rich domain has two defined binding sites and that D301A and D304A mutations may impact the binding of the metal ions. The experiments confirmed that these mutations do actually affect the function of the protein causing a large decrease in the Ca2+-activated Cl− current, fully consistent with our predictions. In addition, other studied mutations (E306A, D312A) did not decrease Ca2+-activated Cl− current in agreement with modeling results