def test_gromacs_merge():
etoh_filename = utils.get_data_filename("chemicals/etoh/etoh.mol2")
benzene_filename = utils.get_data_filename("chemicals/benzene/benzene.mol2")
with utils.enter_temp_directory(): #Prevents creating lots of tleap/antechamber files everywhere
#Generate frcmod files, mol2 files
gaff_mol2_filename1, frcmod_filename1 = amber.run_antechamber( "etoh", etoh_filename, charge_method = None)
gaff_mol2_filename2, frcmod_filename2 = amber.run_antechamber( "benzene", benzene_filename, charge_method = None)
#Set file names
prmtop_filename1 = "./out1.prmtop"
prmtop_filename2 = "./out2.prmtop"
crd_filename1 = "./out1.inpcrd"
crd_filename2 = "./out2.inpcrd"
top_filename1 = "./out1.top"
top_filename2 = "./out2.top"
gro_filename1 = "./out1.gro"
gro_filename2 = "./out2.gro"
#Generate AMBER files
amber.run_tleap( 'etoh', gaff_mol2_filename1, frcmod_filename1, prmtop_filename1, crd_filename1 )
amber.run_tleap( 'benzene', gaff_mol2_filename2, frcmod_filename2, prmtop_filename2, crd_filename2 )
#Convert to GROMACS
utils.convert_via_acpype( "etoh", prmtop_filename1, crd_filename1, out_top = top_filename1, out_gro = gro_filename1 )
utils.convert_via_acpype( "benzene", prmtop_filename2, crd_filename2, out_top = top_filename2, out_gro = gro_filename2 )
#Merge topologies
gromacs.merge_topologies( [ top_filename1, top_filename2], './combined.top', 'combined', molecule_numbers = [1, 5], molecule_names = ['etoh', 'benzene'] )
#Test editing of molecule numbers in topology file
gromacs.change_molecules_section( './combined.top', './edited.top', ['etoh', 'benzene'], [10, 20] )
def test_gromacs_solvate():
etoh_filename = utils.get_data_filename("chemicals/etoh/etoh.mol2")
with utils.enter_temp_directory(): #Prevents creating lots of tleap/antechamber files everywhere
#Generate frcmod files, mol2 files
gaff_mol2_filename, frcmod_filename = amber.run_antechamber( "etoh", etoh_filename, charge_method = None)
#Amber setup
amber.run_tleap( 'etoh', gaff_mol2_filename, frcmod_filename, 'etoh.prmtop', 'etoh.crd' )
#GROMACS conversion
utils.convert_via_acpype( 'etoh', 'etoh.prmtop', 'etoh.crd', 'etoh.top', 'etoh.gro' )
#Solvate
gromacs.do_solvate( 'etoh.top', 'etoh.gro', 'etoh_solvated.top', 'etoh_solvated.gro', 1.2, 'dodecahedron', 'spc216', 'tip3p.itp' )
def test_acpype_conversion():
molecule_name = 'sustiva'
input_filename = utils.get_data_filename("chemicals/sustiva/sustiva.mol2")
with utils.enter_temp_directory(): # Prevents creating tons of GAFF files everywhere.
gaff_mol2_filename, frcmod_filename = amber.run_antechamber(molecule_name, input_filename, charge_method=None)
prmtop, inpcrd = amber.run_tleap(molecule_name, gaff_mol2_filename, frcmod_filename)
out_top, out_gro = utils.convert_via_acpype( molecule_name, prmtop, inpcrd )
def test_acpype_conversion():
molecule_name = 'sustiva'
input_filename = utils.get_data_filename("chemicals/sustiva/sustiva.mol2")
with utils.enter_temp_directory(
): # Prevents creating tons of GAFF files everywhere.
gaff_mol2_filename, frcmod_filename = amber.run_antechamber(
molecule_name, input_filename, charge_method=None)
prmtop, inpcrd = amber.run_tleap(molecule_name, gaff_mol2_filename,
frcmod_filename)
out_top, out_gro = utils.convert_via_acpype(molecule_name, prmtop,
inpcrd)