Fast Atomic Charge Calculation for Implementation into a Polarizable Force Field and Application to an Ion Channel Protein

<div>The side chain torsion angles of TRP-41: <svg height="14.8173pt" id="M85" style="vertical-align:-5.52897pt" version="1.1" viewbox="-0.0498162 -9.28833 89.1676 14.8173" width="89.1676pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M565 434L556 448L482 423L322 189L287 373C280 410 261 448 236 448C207 448 161 415 110 338L131 317C160 355 181 374 193 374C204 374 208 366 217 330C230 279 247 185 259 116L50 -175C18 -223 29 -260 41 -278H51C97 -198 207 -33 270 57L313 -165C328 -235 349 -261 376 -261S445 -244 510 -150L487 -127C444 -178 423 -186 411 -186C402 -186 398 -177 387 -130C371 -62 349 39 332 135C390 217 524 392 565 434Z" id="g113-244"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="582" vert-adv-y="582"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,6.981,3.132)"><path d="M389 0V32C297 38 291 46 291 118V635C234 613 175 595 109 583V556L161 554C203 552 207 547 207 497V118C207 46 201 38 110 32V0H389Z" id="g50-50"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,11.927,0)"><path d="M300 -147C201 -63 143 98 143 270S200 602 300 686L282 710C136 610 70 450 70 271V270C70 89 136 -72 282 -170L300 -147Z" id="g113-41"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="346" vert-adv-y="346"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,16.425,0)"><path d="M719 650H483V622C553 619 576 606 581 560C585 532 588 487 588 394V148H585L167 650H20V622C67 618 88 610 109 585C127 562 129 557 129 484V264C129 173 126 124 123 93C118 44 92 31 35 28V0H272V28C204 32 181 44 176 95C173 124 170 172 170 264V515H172L598 -9H629V394C629 487 631 532 635 563C639 606 663 619 719 622V650Z" id="g121-76"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="743" vert-adv-y="743"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,26.084,0)"><path d="M95 130C70 130 46 113 46 88C46 72 54 64 59 64C93 55 121 33 121 -3C121 -41 93 -68 44 -88L55 -117C117 -98 186 -56 186 22C186 91 131 130 95 130Z" id="g113-45"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="228" vert-adv-y="228"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,31.228,0)"><path d="M614 175C564 76 510 21 408 21C256 21 146 149 146 336C146 488 235 629 402 629C510 629 570 586 597 480L626 488C620 541 614 582 606 638C578 643 510 665 429 665C206 665 44 527 44 316C44 157 153 -15 402 -15C474 -15 558 5 586 11C604 45 629 119 643 165L614 175Z" id="g121-65"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="665" vert-adv-y="665"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,39.873,3.132)"><path d="M558 106L533 129C501 84 476 67 464 67C446 67 435 116 414 239C457 297 508 367 550 440L539 451L484 438C458 390 429 338 404 300H402C375 407 351 451 285 451C168 451 24 312 24 153C24 56 66 -12 132 -12C206 -12 284 64 345 156H347C367 24 386 -12 422 -12C457 -12 503 7 558 106ZM338 206C272 99 215 56 173 56C140 56 116 96 116 170C116 303 194 405 255 405C304 405 324 299 338 206Z" id="g50-223"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="582" vert-adv-y="582"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,45.687,0)"><path d="M95 130C70 130 46 113 46 88C46 72 54 64 59 64C93 55 121 33 121 -3C121 -41 93 -68 44 -88L55 -117C117 -98 186 -56 186 22C186 91 131 130 95 130Z" id="g113-45"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="228" vert-adv-y="228"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,50.83,0)"><path d="M614 175C564 76 510 21 408 21C256 21 146 149 146 336C146 488 235 629 402 629C510 629 570 586 597 480L626 488C620 541 614 582 606 638C578 643 510 665 429 665C206 665 44 527 44 316C44 157 153 -15 402 -15C474 -15 558 5 586 11C604 45 629 119 643 165L614 175Z" id="g121-65"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="665" vert-adv-y="665"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,59.475,3.132)"><path d="M563 583C563 664 501 710 433 710C381 710 334 692 285 648C212 582 178 506 149 350L74 -56C49 -190 36 -222 24 -226L30 -257C52 -256 95 -247 121 -229C131 -222 134 -206 165 -21L219 303C254 510 296 665 397 665C447 665 481 626 481 572C481 512 446 457 390 420C367 405 351 400 323 397L302 349C398 338 466 285 466 197C466 117 419 45 306 45C268 45 224 59 198 74L186 50C197 20 228 -14 274 -14C318 -14 359 3 401 26C478 68 554 144 554 225C554 322 493 367 412 396C490 446 563 500 563 583Z" id="g50-224"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,65.335,0)"><path d="M95 130C70 130 46 113 46 88C46 72 54 64 59 64C93 55 121 33 121 -3C121 -41 93 -68 44 -88L55 -117C117 -98 186 -56 186 22C186 91 131 130 95 130Z" id="g113-45"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="228" vert-adv-y="228"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,70.478,0)"><path d="M614 175C564 76 510 21 408 21C256 21 146 149 146 336C146 488 235 629 402 629C510 629 570 586 597 480L626 488C620 541 614 582 606 638C578 643 510 665 429 665C206 665 44 527 44 316C44 157 153 -15 402 -15C474 -15 558 5 586 11C604 45 629 119 643 165L614 175Z" id="g121-65"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="665" vert-adv-y="665"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,79.123,3.132)"><path d="M493 373C493 420 472 451 445 451C421 451 400 435 400 410C400 404 402 400 406 394S417 371 417 350C417 259 317 128 256 55H254C260 126 256 262 238 340C219 424 194 451 163 451C128 451 78 413 24 329L52 302C89 353 106 370 121 370C129 370 137 363 147 336C186 231 196 73 184 -16C152 -88 120 -192 112 -238L130 -257C157 -253 213 -227 231 -213C231 -187 234 -77 240 -23C251 -6 277 30 309 69C424 208 493 299 493 373Z" id="g50-225"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="517" vert-adv-y="517"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,84.343,0)"><path d="M275 270C275 450 212 609 64 710L45 686C145 604 203 442 203 270S147 -63 45 -147L64 -170C213 -68 275 89 275 270Z" id="g113-42"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="346" vert-adv-y="346"></glyph.data></g></svg> and <svg height="14.8173pt" id="M86" style="vertical-align:-5.52897pt" version="1.1" viewbox="-0.0498162 -9.28833 97.8916 14.8173" width="97.8916pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M565 434L556 448L482 423L322 189L287 373C280 410 261 448 236 448C207 448 161 415 110 338L131 317C160 355 181 374 193 374C204 374 208 366 217 330C230 279 247 185 259 116L50 -175C18 -223 29 -260 41 -278H51C97 -198 207 -33 270 57L313 -165C328 -235 349 -261 376 -261S445 -244 510 -150L487 -127C444 -178 423 -186 411 -186C402 -186 398 -177 387 -130C371 -62 349 39 332 135C390 217 524 392 565 434Z" id="g113-244"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="582" vert-adv-y="582"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,6.981,3.132)"><path d="M414 144C384 79 371 75 317 75H135L276 221C367 316 408 376 408 465C408 570 327 635 237 635C179 635 131 609 100 575L42 494L67 471C94 510 138 565 205 565C277 565 321 517 321 435C321 348 258 270 195 195C146 137 88 81 33 26V0H411C423 44 433 88 446 135L414 144Z" id="g50-51"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,11.927,0)"><path d="M300 -147C201 -63 143 98 143 270S200 602 300 686L282 710C136 610 70 450 70 271V270C70 89 136 -72 282 -170L300 -147Z" id="g113-41"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="346" vert-adv-y="346"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,16.425,0)"><path d="M614 175C564 76 510 21 408 21C256 21 146 149 146 336C146 488 235 629 402 629C510 629 570 586 597 480L626 488C620 541 614 582 606 638C578 643 510 665 429 665C206 665 44 527 44 316C44 157 153 -15 402 -15C474 -15 558 5 586 11C604 45 629 119 643 165L614 175Z" id="g121-65"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="665" vert-adv-y="665"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,25.07,3.132)"><path d="M558 106L533 129C501 84 476 67 464 67C446 67 435 116 414 239C457 297 508 367 550 440L539 451L484 438C458 390 429 338 404 300H402C375 407 351 451 285 451C168 451 24 312 24 153C24 56 66 -12 132 -12C206 -12 284 64 345 156H347C367 24 386 -12 422 -12C457 -12 503 7 558 106ZM338 206C272 99 215 56 173 56C140 56 116 96 116 170C116 303 194 405 255 405C304 405 324 299 338 206Z" id="g50-223"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="582" vert-adv-y="582"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,30.885,0)"><path d="M95 130C70 130 46 113 46 88C46 72 54 64 59 64C93 55 121 33 121 -3C121 -41 93 -68 44 -88L55 -117C117 -98 186 -56 186 22C186 91 131 130 95 130Z" id="g113-45"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="228" vert-adv-y="228"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,36.028,0)"><path d="M614 175C564 76 510 21 408 21C256 21 146 149 146 336C146 488 235 629 402 629C510 629 570 586 597 480L626 488C620 541 614 582 606 638C578 643 510 665 429 665C206 665 44 527 44 316C44 157 153 -15 402 -15C474 -15 558 5 586 11C604 45 629 119 643 165L614 175Z" id="g121-65"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="665" vert-adv-y="665"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,44.673,3.132)"><path d="M563 583C563 664 501 710 433 710C381 710 334 692 285 648C212 582 178 506 149 350L74 -56C49 -190 36 -222 24 -226L30 -257C52 -256 95 -247 121 -229C131 -222 134 -206 165 -21L219 303C254 510 296 665 397 665C447 665 481 626 481 572C481 512 446 457 390 420C367 405 351 400 323 397L302 349C398 338 466 285 466 197C466 117 419 45 306 45C268 45 224 59 198 74L186 50C197 20 228 -14 274 -14C318 -14 359 3 401 26C478 68 554 144 554 225C554 322 493 367 412 396C490 446 563 500 563 583Z" id="g50-224"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,50.533,0)"><path d="M95 130C70 130 46 113 46 88C46 72 54 64 59 64C93 55 121 33 121 -3C121 -41 93 -68 44 -88L55 -117C117 -98 186 -56 186 22C186 91 131 130 95 130Z" id="g113-45"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="228" vert-adv-y="228"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,55.676,0)"><path d="M614 175C564 76 510 21 408 21C256 21 146 149 146 336C146 488 235 629 402 629C510 629 570 586 597 480L626 488C620 541 614 582 606 638C578 643 510 665 429 665C206 665 44 527 44 316C44 157 153 -15 402 -15C474 -15 558 5 586 11C604 45 629 119 643 165L614 175Z" id="g121-65"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="665" vert-adv-y="665"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,64.321,3.132)"><path d="M493 373C493 420 472 451 445 451C421 451 400 435 400 410C400 404 402 400 406 394S417 371 417 350C417 259 317 128 256 55H254C260 126 256 262 238 340C219 424 194 451 163 451C128 451 78 413 24 329L52 302C89 353 106 370 121 370C129 370 137 363 147 336C186 231 196 73 184 -16C152 -88 120 -192 112 -238L130 -257C157 -253 213 -227 231 -213C231 -187 234 -77 240 -23C251 -6 277 30 309 69C424 208 493 299 493 373Z" id="g50-225"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="517" vert-adv-y="517"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,69.541,0)"><path d="M95 130C70 130 46 113 46 88C46 72 54 64 59 64C93 55 121 33 121 -3C121 -41 93 -68 44 -88L55 -117C117 -98 186 -56 186 22C186 91 131 130 95 130Z" id="g113-45"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="228" vert-adv-y="228"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,74.684,0)"><path d="M614 175C564 76 510 21 408 21C256 21 146 149 146 336C146 488 235 629 402 629C510 629 570 586 597 480L626 488C620 541 614 582 606 638C578 643 510 665 429 665C206 665 44 527 44 316C44 157 153 -15 402 -15C474 -15 558 5 586 11C604 45 629 119 643 165L614 175Z" id="g121-65"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="665" vert-adv-y="665"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,83.329,3.132)"><path d="M501 517C491 587 426 710 310 710C240 710 175 662 175 601C175 561 201 517 250 452C219 438 187 423 157 404C88 360 24 280 24 178C24 71 90 -12 195 -12C250 -12 297 6 334 32C412 87 452 170 452 249C452 335 408 396 333 481C259 565 229 601 229 621C229 642 245 652 275 652C359 652 427 583 490 499L501 517ZM364 233C364 142 323 33 222 33C173 33 119 81 119 178C119 265 156 335 193 374C213 395 242 413 272 425C308 384 364 319 364 233Z" id="g50-226"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="525" vert-adv-y="525"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,88.108,3.132)"><path d="M389 0V32C297 38 291 46 291 118V635C234 613 175 595 109 583V556L161 554C203 552 207 547 207 497V118C207 46 201 38 110 32V0H389Z" id="g50-50"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,93.071,0)"><path d="M275 270C275 450 212 609 64 710L45 686C145 604 203 442 203 270S147 -63 45 -147L64 -170C213 -68 275 89 275 270Z" id="g113-42"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="346" vert-adv-y="346"></glyph.data></g></svg>.</div>

Journal of Chemistry

fig5

Figure 5

Figure 5: Fast Atomic Charge Calculation for Implementation into a Polarizable Force Field and Application to an Ion Channel Protein 