Dataset Information

Independent-Trajectories Thermodynamic-Integration Free-Energy Changes for Biomolecular Systems: Determinants of H5N1 Avian Influenza Virus Neuraminidase Inhibition by Peramivir.

ABSTRACT: Free-energy changes are essential physicochemical quantities for understanding most biochemical processes. Yet, the application of accurate thermodynamic-integration (TI) computation to biological and macromolecular systems is limited by finite-sampling artifacts. In this paper, we employ independent-trajectories thermodynamic-integration (IT-TI) computation to estimate improved free-energy changes and their uncertainties for (bio)molecular systems. IT-TI aids sampling statistics of the thermodynamic macrostates for flexible associating partners by ensemble averaging of multiple, independent simulation trajectories. We study peramivir (PVR) inhibition of the H5N1 avian influenza virus neuraminidase flexible receptor (N1). Binding site loops 150 and 119 are highly mobile, as revealed by N1-PVR 20-ns molecular dynamics. Due to such heterogeneous sampling, standard TI binding free-energy estimates span a rather large free-energy range, from a 19% underestimation to a 29% overestimation of the experimental reference value (-62.2 +/- 1.8 kJ mol(-1)). Remarkably, our IT-TI binding free-energy estimate (-61.1 +/- 5.4 kJ mol(-1)) agrees with a 2% relative difference. In addition, IT-TI runs provide a statistics-based free-energy uncertainty for the process of interest. Using approximately 800 ns of overall sampling, we investigate N1-PVR binding determinants by IT-TI alchemical modifications of PVR moieties. These results emphasize the dominant electrostatic contribution, particularly through the N1 E277-PVR guanidinium interaction. Future drug development may be also guided by properly tuning ligand flexibility and hydrophobicity. IT-TI will allow estimation of relative free energies for systems of increasing size, with improved reliability by employing large-scale distributed computing.

SUBMITTER: Lawrenz M

PROVIDER: S-EPMC2669418 | biostudies-literature | 2009 Apr

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Independent-Trajectories Thermodynamic-Integration Free-Energy Changes for Biomolecular Systems: Determinants of H5N1 Avian Influenza Virus Neuraminidase Inhibition by Peramivir.

Lawrenz Morgan M Baron Riccardo R McCammon J Andrew JA

Journal of chemical theory and computation 20090325 4

Free-energy changes are essential physicochemical quantities for understanding most biochemical processes. Yet, the application of accurate thermodynamic-integration (TI) computation to biological and macromolecular systems is limited by finite-sampling artifacts. In this paper, we employ independent-trajectories thermodynamic-integration (IT-TI) computation to estimate improved free-energy changes and their uncertainties for (bio)molecular systems. IT-TI aids sampling statistics of the thermody ...[more]

PMID: 19461872

Similar Datasets

Project description:The H7N9 virus attaches itself to the human cell receptor protein containing the polysaccharide that terminates with sialic acid. The mutation of neuraminidase at residue E119 has been explored experimentally. However, there is no adequate information on the substitution with E119V in peramivir at the intermolecular level. Therefore, a good knowledge of the interatomic interactions is a prerequisite in understanding its transmission mode and subsequent effective inhibitions of the sialic acid receptor cleavage by neuraminidase. Herein, we investigated the mechanism and dynamism on the susceptibility of the E119V mutation on the peramivir-neuraminidase complex relative to the wildtype complex at the intermolecular level. This study aims to investigate the impact of the 119V substitution on the neuraminidase-peramivir complex and unveil the residues responsible for the complex conformations. We employed molecular dynamic (MD) simulations and extensive post-MD analyses in the study. These extensive computational investigations were carried out on the wildtype and the E119V mutant complex of the protein for holistic insights in unveiling the effects of this mutation on the binding affinity and the conformational terrain of peramivir-neuraminidase E119V mutation. The calculated total binding energy (ΔGbind) for the peramivir wildtype is -49.09 ± 0.13 kcal/mol, while the E119V mutant is -58.55 ± 0.15 kcal/mol. The increase in binding energy (9.46 kcal/mol) is consistent with other post-MD analyses results, confirming that E119V substitution confers a higher degree of stability on the protein complex. This study promises to proffer contributory insight and additional knowledge that would enhance future drug designs and help in the fight targeted at controlling the avian influenza H7N9 virus. Therefore, we suggest that experimentalists collaborate with computational chemists for all investigations of this topic, as we have done in our previous studies.

Project description:BACKGROUND: The agro-ecology and poultry husbandry of the south Asian and south-east Asian countries share common features, however, with noticeable differences. Hence, the ecological determinants associated with risk of highly pathogenic avian influenza (HPAI-H5N1) outbreaks are expected to differ between Bangladesh and e.g., Thailand and Vietnam. The primary aim of the current study was to establish ecological determinants associated with the risk of HPAI-H5N1 outbreaks at subdistrict level in Bangladesh. The secondary aim was to explore the performance of two different statistical modeling approaches for unmeasured spatially correlated variation. METHODOLOGY/PRINCIPAL FINDINGS: An ecological study at subdistrict level in Bangladesh was performed with 138 subdistricts with HPAI-H5N1 outbreaks during 2007-2008, and 326 subdistricts with no outbreaks. The association between ecological determinants and HPAI-H5N1 outbreaks was examined using a generalized linear mixed model. Spatial clustering of the ecological data was modeled using 1) an intrinsic conditional autoregressive (ICAR) model at subdistrict level considering their first order neighbors, and 2) a multilevel (ML) model with subdistricts nested within districts. Ecological determinants significantly associated with risk of HPAI-H5N1 outbreaks at subdistrict level were migratory birds' staging areas, river network, household density, literacy rate, poultry density, live bird markets, and highway network. Predictive risk maps were derived based on the resulting models. The resulting models indicate that the ML model absorbed some of the covariate effect of the ICAR model because of the neighbor structure implied in the two different models. CONCLUSIONS/SIGNIFICANCE: The study identified a new set of ecological determinants related to river networks, migratory birds' staging areas and literacy rate in addition to already known risk factors, and clarified that the generalized concept of free grazing duck and duck-rice cultivation interacted ecology are not significant determinants for Bangladesh. These findings will refine current understanding of the HPAI-H5N1 epidemiology in Bangladesh.

Project description:BackgroundOngoing outbreaks of H5N1 highly pathogenic avian influenza (HPAI) viruses and the emergence of the genetic-related hemagglutinin (HA) gene of reassortant H5Nx viruses currently circulating in wild birds and poultries pose a great global public health concern. In this study, we comprehensively analyzed the genetic evolution of Thai H5N1 HA and neuraminidase (NA) genes between 2003 and 2010. The H5N1 Thailand virus clade 2.3.4 was also genetically compared to the currently circulating clade 2.3.4.4 of H5Nx viruses.MethodsFull-length nucleotide sequences of 178 HA and 143 NA genes of H5N1 viruses circulating between 2003 and 2010 were phylogenetically analyzed using maximum likelihood (ML) phylogenetic construction. Bayesian phylogenetic trees were reconstructed using BEAST analysis with a Bayesian Markov chain Monte Carlo (MCMC) approach. The maximum clade credibility (MCC) tree was determined, and the time of the most recent common ancestor (tMRCA) was estimated. The H5N1 HA nucleotide sequences of clade 2.3.4 Thailand viruses were phylogenetically analyzed using ML phylogenetic tree construction and analyzed for nucleotide similarities with various subtypes of reassortant H5Nx HA clade 2.3.4.4.ResultsML phylogenetic analysis revealed two distinct HA clades, clade 1 and clade 2.3.4, and two distinct NA groups within the corresponding H5 clade 1 viruses. Bayesian phylogenetic reconstruction for molecular clock suggested that the Thai H5N1 HA and NA emerged in 2001.87 (95% HPD: 2001.34-2002.49) and 2002.38 (95% HPD: 2001.99-2002.82), respectively, suggesting that the virus existed before it was first reported in 2004. The Thai H5N1 HA clade 2.3.4 was grouped into corresponding clades 2.3.4, 2.3.4.1, 2.3.4.2, and 2.3.4.3, and shared nucleotide similarities to reassortant H5Nx clade 2.3.4.4 ranged from 92.4-96.8%. Phylogenetic analysis revealed monophyletic H5Nx clade 2.3.4.4 evolved from H5N1 clade 2.3.4.ConclusionH5N1 viruses existed, and were presumably introduced and circulated in avian species in Thailand, before they were officially reported in 2004. HA and NA genes continuously evolved during circulation between 2004 and 2010. This study provides a better understanding of genetic evolution with respect to molecular epidemiology. Monitoring and surveillance of emerging variants/reassortants should be continued.

Dataset Information

Independent-Trajectories Thermodynamic-Integration Free-Energy Changes for Biomolecular Systems: Determinants of H5N1 Avian Influenza Virus Neuraminidase Inhibition by Peramivir.

Publications

Independent-Trajectories Thermodynamic-Integration Free-Energy Changes for Biomolecular Systems: Determinants of H5N1 Avian Influenza Virus Neuraminidase Inhibition by Peramivir.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets