Impact of calcium on N1 influenza neuraminidase dynamics and binding free energy

The highly pathogenic influenza strains H5N1 and H1N1 are currently treated with inhibitors of the viral surface protein neuraminidase (N1). Crystal structures of N1 indicate a conserved, high affinity calcium binding site located near the active site. The specific role of this calcium in the enzyme mechanism is unknown, though it has been shown to be important for enzymatic activity and thermostability. We report molecular dynamics (MD) simulations of calcium‐bound and calcium‐free N1 complexes with the inhibitor oseltamivir (marketed as the drug Tamiflu), independently using both the AMBER FF99SB and GROMOS96 force fields, to give structural insight into calcium stabilization of key framework residues. Y347, which demonstrates similar sampling patterns in the simulations of both force fields, is implicated as an important N1 residue that can “clamp” the ligand into a favorable binding pose. Free energy perturbation and thermodynamic integration calculations, using two different force fields, support the importance of Y347 and indicate a +3 to +5 kcal/mol change in the binding free energy of oseltamivir in the absence of calcium. With the important role of structure‐based drug design for neuraminidase inhibitors and the growing literature on emerging strains and subtypes, inclusion of this calcium for active site stability is particularly crucial for computational efforts such as homology modeling, virtual screening, and free energy methods. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Characterization of a clinical polymer‐drug conjugate using multiscale modeling

The molecular conformation of certain therapeutic agents has been shown to affect the ability to gain access to target cells, suggesting potential value in defining conformation of candidate molecules. This study explores how the shape and size of poly‐γ‐glutamyl‐glutamate paclitaxel (PGG‐PTX), an amphiphilic polymer‐drug with potential chemotherapeutic applications, can be systematically controlled by varying hydrophobic and hydrophilic entities. Eighteen different formulations of PGG‐PTX varying in three PTX loading fractions (f_PTX) of 0.18, 0.24, and 0.37 and six spatial arrangements of PTX (‘clusters’, ‘ends’, ‘even’, ‘middle’, ‘random’, and ‘side’) were explored. Molecular dynamics (MD) simulations of all‐atom (AA) models of PGG‐PTX were run until a statistical equilibrium was reached at 100 ns and then continued as coarse‐grained (CG) models until a statistical equilibrium was reached at an effective time of 800 ns. Circular dichroism spectroscopy was used to suggest initial modeling configurations. Results show that a PGG‐PTX molecule has a strong tendency to form coil shapes, regardless of the PTX loading fraction and spatial PTX arrangement, although globular shapes exist at f_PTX = 0.24. Also, less uniform PTX arrangements such as ‘ends’, ‘middle’, and ‘side’ produce coil geometries with more curvature. The prominence of coil shapes over globules suggests that PGG‐PTX may confer a long circulation half‐life and high propensity for accumulation to tumor endothelia. This multiscale modeling approach may be advantageous for the design of cancer therapeutic delivery systems. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 936–951, 2010.     

Germ band retraction as a landmark in glucose metabolism during <Emphasis Type="Italic">Aedes aegypti</Emphasis> embryogenesis

Background  

The mosquito A. aegypti is vector of dengue and other viruses. New methods of vector control are needed and can be achieved by a better understanding of the life cycle of this insect. Embryogenesis is a part of A. aegypty life cycle that is poorly understood. In insects in general and in mosquitoes in particular energetic metabolism is well studied during oogenesis, when the oocyte exhibits fast growth, accumulating carbohydrates, lipids and proteins that will meet the regulatory and metabolic needs of the developing embryo. On the other hand, events related with energetic metabolism during A. aegypti embryogenesis are unknown.     

大鼠心肌缺血再灌注损伤模型的实验研究

目的 制备一种新型的心肌急性缺血再灌注损伤模型,以探讨一种更符合临床实际需求的实验方法.方法 将20只雌性SD(Sprague-Dawley)大鼠随机分成2组(对照组、实验组),采用结扎主动脉根部引起心肌缺血5min再灌注30 min建立心肌急性缺血再灌注模型;通过应用透射电镜观察心肌细胞超微结构的改变,同时检测心肌组织匀浆丙二醛(Maleic Dialdehyde,MDA)含量、超氧化物歧化酶(Superoxide Dismutase,SOD)活力.结果 透射电镜下超微结构显示实验组较对照组明显加重了心肌组织结构和线粒体的损害;实验组心肌组织MDA明显高于对照组(P<0.01),而SOD明显低于对照组(P<0.01).结论 本实验成功建立了方法简便、易于操作、取材范围广泛的心肌缺血再灌注损伤模型,为心肌缺血再灌注损伤研究提供了一种更为可行的模型.     

Bio3d: an R package for the comparative analysis of protein structures

An automated procedure for the analysis of homologous protein structures has been developed. The method facilitates the characterization of internal conformational differences and inter-conformer relationships and provides a framework for the analysis of protein structural evolution. The method is implemented in bio3d, an R package for the exploratory analysis of structure and sequence data. AVAILABILITY: The bio3d package is distributed with full source code as a platform-independent R package under a GPL2 license from: http://mccammon.ucsd.edu/~bgrant/bio3d/     

Dynamical properties of fasciculin-2.

Fasciculin-2 (FAS2) is a potent protein inhibitor of the hydrolytic enzyme acetylcholinesterase. A 2-ns isobaric-isothermal ensemble molecular dynamics simulation of this toxin was performed to examine the dynamic structural properties which may play a role in this inhibition. Conformational fluctuations of the FAS2 protein were examined by a variety of techniques to identify flexible residues and determine their characteristic motion. The tips of the toxin "finger" loops and the turn connecting loops I and II were found to fluctuate, while the rest of the protein remained fairly rigid throughout the simulation. Finally, the structural fluctuations were compared to NMR data of fluctuations on a similar timescale in a related three-finger toxin. The molecular dynamics results were in good qualitative agreement with the experimental measurements. Proteins 1999;36:447-453.     

Ribosome motions modulate electrostatic properties

The electrostatic properties of the 70S ribosome of Thermus thermophilus were studied qualitatively by solving the Poisson-Boltzmann (PB) equation in aqueous solution and with physiological ionic strength. The electrostatic potential was calculated for conformations of the ribosome derived by recent normal mode analysis (Tama, F., et al. Proc Natl Acad Sci USA 2003 100, 9319-9323) of the ratchet-like reorganization that occurs during translocation (Frank, J.; Agrawal, R. K. Nature 2000 406, 318-322). To solve the PB equation, effective parameters (charges and radii), applicable to a highly charged backbone model of the ribosome, were developed. Regions of positive potential were found at the binding site of the elongation factors G and Tu, as well as where the release factors bind. Large positive potential areas are especially pronounced around the L11 and L6 proteins. The region around the L1 protein is also positively charged, supporting the idea that L1 may interact with the E-site tRNA during its release from the ribosome after translocation. Functional rearrangement of the ribosome leads to electrostatic changes which may help the translocation of the tRNAs during the elongation stage.     

Brownian dynamics simulations of ion atmospheres around polyalanine and B-DNA: effects of biomolecular dielectric

We have extended an earlier Brownian dynamics simulation algorithm for simulating the structural dynamics of ions around biomolecules to accommodate dielectric inhomogeneity. The electrostatic environment of a biomolecule immersed in water was obtained by numerically solving the Poisson equation with the biomolecule treated as a low dielectric region and the solvent treated as a high dielectric region. Instead of using the mean-field type approximations of ion interactions as in the Poisson-Boltzmann model, the ions were treated explicitly by allowing them to evolve dynamically under the electrostatic field of the biomolecule. This model thus accounts for ion-ion correlations and the finite-size effects of the ions. For a 13-residue alpha-helical polyalanine and a 12-base-pair bp B-form DNA, we found that the choice of the dielectric constant of the biomolecule has much larger effects on the mean ionic structure around the biomolecule than on the fluctuational and dynamical properties of the ions surrounding the biomolecule.     

The relaxed complex method: Accommodating receptor flexibility for drug design with an improved scoring scheme

An extension of the new computational methodology for drug design, the "relaxed complex" method (J.-H. Lin, A. L. Perryman, J. R. Schames, and J. A. McCammon, Journal of the American Chemical Society, 2002, vol. 24, pp. 5632-5633), which accommodates receptor flexibility, is described. This relaxed complex method recognizes that ligand may bind to conformations that occur only rarely in the dynamics of the receptor. We have shown that the ligand-enzyme binding modes are very sensitive to the enzyme conformations, and our approach is capable of finding the best ligand enzyme complexes. Rapid docking serves as an efficient initial filtering method to screen a myriad of docking modes to a limited set, and it is then followed by more accurate scoring with the MM/PBSA (Molecular Mechanics/Poisson Boltzmann Surface Area) approach to find the best ligand-receptor complexes. The MM/PBSA scorings consistently indicate that the calculated binding modes that are most similar to those observed in the x-ray crystallographic complexes are the ones with the lowest free energies.     

The pH dependence of stability of the activation helix and the catalytic site of GART

We have predicted the free energy of unfolding for the pH-dependent helix-coil transition of the activation helix of GART using continuum electrostatic calculations and structural modeling. We have assigned the contributions of each element of secondary structure and of each ionizable residue, within and in the vicinity of the activation helix, to the stability of several fragments of GART that participate in the formation of the catalytic site. We demonstrate that the interaction of His121-His132 contributes 2.2 kcal/mol to the ionization free energy between pH 0 and approximately 6. We also show that the ionization state of a network of five histidines, His108, His119, His121, His132 and His137, and two aspartic acids Asp141 and Asp144, contributes approximately 12 kcal/mol to the stability of the catalytic site of GART, out of a total stability of 16 kcal/mol of the whole enzyme. These interactions are important for the formation of the catalytic site of GART.