Generating artificial homologous proteins according to the representative family character in molecular mechanics properties - an attempt in validating an underlying rule of protein evolution

The molecular mechanics property is the foundation of many characters of proteins. Based on intramolecular hydrophobic force network, the representative family character underlying a protein’s mechanics property is described by a simple two-letter scheme. The tendency of a sequence to become a member of a protein family is scored according to this mathematical representation. Remote homologs of the WW-domain family could be easily designed using such a mechanistic signature of protein homology. Experimental validation showed that nearly all artificial homologs have the representative folding and bioactivity of their assigned family. Since the molecular mechanics property is the only consideration in this study, the results indicate its possible role in the generation of new members of a protein family during evolution.

Molecular Mechanics Simulation of Recognition Identity of tRNA~(His/GUG)

本文选用经过实验验证的碱基序列 ,用简化的方式 ,构建了被水分子和镁离子修饰的核酸序列的分子模型 ,应用分子力学模拟方法对序列进行能量优化 ,对优化后序列的构象参数、成键状况和能量数据等进行了分析。对tRNAHHis GUG的识别特性作了初步的探索 ,得到了和实验结果相近的结论。此外 ,还从能力学的角度讨论了溶剂 -溶质 -溶剂相互作用形成的网状氢键网络对核酸结构稳定性的影响 ,探讨了非Crick_WatsonGU、UU配对的能力学特征并存在于被水分子和镁离子修饰的核酸序列中的GU、UU配对情况。

X-ray crystal structure and molecular mechanics calculations of (2,6-iso-dipropyl-phenylamide) dimethyl (tetra-methylcyclopenta-dienyl) silane titanium dichloride

Crystal and molecular structure of (2.6-dipropylphenylamide) dimethyl (tetra-methyl cyclopentadienyl) silane titanium dichloride (I) was fully characterized by X-ray diffraction. The crystal is obtained from a mixture of ether/hexane as orthorhombic. with a = 12.658 (3) Angstrom. b = 16.62 (3) Angstrom. c = 11.760 (2) Angstrom. V = 2474.2 (9) Angstrom(3). Z = 4, space group Pnma. R = 0.0399; Componud I compose of the pi-bounded ring with its dimethylsilyl-dipropyl phenyl amido group and the two terminal chloride atoms coordinated to central metal to form a so-called constrained geometry catalyst (CGC) structure. The result of molecular mechanics (MM) calculations on compound I shows that bond lengths and bond angles from the MM calculation are comparable to the data obtained from the X-ray diffraction study. The relation of the structure of CGCs and their catalytic activity by MM calculations is also discussed.

A study of a molecular mechanics field used in simulating enantiorecognition

A molecular mechanics field, Alchemy II, was utilized to model the chiral recognition between S-N-acetyl-alpha-methyl-alpha-naphthylamine and (R, S)-N-(3, 5-dinitrophrnyl)-alpha-methyl-benzeneacetamide and between beta-cyclodextrin and (R, S)-fenoprofen. Some preliminary results have: been obtained to sustain the three-point action models and the induce-fit action in enantiorecognition.

QM/MM and Classical Molecular Dynamics Simulation of Histidine-Tagged Peptide Immobilization on Nickel Surface

The hybrid quantum mechanics (QM) and molecular mechanics (MM) method is employed to simulate the His-tagged peptide adsorption to ionized region of nickel surface. Based on the previous experiments, the peptide interaction with one Ni ion is considered. In the QM/MM calculation, the imidazoles on the side chain of the peptide and the metal ion with several neighboring water molecules are treated as QM part calculated by “GAMESS”, and the rest atoms are treated as MM part calculated by “TINKER”. The integrated molecular orbital/molecular mechanics (IMOMM) method is used to deal with theQMpart with the transitional metal. By using the QM/MM method, we optimize the structure of the synthetic peptide chelating with a Ni ion. Different chelate structures are considered. The geometry parameters of the QM subsystem we obtained by QM/MM calculation are consistent with the available experimental results. We also perform a classical molecular dynamics (MD) simulation with the experimental parameters for the synthetic peptide adsorption on a neutral Ni(1 0 0) surface. We find that half of the His-tags are almost parallel with the substrate, which enhance the binding strength. Peeling of the peptide from the Ni substrate is simulated in the aqueous solvent and in vacuum, respectively. The critical peeling forces in the two environments are obtained. The results show that the imidazole rings are attached to the substrate more tightly than other bases in this peptide.

QM/MM and Classical Molecular Dynamics Simulation of Histidine-Tagged Peptide Immobilization on Nickel Surface

The hybrid quantum mechanics (QM) and molecular mechanics (MM) method is employed to simulate the His-tagged peptide adsorption to ionized region of nickel surface. Based on the previous experiments, the peptide interaction with one Ni ion is considered. In the QM/MM calculation, the imidazoles on the side chain of the peptide and the metal ion with several neighboring water molecules are treated as QM part calculated by "GAMESS", and the rest atoms are treated as MM part calculated by "TINKER". The integrated molecular orbital/molecular mechanics (IMOMM) method is used to deal with the QM part with the transitional metal. By using the QM/MM method, we optimize the structure of the synthetic peptide chelating with a Ni ion. Different chelate structures are considered. The geometry parameters of the QM subsystem we obtained by QM/MM calculation are consistent with the available experimental results. We also perform a classical molecular dynamics (MD) simulation with the experimental parameters for the synthetic peptide adsorption on a neutral Ni(100) surface. We find that half of the His-tags are almost parallel with the substrate, which enhance the binding strength. Peeling of the peptide from the Ni substrate is simulated in the aqueous solvent and in vacuum, respectively. The critical peeling forces in the two environments are obtained. The results show that the in-tidazole rings are attached to the substrate more tightly than other bases in this peptide.

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.

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.

Theoretical studies on conformation comparison of braid-like and triplex DNA

Based on the experimental data of scanning tunneling microscopy (STM), models of three-stranded braid-like DNAs composed by three kinds of base triplets AAA, TAT and GCA were constructed. We investigated the braid-like DNAs and their comparative tripler DNAs using a molecular mechanics method. The three strands of braid-like DNAs are proven equivalent, while those of tripler DNAs are not. The conformational energies for braid-like DNAs were found to be higher than that for tripler DNAs. Each period in one strand of braid-like DNA has 18 nucleotides, half of which are right-handed, while the other half are left-handed. Additional discussions concerning sugar puckering modes and the H-bonds are also included. (C) 1999 Elsevier Science B.V. All rights reserved.

Interface energy and its influence on interface fracture between metal and ceramic thin films in nanoscale

A theoretical model about the size-dependent interface energy between two thin films with different materials is developed by considering the chemical bonding contribution based on the thermodynamic expressions and the structure strain contribution based on the mechanical characteristics. The interface energy decreases with reducing thickness of thin films, and is determined by such available thermodynamic and mechanical parameters as the melting entropy, the melting enthalpy, the shear modulus of two materials, etc. The predicted interface energies of some metal/MgO and metal/Al2O3 interfaces based on the model are consistent with the results based on the molecular mechanics calculation. Furthermore, the interface fracture properties of Ag/MgO and Ni/Al2O3 based on the atomistic simulation are further compared with each other. The fracture strength and the toughness of the interface with the smaller structure interface energy are both found to be lower. The intrinsic relations among the interface energy, the interface strength, and the fracture toughness are discussed by introducing the related interface potential and the interface stress. The microscopic interface fracture toughness is found to equal the structure interface energy in nanoscale, and the microscopic fracture strength is proportional to the fracture toughness. (C) 2010 American Institute of Physics. [doi:10.1063/1.3501090]

Study of adsorption of methylene blue and new methylene blue in liquid-solid interface by slab optical waveguide spectroscopy

A slab optical waveguide (SOWG) has been used for study of adsorption of both methylene blue (MB) and new methylene blue (NMB) in liquid-solid interface. Adsorption characteristics of MB and NMB on both bare SOWG and silanized SOWG by octadecyltrichlorosilane (ODS) were compared. Effect of pH on adsorption on MB and NMB was investigated. Binding rate constant analysis showed that both MB and NMB on bare SOWG demonstrates larger association constants than those on ODS-SOWG. Interactions of NIB and NMB on bare SOWG and ODS-SOWG were analyzed by molecular mechanics calculation method. The binding energy change was in the following order: ENMB-bare > EMB-bare > ENMB-ODS > EMB-ODS. (c) 2004 Elsevier B.V. All rights reserved.

Effect Of Water On Brittle Fracture Of Sio2 By Molecular Dynamics Study

The influence of water on the brittle behavior of beta-cristobalite is studied by means of molecular dynamics (MD) simulation With the TTAM potential. Crack extension of mode 1 type is observed as the crack opening is filled LIP With water. The critical stress intensity factor K-lc(MD) is used to characterize the crack extension of MD simulation. The surface energy of SiO2 covered with layers of water is calculated at temperature of 300 K. Based oil the Griffith fracture criterion, the critical stress intensity factor K-lc(Griffith) is calculated, and it is in good agreement with that of MD simulation. (C) 2008 Elsevier B.V. All rights reserved.

Single-walled carbon nanotubes binding to human telomeric i-motif DNA under molecular-crowding conditions: More water molecules released

The natural occurrence of the human telomeric G-quadruplex or i-motif in vivo has not been demonstrated and the biological effects of the induction of these structures need to be clarified. Intracellular environments are highly crowded with various biomolecules and in vitro studies under molecular-crowding conditions will provide important information on how biomolecules behave in cells. Here we report that cell-mimic crowding can increase i-motif stability at acid pH and cause dehydration.