CN100468427C - Method for simulating protein interaction using computer - Google Patents

Method for simulating protein interaction using computer Download PDF

Info

Publication number
CN100468427C
CN100468427C CNB2007100154936A CN200710015493A CN100468427C CN 100468427 C CN100468427 C CN 100468427C CN B2007100154936 A CNB2007100154936 A CN B2007100154936A CN 200710015493 A CN200710015493 A CN 200710015493A CN 100468427 C CN100468427 C CN 100468427C
Authority
CN
China
Prior art keywords
dcd
atom
force
called
movement locus
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CNB2007100154936A
Other languages
Chinese (zh)
Other versions
CN101051335A (en
Inventor
时永香
连鹏
张楠
刘洁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University
Original Assignee
Shandong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shandong University filed Critical Shandong University
Priority to CNB2007100154936A priority Critical patent/CN100468427C/en
Publication of CN101051335A publication Critical patent/CN101051335A/en
Application granted granted Critical
Publication of CN100468427C publication Critical patent/CN100468427C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Peptides Or Proteins (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

A method of utilizing computer to simulate interaction of protein includes preparing protein structure in PDB databank to be proper study system, carrying out energy minimization treatment on said system, heating said system up to set temperature then carrying out multi-step balance treatment on said system, executing molecular dynamic simulation on said system and studying protein amino acid residual base acted as key function in atom movement race and confirmation variation as well as variation process under set conditions.

Description

A kind of method of utilizing the computer simulation protein interaction
Technical field
The present invention relates to a kind of method of utilizing the computer simulation protein interaction, relate in particular to interactional method between the transcription factor of a kind of NAMD_2.6 of utilization and VMD-1.8.5 software study protein.
Background technology
Protein interaction is the committed step of vital movement process, as the transduction of synthetic, the molecular recognition of coordinated regulation, Protein Folding and the enzyme of many transcription factors and signaling molecule, immune response etc.Conventional experimental technique adopts the assaying reaction equilibrium constant to detect bond strength between the protein, also can pass through atomic force microscope, and unimolecule is operated, and protein is separated from each other provides information for dissociation process.But be the assaying reaction equilibrium constant or operate the non-bonded interaction that all can not well disclose between protein, and the method for molecular dynamics simulation can remedy above-mentioned deficiency with atomic force microscope.Moreover, can also obtain the detailed process (movement locus of each atom) of protein interaction and the many important dynamics parameters in this process by molecular dynamics simulation.
From McCammon J.A. in 1977 etc. so far first with the method research protein of molecular dynamics, the computing method of molecular dynamics simulation develop rapidly, the yardstick of the system of simulation from the yardstick of hundreds of initial several psecs of atom to more than 260 ten thousand an atom hundreds of nanosecond of today, the object of simulation comprises protein, DNA, RNA etc.But, the process of these molecular dynamics simulations is mostly finished at mainframe computer or HLRS, and utilize NAMD_2.6_Linux-i686 software for calculation and VMD-1.8.5 analysis software on PC, to carry out molecular dynamics simulation, do not appear in the newspapers both at home and abroad at present.
Summary of the invention
At the deficiencies in the prior art, the problem to be solved in the present invention provides a kind of NAMD_2.6_Linux-i686 of utilization software for calculation and VMD-1.8.5 analysis software (http://www.ks.uiuc.edu/Research/namd/) in the interactional method of ordinary PC patrix albuminoid.
Method of the present invention, carry out as follows:
(1) obtains to contain the protein structure of one or more peptide chain from the PDB database;
(2) utilize VMD-1.8.5 software to extract the target peptide chain of desire research, it is dissolved in the rectangular parallelepiped water tank (Water Box) that fills 0.9% NaCl solution, and makes its minor increment be apart from each wall of water tank
Figure C200710015493D0003164234QIETU
, obtain peptide-salt solusion system, this system is made the structure that contains the non-hydrogen atom spatial position data, called after ionized.pdb makes the structure that contains hydrogen atom and hydrogen bonding parameter with this system simultaneously, called after ionized.psf;
(3) utilize NAMD_2.6_Linux-i686 software, adopt the Charmm27 field of force, with Steepest Descent method, under periodic boundary condition (Periodic Boundary Conditions PBC) above-mentioned peptide-salt solusion system being carried out T.T. length is that the energy minimization (Minimize) of 20 psecs (ps) is handled, and to calculate be that the centre of sphere, radius are with each atom once to per 2 femtoseconds (fs) in this process
Figure C200710015493D00031
Diameter of Spherical Volume in the Van der Waals force of other this atom of atom pair, it is that the centre of sphere, radius are with each atom once that every 4fs calculates
Figure C200710015493D00032
Diameter of Spherical Volume in the electric field force of other this atom of atom pair, calculate the movement locus that obtains each atom after finishing, called after minimize.dcd;
(4) then, fixing all atoms in the peptide chain, adopt with (3) in the identical field of force, periodic boundary condition, Van der Waals force and electric field force computing method, use the time of Langevin dynamics (Langevin Dynamics) temperature control method heating 20ps, system heats up gradually and is stabilized to 280-320K as a result, obtain the movement locus of each atom in this process, called after heat.dcd;
(5) use 100kcal/mol successively
Figure C200710015493D00041
, 50kcal/mol
Figure C200710015493D00042
, 20kcal/mol
Figure C200710015493D00043
And 10kcal/mol
Figure C200710015493D00044
Harmonious effect of contraction (Harmonic Constraint) retrain above-mentioned peptide chain, adopt with (3) in the identical field of force, periodic boundary condition, Van der Waals force and electric field force computing method, use Langevin dynamics (Langevin Dynamics) temperature control method and Nose-Hoover Langevin piston pressure control method that temperature is controlled at 300K, pressure is controlled at 1atm, difference equilibrium system 20ps, obtain the movement locus of all atoms in each step successively respectively, distinguish called after equ100.dcd, equ50.dcd, equ20.dcd and equ10.dcd successively;
(6) harmony to peptide chain retrains in the removal (5), adopt the field of force, periodic boundary condition, Van der Waals force and electric field force computing method, the temperature-control pressure-control method identical with (5), PME grid (the Particle Mesh Ewald grid) length of side is set to 1 ± 0.1 , utilizing PME long-range electric field force computing method, all hydrogen-oxygen bond distance in the fixed system carries out the molecular dynamics simulation of 2 nanoseconds (ns) to whole system, obtains the movement locus of each atom in this process, called after simulate.dcd;
(7) result of calculation minimize.dcd, heat.dcd, equ100.dcd, equ50.dcd, equ20.dcd, equ10.dcd and simulate.dcd in above-mentioned (3)-(6) are written into the movement locus of all atoms of VMD-1.8.5 software observes;
(8) the conformation change process by above-mentioned movement locus analysing protein;
(9) find out the amino acid residue that plays a crucial role in the protein conformation change procedure.
In the above-mentioned method of utilizing the computer simulation protein interaction: the system equilibrium temperature described in the step (4) is preferably 310K.
The simulation that utilizes method of the present invention to carry out protein interaction is compared with classic method has obvious superiority:
(1) to the computational accuracy height of Van der Waals force and electric field force, computer capacity is respectively 10 With 12 , and classic method mostly uses 9
(2) do not need by means of mainframe computer and HLRS, equipment requirements is low, can carry out on common PC;
(3) cpu busy percentage height, utilization factor can reach more than 95% under the situation that does not have interference;
(4) short during used machine, get final product for system about week age on the PC of double-core Pentium4 3.0G Hz processor, 1G internal memory, Redfleg Desktop 5.0 operating systems of 30,000 to 40,000 atoms;
(5) result accurately and reliably;
(6) be convenient in the life science relevant with the protein dynamics behavior and physical chemistry and field of medicaments widespread use with molecular recognition.
Description of drawings
The initial state of Fig. 1 HMG domain and POUs domain and final state conformation
A is the relative position and the conformation of the POUs domain of the HMG domain of Sox and Oct4 before the simulation among the figure, and B is both relative position and conformation (HMG domain in the above, the POUs domain below) after the simulation.
Embodiment
The invention will be further described below in conjunction with embodiment:
Embodiment 1:
(1) obtaining code from PDB database (http://www.rcsb.org/pdb/home/home.do) is the protein structure called after 1o4x.pdb of 1o4x;
(2) utilize VMD-1.8.5 software to extract the target peptide chain of desire research, it is dissolved in the rectangular parallelepiped water tank (Water Box) that fills 0.9% NaCl solution, and makes its minor increment be apart from each wall of water tank
Figure C200710015493D0005164350QIETU
, obtain peptide-salt solusion system, this system is made the structure that contains the non-hydrogen atom spatial position data, called after 1o4x_ionized.pdb makes the structure that contains hydrogen atom and hydrogen bonding parameter with this system simultaneously, called after 1o4x_ionized.psf;
(3) utilize NAMD_2.6_Linux-i686 software, adopt Charmm 27 field of forces, with Steepest Descent method, under periodic boundary condition (Periodic Boundary Conditions PBC) above-mentioned peptide-salt solusion system being carried out T.T. length is that the energy minimization (Minimize) of 20 psecs (ps) is handled, and to calculate be that the centre of sphere, radius are with each atom once to per 2 femtoseconds (fs) in this process
Figure C200710015493D0005164408QIETU
Diameter of Spherical Volume in the Van der Waals force of other this atom of atom pair, it is that the centre of sphere, radius are with each atom once that every 4fs calculates
Figure C200710015493D0005164418QIETU
Diameter of Spherical Volume in the electric field force of other this atom of atom pair, calculate the movement locus that obtains each atom after finishing, called after 1o4x_minimize.dcd;
(4) then, fixing all atoms in the peptide chain, adopt with (3) in the identical field of force, periodic boundary condition, Van der Waals force and electric field force computing method, use the time of Langevin dynamics (Langevin Dynamics) temperature control method heating 20ps, system heats up gradually and is stabilized to 310K as a result, obtain the movement locus of each atom in this process, called after 1o4x_heat.dcd;
(5) use 100kcal/mol successively
Figure C200710015493D0005164430QIETU
, 50kcal/mol , 20kcal/mol
Figure C200710015493D00053
And 10kcal/mol
Figure C200710015493D00054
Harmonious effect of contraction (Harmonic Constraint) retrain above-mentioned peptide chain, adopt with (3) in the identical field of force, periodic boundary condition, Van der Waals force and electric field force computing method, use Langevin dynamics (Langevin Dynamics) temperature control method and Nose-Hoover Langevin piston pressure control method that temperature is controlled at 300K, pressure is controlled at 1atm, difference equilibrium system 20ps, obtain the movement locus of all atoms in each step successively respectively, distinguish called after 1o4x_equ100.dcd successively, 1o4x_equ50.dcd, 1o4x_equ20.dcd and 1o4x_equ10.dcd;
(6) harmony to peptide chain retrains in the removal (5), adopt the field of force, periodic boundary condition, Van der Waals force and electric field force computing method, the temperature-control pressure-control method identical with (5), PME grid (the Particle Mesh Ewald grid) length of side is set to 1 ± 0.1
Figure C200710015493D0005164451QIETU
, utilizing PME long-range electric field force computing method, all hydrogen-oxygen bond distance in the fixed system carries out the molecular dynamics simulation of 2 nanoseconds (ns) to whole system, obtains the movement locus of each atom in this process, called after 1o4x_simulate.dcd;
(7) result of calculation 1o4x_minimize.dcd, 1o4x_heat.dcd, 1o4x_equ100.dcd, 1o4x_equ50.dcd, 1o4x_equ20.dcd, 1o4x_equ10.dcd and 1o4x_simulate.dcd in above-mentioned (3)-(6) are written into the movement locus of all atoms of VMD-1.8.5 software observes;
(8) by the conformation change process of above-mentioned movement locus analysing protein, initial state and final state conformation are seen Fig. 1;
(9) find the ASP274 of the amino acid residue HMG that plays a crucial role in the protein conformation change procedure and the LYS18 of POUs.
Embodiment 2:
(1) obtaining code from PDB database (http://www.rcsb.org/pdb/home/home.do) is the protein structure called after 1JLU.pdb of 1JLU;
(2) utilize VMD-1.8.5 software to extract the target peptide chain of desire research, it is dissolved in the rectangular parallelepiped water tank (Water Box) that fills 0.9% NaCl solution, and makes its minor increment be apart from each wall of water tank
Figure C200710015493D00055
Obtain peptide-salt solusion system, this system is made the structure that contains the non-hydrogen atom spatial position data, called after 1JLU_ionized.pdb makes the structure that contains hydrogen atom and hydrogen bonding parameter with this system simultaneously, called after 1JLU_ionized.psf;
(3) utilize NAMD_2.6_Linux-i686 software, adopt Charmm 27 field of forces, with Steepest Descent method, under periodic boundary condition (Periodic Boundary Conditions PBC) above-mentioned peptide-salt solusion system being carried out T.T. length is that the energy minimization (Minimize) of 20 psecs (ps) is handled, and to calculate be that the centre of sphere, radius are with each atom once to per 2 femtoseconds (fs) in this process
Figure C200710015493D0006164528QIETU
Diameter of Spherical Volume in the Van der Waals force of other this atom of atom pair, it is that the centre of sphere, radius are with each atom once that every 4fs calculates
Figure C200710015493D0006164536QIETU
Diameter of Spherical Volume in the electric field force of other this atom of atom pair, calculate the movement locus that obtains each atom after finishing, called after 1JLU_minimize.dcd;
(4) then, fixing all atoms in the peptide chain, adopt with (3) in the identical field of force, periodic boundary condition, Van der Waals force and electric field force computing method, use the time of Langevin dynamics (Langevin Dynamics) temperature control method heating 20ps, system heats up gradually and is stabilized to 320K as a result, obtain the movement locus of each atom in this process, called after 1JLU_heat.dcd;
(5) use 100kcal/mol successively
Figure C200710015493D00061
, 50kcal/mol
Figure C200710015493D00062
, 20kcal/mol
Figure C200710015493D00063
And 10kcal/mol Harmonious effect of contraction (Harmonic Constraint) retrain above-mentioned peptide chain, adopt with (3) in the identical field of force, periodic boundary condition, Van der Waals force and electric field force computing method, use Langevin dynamics (Langevin Dynamics) temperature control method and Nose-Hoover Langevin piston pressure control method that temperature is controlled at 300K, pressure is controlled at 1atm, difference equilibrium system 20ps, obtain the movement locus of all atoms in each step successively respectively, distinguish called after 1JLU_equ100.dcd successively, 1JLU_equ50.dcd, 1JLU_equ20.dcd and 1JLU_equ10.dcd;
(6) harmony to peptide chain retrains in the removal (5), adopt the field of force, periodic boundary condition, Van der Waals force and electric field force computing method, the temperature-control pressure-control method identical with (5), PME grid (the Particle Mesh Ewald grid) length of side is set to 1 ± 0.1
Figure C200710015493D0006164604QIETU
, utilizing PME long-range electric field force computing method, all hydrogen-oxygen bond distance in the fixed system carries out the molecular dynamics simulation of 2 nanoseconds (ns) to whole system, obtains the movement locus of each atom in this process, called after 1JLU_simulate.dcd;
(7) result of calculation 1JLU_minimize.dcd, 1JLU_heat.dcd, 1JLU_equ100.dcd, 1JLU_equ50.dcd, 1JLU_equ20.dcd, 1JLU_equ10.dcd and 1JLU_simulate.dcd in above-mentioned (3)-(6) are written into the movement locus of all atoms of VMD-1.8.5 software observes;
(8) the conformation change process by above-mentioned movement locus analysing protein;
(9) find the amino acid residue that plays a crucial role in the protein conformation change procedure.
Embodiment 3:
(1) obtaining code from PDB database (http://www.rcsb.org/pdb/home/home.do) is the protein structure called after 1HN2.pdb of 1HN2;
(2) utilize VMD-1.8.5 software to extract the target peptide chain of desire research, it is dissolved in the rectangular parallelepiped water tank (Water Box) that fills 0.9% NaCl solution, and to make its minor increment apart from each wall of water tank be 10
Figure C200710015493D0006164629QIETU
, obtain peptide-salt solusion system, this system is made the structure that contains the non-hydrogen atom spatial position data, called after 1HN2_ionized.pdb makes the structure that contains hydrogen atom and hydrogen bonding parameter with this system simultaneously, called after 1HN2_ionized.psf;
(3) utilize NAMD_2.6_Linux-i686 software, adopt Charmm 27 field of forces, with Steepest Descent method, under periodic boundary condition (Periodic Boundary Conditions PBC) above-mentioned peptide-salt solusion system being carried out T.T. length is that the energy minimization (Minimize) of 20 psecs (ps) is handled, and to calculate be that the centre of sphere, radius are with each atom once to per 2 femtoseconds (fs) in this process
Figure C200710015493D00065
Diameter of Spherical Volume in the Van der Waals force of other this atom of atom pair, it is that the centre of sphere, radius are with each atom once that every 4fs calculates
Figure C200710015493D0007164701QIETU
Diameter of Spherical Volume in the electric field force of other this atom of atom pair, calculate the movement locus that obtains each atom after finishing, called after 1HN2_minimize.dcd;
(4) then, fixing all atoms in the peptide chain, adopt with (3) in the identical field of force, periodic boundary condition, Van der Waals force and electric field force computing method, use the time of Langevin dynamics (Langevin Dynamics) temperature control method heating 20ps, system heats up gradually and is stabilized to 280K as a result, obtain the movement locus of each atom in this process, called after 1HN2_heat.dcd;
(5) use 100kcal/mol successively
Figure C200710015493D0007164717QIETU
, 50kcal/mol
Figure C200710015493D0007164728QIETU
, 20kcal/m0l
Figure C200710015493D0007164736QIETU
And 10kcal/mol
Figure C200710015493D0007164744QIETU
Harmonious effect of contraction (Harmonic Constraint) constraint peptide chain, adopt with (3) in the identical field of force, periodic boundary condition, Van der Waals force and electric field force computing method, use Langevin dynamics (Langevin Dynamics) temperature control method and Nose-Hoover Langevin piston pressure control method that temperature is controlled at 300K, pressure is controlled at 1atm, difference equilibrium system 20ps, obtain the movement locus of all atoms in each step successively respectively, distinguish called after 1HN2_equ100.dcd successively, 1HN2_equ50.dcd, 1HN2_equ20.dcd and 1HN2_equ10.dcd;
(6) harmony to peptide chain retrains in the removal (5), adopt the field of force, periodic boundary condition, Van der Waals force and electric field force computing method, the temperature-control pressure-control method identical with (5), PME grid (the Particle Mesh Ewald grid) length of side is set to 1 ± 0.1
Figure C200710015493D0007164811QIETU
, utilizing PME long-range electric field force computing method, all hydrogen-oxygen bond distance in the fixed system carries out the molecular dynamics simulation of 2 nanoseconds (ns) to whole system, obtains the movement locus of each atom in this process, called after 1HN2_simulate.dcd;
(7) result of calculation 1HN2_minimize.dcd, 1HN2_heat.dcd, 1HN2_equ100.dcd, 1HN2_equ50.dcd, 1HN2_equ20.dcd, 1HN2_equ10.dcd and 1HN2_simulate.dcd in above-mentioned (3)-(6) are written into the movement locus of all atoms of VMD-1.8.5 software observes;
(8) the conformation change process by above-mentioned movement locus analysing protein;
(9) find out the amino acid residue that plays a crucial role in the protein conformation change procedure.

Claims (2)

1. method of utilizing the computer simulation protein interaction, carry out as follows:
(1) obtains to contain the protein structure of one or more peptide chain from Protein Data Bank database;
(2) utilize VMD-1.8.5 software to extract the target peptide chain of desire research, it is dissolved in the rectangular parallelepiped water tank that fills 0.9%NaCl solution, and makes its minor increment be apart from each wall of water tank
Figure C200710015493C0002104134QIETU
Obtain peptide-salt solusion system, this system is made the structure that contains the non-hydrogen atom spatial position data, called after ionized.pdb makes the structure that contains hydrogen atom and hydrogen bonding parameter with this system simultaneously, called after ionized.psf;
(3) utilize NAMD_2.6_Linux-i686 software, adopt the Charmm27 field of force, with Steepest Descent method, under periodic boundary condition above-mentioned peptide-salt solusion system being carried out T.T. length is that the energy minimization of 20 psecs is handled, and to calculate be that the centre of sphere, radius are with each atom once to per 2 femtoseconds in this process
Figure C200710015493C0002104449QIETU
Diameter of Spherical Volume in the Van der Waals force of other this atom of atom pair, it is that the centre of sphere, radius are with each atom once that per 4 femtoseconds calculate
Figure C200710015493C0002104455QIETU
Diameter of Spherical Volume in the electric field force of other this atom of atom pair, calculate the movement locus that obtains each atom after finishing, called after minimize.dcd;
(4) then, fixing all atoms in the peptide chain, adopt with (3) in the identical field of force, periodic boundary condition, Van der Waals force and electric field force computing method, the time of using Langevin dynamics temperature control method to heat 20 psecs, system heats up gradually and is stabilized to 280-320K as a result, obtain the movement locus of each atom in this process, called after heat.dcd;
(5) use 100kcal/mol successively
Figure C200710015493C0002104503QIETU
, 50kcal/mol
Figure C200710015493C0002104508QIETU
, 20kcal/mol
Figure C200710015493C0002104511QIETU
And 10kcal/mol
Figure C200710015493C0002104515QIETU
Harmonious effect of contraction retrain above-mentioned peptide chain, adopt with (3) in the identical field of force, periodic boundary condition, Van der Waals force and electric field force computing method, use Langevin dynamics temperature control method and Nose-Hoover Langevin piston pressure control method that temperature is controlled at 300K, pressure is controlled at 1atm, difference equilibrium system 20 psecs, obtain the movement locus of all atoms in each step successively respectively, distinguish called after equ100.dcd, equ50.dcd, equ20.dcd and equ10.dcd successively;
(6) harmony to peptide chain retrains in the removal (5), adopts the field of force, periodic boundary condition, Van der Waals force and electric field force computing method, the temperature-control pressure-control method identical with (5), and the Particle Mesh Ewald grid length of side is set to 1 ± 0.1 Utilize PME long-range electric field force computing method, all hydrogen-oxygen bond distance in the fixed system carries out the molecular dynamics simulations of 2 nanoseconds to whole system, obtains the movement locus of each atom in this process, called after simulate.dcd;
(7) result of calculation minimize.dcd, heat.dcd, equ100.dcd, equ50.dcd, equ20.dcd, equ10.dcd and simulate.dcd in above-mentioned (3)-(6) are written into the movement locus of all atoms of VMD-1.8.5 software observes;
(8) the conformation change process by above-mentioned movement locus analysing protein;
(9) find out the amino acid residue that plays a crucial role in the protein conformation change procedure.
2. the method for utilizing the computer simulation protein interaction as claimed in claim 1 is characterized in that: the system equilibrium temperature described in the step (4) is 310K.
CNB2007100154936A 2007-05-11 2007-05-11 Method for simulating protein interaction using computer Expired - Fee Related CN100468427C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2007100154936A CN100468427C (en) 2007-05-11 2007-05-11 Method for simulating protein interaction using computer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2007100154936A CN100468427C (en) 2007-05-11 2007-05-11 Method for simulating protein interaction using computer

Publications (2)

Publication Number Publication Date
CN101051335A CN101051335A (en) 2007-10-10
CN100468427C true CN100468427C (en) 2009-03-11

Family

ID=38782745

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2007100154936A Expired - Fee Related CN100468427C (en) 2007-05-11 2007-05-11 Method for simulating protein interaction using computer

Country Status (1)

Country Link
CN (1) CN100468427C (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102446235A (en) * 2010-10-11 2012-05-09 中国石油化工股份有限公司 Method for simulating and calculating interaction parameters among chemical components by using computer
CN105893759B (en) * 2016-04-01 2018-08-24 南京大学 A kind of thyroid hormone replacement therapy virtual screening and its active quantitative calculation method of interference being total to regulatory factor based on nuclear receptor
CN107729717B (en) * 2017-11-03 2019-09-27 四川大学 A kind of method that computer simulation obtains g protein coupled receptor intermediate structure
CN109994150B (en) * 2019-03-12 2022-11-01 华东师范大学 Method for dominantly expressing degree of interaction satisfaction between protein pocket surface layer atoms and ligands
CN112116948B (en) * 2020-09-25 2022-06-03 山东大学 Simulation analysis method of DNA polyhedron with special branch number
CN113009828B (en) * 2021-02-23 2022-07-22 桂林电子科技大学 Anti-interference calibration method for dynamic parameters of complex mechanical system

Also Published As

Publication number Publication date
CN101051335A (en) 2007-10-10

Similar Documents

Publication Publication Date Title
CN100468427C (en) Method for simulating protein interaction using computer
Kumar et al. Hierarchical virtual screening approaches in small molecule drug discovery
Gervasio et al. Flexible docking in solution using metadynamics
Kirchmair et al. Comparative performance assessment of the conformational model generators omega and catalyst: a large-scale survey on the retrieval of protein-bound ligand conformations
Rao et al. Modern drug discovery process: An in silico approach
Kalyaanamoorthy et al. Structure-based drug design to augment hit discovery
Huang et al. PhDD: a new pharmacophore-based de novo design method of drug-like molecules combined with assessment of synthetic accessibility
Muhammed et al. Pharmacophore modeling in drug discovery: methodology and current status
Rahman et al. Use of computer in drug design and drug discovery: A review
CN100489874C (en) Method for calculating force between protein and DNA by computer simulation
Burger et al. Quantum mechanics/molecular mechanics strategies for docking pose refinement: distinguishing between binders and decoys in cytochrome C peroxidase
Mao et al. Transformer-based molecular generative model for antiviral drug design
Paul et al. Determination of the reaction coordinate for a key conformational fluctuation in human carbonic anhydrase II
Fradelos et al. The importance of going beyond Coulombic potential in embedding calculations for molecular properties: The case of Iso-G for biliverdin in protein-like environment
Liu et al. Bridging protein local structures and protein functions
Gerogiokas et al. Assessment of hydration thermodynamics at protein interfaces with grid cell theory
Yeh et al. Calculation of protein heat capacity from replica-exchange molecular dynamics simulations with different implicit solvent models
Jung et al. Thermodynamic description of interfaces applying the 2PT method on ReaxFF molecular dynamics simulations
Saha et al. Binding activity prediction of cyclin-dependent inhibitors
Ramazi et al. Computational prediction of proteins sumoylation: a review on the methods and databases
Gao et al. Molecular dynamics simulations of the Escherichia coli HPPK apo-enzyme reveal a network of conformational transitions
Lawrenz et al. A computational physics-based approach to predict unbound brain-to-plasma partition coefficient, Kp, uu
Niedbala et al. Comparative molecular surface analysis (CoMSA) for virtual combinatorial library screening of styrylquinoline HIV-1 blocking agents
Yanagisawa et al. Effective Protein–Ligand Docking Strategy via Fragment Reuse and a Proof-of-Concept Implementation
Mukherjee et al. Profile-QSAR and Surrogate AutoShim protein-family modeling of proteases

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20090311

Termination date: 20100511