Brownian dynamics simulations of protein folding: Access to milliseconds time scale and beyond

Protein folding occurs on a time scale ranging from milliseconds to minutes for a majority of proteins. Computer simulation of protein folding, from a random configuration to the native structure, is nontrivial owing to the large disparity between the simulation and folding time scales. As an effort to overcome this limitation, simple models with idealized protein subdomains, e.g., the diffusion–collision model of Karplus and Weaver, have gained some popularity. We present here new results for the folding of a four-helix bundle within the framework of the diffusion–collision model. Even with such simplifying assumptions, a direct application of standard Brownian dynamics methods would consume 10,000 processor-years on current supercomputers. We circumvent this difficulty by invoking a special Brownian dynamics simulation. The method features the calculation of the mean passage time of an event from the flux overpopulation method and the sampling of events that lead to productive collisions even if their probability is extremely small (because of large free-energy barriers that separate them from the higher probability events). Using these developments, we demonstrate that a coarse-grained model of the four-helix bundle can be simulated in several days on current supercomputers. Furthermore, such simulations yield folding times that are in the range of time scales observed in experiments.

Tropomyosin implicated in host protective responses to microfilariae in onchocerciasis

A cDNA from adult female Onchocerca volvulus encoding the C-terminal portion of a tropomyosin isoform (termed MOv-14) has been shown previously to confer protective immunity in rodent models of onchocerciasis. The full-length sequence (designated Ov-tmy-1) obtained by PCR amplification, codes for a protein of 33 kDa and shares 91% identity with tropomyosins from other nematodes, falling to 57% identity with human α-tropomyosin. Ov-TMY-1 migrates with an apparent molecular mass of 42 kDa on SDS/PAGE and is present in all life-cycle stages, as determined by immunoblotting. Immunogold electron microscopy identified antigenic sites within muscle blocks and the cuticle of microfilariae and infective larvae. Anti-MOv14 antibodies were abundant in mice exhibiting serum-transferable protection against microfilariae conferred by vaccination with a PBS-soluble parasite extract. In contrast, little or no MOv14-specific antibody was present in mice inoculated with live microfilariae, in which resistance is mediated by antibody-independent mechanisms. In human infections, there was an inverse correlation between anti-tropomyosin IgG levels and densities of microfilariae in the skin. Seropositivity varied with the relative endemicity of infection. An immunodominant B cell epitope within Ov-TMY-1 (AQLLAEEADRKYD) was mapped to the N terminus of the MOv14 protein by using sera from protectively vaccinated mice. Intriguingly, the sequence coincides with an IgE-binding epitope within shrimp tropomyosin, believed to be responsible for hypersensitivity in individuals exhibiting allergy to shellfish. IgG and IgE antibodies reacting with the O. volvulus epitope were detected in human infections. It is concluded that antibody responses to tropomyosin may be important in limiting microfilarial densities in a proportion of individuals with onchocerciasis and have the potential to mediate hypersensitivity reactions to dead microfilariae, raising the possibility of a link with the immunopathology of infection.

Transgenic DNA integrated into the oat genome is frequently interspersed by host DNA

Integration of transgenic DNA into the plant genome was investigated in 13 transgenic oat (Avena sativa L.) lines produced using microprojectile bombardment with one or two cotransformed plasmids. In all transformation events, the transgenic DNA integrated into the plant genome consisted of intact transgene copies that were accompanied by multiple, rearranged, and/or truncated transgene fragments. All fragments of transgenic DNA cosegregated, indicating that they were integrated at single gene loci. Analysis of the structure of the transgenic loci indicated that the transgenic DNA was interspersed by the host genomic DNA. The number of insertions of transgenic DNA within the transgene loci varied from 2 to 12 among the 13 lines. Restriction endonucleases that do not cleave the introduced plasmids produced restriction fragments ranging from 3.6 to about 60 kb in length hybridizing to a probe comprising the introduced plasmids. Although the size of the interspersing host DNA within the transgene locus is unknown, the sizes of the transgene-hybridizing restriction fragments indicated that the entire transgene locus must be at least from 35–280 kb. The observation that all transgenic lines analyzed exhibited genomic interspersion of multiple clustered transgenes suggests a predominating integration mechanism. We propose that transgene integration at multiple clustered DNA replication forks could account for the observed interspersion of transgenic DNA with host genomic DNA within transgenic loci.

Crystal structure of Taq DNA polymerase in complex with an inhibitory Fab: The Fab is directed against an intermediate in the helix-coil dynamics of the enzyme

We report the crystal structure of Thermus aquaticus DNA polymerase I in complex with an inhibitory Fab, TP7, directed against the native enzyme. Some of the residues present in a helical conformation in the native enzyme have adopted a γ turn conformation in the complex. Taken together, structural information that describes alteration of helical structure and solution studies that demonstrate the ability of TP7 to inhibit 100% of the polymerase activity of the enzyme suggest that the change in conformation is probably caused by trapping of an intermediate in the helix-coil dynamics of this helix by the Fab. Antibodies directed against modified helices in proteins have long been anticipated. The present structure provides direct crystallographic evidence. The Fab binds within the DNA binding cleft of the polymerase domain, interacting with several residues that are used by the enzyme in binding the primer:template complex. This result unequivocally corroborates inferences drawn from binding experiments and modeling calculations that the inhibitory activity of this Fab is directly attributable to its interference with DNA binding by the polymerase domain of the enzyme. The combination of interactions made by the Fab residues in both the polymerase and the vestigial editing nuclease domain of the enzyme reveal the structural basis of its preference for binding to DNA polymerases of the Thermus species. The orientation of the structure-specific nuclease domain with respect to the polymerase domain is significantly different from that seen in other structures of this polymerase. This reorientation does not appear to be antibody-induced and implies remarkably high relative mobility between these two domains.

Inactivation of the mitogen-activated protein kinase Mps1 from the rice blast fungus prevents penetration of host cells but allows activation of plant defense responses

The rice blast fungus, Magnaporthe grisea, generates enormous turgor pressure within a specialized cell called the appressorium to breach the surface of host plant cells. Here, we show that a mitogen-activated protein kinase, Mps1, is essential for appressorium penetration. Mps1 is 85% similar to yeast Slt2 mitogen-activated protein kinase and can rescue the thermosensitive growth of slt2 null mutants. The mps1–1Δ mutants of M. grisea have some phenotypes in common with slt2 mutants of yeast, including sensitivity to cell-wall-digesting enzymes, but display additional phenotypes, including reduced sporulation and fertility. Interestingly, mps1–1Δ mutants are completely nonpathogenic because of the inability of appressoria to penetrate plant cell surfaces, suggesting that penetration requires remodeling of the appressorium wall through an Mps1-dependent signaling pathway. Although mps1–1Δ mutants are unable to cause disease, they are able to trigger early plant-cell defense responses, including the accumulation of autofluorescent compounds and the rearrangement of the actin cytoskeleton. We conclude that MPS1 is essential for pathogen penetration; however, penetration is not required for induction of some plant defense responses.

A novel alternate secretory pathway for the export of Plasmodium proteins into the host erythrocyte

The malarial parasite dramatically alters its host cell by exporting and targeting proteins to specific locations within the erythrocyte. Little is known about the mechanisms by which the parasite is able to carry out this extraparasite transport. The fungal metabolite brefeldin A (BFA) has been used to study the secretory pathway in eukaryotes. BFA treatment of infected erythrocytes inhibits protein export and results in the accumulation of exported Plasmodium proteins into a compartment that is at the parasite periphery. Parasite proteins that are normally localized to the erythrocyte membrane, to nonmembrane bound inclusions in the erythrocyte cytoplasm, or to the parasitophorous vacuolar membrane accumulate in this BFA-induced compartment. A single BFA-induced compartment is detected per parasite and the various exported proteins colocalize to this compartment regardless of their final destinations. Parasite membrane proteins do not accumulate in this novel compartment, but accumulate in the endoplasmic reticulum (ER), suggesting that the parasite has two secretory pathways. This alternate secretory pathway is established immediately after merozoite invasion and at least some dense granule proteins also use the alternate pathway. The BFA-induced compartment exhibits properties that are similar to the ER, but it is clearly distinct from the ER. We propose to call this new organelle the secondary ER of apicomplexa. This ER-like organelle is an early, if not the first, step in the export of Plasmodium proteins into the host erythrocyte.

Molecular dynamics of MHC genesis unraveled by sequence analysis of the 1,796,938-bp HLA class I region

The intensely studied MHC has become the paradigm for understanding the architectural evolution of vertebrate multigene families. The 4-Mb human MHC (also known as the HLA complex) encodes genes critically involved in the immune response, graft rejection, and disease susceptibility. Here we report the continuous 1,796,938-bp genomic sequence of the HLA class I region, linking genes between MICB and HLA-F. A total of 127 genes or potentially coding sequences were recognized within the analyzed sequence, establishing a high gene density of one per every 14.1 kb. The identification of 758 microsatellite provides tools for high-resolution mapping of HLA class I-associated disease genes. Most importantly, we establish that the repeated duplication and subsequent diversification of a minimal building block, MIC-HCGIX-3.8–1-P5-HCGIV-HLA class I-HCGII, engendered the present-day MHC. That the currently nonessential HLA-F and MICE genes have acted as progenitors to today’s immune-competent HLA-ABC and MICA/B genes provides experimental evidence for evolution by “birth and death,” which has general relevance to our understanding of the evolutionary forces driving vertebrate multigene families.

Mapping regions containing binding residues within functional domains of Plasmodium vivax and Plasmodium knowlesi erythrocyte-binding proteins

Invasion of erythrocytes by malaria parasites is mediated by specific molecular interactions. Whereas Plasmodium vivax and Plasmodium knowlesi use the Duffy blood group antigen, Plasmodium falciparum uses sialic acid residues of glycophorin A as receptors to invade human erythrocytes. P. knowlesi uses the Duffy antigen as well as other receptors to invade rhesus erythrocytes by multiple pathways. Parasite ligands that bind these receptors belong to a family of erythrocyte-binding proteins (EBP). The EBP family includes the P. vivax and P. knowlesi Duffy-binding proteins, P. knowlesi β and γ proteins, which bind alternate receptors on rhesus erythrocytes, and P. falciparum erythrocyte-binding antigen (EBA-175), which binds sialic acid residues of human glycophorin A. Binding domains of each EBP lie in a conserved N-terminal cysteine-rich region, region II, which contains around 330 amino acids with 12 to 14 conserved cysteines. Regions containing binding residues have now been mapped within P. vivax and P. knowlesi β region II. Chimeric domains containing P. vivax region II sequences fused to P. knowlesi β region II sequences were expressed on the surface of COS cells and tested for binding to erythrocytes. Binding residues of P. vivax region II lie in a 170-aa stretch between cysteines 4 and 7, and binding residues of P. knowlesi β region II lie in a 53-aa stretch between cysteines 4 and 5. Mapping regions responsible for receptor recognition is an important step toward understanding the structural basis for the interaction of these parasite ligands with host receptors.

Antigen-driven CD4+ T cell and HIV-1 dynamics: Residual viral replication under highly active antiretroviral therapy

Antigen-induced stimulation of the immune system can generate heterogeneity in CD4+ T cell division rates capable of explaining the temporal patterns seen in the decay of HIV-1 plasma RNA levels during highly active antiretroviral therapy. Posttreatment increases in peripheral CD4+ T cell counts are consistent with a mathematical model in which host cell redistribution between lymph nodes and peripheral blood is a function of viral burden. Model fits to patient data suggest that, although therapy reduces HIV replication below replacement levels, substantial residual replication continues. This residual replication has important consequences for long-term therapy and the evolution of drug resistance and represents a challenge for future treatment strategies.

The replication initiation protein of the broad-host-range plasmid RK2 is activated by the ClpX chaperone

Initiation and control of replication of the broad-host-range plasmid RK2 requires two plasmid-encoded elements, the replication origin (oriV) and the initiation protein TrfA. Purified TrfA is largely in the form of a dimer; however, only the monomeric form of the protein can bind specifically to the direct repeats (iterons) at the RK2 origin. The largely dimeric form of wild-type TrfA is inactive in the initiation of replication of RK2 in an in vitro replication system reconstituted from purified components. However, preincubation of the TrfA protein with the ClpX molecular chaperone isolated from Escherichia coli activates the initiator protein for replication in the purified system. We further observed that ClpX, in an ATP-dependent reaction, greatly increases the proportion of TrfA monomers and, therefore, the ability of this protein to bind to iterons localized within RK2 origin. Finally, a copy-up mutant of the TrfA protein which is largely in the monomer form is active in the reconstituted in vitro replication system, and its activity is not affected by ClpX.

In Vivo Observations of Myosin II Dynamics Support a Role in Rear RetractionV⃞

To investigate myosin II function in cell movement within a cell mass, we imaged green fluorescent protein-myosin heavy chain (GFP-MHC) cells moving within the tight mound of Dictyostelium discoideum. In the posterior cortex of cells undergoing rotational motion around the center of the mound, GFP-MHC cyclically formed a “C,” which converted to a spot as the cell retracted its rear. Consistent with an important role for myosin in rotation, cells failed to rotate when they lacked the myosin II heavy chain (MHC−) or when they contained predominantly monomeric myosin II (3xAsp). In cells lacking the myosin II regulatory light chain (RLC−), rotation was impaired and eventually ceased. These rotational defects reflect a mechanical problem in the 3xAsp and RLC− cells, because these mutants exhibited proper rotational guidance cues. MHC− cells exhibited disorganized and erratic rotational guidance cues, suggesting a requirement for the MHC in organizing these signals. However, the MHC− cells also exhibited mechanical defects in rotation, because they still moved aberrantly when seeded into wild-type mounds with proper rotational guidance cues. The mechanical defects in rotation may be mediated by the C-to-spot, because RLC− cells exhibited a defective C-to-spot, including a slower C-to-spot transition, consistent with this mutant’s slower rotational velocity.

Transmembrane Insertion of the Toxoplasma gondii GRA5 Protein Occurs after Soluble Secretion into the Host Cell

The intracellular parasite Toxoplasma gondii resides within a specialized compartment, the parasitophorous vacuole (PV), that resists fusion with host cell endocytic and lysosomal compartments. The PV is extensively modified by secretion of parasite proteins, including the dense granule protein GRA5 that is specifically targeted to the delimiting membrane of the PV (PVM). We show here that GRA5 is present both in a soluble form and in hydrophobic aggregates. GRA5 is secreted as a soluble form into the PV after which it becomes stably associated with the PVM. Topological studies demonstrated that GRA5 was inserted into the PVM as a transmembrane protein with its N-terminal domain extending into the cytoplasm and its C terminus in the vacuole lumen. Deletion of 8 of the 18 hydrophobic amino acids of the single predicted transmembrane domain resulted in the failure of GRA5 to associate with the PVM; yet it remained correctly packaged in the dense granules and was secreted as a soluble protein into the PV. Collectively, these studies demonstrate that the secretory pathway in Toxoplasma is unusual in two regards; it allows soluble export of proteins containing typical transmembrane domains and provides a mechanism for their insertion into a host cell membrane after secretion from the parasite.

The Dynamics of Golgi Protein Traffic Visualized in Living Yeast Cells

We describe for the first time the visualization of Golgi membranes in living yeast cells, using green fluorescent protein (GFP) chimeras. Late and early Golgi markers are present in distinct sets of scattered, moving cisternae. The immediate effects of temperature-sensitive mutations on the distribution of these markers give clues to the transport processes occurring. We show that the late Golgi marker GFP-Sft2p and the glycosyltransferases, Anp1p and Mnn1p, disperse into vesicle-like structures within minutes of a temperature shift in sec18, sft1, and sed5 cells, but not in sec14 cells. This is consistent with retrograde vesicular traffic, mediated by the vesicle SNARE Sft1p, to early cisternae containing the target SNARE Sed5p. Strikingly, Sed5p itself moves rapidly to the endoplasmic reticulum (ER) in sec12 cells, implying that it cycles through the ER. Electron microscopy shows that Golgi membranes vesiculate in sec18 cells within 10 min of a temperature shift. These results emphasize the dynamic nature of Golgi cisternae and satisfy the kinetic requirements of a cisternal maturation model in which all resident proteins must undergo retrograde vesicular transport, either within the Golgi complex or from there to the ER, as anterograde cargo advances.

Density-dependent decline of host abundance resulting from a new infectious disease

Although many new diseases have emerged within the past 2 decades [Cohen, M. L. (1998) Brit. Med. Bull. 54, 523–532], attributing low numbers of animal hosts to the existence of even a new pathogen is problematic. This is because very rarely does one have data on host abundance before and after the epizootic as well as detailed descriptions of pathogen prevalence [Dobson, A. P. & Hudson, P. J. (1985) in Ecology of Infectious Diseases in Natural Populations, eds. Grenfell, B. T. & Dobson, A. P. (Cambridge Univ. Press, Cambridge, U.K.), pp. 52–89]. Month by month we tracked the spread of the epizootic of an apparently novel strain of a widespread poultry pathogen, Mycoplasma gallisepticum, through a previously unknown host, the house finch, whose abundance has been monitored over past decades. Here we are able to demonstrate a causal relationship between high disease prevalence and declining house finch abundance throughout the eastern half of North America because the epizootic reached different parts of the house finch range at different times. Three years after the epizootic arrived, house finch abundance stabilized at similar levels, although house finch abundance had been high and stable in some areas but low and rapidly increasing in others. This result, not previously documented in wild populations, is as expected from theory if transmission of the disease was density dependent.

N-terminal transcriptional activation domain of LZIP comprises two LxxLL motifs and the Host Cell Factor-1 binding motif

Host Cell Factor-1 (HCF-1, C1) was first identified as a cellular target for the herpes simplex virus transcriptional activator VP16. Association between HCF and VP16 leads to the assembly of a multiprotein enhancer complex that stimulates viral immediate-early gene transcription. HCF-1 is expressed in all cells and is required for progression through G1 phase of the cell cycle. In addition to VP16, HCF-1 associates with a cellular bZIP protein known as LZIP (or Luman). Both LZIP and VP16 contain a four-amino acid HCF-binding motif, recognized by the N-terminal β-propeller domain of HCF-1. Herein, we show that the N-terminal 92 amino acids of LZIP contain a potent transcriptional activation domain composed of three elements: the HCF-binding motif and two LxxLL motifs. LxxLL motifs are found in a number of transcriptional coactivators and mediate protein–protein interactions, notably recognition of the nuclear hormone receptors. LZIP is an example of a sequence-specific DNA-binding protein that uses LxxLL motifs within its activation domain to stimulate transcription. The LxxLL motifs are not required for association with the HCF-1 β-propeller and instead interact with other regions in HCF-1 or recruit additional cofactors.