The micro-capsules used for drug delivery are fabricated using polylactic acid（PLA）,which is a biomedical material approved by the FDA.A coarse-grained model of long-chain PLA was built,and molecular dynamics（MD）s...The micro-capsules used for drug delivery are fabricated using polylactic acid（PLA）,which is a biomedical material approved by the FDA.A coarse-grained model of long-chain PLA was built,and molecular dynamics（MD）simulations of the model were performed using a MARTINI force field.Based on the nonlocal theory,the formula for the initial elastic modulus of polymers considering the nonlocal effect was derived,and the scaling law of internal characteristic length of polymers was proposed,which was used to adjust the cut-off radius in the MD simulations of PLA.The results show that the elastic modulus should be computed using nonlinear regression.The nonlocal effect has a certain influence on the simulation results of PLA.According to the scaling law,the cut-off radius was determined and applied to the MD simulations,the results of which reflect the influence of the molecular weight change on the elastic moduli of PLA,and are in agreement with the experimental outcome.展开更多
Nuclear magnetic resonance spectroscopy offers a powerful method for validation of molecular dynamics simulations as it provides information on the molecular structure and dynamics in solution. We performed 10 ns MD s...Nuclear magnetic resonance spectroscopy offers a powerful method for validation of molecular dynamics simulations as it provides information on the molecular structure and dynamics in solution. We performed 10 ns MD simulations using the CHARMM27 force field of four palindromic oligonucleotides and compared the results with experimental NOESY data using the full relaxation matrix formalism. The correlation coefficients between theoretical and experimental data for the four molecular species under study ranged from 0.82 to 0.98 confirming the high quality of the selected force field and providing a valid basis for the identification of force field imperfections. Hence, we observed an unsatisfactory treatment of deoxyribose conformational equilibrium, which resulted in the overrepresentation of the energetically favorable C3'-endo conformation in the MD trajectory. Our developed approach for force field validation based on NMR NOESY spectral data is applicable to a wide range of molecular systems and appropriate force fields.展开更多
The DNA conformational transition depends on both the DNA sequences and environment such as solvent as well as electrolyte in the solution. This paper uses the AMBER8 package to investigate the electrolyte concentrati...The DNA conformational transition depends on both the DNA sequences and environment such as solvent as well as electrolyte in the solution. This paper uses the AMBER8 package to investigate the electrolyte concentration influence on the dynamics of the A→B conformational transition of DNA duplex d（CGCGAATTCGCG）2. The results from the restrained molecular dynamics （MD） simulations indicate that the total energies of the systems for A-DNA are always higher than those for B-DNA, and that the A→B conformational transition in aqueous NaCl solution is a downhill process. The results from the unrestrained MD simulations, as judged by the average distance between the C5＇ atoms （average helical rise per ten base pair）, show that the concentrated NaC1 solution slows down the A→B conformational transition. This observation can be well understood by analyses of the difference between the counterion distributions around A-DNA and B-DNA.展开更多
The contribution of water molecules in molecular dynamics simulation (MDS) is unquestionably high, particularly for enzymatic interaction which occurred in the cytoplasmic environment. The addition of water molecules ...The contribution of water molecules in molecular dynamics simulation (MDS) is unquestionably high, particularly for enzymatic interaction which occurred in the cytoplasmic environment. The addition of water molecules to the system will surely influence different direct interaction between active site residues and substrate. We try to theoretically investigate to what extent the pathogenicity characterization will varies in different neuraminidase-sialic acid complex systems. The heating dynamics simulations were produced with and without TIP3P water molecules. The periodic boundary system was made for explicitly added TIP3P water molecules and generalized born molecular volume (GBMV) energy contribution was added for implicit solvent system. Both complexes, neuraminidasesialic acid of A/Tokyo/3/67 and A/Pennsylvania/10218/ 84, which have a different pathogenicity levels were minimized and simulated. The result shows more residues produced hydrogen bonds with substrate when water molecules were not added to the system. The binding free energies also show differences. Overall, even the values of energy differences, but an implicit solvent provides the similar result (HPAI complex has higher activity than LPAI for both systems) in characterization of pathogenic virus neuraminidase activity.展开更多
The molecular mechanics/Poisson-Boltzmann surface area(MM/PBSA) method has been widely used in predicting the binding affinity among ligands,proteins,and nucleic acids.However,the accuracy of the predicted binding ene...The molecular mechanics/Poisson-Boltzmann surface area(MM/PBSA) method has been widely used in predicting the binding affinity among ligands,proteins,and nucleic acids.However,the accuracy of the predicted binding energy by the standard MM/PBSA is not always good,especially in highly charged systems.In this work,we take the protein-nucleic acid complexes as an example,and showed that the use of screening electrostatic energy(instead of Coulomb electrostatic energy) in molecular mechanics can greatly improve the performance of MM/PBSA.In particular,the Pearson correlation coefficient of dataset Ⅱ in the modified MM/PBSA(i.e.,screening MM/PBSA) is about 0.52,much better than that(<0.33)in the standard MM/PBSA.Further,we also evaluate the effect of solute dielectric constant and salt concentration on the performance of the screening MM/PBSA.The present study highlights the potential power of the screening MM/PBSA for predicting the binding energy in highly charged bio-systems.展开更多
基金Project supported by the National Natural Science Foundation of China(no.11272360)the Natural Science Foundation of Guangdong Province(no.2014A030313793)+1 种基金the Special Program for Applied Research on Super Computation of the NSFC-Guangdong Joint Fund(the second phase)National Supercomputer Center in Guangzhou
文摘The micro-capsules used for drug delivery are fabricated using polylactic acid（PLA）,which is a biomedical material approved by the FDA.A coarse-grained model of long-chain PLA was built,and molecular dynamics（MD）simulations of the model were performed using a MARTINI force field.Based on the nonlocal theory,the formula for the initial elastic modulus of polymers considering the nonlocal effect was derived,and the scaling law of internal characteristic length of polymers was proposed,which was used to adjust the cut-off radius in the MD simulations of PLA.The results show that the elastic modulus should be computed using nonlinear regression.The nonlocal effect has a certain influence on the simulation results of PLA.According to the scaling law,the cut-off radius was determined and applied to the MD simulations,the results of which reflect the influence of the molecular weight change on the elastic moduli of PLA,and are in agreement with the experimental outcome.
文摘Nuclear magnetic resonance spectroscopy offers a powerful method for validation of molecular dynamics simulations as it provides information on the molecular structure and dynamics in solution. We performed 10 ns MD simulations using the CHARMM27 force field of four palindromic oligonucleotides and compared the results with experimental NOESY data using the full relaxation matrix formalism. The correlation coefficients between theoretical and experimental data for the four molecular species under study ranged from 0.82 to 0.98 confirming the high quality of the selected force field and providing a valid basis for the identification of force field imperfections. Hence, we observed an unsatisfactory treatment of deoxyribose conformational equilibrium, which resulted in the overrepresentation of the energetically favorable C3'-endo conformation in the MD trajectory. Our developed approach for force field validation based on NMR NOESY spectral data is applicable to a wide range of molecular systems and appropriate force fields.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.20876083 and 20736003)Specialized Research Fund for the Doctorial Program of Higher Education of China(Grant No.2007003009)
文摘The DNA conformational transition depends on both the DNA sequences and environment such as solvent as well as electrolyte in the solution. This paper uses the AMBER8 package to investigate the electrolyte concentration influence on the dynamics of the A→B conformational transition of DNA duplex d（CGCGAATTCGCG）2. The results from the restrained molecular dynamics （MD） simulations indicate that the total energies of the systems for A-DNA are always higher than those for B-DNA, and that the A→B conformational transition in aqueous NaCl solution is a downhill process. The results from the unrestrained MD simulations, as judged by the average distance between the C5＇ atoms （average helical rise per ten base pair）, show that the concentrated NaC1 solution slows down the A→B conformational transition. This observation can be well understood by analyses of the difference between the counterion distributions around A-DNA and B-DNA.
文摘The contribution of water molecules in molecular dynamics simulation (MDS) is unquestionably high, particularly for enzymatic interaction which occurred in the cytoplasmic environment. The addition of water molecules to the system will surely influence different direct interaction between active site residues and substrate. We try to theoretically investigate to what extent the pathogenicity characterization will varies in different neuraminidase-sialic acid complex systems. The heating dynamics simulations were produced with and without TIP3P water molecules. The periodic boundary system was made for explicitly added TIP3P water molecules and generalized born molecular volume (GBMV) energy contribution was added for implicit solvent system. Both complexes, neuraminidasesialic acid of A/Tokyo/3/67 and A/Pennsylvania/10218/ 84, which have a different pathogenicity levels were minimized and simulated. The result shows more residues produced hydrogen bonds with substrate when water molecules were not added to the system. The binding free energies also show differences. Overall, even the values of energy differences, but an implicit solvent provides the similar result (HPAI complex has higher activity than LPAI for both systems) in characterization of pathogenic virus neuraminidase activity.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11874045 and 11774147)。
文摘The molecular mechanics/Poisson-Boltzmann surface area(MM/PBSA) method has been widely used in predicting the binding affinity among ligands,proteins,and nucleic acids.However,the accuracy of the predicted binding energy by the standard MM/PBSA is not always good,especially in highly charged systems.In this work,we take the protein-nucleic acid complexes as an example,and showed that the use of screening electrostatic energy(instead of Coulomb electrostatic energy) in molecular mechanics can greatly improve the performance of MM/PBSA.In particular,the Pearson correlation coefficient of dataset Ⅱ in the modified MM/PBSA(i.e.,screening MM/PBSA) is about 0.52,much better than that(<0.33)in the standard MM/PBSA.Further,we also evaluate the effect of solute dielectric constant and salt concentration on the performance of the screening MM/PBSA.The present study highlights the potential power of the screening MM/PBSA for predicting the binding energy in highly charged bio-systems.
