Electrostatic Interactions Between Glycosaminoglycan Molecules

The electrostatic interactions between nearest-neighbouring chondroitin sulfate glycosaminoglycan (CS-GAG) molecular chains are obtained on the bottle brush conformation of proteoglycan aggrecan based on an asymptotic solution of the Poisson-Boltzmann equation the CS-GAGs satisfy under the physiological conditions of articular cartilage. The present results show that the interactions are associated intimately with the minimum separation distance and mutual angle between the molecular chains themselves. Further analysis indicates that the electrostatic interactions are not only expressed to be purely exponential in separation distance and decrease with the increasing mutual angle but also dependent sensitively on the saline concentration in the electrolyte solution within the tissue, which is in agreement with the existed relevant conclusions.

Numerical Analysis on Intermittent Flow in a Longitudinal Section of Bingham Fluid along a Slope

The transition process of intermittent flow in a longitudinal section of Bingham fluid from initial distribution to fully developed state was numerically investigated in this paper. The influences of slope dimensionless runoff Q* and viscosity μ0* on the dimensionless surge speed U* were also presented in a wide range of parameters. By one typical example, the intermittent flow possessed wave characteristics and showed a supercritical flow conformation for a fully developed flow. The distributions of gravity and bed drag along the flow path and the velocity distribution of flow field were also analyzed.

Forced Dissociation of Selectin-ligand Complexes Using Steered Molecular Dynamics Simulation

Selectin-ligand interactions are crucial to such biological processes as inflammatory cascade or tumor metastasis. How transient formation and dissociation of selectin-ligand bonds in blood flow are coupled to molecular conformation at atomic level, however, has not been well understood. In this study, steered molecular dynamics (SMD) simulations were used to elucidate the intramolecular and intermolecular conformational evolutions involved in forced dissociation of three selectin-ligand systems: the construct consisting of P-selectin lectin (Lec) and epidermal growth factor (EGF)-like domains (P-LE) interacting with synthesized sulfoglycopeptide or SGP-3, P-LE with sialyl Lewis X (sLeX), and E-LE with sLeX. SMD simulations were based on newly built-up force field parameters including carbohydrate units and sulfated tyrosine(s) using an analogy approach. The simulations demonstrated that the complex dissociation was coupled to the molecular extension. While the intramolecular unraveling in P-LESGP-3 system mainly resulted from the destroy of the two anti-parallel sheets of EGF domain and the breakage of hydrogen-bond cluster at the Lec-EGF interface, the intermolecular dissociation was mainly determined by separation of fucose (FUC) from Ca2+ ion in all three systems. Conformational changes during forced dissociations depended on pulling velocities and forces, as well as on how the force was applied. This work provides an insight into better understanding of conformational changes and adhesive functionality of selectin-ligand interactions under external forces.

Simulated pathogenic conformational switch regions matched well with the biochemical findings

Pathogenic conformational conversion is a general causation of many disease, such as transmissible spon- giform encephalopathy (TSE) caused by misfolding of prion, sickle cell anemia, and etc. In such structural changes, misfolding occurs in regions important for the stability of native structure firstly. This destabi- lizes the normal conformation and leads to subsequent errors in folding pathway. Sites involved in the first stage can be deemed switch regions of the protein, and are vital for conformational conversion. Namely it could be a switch of disease at residue level. Here we report an algorithm that can identify such sites computationally with an accuracy of 93%, by calculating the probability of the native structure of a short segment jumping to a mistake one. Knowledge of such switch sites could be used to target clinical therapy, study physiological and pathologic mechanism of protein, and etc.

Switch Region for Pathogenic Structural Change in Conformational Disease and Its Prediction

Many diseases are believed to be related to abnormal protein folding. In the first step of such pathogenic structural changes, misfolding occurs in regions important for the stability of the native structure. This destabilizes the normal protein conformation, while exposing the previously hidden aggregation-prone regions, leading to subsequent errors in the folding pathway. Sites involved in this first stage can be deemed switch regions of the protein, and can represent perfect binding targets for drugs to block the abnormal folding pathway and prevent pathogenic conformational changes. In this study, a prediction algorithm for the switch regions responsible for the start of pathogenic structural changes is introduced. With an accuracy of 94%, this algorithm can successfully find short segments covering sites significant in triggering conformational diseases (CDs) and is the first that can predict switch regions for various CDs. To illustrate its effectiveness in dealing with urgent public health problems, the reason of the increased pathogenicity of H5N1 influenza virus is analyzed; the mechanisms of the pandemic swine-origin 2009 A(H1N1) influenza virus in overcoming species barriers and in infecting large number of potential patients are also suggested. It is shown that the algorithm is a potential tool useful in the study of the pathology of CDs because: (1) it can identify the origin of pathogenic structural conversion with high sensitivity and specificity, and (2) it provides an ideal target for clinical treatment.

Molecular Modeling and Affinity Determination of scFv Antibody: Proper Linker Peptide Enhances Its Activity

One of existing strategies to engineer active antibody is to link VH and VL domains via a linker peptide. How the composition, length, and conformation of the linker affect antibody activity, however, remains poorly understood. In this study, a dual approach that coordinates molecule modeling, biological measurements, and affinity evaluation was developed to quantify the binding activity of a novel stable miniaturized anti-CD20 antibody or singlechain fragment variable (scFv) with a linker peptide. Upon computer-guided homology modeling, distance geometry analysis, and molecular superimposition and optimization, three new linker peptides PT1, PT2, and PT3 with respective 7, 10, and 15 residues were proposed and three engineered antibodies were then constructed by linking the cloned VH and VL domains and fusing to a derivative of human IgG1. The binding stability and activity of scFv-Fc chimera to CD20 antigen was quantified using a micropipette adhesion frequency assay and a Scatchard analysis. Our data indicated that the binding affinity was similar for the chimera with PT2 or PT3 and ~24-fold higher than that for the chimera with PT1, supporting theoretical predictions in molecular modeling. These results further the understanding in the impact of linker peptide on antibody structure and activity.

Detection of constitutive heterodimerization of the integrin Mac-1 subunits by fluorescence resonance energy transfer in living cells

Macrophage differentiation antigen associated with complement three receptor function (Mac-1) belongs to beta(2) subfamily of integrins that mediate important cell-cell and cell-extracellular matrix interactions. Biochemical studies have indicated that Mac-1 is a constitutive heterodimer in vitro. Here, we detected the heterodimerization of Mac-1 subunits in living cells by means of two fluorescence resonance energy transfer (FRET) techniques (fluorescence microscopy and fluorescence spectroscopy) and our results demonstrated that there is constitutive heterodimerization of the Mac-1 subunits and this constitutive heterodimerization of the Mac-1 subunits is cell-type independent. Through FRET imaging, we found that heterodimers of Mac-1 mainly localized in plasma membrane, perinuclear, and Golgi area in living cells. Furthermore, through analysis of the estimated physical distances between cyan fluorescent protein (CFP) and yellow fluorescent protein (YFP) fused to Mac-1 subunits, we suggested that the conformation of Mac-1 subunits is not affected by the fusion of CFP or YFP and inferred that Mac-1 subunits take different conformation when expressed in Chinese hamster ovary (CHO) and human embryonic kidney (HEK) 293T cells, respectively. (c) 2006 Elsevier Inc. All rights reserved.

魔芋葡甘聚糖结构及其应用研究

本文以葡甘聚糖为试材，运用分子模拟同仪器分析相结合的手段，预测了葡甘聚糖分子链的高级结构，分析了无机分子对其结构、性能的影响，探讨了葡甘聚糖与卡拉胶微观作用机理。 主要研究内容与结果如下： 1. 葡甘聚糖单链高级结构的预测 利用Hyperchem7.0、VM2.0分子结构计算软件采用分子动力学和分子力学的方法，以真空中葡甘聚糖单链为研究模型，研究了聚合度、取代基对动态构象的影响及影响链构象的作用力。首次提出了KGM链的动态模型，得到了以下结果：聚合度影响其链形态和稳定性，对于高聚合度的魔芋葡甘聚糖来说，其链呈现无规卷曲状态且稳定性下降，在整个动态运动过程中KGM链脱乙酰基前后都呈现无规卷曲状态，而且其伸展和卷曲的变化是周期性的，表现出了很好的柔性，说明乙酰基不是影响其链形态的主要因素，二面角能和静电作用是真空中影响单链构象的主要的键合作用力和非键合作用力，但是乙酰基对氢键作用的影响较大。 2. 无机分子对葡甘聚糖溶液体系结构性能影响的研究 利用Hyperchem7.0分子结构计算软件，采用分子动力学及红外光谱、核磁共振等技术，对无机分子对葡甘聚糖体系的影响进行分析，很好的解释了性能变化的结构原因，结合以往的研究及参考文献得出以下结论：KGM在碱性条件下由于化学作用乙酰基的脱除分子间氢键作用的加强提高了凝胶强度，分子间氢键的主要作用位点是葡萄糖 的O(6)与甘露糖的O(2)之间;硼与KGM形成的分子内和分子间配合作用及分子间作用力氢键的增强是KGM特性粘度和致密性提高的主要原因，分子间型配位反应发生在葡萄糖和甘露糖两个糖环之间的几率最大;加入尿素后，表现为宏观性能的下降，葡甘聚糖氢键网络被破坏，氢键的作用位点由甘露糖的O(2)、O（3）变为O(4)，葡萄糖的O（3）、O(6)变为O（1）、O（2）。 3. 葡甘聚糖与卡拉胶共混作用的研究 利用Hyperchem7.0分子结构计算软件运用分子动力学方法、DSC、红外光谱技术，研究了葡甘聚糖与卡拉胶的微观结构及作用过程，揭示了性能变化的结构原因和分子之间的作用位点。得出了以下结论：葡甘聚糖同卡拉胶共混后通过分子间氢键作用形成了强度高、弹性好的热可逆凝胶。其凝胶强度与单一胶相比较，凝胶特性得到了很大的改善。通过红外光谱可以发现形成复合溶胶后化学基团没有发生本质上的改变，但是氢键缔合作用增强;通过DSC分析可以发现仅出现1个吸热峰，两种生物大分子达到了相容的结果，经过分子动力学模拟表明，与单一体系比较，葡甘聚糖与卡拉胶共混时稳定性提高，分子间氢键作用力明显增强，主要的作用位点是葡甘聚糖的上甘露糖的O(2)、O(4)、O(6)、乙酰基位置及卡拉胶上糖环上的 O(6)、硫酸基。

Molecular dynamics research of G249 and S249 substitutions of p53 protein.

,The molecular dynamics research of the core domain of p53 protein crystal structure shows that besides the stability in biochemistry this domain also shows a high stability in molecular mechanics. Based on that work, the residue R249 was substituted with amino acids Gly and Ser respectively, and molecular dynamics researches were performed separately. The results show that these substitutions cause a relax tendency between loop2 and 3 domains, leading to an alteration of the whole conformation of p53 core domain and ruining its stability. The results visually explains the mechanism of p53 changes in immunological and biochemical reactions, which are caused by 249 residue substitutions from 3-D structure variations.

Novel cathelicidin-derived antimicrobial peptides from Equus asinus

In the present study, EA-CATH1 and EA-CATH2 were identified from a constructed lung cDNA library of donkey (Equus asinus) as members of cathelicidin-derived antimicrobial peptides, using a nested PCR-based cloning strategy. Composed of 25 and 26 residues, respectively, EA-CATH1 and EA-CATH2 are smaller than most other cathelicidins and have no sequence homology to other cathelicidins identified to date. Chemically synthesized EA-CATH1 exerted potent antimicrobial activity against most of the 32 strains of bacteria and fungi tested, especially the clinically isolated drug-resistant strains, and minimal inhibitory concentration values against Gram-positive bacteria were mostly in the range of 0.3-2.4 mu g center dot mL-1. EA-CATH1 showed an extraordinary serum stability and no haemolytic activity against human erythrocytes in a dose up to 20 mu g center dot mL-1. CD spectra showed that EA-CATH1 mainly adopts an alpha-helical conformation in a 50% trifluoroethanol/water solution, but a random coil in aqueous solution. Scanning electron microscope observations of Staphylococcus aureus (ATCC2592) treated with EA-CATH1 demonstrated that EA-CATH could cause rapid disruption of the bacterial membrane, and in turn lead to cell lysis. This might explain the much faster killing kinetics of EA-CATH1 than conventional antibiotics revealed by killing kinetics data. In the presence of CaCl2, EA-CATH1 exerted haemagglutination activity, which might potentiate an inhibition against the bacterial polyprotein interaction with the host erythrocyte surface, thereby possibly restricting bacterial colonization and spread.

The three dimensional structure and biological domains of the p53 C-terminus.

It was expected that there are a coil (289 similar to 325) and two a helix (alpha(1)368 similar to 373, alpha(2)381 similar to 388) structures in p53 protein C-terminal region based on its mRNA secondary structure template and Chou-Fasman's protein secondary structure principle of prediction. The result was conformed by the other four methods of protein secondary structure prediction that are based on the multiple sequence alignment (accuracy = 73.20%). Combine with the 31 amino acids crystal structure of the oligomerization, the three dimensional conformation of p53 C-terminal 108 residues was built using the SGI INDIGO(2) computer. This structure further expounds the relationship among those biological function domains of p53 C- terminus at three-dimensional level.

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.

Dynamics of RNA with A6 bulge loop

In order to understand the behavior of RNAs with large bulges In solution, molecular dynamics was performed on the RNA molecule in water with A6 bulge. The result of simulation showed that nonstacked conformation Is the main conformation in large bulges, and the backbone of large bulge is of great conformational flexibility, but bulges-induced bends are relatively rigid. The fluctuation in bulge has little influence on the bend angle of RNAs.