Trimeresurus stejnegeri snake venom plasminogen activator - Site-directed mutagenesis and molecular modeling

The specific plasminogen activator from Trimeresurus stejnegeri venom (TSV-PA) is a serine proteinase presenting 23% sequence identity with the proteinase domain of tissue type plasminogen activator, and 63% with batroxobin, a fibrinogen clotting enzyme from Bothrops atrox venom that does not activate plasminogen. TSV-PA contains six disulfide bonds and has been successfully overexpressed in Escherichia coli (Zhang, Y., Wisner, A., Xiong, Y. L,, and Bon, C, (1995) J. Biol. Chem. 270, 10246-10255), To identify the functional domains of TSV-PA, we focused on three short peptide fragments of TSV-PA showing important sequence differences with batroxobin and other venom serine proteinases. Molecular modeling shows that these sequences are located in surface loop regions, one of which is next to the catalytic site, When these sequences were replaced in TSV-PA by the equivalent batroxobin residues none generated either fibrinogen-clotting or direct fibrinogenolytic activity, Two of the replacements had little effect in general and are not critical to the specificity of TSV-PA for plasminogen. Nevertheless, the third replacement, produced by the conversion of the sequence DDE 96a-98 to NVI, significantly increased the K-m for some tripeptide chromogenic substrates and resulted in undetectable plasminogen activation, indicating the key role that the sequence plays in substrate recognition by the enzyme.

A novel short neurotoxin, cobrotoxin c, from monocellate cobra (Naja kaouthia) venom: isolation and purification, primary and secondary structure determination, and tertiary structure modeling

A novel short neurotoxin, cobrotoxin c (CBT C) was isolated from the venom of monocellate cobra (Naja kaouthia) using a combination of ion-exchange chromatography and FPLC. Its primary structure was determined by Edman degradation. CBT C is composed of 61 amino acid residues. It differs from cobrotoxin b (CBT B) by only two amino acid substitutions, Thr/Ala11 and Arg/Thr56, which are not located on the functionally important regions by sequence similarity. However, the LD50 is 0.08 mg/g to mice, i.e. approximately five-fold higher than for CBT B. Strikingly, a structure-function relationship analysis suggests the existence of a functionally important domain on the outside of Loop III of CBT C. The functionally important basic residues on the outside of Loop III might have a pairwise interaction with alpha subunit, instead of gamma or delta subunits of the nicotinic acetylcholine receptor (nAChR). (C) 2002 Elsevier Science Inc. All rights reserved.

Structural analysis of human CCR2b and primate CCR2b by molecular modeling and molecular dynamics simulation

CCR2b, a chemokine receptor for MCP-1, -2, -3, -4, plays an important role in a variety of diseases involving infection, inflammation, and/or injury, as well as being a coreceptor for HIV-1 infection. Two models of human CCR2b (hCCR2b) were generated by h

Molecular modeling on some human interleukins sharing gamma c of interleukin-2 receptor, and structure-function relationships

Five models for human interleukin-7 (HIL-7), HIL-9, HIL-13, HIL-15 and HIL-17 have been generated by SYBYL software package. The primary models were optimized using molecular dynamics and molecular mechanics methods. The final models were optimized using a steepest descent algorithm and a subsequent conjugate gradient method. The complexes with these interleukins and the common gamma chain of interleukin-2 receptor (IL-2R) were constructed and subjected to energy minimization. We found residues, such as Gln127 and Tyr103, of the common gamma chain of IL-2R are very important. Other residues, e.g. Lys70, Asn128 and Glu162, are also significant. Four hydrophobic grooves and two hydrophilic sites converge at the active site triad of the gamma chain. The binding sites of these interleukins interaction with the common gamma chain exist in the first helical and/or the fourth helical domains. Therefore, we conclude that these interleukins binds to the common gamma chain of IL-2R by the first and the fourth helix domain. Especially at the binding sites of some residues (lysine, arginine, asparagine, glutamic acid and aspartic acid), with a discontinuous region of the common gamma chain of IL-2R, termed the interleukins binding sites (103-210). The study of these sites can be important for the development of new drugs. (C) 2000 Elsevier Science B.V. All rights reserved.

Stability and molecular modeling of triplex DNA inhibiting DNA binding protein binding to the core promoter of hepatitis B virus.

Two three-dimensional structure models of the 21nt oligodeoxyribonucleotides, CPI (G3TG-2TGT2G5TG2TGT) and CP3 (TGTG2TGST2GTG2TG3), were constructed by InsightII (MSI) software in IRIS Indigo2 (SGI) workstation using the crystal structure of TAT tripler formation as the template. The initial structures subsequently were minimized by molecular mechanics. The final structures were believed as the dominant conformation. The results showed that the energy of CP1 is lower than that of CP3, and the former is more stable than the latter. Moreover, the results further proved that the 21nt oligodeoxyribo-nucleotide CP1 stably combines with the core promoter (Cp) fragment of hepatitis B virus (HBV) to form a tripler DNA, and CP1 specifically inhibits a specific cellular factor (DNA binding protein) binding to Cp fragment. These results indicated that specific repression of gene transcription of HBV DNA might be possible by tripler-formation DNA.

Molecular modeling on human CCR5 receptors and complex with CD4 antigens and HIV-1 envelope glycoprotein gp120

AIM: To investigate the interaction between human CCR5 receptors (CCR5) and HIV-1 envelope glycoprotein gp120 (HIV-1 gp120) and HIV-1 receptor CD4 antigens (CD4). METHODS: The structurally con served regions (SCR) of human CCR5 was built by the SYBYL/Biopolymer module using the corresponding transmembrane (TM) domain of bacteriorhodopsin (bR) as the template. The coordinates for amino-ter minal residue sequence, and carboxyl-terminal residue sequence, extracellular and cytoplasmic loops were generated using LOOP SEARCH algorithm. Subsequently the structural model was merged into the complex with HIV-1 gp120 and CD4. RESULTS: Human CCR5 interacted with both an HIV-1 gp120 and CD4. The N-terminal residues (especially Met1 and Gln4) of human CCR5, contacted with CD4 residues, mainly 7Nith one span (56 - 59) of CD4 in electrostatic interaction and hydrogen-bonds. The binding sites of human CCR5 were buried in a hydrophobic center surrounded by a highly basic periphery. On the other hand, direct interatomic contacts were made between ? CCR5 residues and 6 gp120 amino-acid residues, which included van der Waals contacts, hydrophobic interaction, and hydrogen bonds. CONCLUSION: The interaction model should be helpful for rational design of novel anti-HIV drugs.

Structural and functional characterization of the human CCR5 receptor in complex with HIV gp120 envelope glycoprotein and CD4 receptor by molecular modeling studies

The entry of human immunodeficiency virus (HIV) into cells depends on a sequential interaction of the gp120 envelope glycoprotein with the cellular receptors CD4 and members of the chemokine receptor family. The CC chemokine receptor CCR5 is such a receptor for several chemokines and a major coreceptor for the entry of R5 HIV type-1 (HIV-1) into cells. Although many studies focus on the interaction of CCR5 with HIV-1, the corresponding interaction sites in CCR5 and gp120 have not been matched. Here we used an approach combining protein structure modeling, docking and molecular dynamics simulation to build a series of structural models of the CCR5 in complexes with gp120 and CD4. Interactions such as hydrogen bonds, salt bridges and van der Waals contacts between CCR5 and gp120 were investigated. Three snapshots of CCR5-gp120-CD4 models revealed that the initial interactions of CCR5 with gp120 are involved in the negatively charged N-terminus (Nt) region of CCR5 and positively charged bridging sheet region of gp120. Further interactions occurred between extracellular loop2 (ECL2) of CCR5 and the base of V3 loop regions of gp120. These interactions may induce the conformational changes in gp120 and lead to the final entry of HIV into the cell. These results not only strongly support the two-step gp120-CCR5 binding mechanism, but also rationalize extensive biological data about the role of CCR5 in HIV-1 gp120 binding and entry, and may guide efforts to design novel inhibitors.

Characterize dynamic conformational space of human CCR5 extracellular domain by molecular modeling and molecular dynamics simulation

The chemokine receptor CCR5 is the receptor for several chemokines and major coreceptor for R5 human immunodeficiency virus type-1 strains entry into cell. Three-dimensional models of CCR5 were built by using homology modeling approach and 1 ns molecular dynamics (MD) simulation, because studies of site-directed mutagenesis and chimeric receptors have indicated that the N-terminus (Nt) and extracellular loops (ECLs) of CCR5 are important for ligands binding and viral fusion and entry, special attention was focused on disulfide bond function, conformational flexibility, hydrogen bonding, electrostatic interactions, and solvent-accessible surface area of Nt and ECLs of this protein part. We found that the extracellular segments of CCR5 formed a well-packet globular domain with complex interactions occurred between them in a majority of time of MID simulation, but Nt region could protrude from this domain sometimes. The disulfide bond Cys20-Cys269 is essential in controlling specific orientation of Nt region and maintaining conformational integrity of extracellular domain. RMS comparison analysis between conformers revealed the ECL1 of CCR5 stays relative rigid, whereas the ECL2 and Nt are rather flexible. Solvent-accessible surface area calculations indicated that the charged residues within Nt and ECL2 are often exposed to solvent. Integrating these results with available experimental data, a two-step gp120-CCR5 binding mechanism was proposed. The dynamic interaction of CCR5 extracellular domain with gp120 was emphasized. (C) 2004 Elsevier B.V. All rights reserved.

Cooperation of Mtmr8 with PI3K Regulates Actin Filament Modeling and Muscle Development in Zebrafish

Background: It has been shown that mutations in at least four myotubularin family genes (MTM1, MTMR1, 2 and 13) are causative for human neuromuscular disorders. However, the pathway and regulative mechanism remain unknown. Methodology/Principal Findings: Here, we reported a new role for Mtmr8 in neuromuscular development of zebrafish. Firstly, we cloned and characterized zebrafish Mtmr8, and revealed the expression pattern predominantly in the eye field and somites during early somitogenesis. Using morpholino knockdown, then, we observed that loss-of-function of Mtmr8 led to defects in somitogenesis. Subsequently, the possible underlying mechanism and signal pathway were examined. We first checked the Akt phosphorylation, and observed an increase of Akt phosphorylation in the morphant embryos. Furthermore, we studied the PH/G domain function within Mtmr8. Although the PH/G domain deletion by itself did not result in embryonic defect, addition of PI3K inhibitor LY294002 did give a defective phenotype in the PH/G deletion morphants, indicating that the PH/G domain was essential for Mtmr8's function. Moreover, we investigated the cooperation of Mtmr8 with PI3K in actin filament modeling and muscle development, and found that both Mtmr8-MO1 and Mtmr8-MO2+LY294002 led to the disorganization of the actin cytoskeleton. In addition, we revealed a possible participation of Mtmr8 in the Hedgehog pathway, and cell transplantation experiments showed that Mtmr8 worked in a non-cell autonomous manner in actin modeling. Conclusion/Significance: The above data indicate that a conserved functional cooperation of Mtmr8 with PI3K regulates actin filament modeling and muscle development in zebrafish, and reveal a possible participation of Mtmr8 in the Hedgehog pathway. Therefore, this work provides a new clue to study the physiological function of MTM family members.

Modeling phytoplankton dynamics in the River Darling (Australia) using the radial basis function neural network

A radial basis function neural network was employed to model the abundance of cyanobacteria. The trained network could predict the populations of two bloom forming algal taxa with high accuracy, Nostocales spp. and Anabaena spp., in the River Darling, Australia. To elucidate the population dynamics for both Nostocales spp. and Anabaena spp., sensitivity analysis was performed with the following results. Total Kjeldahl nitrogen had a very strong influence on the abundance of the two algal taxa, electrical conductivity had a very strong negative relationship with the population of the two algal species, and flow was identified as one dominant factor influencing algal blooms after a scatter plot revealed that high flow could significantly reduce the algal biomass for both Nostocales spp. and Anabaena spp. Other variables such as turbidity, color, and pH were less important in determining the abundance and succession of the algal blooms.