X-ray structure and conformational dynamics of the HIV-1 protease in complex with the inhibitor SDZ283-910: agreement of time-resolved spectroscopy and molecular dynamics simulations

Based on the X-ray structure of the human immunodeficiency virus type-1 (HIV-1) protease in complex with the statine-derived inhibitor SDZ283-910, a 542 ps molecular dynamics trajectory was computed. For comparison with the 805 ps trajectory obtained for the uncomplexed enzyme, the theoretical fluorescence anisotropy decay of the unliganded protease and the inhibitor complex was calculated from the trajectories of the Trp6A/Trp6B and Trp42A/Trp42B transition dipole moments. This enabled us to directly compare the simulated data with the experimental picosecond time-resolved fluorescence data. Fitting both experimental and simulated data to the Kohlrausch-Williams-Watts (KWW) function exp(-t/tauk)beta revealed a very good agreement for the uncomplexed protease as well as for the SDZ283-910 complex. Binding of the inhibitor induced a faster decay of both the experimental and the computed protease fluorescence anisotropy decay. By this integrative approach, the atomic detail of inhibitor-induced changes in the conformational dynamics of the HIV-1 protease was experimentally verified and will be used for further inhibitor optimisation.     

Identification of the glycosaminoglycan binding site of the CC chemokine, MCP-1: implications for structure and function in vivo

In a recent study, we demonstrated that glycosaminoglycan (GAG) binding and oligomerization are essential for the in vivo function of the chemokines MCP-1/CCL2, RANTES/CCL5, and MIP-1beta/CCL4 (1). Binding to the GAG chains of cell surface proteoglycans is thought to facilitate the formation of high localized concentrations of chemokines, which in turn provide directional signals for leukocyte migration. To understand the molecular details of the chemokine-GAG interaction, in the present study we identified the GAG binding epitopes of MCP-1/CCL2 by characterizing a panel of surface alanine mutants in a series of heparin-binding assays. Using sedimentation equilibrium and cross-linking methods, we also observed that addition of heparin octasaccharide induces tetramer formation of MCP-1/CCL2. Although MCP-1/CCL2 forms a dimer in solution, both a dimer and tetramer have been observed by x-ray crystallography, providing a glimpse of the putative heparin-bound state. When the GAG binding residues are mapped onto the surface of the tetramer, the pattern that emerges is a continuous ring of basic residues encircling the tetramer, creating a positively charged surface well suited for binding GAGs. The structure also suggests several possible functional roles for GAG-induced oligomerization beyond retention of chemokines at the site of production.     

A structural and dynamic model for the interaction of interleukin-8 and glycosaminoglycans: support from isothermal fluorescence titrations

Binding of interleukin-8 (IL-8) to glycosaminoglycans (GAGs) on the surface of endothelial cells is crucial for the recruitment of neutrophils to an inflammatory site. Deriving structural knowledge about this interaction from in silico docking experiments has proved difficult because of the high flexibility and the size of GAGs. Therefore, we developed a docking method that takes into account ligand and protein flexibility by running approximately 15,000 molecular dynamics simulations of the docking event with different initial orientations of the binding partners. The method was shown to successfully reproduce the residues of basic fibroblast growth factor involved in GAG binding. Docking of a heparin hexasaccharide to IL-8 gave an interaction interface involving the basic residues His18, Lys20, Arg60, Lys64, Lys67, and Arg68. By subjecting IL-8 single-site mutants, in which these amino acids were replaced by alanine, to isothermal fluorescence titrations, the affinities for heparin were determined to be wtIL-8 > IL-8(H18A) > IL-8(R68A) > IL-8(K67A) > IL-8(K20A) > IL-8(R60A) > IL-8(K64A). A comparison with the binding energies calculated from the model revealed high values for wtIL-8 and the H18A mutant and significantly lower but similar energies for the remaining mutants. Connecting the two fully sulfated hexasaccharides bound to each of the two IL-8 monomers in the dimeric chemokine by an N-acetylated dodecasaccharide gave a complex structure in which the GAG molecule aligned in a parallel fashion to the N-terminal alpha-helices of IL-8 like a horseshoe. A 5-ns molecular dynamics simulation of this complex confirmed its structural stability and revealed a reorientation in both binding sites where a disaccharide became the central binding unit. Isothermal fluorescence titration experiments using differently sulfated heparin disaccharides confirmed that a single disaccharide can indeed bind IL-8 with high affinity.     

Rearrangements of integral membrane components during in vitro aging of sheep erythrocyte membranes

In vitro aged sheep erythrocytes and sheep erythrocyte ghosts spontaneously release vesicles that consist of long protrusions affixed to flattened headlike structures. The intramembranous particles seen on the protoplasmic face of freeze fracture electron micrographs of vesicle protrusions are arranged in paired particle rows. On the equivalent fracture face of headlike structures, the particle density is low; if particles are present, they are clustered along the rim of the flattened headlike structure and at the junction with the protrusion. The released vesicles are depleted of the intramembranous particles seen on the exoplasmic face of ghost but retain almost exclusively particles of the protoplasmic face. Correspondingly, the exoplasmic face of ghosts that have released vesicles reveals a 28 percent higher density of intramembranous particles than that of fresh ghosts. Purified vesicles are depleted of spectrin but retain integral membrane proteins, with one of an apparent mol wt of 160,000 accounting for nearly 50 percent of the total protein (Lutz, H.U.,R. Barber, and R.F. McGuire. 1976. J. Biol. Chem. 251:3500-3510). When vesicles are modified with the cleavable cross-linking reagent [(35)S]dithiobis (succinimidyl propionate)at 0 degrees C, the 160,000 mol wt protein is rapidly converted to disulfide-linked dimers and higher oligomers. Exposure of intact ghosts to the reagent in the same way fails to yield equivalent polymers. A comparison of the morphological and biochemical aspects of ghosts and vesicles suggest that a marked rearrangement of membrane proteins accompanies the supramolecular redistribution of intramembranous particles during spontaneous vesiculation. The results also suggest that the paired particles of the protoplasmic face of vesicle protrusions are arranged in paired helices and contain the 160,000 mol wt protein as dimers.     

A proteomic snapshot of the human heat shock protein 90 interactome

Heat shock protein 90 (Hsp90) is a molecular chaperone which modulates several signalling pathways within a cell. By applying co-immunoprecipitation with endogeneous Hsp90, we were able to identify 39 novel protein interaction partners of this chaperone in human embryonic kidney cells (HEK293). Interestingly, levels of DNA-activated protein kinase catalytic subunit, an Hsp90 interaction partner found in this study, were found to be sensitive to Hsp90 inhibitor treatment only in HeLa cells but not in HEK293 cells referring to the tumorgenicity of this chaperone.     

The influence of ATP on the association and unfolding of the tyrosine repressor ligand response domain of Haemophilus influenzae

The secondary structure of the ligand response domain of the Haemophilus influenzae tyrosine repressor, TyrR(lrd), was investigated using CD spectroscopy which revealed 42.5% alpha-helix, 17.6% beta-sheet, and 39.9% loops. Quaternary structure analysis by fluorescence anisotropy showed that TyrR(lrd) is monomeric at a concentration of 100 nM to 2 microM but that the protein readily dimerizes in the presence of its natural ligand ATP. Equilibrium unfolding studies of TyrR(lrd) using guanidinium hydrochloride suggested a two-state model with no detectable stable intermediates. The unfolding transition monitored by CD spectroscopy was responsive to tyrosine and ATP resulting in a shift to higher denaturant concentrations in the presence of these ligands. Differential scanning calorimetry yielded melting temperatures, T(m), of 51.15 and 58.07 degrees C for the unliganded and for the ATP-liganded protein, respectively. ATP is thus proposed to be a major structural cofactor for the molecular architecture of TyrR(lrd).     

Spatial phylogenetics of the native woody plant species in Hainan,China

To better identify biodiversity hotspots for conservation on Hainan Island, a tropical island in southern China, we assessed spatial variation in phylogenetic diversity and species richness using 18,976 georeferenced specimen records and a newly reconstructed molecular phylogeny of 957 native woody plants. Within this framework, we delineated bioregions based on vegetation composition and mapped areas of neoendemism and paleoendemism to identify areas of priority for conservation. Our results reveal that the southwest of Hainan is the most important hot spot for endemism and plant diversity followed by the southeast area. The distribution of endemic species showed a scattered, rather than clustered, pattern on the island. Based on phylogenetic range‐weighted turnover metrics, we delineated three major vegetational zones in Hainan. These largely correspond to natural secondary growth and managed forests (e.g., rubber and timber forests) in central Hainan, old‐growth forests and natural secondary growth forest at the margins of Hainan, and nature reserves on the island (e.g., Jianfeng and Diaoluo National Nature Reserves). Our study helps to elucidate potential botanical conservation priorities for Hainan within an evolutionary, phylogenetic framework.     

A network-based approach to classify the three domains of life

Background

Chromosomal orthologs can reveal the shared ancestral gene set and their evolutionary trends. Additionally, physico-chemical properties of encoded proteins could provide information about functional adaptation and ecological niche requirements.

Results

We analyzed 7080 genes (five groups of 1416 orthologs each) from Rhizobiales species (S. meliloti, R. etli, and M. loti, plant symbionts; A. tumefaciens, a plant pathogen; and B. melitensis, an animal pathogen). We evaluated their phylogenetic relationships and observed three main topologies. The first, with closer association of R. etli to A. tumefaciens; the second with R. etli closer to S. meliloti; and the third with A. tumefaciens and S. meliloti as the closest pair. This was not unusual, given the close relatedness of these three species. We calculated the synonymous (dS) and nonsynonymous (dN) substitution rates of these orthologs, and found that informational and metabolic functions showed relatively low dN rates; in contrast, genes from hypothetical functions and cellular processes showed high dN rates. An alternative measure of sequence variability, percentage of changes by species, was used to evaluate the most specific proportion of amino acid residues from alignments. When dN was compared with that measure a high correlation was obtained, revealing that much of evolutive information was extracted with the percentage of changes by species at the amino acid level. By analyzing the sequence variability of orthologs with a set of five properties (polarity, electrostatic charge, formation of secondary structures, molecular volume, and amino acid composition), we found that physico-chemical characteristics of proteins correlated with specific functional roles, and association of species did not follow their typical phylogeny, probably reflecting more adaptation to their life styles and niche preferences. In addition, orthologs with low dN rates had residues with more positive values of polarity, volume and electrostatic charge.

Conclusions

These findings revealed that even when orthologs perform the same function in each genomic background, their sequences reveal important evolutionary tendencies and differences related to adaptation. This article was reviewed by: Dr. Purificación López-García, Prof. Jeffrey Townsend (nominated by Dr. J. Peter Gogarten), and Ms. Olga Kamneva.     

Hochstetter's frog (Leiopelma hochstetteri) egg,mobile larvae and froglet development

Egg strings and larvae of Hochstetter's frog (Leiopelma hochstetteri) were located at three widely separated North Island sites: in seeps at Brynderwyns in December 2004, in an open pool at Wharerino in March 2009, and in an underground pool near the Kaipawa Track, Coromandel, in late May 2009. Ten egg strings were also laid by captive frogs in water courses at Hamilton Zoo in April 2009. All egg strings held from 11 to 13 eggs. The egg strings laid in the Brynderwyns were regularly observed until metamorphosis was completed in March 2005. Twenty-four swimming larvae emerged from 25 capsules at c. 40 days after discovery, and at least 14 froglets were produced at c. 90 days. All of them developed in darkness, in a 120 ml pool <30 mm deep. The emerged froglets ranged from 9.8 to 10.8 mm snout-vent length. The detection of eggs, larvae and <11 mm froglets indicates that the egg laying period is at least from late September to May.