Research Article  Open Access
Fast Atomic Charge Calculation for Implementation into a Polarizable Force Field and Application to an Ion Channel Protein
Abstract
Polarization of atoms plays a substantial role in molecular interactions. Class I and II force fields mostly calculate with fixed atomic charges which can cause inadequate descriptions for highly charged molecules, for example, ion channels or metalloproteins. Changes in charge distributions can be included into molecular mechanics calculations by various methods. Here, we present a very fast computational quantum mechanical method, the Bond Polarization Theory (BPT). Atomic charges are obtained via a charge calculation method that depend on the 3D structure of the system in a similar way as atomic charges of ab initio calculations. Different methods of population analysis and charge calculation methods and their dependence on the basis set were investigated. A refined parameterization yielded excellent correlation of . The method was implemented in the force field COSMOSNMR and applied to the histidinetryptophancomplex of the transmembrane domain of the M2 protein channel of influenza A virus. Our calculations show that moderate changes of side chain torsion angle and small variations of of Trp41 are necessary to switch from the inactivated into the activated state; and a rough twoside jump model of His37 is supported for proton gating in accordance with a flipping mechanism.
1. Introduction
Substantial effort has been put in the development of polarizable force fields to understand the function of biomolecules in their realistic environment, that is, biomembranes or aqueous solutions. Standard Class I and II molecular mechanics force fields, used for chemical and biological applications, describe electrostatic interactions in terms of fixed, in most cases atomcentred, partial charges. However, polarization plays a substantial role in shifting charges within real molecular systems. The interaction energies of molecular systems are considerably influenced by polarization due to molecular rearrangements that cause charges to move in the region of interaction. The neglect of polarizations is one of the most severe drawbacks of force fields using fixed atomic charges; and many attempts have been made to develop polarizable force fields (for a review see, e.g., Cieplak et al. [1] and references therein).
Partial atomic charges in molecules are not directly measurable and they are not observable in the context of quantum theory. They give a quantified description of the electron distribution in molecules in a simplified view. Atomic charges can only be determined indirectly by measuring strongly charge dependent observables, for example, the dipole moment, the electric field gradient, or the magnetic shielding. Charges induce correct tendencies and magnitudes of interactions utilized for the calculation of electrostatic interactions by applying Coulomb’s law. Therefore, the concept of partial atomic charges is extremely helpful in many fields of biology, chemistry, and physics.
Widely used force fields (and collections) like AMBER [2, 3], CHARMM [4–6], GROMACS [7–10], GROMOS [11–13], and MMX [14] use fixed atomic charges which result in inadequate descriptions for some applications, especially in case of highly charged molecules. Efforts were made to improve existing force fields with applications tested by the developers, for example, including fluctuating charges in OPLS/PROSA [15], applying induced dipole models in AMBER [16], and introducing multipole expansions as well as additional potentials in AMOEBA [17] or utilizing Drude oscillator models in CHARMM [18].
Methods, based on the combination of quantum mechanics (QM) and molecular mechanics (MM), received more attention due to the development of increasingly powerful computation options and fast data processing. Although QM/MM methods always bear the problem of truncation effects caused by the discontinuity at the borderline between the QM and MM treated section within a molecular system, QM/MM methods have been a valuable option to model chemical reactions (breaking and forming of chemical bonds as well as states far away from equilibrium) [19]. The results allow valuable insights into, for example, proton transfer reactions and transition states in biological interesting molecular systems due to the consideration of improved coupling potentials and a critical data analysis. However, quantitative comparisons to experimental and calculated data require careful analysis; see [20]. An excellent article regarding the capabilities and limitations to treat complex biomolecular processes was published by Riccardi et al. [21]. Emphasis was given to the prevention of misinterpretations when only a minimum QM region was defined; moreover, consistent treatment of the electrostatics and sufficient sampling are crucial for the success of the method.
In our approach, partial atomic charges are calculated using the Bond Polarization Theory (BPT) [22–25]. Two parameters per bond type are required to calculate partial charges from bond orbitals. These two parameters linearly enter into the BPT formula and therefore can be determined from a least square solution of a system of equations, provided partial charges are known from ab initio calculations. These two parameters have a physical meaning: One parameter is the bond partial charge of a nonpolarized bond and the other parameter is the change of the charge due to polarization. The polarization energy for each bond is calculated within the BPT framework. BPT charge calculations are easily integrated into a force field by constructing a bond orbital for each bond defined in a molecule. The capability to describe polarization within a force field model potentially serves as key to a better understanding of structural based function of biological systems [23, 26]. The concept of polarization has been also applied to fast calculation of chemical shift tensors, for nuclear magnetic resonance spectroscopic investigations, and for direct structure elucidation using them [25, 27–33].
The semiempirical treatment of charge distributions integrated into a force field improves the description of short and longrange electrostatics within a large molecular system and this should bridge the gap between QM/MM methods and force field calculations with fixed charges. Previously, we developed a force field (COSMOSNMR) [34] with polarisable BPTderived charges to describe intermolecular interactions for most organic compounds: Molecules consisting of C, O, N, and H atoms [26]. Although the BPT is not able to describe chemical reactions, the COSMOSNMR force field provides a novel approach to treat biomolecular systems inclusive polarization effects with consistent electrostatics and allows sufficient sampling due to the low computational effort; for benchmark simulations see Schneider et al. [35].
The interaction of Valinomycin with potassium ions was studied earlier [26] and offers an explanation regarding the selectivity of the ion transporter protein for K^{+} versus Na^{+}. The calculations yielded interaction energies of the cyclic depsipeptide complex that are in the order of the hydration energy; that is, Valinomycin is able to strip off the hydration sphere of a K^{+} ion by forming an interaction complex lower in energy compared to the ion in its hydration shell. The K^{+} selectivity of Valinomycin in nature provides qualitative proof of our quantitative model. The calculations pointed to polarization effects that are mainly attributed to carbonyls of the backbone of Valinomycin, which matches experimental results in related systems [36].
In this work, one of our objectives is the extension of our semiempirical charge calculation method to F, Si, P, S, and Cl containing compounds by a reparameterization. Also, some criticism of the BPT charge calculation was voiced due to the choice of a minimum basis set (STO3G) and population analysis (Mulliken) to obtain reference atomic charges for the parameterization. Thus, several alternative parameterizations of the BPT charge calculation were tested for best performance based on ab initio calculations using different definitions of partial charges and different sized basis sets as reference data.
The new method was applied to the molecular model of a proton conducting ion channel where charges and polarizations are of concern for the function of the gating mechanism. The influenza A virus contains a tetrameric bundle of proteins in its envelope that form a membrane spanning proton channel, the M2 channel. The channel conducts protons essential for viral replication [37–40]. After the virus enters a cell and becomes embedded in an acidic compartment (pH 56), protons are imported into the virion balancing the pHgradient across the membrane between lumen and cytoplasm. The resulting change in proteinprotein and proteinlipid interaction triggers the uncoating of the virus. The channel formation can be inhibited by the antiviral drugs Amantadine [41] and Rimantadine [42]. Nevertheless, the virus becomes increasingly drug resistant [43] due to mutations. Therefore there exists a priority to clearly understand the structural basis for the opening and closing mechanism of the channel in order to find new drugs against the virus and at the same time addressing drug resistant viral strains [44].
A wealth of information is available on the structure and structural details concerning the function of the M2 channel: structures from experiments have been discussed [44–58]; several solid state NMR structures of the M2 transmembrane domain (TMD) [44, 50, 52–54, 59, 60] and solution state NMR structures of M2 in micelles [46, 61, 62] as well as Xray crystal structures of M2 TMD in octylglucoside [47, 63] have been published, and there are also a broad variety of calculated structural models available [49, 57, 58]. These proteins highlight different features in structure and function which seem to be strongly dependent on the membrane environment, pH adjustment, and temperature [56]. Also, a heterogeneous model was reported showing a helixkink angle distribution [50].
One feature these structures have in common is the core of the channel: It comprises four Trp41 (indole rings) and four His37 (imidazoles) amino acids. The latter groups are either protonated or deprotonated at physiologic range between pH 5 and pH 8 [64] which is essential for the channel function/gating mechanism. His37 initiates the opening of the channel [65, 66] when it becomes positively charged through protonation. His37 residues appear to be physically involved in proton conducting, while Trp41 groups mechanically open and close the channel [67].
The timelimiting step of the channel opening is the protonation/deprotonation of His37 (imidazole side chain) proceeding on a mstosec timescale [68, 69]. Khurana et al. [49] performed MD simulations on four structures/channel models in different protonation states to show the relaxation of the His37 groups. The authors propose the channel action as a peristaltic water pump with a “sphincter” formed by the Val27 residues. There was no explanation how the protonation state of His37 is “recognized” by hydrophobic Val27 residues.
Interestingly, the M2 channel has a surprising low conductivity. Two different mechanisms of proton conduction through the channel have been suggested: (1) A continuous “waterwire” enabling proton “hopping” [70] which is very fast [57, 71] or (2) a shuttling/sphincter mechanism that directly involves the protonation/deprotonation of His37. The latter mechanism is more consistent with the relative slow conductance rate of 10^{4} ions per second [49, 69, 72].
There is evidence for a gating function of tryptophan rings [32, 73] from NMR investigations on membrane oriented channels with ^{15}N and ^{19}F labeled Trp41 groups. The NMR investigations provided fluorfluor distances for the opened and closed state. Cationπ interactions exist between the protonated imidazole ring of His37 and the indole ring of Trp41 [66]. Fluorescence investigations on M2 also demonstrated a pH dependence of the Trp41 and His37 interactions [74].
A twosite Trp41 exchange model has been introduced to explain the dynamics of the system. The coupling of Trp41 has been measured to be 2.6 Hz. However, Pérez et al. reported an error on the Karplus relation for calculating of 0.8 Hz [75], which renders any conclusion concerning the conformational state uncertain. The relative large error could be an indication for the simultaneous occurrence of mixed protonation states within a sample, causing different conformational states of the side chain residues His37 and Trp41 at the active site, the “inner” gate. Therefore, the precise conformational states of the channel core are a key to understand the functional mechanics of the M2 channel.
In this paper, we provide evidence for a plausible opening mechanism for the M2 channel: The tryptophan gate of the channel is pushed open by the histidines as consequence of their conformational change triggered by protonation. A conformational search of the gating model at different states offers the energetic rationale for a detailed opening mechanism by including all mutual polarizations in the calculation of the electrostatic energy.
2. Materials and Methods
2.1. Bond Polarization Theory for Partial Atomic Charges
The BPT (Bond Polarization Theory) is the starting point for the calculation of atomic charges, describing the dependence of expectation values from polarization [76]. The definition of a local charge operator (site ) is a prerequisite to calculate atomic charges :The expectation value is expressed as [23, 24, 26]The general atomic charge operator is defined by the difference between the nuclear charge and the electron population. The latter depends strongly on the kind of atomic orbitals, , associated with atom . The molecular wave function is build up as a configuration interaction (CI) series of Slater determinants of bond orbitals (BOs). The BOs are linear combinations of hybrid atomic orbitals (AOs) and . The first sum in (1b) counts all bonds of atom denoted by . The variable refers to the occupation number of th bonds; and is the corresponding bond partial charge. The parameter quantifies the change of the charge due to bond polarization. The outer sum adds the contribution of each bond of to the total charge, consisting of the bond partial charge and the correction due to polarization caused by the surrounding atoms. Therefore, the inner sum, accounting for all polarizing atoms (denoted as system C), has to exclude atom itself. The bond partial charges and their variation caused by bond polarization are the parameters that have to be determined by a parametrization procedure. This is achieved with known charges . Two parameters need to be determined for each type of bond/respective combination of elements . A sufficient large set of charges is necessary to obtain an overdetermined system of equations (see (1b)). The least square solution of this set of equations will provide the parameters and .
The occupation number of bonds is introduced to treat molecular systems with conjugated bonds between two atoms and . An empirical relationship between the contraction of an ideal single bond and the valence of the bond was applied to estimate the bond occupation numbers (see (1b)) between two atoms and via the corrected/real bond length ():The ideal single bond lengths are obtained from a parameterization of O’Keeffe and Brese [77].
BPT charges are calculated from all other atomic charges of the molecular system (see (1b)). The resulting system of equations has the dimension of the number of atoms in the molecular system.
2.2. Parameterization
A set of 163 molecules consisting of H, C, N, O, F, Si, P, S, and Cl atoms (with 25 bond types: HC, HN, HO, CF, CN, C=N, CC, C=C, CO, C=O, ClC, PO, P=O, SiH, SiC, SiO, SiCl, SH, SC, S=C, SO, S=O, SS, ZnO, and ZnN) was selected to calibrate the BPT charge parameters (22 numbers) and (25 parameters) with 25 valences and global scaling factor for the electrostatic energy.
Based on experimentally available structural data [78], the initial molecular geometries were optimized with GAUSSIAN 98 [79] using hybrid functional option (DFT/B3LYP) for a 631G(d,p) basis set. The distances were checked against the original experimental data from literature [79, 80]. No major deviations in bond distances (≥0.02 Å) and bond angles (≥2°) were observed.
The BPT method was parameterized using three different atomic charge models: Mulliken population analysis (MPA), charges derived from electrostatic potentials (ESP), and charges from a natural population analysis (NPA). The ESP method does not fit into the context of methods for population analysis since these charges are fitted to represent the electrostatic potential at the molecular surface. It is possible to derive the partial charge as expectation value of an operator (see (1a)) mathematically correct only for the MPA and NPA method. The basis set dependence of the charges was investigated using 11 different basis sets for each charge model (MPA, ESP, and NPA) (see Supplementary Material available online at http://dx.doi.org/10.1155/2015/908204). Table 1 shows the correlation coefficients of the parameterization of the BPT charge calculation of six selected charge calculation models; only results of the best and worst performing basis sets are listed. The overall best performance has been obtained with NPA/631G(d,p).

For a special application, zinc was introduced as new element to prove the general applicability, the quality of the correlation, and the stability of the results. The new training set for the parameterization consisted of the same 163 molecules that was used earlier plus 12 additional zinc compounds. Again, the NPA/631G(d,p) basis set yielded the highest correlation coefficient compared to different charge models (MPA/STO6G and ESP/321G) that showed also high correlation coefficients (see Table 1) for the 175 compounds. Figure 1 depicts the correlation of the BPT charge calculation against NPA charges from DFT calculations with the 631G(d,p) basis set for 175 compounds of the calibration set. The ESP method did not perform such well as the initial parameterization.
A correlation coefficient of 0.9961 was obtained between BPT/NPA and DFT/NPA, and the standard deviation was 0.05 e. The inclusion of rather exotic zinc compounds did not drastically change correlation coefficient and kept the BPT parameters nearly unchanged. This demonstrates that (1) the BPT is able to reproduce atomic charges from ab initio calculations in a very satisfactory manner and (2) the interpretation of the BPT parameters regarding the polarization is physically correct.
One major point to be addressed is the basis set dependence of the calculated atomic charges. The correlations of BPT/NPA parameters for different compact basis sets (321G and 631G(d,p)) are shown in Figure 2. The correlation of the bond partial charges yields a correlation coefficient of (SD = 0.03 e), while the polarization parameters show a value of 0.9969 and a standard deviation of 0.26 e/H. This indicates that sufficient parameter stability is maintained regarding the usage of different compact basis sets (no polarization functions are used); no basis set dependence of the calculated atomic charges was found in our tests.
(a)
(b)
In a next step, the BPT/NPA method was compared with DFT/NPA results of a pseudopeptide zinc complex that was not included in the training set, consisting of 64 atoms (H, C, N, O, and Zn). The correlation of the calculated atomic charges for the zinc complex (correlation coefficient , standard deviation SD = 0.07 e) is shown in Figure 3.
A detailed comparison of all correlation coefficients regarding different charge models, basis set sizes, and numbers of the corresponding BPT parameters is given in Supplementary Material.
One major focus of this work was to improve the description of intermolecular interactions for biomolecular applications. Intermolecular interactions are mainly described by the sum of Coulomb and Van der Waals energy terms within the COSMOSNMR force field with fluctuating charges (for details see [26]). Since interacting electron distributions are replaced by potential functions, this simple approach will lead to a systematic bias in interaction energies.
A scaling factor was introduced into the potential energy term in order to finetune the interaction energies. This scaling factor was determined iteratively by comparing interaction energies calculated with the COSMOSNMR force field with interaction energies from the web service based force field evaluation suite by Halgren [81–83]. Among other types of data, this data base contains ab initio calculated intermolecular interaction energies of small dimeric molecules. 66 hydrogenbonded dimers were optimized with the COSMOSNMR force field [25]. The intermolecular interaction energies were calculated as the difference between the energy of the dimer and the sum of the energies of the two monomers. The optimized structures were compared to the quantum mechanically optimized geometries from Halgren [81–83] by calculating the RMS deviations for the whole structures as well as for the hydrogen bond structures (distances and angles). A scaling factor was calculated as arithmetic mean of the scaling factors of each interaction pair. The whole procedure was repeated until the difference of the scaling factor per iteration was smaller than 10^{−5}, which is 10^{3} times below the error limit of the calculated BPT charges. The calculated scaling factor for the Coulomb energy of the COSMOSNMR force field is 0.87245. A comparison of the mean deviations for the interaction energies and hydrogen bonding geometries from several force fields calculations including our own developed COSMOS force field is given in Table 2.
 
Each of the force fields handles all 66 hydrogenbonded dimers except for (i) CVFF, for which the Cerius2 implementation could not perform the atom typing for imidazole, (ii) MM3, which lacks parameters for 12 of the dimers, and (iii) CHARMM 22, which lacked parameters for about half the systems even with extensions made locally. Some problems occurred with the interaction energies of two cyclic dimers that were not included in the comparison (oxalic acid/water, cyclic ammonia dimer) and three charged dimers (formamidinium/water, formaldehyde iminium/water, and bidentate/water) that show larger deviations of the interaction energies (≥2 kcal/mol). 
The COSMOSNMR force field shows slight overestimation of the interaction energies compared to the quantum mechanical data, whereas the hydrogen bond geometries show relative small deviations from the quantum mechanical calculated geometries; the distances are too short (0.04 Å) and the angles are larger (10.6°) on average than the HF/ calculated data.
2.3. Program Implementation
The routines for the BPT charge calculations were integrated into the COSMOSNMR force field to enable the recalculation of the charges at every step of a molecular mechanics calculation. A graphical user interface for Windows (Win32) was implemented that integrates COSMOSNMR into the general modelling and graphics program COSMOS [34].
All calculations can be performed using the COSMOSbackend program (C^{++}) that was compiled for several operating systems as, for instance, UNIX, LINUX, AIX, and Win32. The COSMOSbackend is controlled by command line parameters and a project file . The project file (ACII) can be edited or alternatively generated using the GUI program COSMOS. A MPI parallelized version is available [35] for MD calculations.
The most time consuming part of the charge calculation is the evaluation of the integrals over bond orbitals describing the bond polarization energies (see (1a) and (1b)). The expressions for the integrals that are needed as coefficients of the system of equations are given in the supporting material. These integrals are used to build up a matrix of the dimension equal to the number of atoms . The computational time for setting up all charge equations is proportional to the number of bonds multiplied by . Charge calculation within the BPT framework means solving this set of linear equations for which the number of floating point operations is proportional to . A time profile for the routines was generated and showed that the solution of the equation is about ten times faster than the calculation of all integrals. Highly optimized standard packages (LAPACK, see, e.g., [35]) are used to speed up the matrix inversion.
3. Application and Results
The membrane protein M2 of influenza A virus consists of 97 residues including a 24residue Nterminal and a 54residue Cterminal segment. The transmembrane domain (TMD) consists of 19 residues that form an αhelical secondary structure. The αhelices associate with a tetrameric channel forming a pore in a lipid bilayer as illustrated in Figures 4 and 6 [74, 84, 85].
(a)
(b)
Investigations under different conditions indicate that the tilt angle of the helix backbone orientation varies between 15° and 45° [46, 47, 86–93]. Channels without Amantadine show much less variation of the tilt angle. A tilt angle between 32° and 38° was detected for the channel in lipid DMPC. The tertiary structure seems not to change significantly during the switch between an open and closed state. Thus, we use the same backbone structure for our calculations (DMPC/water box) at different pH. A 25residue peptide (SSDPLVVAASIIGILH^{37}LILW^{41}ILDRL) containing the hydrophobic TMD (P^{25}L^{46}) was investigated. This TMD model and its equivalent mutants are widely used for structural investigations [52]. Here, we focused on the pHdependent conformational change of the His37 and Trp41 side chains inside the channel pore.
The COSMOSNMR force field is especially suited for calculations of charged systems since all mutual polarizations are included into the electrostatic energy. The PDB structure 1NYJ [52] was used as initial model of M2 TMD. The protonation states determined by Hu et al. [64] were applied for the His37tetramer. The charged groups were modeled according to the pH values for the opened and closed channel configuration. Figure 5 illustrates Trp41; the side chain torsion angles of His37 are named accordingly.
(a)
(b)
(c)
COSMOSNMR force field energies were calculated for different His37 and Trp41 side chain conformations at pH 5 and mapped onto the total force field energy of the channel.
Four His37 have an average total charge of +1.3 at pH 8.0 and +3.0 at pH 5.0. The charged tetrameric M2 TMD structure was placed into an equilibrated box with a size of 60^{3} Ǻ^{3}, containing 128 DMPC lipid and 3000 water molecules; see Figure 6(c). All charges and energies, including the coordinatedependent electrostatic interactions, were calculated for conformations of the four His37 and Trp41 residues under these conditions.
Figures 7 and 8 show conformational energy maps at pH 8 (closed channel) and at pH 5 (open channel), respectively. Panels (a), (d), and (g) picture the calculated total energies as a function of the His37 (a) and Trp41 ((d) and (g)) side chain torsion angles and , as well as and in different configurations. The dark areas represent allowed rotamers within the channel in a DMPC/water bilayer environment. Figures 7(b), 7(e), 7(h), 8(b), 8(e) and 8(h) depict the calculated Van der Waals and (c), (f), and (i) the electrostatic energy contour plots for the different combinations of side chain torsion angles, and .
(a) His37
(b) His37
(c) His37
(d) Trp41
(e) Trp41
(f) Trp41
(g) Trp41
(h) Trp41
(i) Trp41
(a) His37
(b) His37
(c) His37
(d) Trp41
(e) Trp41
(f) Trp41
(g) Trp41
(h) Trp41
(i) Trp41
The calculations allow only an estimate of forbidden conformational areas for the side chains of the gating residues His37 and Trp41 and provide reasons for the conformational restrictions that are concluded from the energy decomposition analysis, since the structures are not geometry optimized.
An analysis of the calculations for the inactivated/closed state at pH 8 yields the following results: (1) Dihedral angels of Trp41 were fixed according to experimental results (, ) = (−100 ± 10°, +110 ± 10°) by Witter et al. [32]: (i) There is a broad forbidden area (between −110° and 60°) for the angles of His37 due to Coulomb interaction (see Figure 7(c)). (ii) Van der Waals interactions cause two main allowed conformational areas (see Figure 7(b)). (iii) From these two conformational areas each can be split into two rough torsion angles combinations of His37: (, ) = (−175°, +50°)/(+120°, +110°) and (−175°, −150°)/(+110°, −65°) which can be considered as rough energetically allowed/preferred conformations. The angle distribution around these values is between 20° and 50°. Transitions between these two regions possibly take place by 180° flips of which correspond to flipping mechanism [72]. (2) Two representing angles: (, ) = (+120°, +110°) and (−180°, −180°) of His37 were used for the analysis of the torsion angles of the mechanically gating Trp41 side chains. The resulting angle distributions for Trp41 are shown in Figures 7(d)–7(i). (i) No major restrictions of the allowed conformational areas are attributed to the Coulomb interactions of Trp41 (see widespread dark patches, Figures 7(f) and 7(i)). (ii) Van der Waals interactions cause major forbidden areas in the conformational space due to the bulky side chains of Trp41 residues. (iii) Torsion angles of (, ) = (−115° ± 50°, +130° ± 30°) are identified as energetically lowest conformations, which are close to the initial experimental values for the Trp41 torsion angle calculation; see above: (, ) = (−100° ± 10°, +110° ± 10°) [32]. Nevertheless view values with higher energies might be possible, for example, (−100°, −100°), depending on His37 conformation. The resulting structure of a closed channel configuration is pictured in Figure 9.
(a)
(b)
A similar analysis of the calculations for the activated/open state at pH 5 yields the following results: The panels in Figure 8 show similar energy maps as panels in Figure 7 that lead to corresponding conclusions regarding forbidden conformational areas for the side chains of the gating residues His37 and Trp41. 1. His37 was analyzed for fixed dihedral angels of Trp41 according to experimental results (, ) = (−50° ± 10°, +115° ± 10°) [32]; see Figures 8(a), 8(b), and 8(c). (i) Again, the Coulomb interactions cause a large forbidden conformational area for the angles of His37 (panel (c)). (ii) There are mostly only two areas for the torsion angles of His37, (, ) = (+95°, +130°) and (+100°, −50°), that keep the side chain locked in the open position of the channel due to Van der Waals interaction of the protonated His37 side chains. The angle distribution is roughly 20°. The conformational difference is just about a 180° flip. These two dihedral angles of His37, (+95°, +130°) and (+100°, −50°), were used for the conformational analysis of the side chain angles of Trp41; see Figures 8(d)–8(i). (i) Again, restrictions of the conformational space are mainly caused by the bulky side chains of Trp41 due to Van der Waals interaction that lock the side chain. The position with the lowest energy can be found around the initial experimental open state position (, ) = (−70° ± 20, +110° ± 20). Still other conformations at higher energies are valid. During proton gating molecules of water (not considered in these calculations) interacting with His37 side chains will restrict Trp41 to such open configuration, presented in Figure 9.
Figures 7(a), 7(d), 7(g), 8(a), 8(d), and 8(g) show the sum of the electrostatic and Van der Waals energies. The total energy is dominated by Van der Waals interactions that are altered by the electrostatic contributions.
Four Trp41 and four His37 residues participate actively in the gating of the channel structure. The function of the Trp41 residues is a mechanical closing and opening mechanism of the channel. A change in pH causes the protonation/deprotonation of His37 and causes a structural response of His37 and Trp41 driven by electrostatic and Van der Waals interactions between these residues. We conclude from our calculations that the bulky Trp41 side chains enhance the steric effect of the His37 residues that are directly involved in the H^{+} gating. In the open state the conformational space of His37 is more restricted such that they are locked for gating. It also implies the possibility of a twoside jump model which supports the function of the imidazole rings flipping around by almost 180°, thereby able to perform the “turning step” of passing the protons through the channel. At the closed state the His37 side chains are locked due to hydrogen bridges and are not accessible for H^{+} gating. When the channel opens, the Trp41 (folded back into the sides of the channel) are exposing the His side chains that only now are able to bind the hydrogen ions for passing through the open channel. That is in accordance with the flipping mechanism [72].
4. Discussion
BPT procedures for atomic charge calculation were integrated into traditional molecular mechanic force fields. Thus, a more detailed description of Coulomb interactions becomes feasible within the COSMOSNMR force field without defining special QM regions. A careful parameterization of the BPT procedure allows the calculation of atomic charges virtually on ab initio level for very large molecular systems. The charges depend on all mutual polarizations and provide with a correlation coefficient of 0.9967 improved accesses to electrostatic interaction energies.
The COSMOSNMR force field was applied to the proton channel of the influenza A virus. A detailed structurefunction relation of the channel and its gating residues was obtained and attributed to Van der Waals and Coulomb interactions of the gating residues His37 and Trp41. Experimentally derived torsion angles of Trp41 at pH 8 and pH 5 were theoretically approved. The presented states of the virus channel are not single valued or very fixed but represent preferred conformational distributions. His37 exists in mainly two different conformational distributions which can be transformed by 180° flips around . His37 are pointed to the centre of the channel at the open state. At the closed state they are locked via hydrogen bridges to the channel’s side. To change Trp41 from the inactivated into the activated state, only moderate changes of torsion angle are necessary. With our analysis we confirm former results obtained with ^{15}N, ^{19}FNMR [32, 73] and the gating function of His37 for proton conduction related to a His37 flipping [72].
Supporting Information
All coordinates of used structures within this publication are available from the authors upon request. The COSMOSbackend can be also obtained from the authors. The program is available without charge for scientific applications (request to Ulrich Sternberg: postmaster@COSMOSsoftware.de).
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors gratefully appreciate European Social Fund, Estonian Science Foundation (ESF), Tallinn University of Technology, for funding project MTT68. They highly value the financial and infrastructural support from Karlsruhe Institute of Technology (KIT). Also, parts of this work were funded by project MO 923/21; the authors thank the German Science Foundation (DFG) for their financial support.
Supplementary Materials
Details related to the Bond Polarization Theory, charge model evaluation, charge calculations, the full spectrum of parameterizations with 11 different basis sets and complete set of tables with all parameters are provided.
References
 P. Cieplak, F.Y. Dupradeau, Y. Duan, and J. Wang, “Polarization effects in molecular mechanical force fields,” Journal of Physics Condensed Matter, vol. 21, no. 33, Article ID 333102, 2009. View at: Publisher Site  Google Scholar
 P. K. Weiner and P. A. Kollman, “AMBER: assisted model building with energy refinement. A general program for modeling molecules and their interactions,” Journal of Computational Chemistry, vol. 2, no. 3, pp. 287–303, 1981. View at: Publisher Site  Google Scholar
 W. D. Cornell, P. Cieplak, C. I. Bayly et al., “A second generation force field for the simulation of proteins, nucleic acids, and organic molecules,” Journal of the American Chemical Society, vol. 117, no. 19, pp. 5179–5197, 1995. View at: Publisher Site  Google Scholar
 B. R. Brooks, R. E. Bruccoleri, B. D. Olafson, D. J. States, S. Swaminathan, and M. Karplus, “CHARMM: a program for macromolecular energy, minimization, and dynamics calculations,” Journal of Computational Chemistry, vol. 4, no. 2, pp. 187–217, 1983. View at: Publisher Site  Google Scholar
 F. A. Momany and R. Rone, “Validation of the general purpose QUANTA 3.2/CHARMm force field,” Journal of Computational Chemistry, vol. 13, no. 7, pp. 888–900, 1992. View at: Publisher Site  Google Scholar
 F. A. Momany, R. Rone, and H. Kunz, “Quanta3.3/Charm22 parameters,” Abstracts of Papers of the American Chemical Society, vol. 204, p. 42, 1992. View at: Google Scholar
 H. J. C. Berendsen, D. van der Spoel, and R. van Drunen, “GROMACS: a messagepassing parallel molecular dynamics implementation,” Computer Physics Communications, vol. 91, no. 1–3, pp. 43–56, 1995. View at: Publisher Site  Google Scholar
 E. Lindahl, B. Hess, and D. van der Spoel, “GROMACS 3.0: a package for molecular simulation and trajectory analysis,” Journal of Molecular Modeling, vol. 7, no. 8, pp. 306–317, 2001. View at: Publisher Site  Google Scholar
 D. van der Spoel, D. A. R. van Buuren, E. Apol et al., Gromacs User Manual, Version 3.1.1, Department of Biophysical Chemistry, University of Groningen, Groningen, The Netherlands, 2002.
 W. F. Van Gunsteren and H. J. C. Berendsen, “Computer simulation of molecular dynamics: methodology, applications, and perspectives in chemistry,” Angewandte Chemie International Edition in English, vol. 29, no. 9, pp. 992–1023, 1990. View at: Publisher Site  Google Scholar
 M. Christen, P. H. Hünenberger, D. Bakowies et al., “The GROMOS software for biomolecular simulation: GROMOS05,” Journal of Computational Chemistry, vol. 26, no. 16, pp. 1719–1751, 2005. View at: Publisher Site  Google Scholar
 W. R. P. Scott, P. H. Hünenberger, I. G. Tironi et al., “The GROMOS biomolecular simulation program package,” Journal of Physical Chemistry A, vol. 103, no. 19, pp. 3596–3607, 1999. View at: Publisher Site  Google Scholar
 W. F. Van Gunsteren, S. R. Billeter, A. A. Eising et al., Biomolecular Simulation: The GROMOS96 Manual and User Guide, Vdf Hochschulverlag AG an der ETH Zürich, Zürich, Switzerland, 1996.
 N. L. Allinger, K. Chen, and J.H. Lii, “An improved force field (MM4) for saturated hydrocarbons,” Journal of Computational Chemistry, vol. 17, no. 56, pp. 642–668, 1996. View at: Publisher Site  Google Scholar
 J. L. Banks, G. A. Kaminski, R. Zhou, D. T. Mainz, B. J. Berne, and R. A. Friesner, “Parametrizing a polarizable force field from ab initio data. I. The fluctuating point charge model,” The Journal of Chemical Physics, vol. 110, no. 2, article 741, 1999. View at: Publisher Site  Google Scholar
 P. Cieplak, J. Caldwell, and P. Kollman, “Molecular mechanical models for organic and biological systems going beyond the atom centered two body additive approximation: aqueous solution free energies of methanol and Nmethyl acetamide, nucleic acid base, and amide hydrogen bonding and chloroform/water partition coefficients of the nucleic acid bases,” Journal of Computational Chemistry, vol. 22, no. 10, pp. 1048–1057, 2001. View at: Publisher Site  Google Scholar
 P. Y. Ren and J. W. Ponder, “Polarizable atomic multipole water model for molecular mechanics simulation,” Journal of Physical Chemistry B, vol. 107, no. 24, pp. 5933–5947, 2003. View at: Google Scholar
 S. Patel, A. D. Mackerell Jr., and C. L. Brooks III, “CHARMM fluctuating charge force field for proteins: II protein/solvent properties from molecular dynamics simulations using a nonadditive electrostatic model,” Journal of Computational Chemistry, vol. 25, no. 12, pp. 1504–1514, 2004. View at: Publisher Site  Google Scholar
 J. L. Gao and D. G. Truhlar, “Quantum mechanical methods for enzyme kinetics,” Annual Review of Physical Chemistry, vol. 53, pp. 467–505, 2002. View at: Publisher Site  Google Scholar
 R. A. Friesner and V. Guallar, “Ab initio quantum chemical and mixed quantum mechanics/molecular mechanics (QM/MM) methods for studying enzymatic catalysis,” Annual Review of Physical Chemistry, vol. 56, pp. 389–427, 2005. View at: Publisher Site  Google Scholar
 D. Riccardi, P. Schaefer, Y. Yang et al., “Development of effective quantum mechanical/molecular mechanical (QM/MM) methods for complex biological processes,” Journal of Physical Chemistry B, vol. 110, no. 13, pp. 6458–6469, 2006. View at: Publisher Site  Google Scholar
 U. Sternberg, “Theory of the influence of the second coordination sphere on the chemical shift,” Molecular Physics, vol. 63, no. 2, p. 249, 1988. View at: Google Scholar
 U. Sternberg, F. T. Koch, and M. Möllhoff, “New approach to the semiempirical calculation of atomic charges for polypeptides and large molecular systems,” Journal of Computational Chemistry, vol. 15, no. 5, pp. 524–531, 1994. View at: Publisher Site  Google Scholar
 U. Sternberg, F.T. Kock, M. Bräuer, M. Kunert, and E. Anders, “Molecular mechanics for zinc complexes with fluctuating atomic charges,” Journal of Molecular Modeling, vol. 7, no. 4, pp. 54–64, 2001. View at: Publisher Site  Google Scholar
 R. Witter, “Structure elucidation with the COSMOSNMR force field,” in PhysikalischAstronomischeFakultät, FriedrichSchillerUniversität, Jena, Germany, 2003. View at: Google Scholar
 M. Möllhoff and U. Sternberg, “Molecular mechanics with fluctuating atomic charges—a new force field with a semiempirical charge calculation,” Journal of Molecular Modeling, vol. 7, no. 4, pp. 90–102, 2001. View at: Publisher Site  Google Scholar
 R. Witter, W. Priess, and U. Sternberg, “Chemical shift driven geometry optimization,” Journal of Computational Chemistry, vol. 23, no. 2, pp. 298–305, 2002. View at: Publisher Site  Google Scholar
 R. Witter, L. Seyfart, G. Greiner et al., “Structure determination of a pseudotripeptide zinc complex with the COSMOSNMR force field and DFT methods,” Journal of Biomolecular NMR, vol. 24, no. 4, pp. 277–289, 2002. View at: Publisher Site  Google Scholar
 U. Sternberg, R. Witter, and A. S. Ulrich, “3D structure elucidation using NMR chemical shifts,” Annual Reports on NMR Spectroscopy, vol. 52, pp. 53–104, 2004. View at: Publisher Site  Google Scholar
 R. Witter, U. Sternberg, S. Hesse, T. Kondo, F.T. Koch, and A. S. Ulrich, “^{13}C chemical shift constrained crystal structure refinement of cellulose I_{α} and its verification by NMR anisotropy experiments,” Macromolecules, vol. 39, no. 18, pp. 6125–6132, 2006. View at: Publisher Site  Google Scholar
 U. Sternberg, R. Witter, and A. S. Ulrich, “Allatom molecular dynamics simulations using orientational constraints from anisotropic NMR samples,” Journal of Biomolecular NMR, vol. 38, no. 1, pp. 23–39, 2007. View at: Publisher Site  Google Scholar
 R. Witter, F. Nozirov, U. Sternberg, T. A. Cross, A. S. Ulrich, and R. Fu, “Solidstate ^{19}F NMR spectroscopy reveals that Trp_{41} participates in the gating mechanism of the M2 proton channel of influenza A virus,” Journal of the American Chemical Society, vol. 130, no. 3, pp. 918–924, 2008. View at: Publisher Site  Google Scholar
 U. Sternberg, M. Klipfel, S. L. Grage, R. Witter, and A. S. Ulrich, “Calculation of fluorine chemical shift tensors for the interpretation of oriented ^{19}FNMR spectra of gramicidin A in membranes,” Physical Chemistry Chemical Physics, vol. 11, no. 32, pp. 7048–7060, 2009. View at: Publisher Site  Google Scholar
 U. Sternberg, F.T. Koch, and P. Losso, COSMOS: Computer Simulation of Molecular Structures, 2006, http://www.cosmossoftware.de.
 O. Schneider, U. Sternberg, and A. Vondrous, Parallel Molecular Dynamics Simulation Using a Hybrid Force Field with Quantum Chemical Charge Calculations, Nova Science, New York, NY, USA, 2009.
 D. A. Doyle, J. M. Cabral, R. A. Pfuetzner et al., “The structure of the potassium channel: molecular basis of K^{+} conduction and selectivity,” Science, vol. 280, no. 5360, pp. 69–77, 1998. View at: Publisher Site  Google Scholar
 N. Agmon, “The Grotthuss mechanism,” Chemical Physics Letters, vol. 244, no. 56, pp. 456–462, 1995. View at: Publisher Site  Google Scholar
 L. H. Pinto, L. J. Holsinger, and R. A. Lamb, “Influenzavirus M2 protein has ion channel activity,” Cell, vol. 69, no. 3, pp. 517–528, 1992. View at: Publisher Site  Google Scholar
 L. H. Pinto, L. J. Holsinger, and R. A. Lamb, “Influenza virus M2 protein has ion channel activity,” Cell, vol. 69, no. 3, pp. 517–528, 1992. View at: Publisher Site  Google Scholar
 A. Okada, T. Miura, and H. Takeuchi, “Protonation of histidine and histidine−tryptophan interaction in the activation of the M2 ion channel from influenza A virus,” Biochemistry, vol. 40, no. 20, pp. 6053–6060, 2001. View at: Publisher Site  Google Scholar
 C. Wang, K. Takeuchi, L. H. Pinto, and R. A. Lamb, “Ion channel activity of influenza A virus M2 protein: characterization of the amantadine block,” Journal of Virology, vol. 67, no. 9, pp. 5585–5594, 1993. View at: Google Scholar
 R. Dolin, R. C. Reichman, H. P. Madore, R. Maynard, P. N. Linton, and J. WebberJones, “A controlled trial of amantadine and rimantadine in the prophylaxis of influenza A infection,” The New England Journal of Medicine, vol. 307, no. 10, pp. 580–584, 1982. View at: Publisher Site  Google Scholar
 T. Jefferson, V. Demicheli, D. Rivetti, M. Jones, C. Di Pietrantonj, and A. Rivetti, “Antivirals for influenza in healthy adults: systematic review,” The Lancet, vol. 367, no. 9507, pp. 303–313, 2006. View at: Publisher Site  Google Scholar
 S. D. Cady, K. SchmidtRohr, J. Wang, C. S. Soto, W. F. Degrado, and M. Hong, “Structure of the amantadine binding site of influenza: M2 proton channels in lipid bilayers,” Nature, vol. 463, no. 7281, pp. 689–692, 2010. View at: Publisher Site  Google Scholar
 C. Miller, “Ion channels: coughing up flu's proton channels,” Nature, vol. 451, no. 7178, pp. 532–533, 2008. View at: Publisher Site  Google Scholar
 J. R. Schnell and J. J. Chou, “Structure and mechanism of the M2 proton channel of influenza A virus,” Nature, vol. 451, no. 7178, pp. 591–595, 2008. View at: Publisher Site  Google Scholar
 A. L. Stouffer, R. Acharya, D. Salom et al., “Structural basis for the function and inhibition of an influenza virus proton channel,” Nature, vol. 451, pp. 596–599, 2008. View at: Publisher Site  Google Scholar
 H. Leonov and I. T. Arkin, “Structure and dynamics of the influenza A M2 Channel: a comparison of three structures,” Journal of Molecular Modeling, vol. 15, no. 11, pp. 1317–1328, 2009. View at: Publisher Site  Google Scholar
 E. Khurana, M. Dal Peraro, R. DeVane, S. Vemparala, W. F. DeGrado, and M. L. Klein, “Molecular dynamics calculations suggest a conduction mechanism for the M2 proton channel from influenza A virus,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 4, pp. 1069–1074, 2009. View at: Publisher Site  Google Scholar
 M. G. Yi, T. A. Cross, and H.X. Zhou, “Conformational heterogeneity of the M2 proton channel and a structural model for channel activation,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 32, pp. 13311–13316, 2009. View at: Publisher Site  Google Scholar
 Y. Ohigashi, C. L. Ma, X. H. Jing, V. Balannick, L. H. Pinto, and R. A. Lamb, “An amantadinesensitive chimeric BM2 ion channel of influenza B virus has implications for the mechanism of drug inhibition,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 44, pp. 18775–18779, 2009. View at: Publisher Site  Google Scholar
 K. Nishimura, S. G. Kim, L. Zhang, and T. A. Cross, “The closed state of a H^{+} channel helical bundle combining precise orientational and distance restraints from solid state NMR,” Biochemistry, vol. 41, no. 44, pp. 13170–13177, 2002. View at: Publisher Site  Google Scholar
 J. Hu, T. Asbury, S. Achuthan et al., “Backbone structure of the amantadineblocked transmembrane domain M2 proton channel from influenza A virus,” Biophysical Journal, vol. 92, no. 12, pp. 4335–4343, 2007. View at: Publisher Site  Google Scholar
 S. D. Cady, T. V. Mishanina, and M. Hong, “Structure of amantadinebound M2 transmembrane peptide of influenza A in lipid bilayers from magicanglespinning solidstate NMR: the role of Ser31 in amantadine binding,” Journal of Molecular Biology, vol. 385, no. 4, pp. 1127–1141, 2009. View at: Publisher Site  Google Scholar
 X. H. Jing, C. L. Ma, Y. Ohigashi et al., “Functional studies indicate amantadine binds to the pore of the influenza A virus M2 protonselective ion channel,” Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 31, pp. 10967–10972, 2008. View at: Publisher Site  Google Scholar
 T. A. Cross, “Flu BM2 structure and function,” Nature Structural & Molecular Biology, vol. 16, no. 12, pp. 1207–1209, 2009. View at: Publisher Site  Google Scholar
 I. Kass and I. T. Arkin, “How pH opens a H^{+} channel: the gating mechanism of influenza A M2,” Structure, vol. 13, no. 12, pp. 1789–1798, 2005. View at: Publisher Site  Google Scholar
 E. Khurana, M. Dal Peraro, R. DeVane, S. Vemparala, W. F. DeGrado, and M. L. Klein, “Molecular dynamics calculations suggest a conduction mechanism for the M2 proton channel from influenza A virus,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 4, pp. 1069–1074, 2009. View at: Publisher Site  Google Scholar
 M. Sharma, M. G. Yi, H. Dong et al., “Insight into the mechanism of the influenza A proton channel from a structure in a lipid bilayer,” Science, vol. 330, no. 6003, pp. 509–512, 2010. View at: Publisher Site  Google Scholar
 J. Wang, S. Kim, F. Kovacs, and T. A. Cross, “Structure of the transmembrane region of the M2 protein H^{+} channel,” Protein Science, vol. 10, no. 11, pp. 2241–2250, 2001. View at: Publisher Site  Google Scholar
 R. M. Pielak, J. R. Schnell, and J. J. Chou, “Mechanism of drug inhibition and drug resistance of influenza A M2 channel,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 18, pp. 7379–7384, 2009. View at: Publisher Site  Google Scholar
 R. M. Pielak and J. J. Chou, “Solution NMR structure of the V27A drug resistant mutant of influenza A M2 channel,” Biochemical and Biophysical Research Communications, vol. 401, no. 1, pp. 58–63, 2010. View at: Publisher Site  Google Scholar
 R. Acharya, V. Carnevale, G. Fiorin et al., “Structure and mechanism of proton transport through the transmembrane tetrameric M2 protein bundle of the influenza A virus,” Proceedings of the National Academy of Sciences of the United States of America, vol. 107, no. 34, pp. 15075–15080, 2010. View at: Publisher Site  Google Scholar
 J. Hu, R. Fu, K. Nishimura et al., “Histidines, heart of the hydrogen ion channel from influenza A virus: toward an understanding of conductance and proton selectivity,” Proceedings of the National Academy of Sciences of the United States of America, vol. 103, no. 18, pp. 6865–6870, 2006. View at: Publisher Site  Google Scholar
 C. Wang, R. A. Lamb, and L. H. Pinto, “Activation of the M2 ion channel of influenza virus: a role for the transmembrane domain histidine residue,” Biophysical Journal, vol. 69, no. 4, pp. 1363–1371, 1995. View at: Publisher Site  Google Scholar
 A. Okada, T. Miura, and H. Takeuchi, “Protonation of histidine and histidinetryptophan interaction in the activation of the M2 ion channel from influenza A virus,” Biochemistry, vol. 40, no. 20, pp. 6053–6060, 2001. View at: Publisher Site  Google Scholar
 Y. J. Tang, F. Zaitseva, R. A. Lamb, and L. H. Pinto, “The gate of the influenza virus M2 proton channel is formed by a single tryptophan residue,” The Journal of Biological Chemistry, vol. 277, no. 42, pp. 39880–39886, 2002. View at: Publisher Site  Google Scholar
 L. H. Pinto and R. A. Lamb, “Influenza virus proton channels,” Photochemical and Photobiological Sciences, vol. 5, no. 6, pp. 629–632, 2006. View at: Publisher Site  Google Scholar
 J. A. Mould, H.C. Li, C. S. Dudlak et al., “Mechanism for proton conduction of the M2 ion channel of influenza A virus,” The Journal of Biological Chemistry, vol. 275, no. 12, pp. 8592–8599, 2000. View at: Publisher Site  Google Scholar
 M. L. Brewer, U. W. Schmitt, and G. A. Voth, “The formation and dynamics of proton wires in channel environments,” Biophysical Journal, vol. 80, no. 4, pp. 1691–1702, 2001. View at: Publisher Site  Google Scholar
 B. Martinac, Y. Saimi, and C. Kung, “Ion channels in microbes,” Physiological Reviews, vol. 88, no. 4, pp. 1449–1490, 2008. View at: Publisher Site  Google Scholar
 C. Y. Wei and A. Pohorille, “Activation and proton transport mechanism in influenza A M2 channel,” Biophysical Journal, vol. 105, no. 9, pp. 2036–2045, 2013. View at: Publisher Site  Google Scholar
 R. Witter, R. Fu, C. Li, T. A. Cross, U. Sternberg, and A. S. Ulrich, “Solidstate 15NNMR investigations on the pH dependent conformation of Trp41 of the M2 proton channel of influenza A virus,” WSEAS Transactions on Biology and Biomedicine, pp. 277–283, 2011. View at: Google Scholar
 P. E. Czabotar, S. R. Martin, and A. J. Hay, “Studies of structural changes in the M2 proton channel of influenza a virus by tryptophan fluorescence,” Virus Research, vol. 99, no. 1, pp. 57–61, 2004. View at: Publisher Site  Google Scholar
 C. Pérez, F. Löhr, H. Rüterjans, and J. M. Schmidt, “Selfconsistent Karplus parametrization of ^{3}J couplings depending on the polypeptide sidechain torsion Χ_{1},” Journal of the American Chemical Society, vol. 123, no. 29, pp. 7081–7093, 2001. View at: Publisher Site  Google Scholar
 I. Jakovkin, M. Klipfel, C. MuhleGoll, A. S. Ulrich, B. Luy, and U. Sternberg, “Rapid calculation of protein chemical shifts using bond polarization theory and its application to protein structure refinement,” Physical Chemistry Chemical Physics, vol. 14, pp. 12263–12276, 2012. View at: Google Scholar
 M. O'Keeffe and N. E. Brese, “Atom sizes and bond lengths in molecules and crystals,” Journal of the American Chemical Society, vol. 113, no. 9, pp. 3226–3229, 1991. View at: Publisher Site  Google Scholar
 F. H. Allen, O. Kennard, D. G. Watson, L. Brammer, A. G. Orpen, and R. Taylor, “Tables of bond lengths determined by Xray and neutron diffraction. Part 1. Bond lengths in organic compounds,” Journal of the Chemical Society, Perkin Transactions, vol. 2, supplement, no. 12, pp. S1–S19, 1987. View at: Publisher Site  Google Scholar
 M. J. Frisch, G. W. Trucks, H. B. Schlegel et al., Gaussian 98 Revision A.5, 1998.
 H. J. M. Bowen and L. E. Sutton, Tables of Interatomic Distances and Configuration in Molecules and Ions, Chemical Society, Great Britain, UK, 1958.
 T. A. Halgren, “MMFF VI. MMFF94s option for energy minimization studies,” Journal of Computational Chemistry, vol. 20, no. 7, pp. 720–729, 1999. View at: Publisher Site  Google Scholar
 T. A. Halgren, “MMFF VII. Characterization of MMFF94, MMFF94s, and other widely available force fields for conformational energies and for intermolecularinteraction energies and geometries,” Journal of Computational Chemistry, vol. 20, no. 7, pp. 730–748, 1999. View at: Publisher Site  Google Scholar
 T. A. Halgren, “Force Field Evaluation Suite,” T. A. Halgren, Ed., Jersey City, 1999, http://server.ccl.net/cca/data/ff_evaluation_suite/. View at: Google Scholar
 L. J. Holsinger, D. Nichani, L. H. Pinto, and R. A. Lamb, “Influenza A virus M2 ion channel protein: a structurefunction analysis,” Journal of Virology, vol. 68, pp. 1551–1563, 1994. View at: Google Scholar
 I. V. Chizhmakov, F. M. Geraghty, D. C. Ogden, A. Hayhurst, M. Antoniou, and A. J. Hay, “Selective proton permeability and pH regulation of the influenza virus M2 channel expressed in mouse erythroleukaemia cells,” Journal of Physiology, vol. 494, no. 2, pp. 329–336, 1996. View at: Publisher Site  Google Scholar
 C. Li, H. Qin, F. P. Gao, and T. A. Cross, “Solidstate NMR characterization of conformational plasticity within the transmembrane domain of the influenza A M2 proton channel,” Biochimica et Biophysica Acta: Biomembranes, vol. 1768, no. 12, pp. 3162–3170, 2007. View at: Publisher Site  Google Scholar
 W. B. Luo, S. D. Cady, and M. Hong, “Immobilization of the influenza A M2 transmembrane peptide in virus envelopemimetic lipid membranes: a solidstate NMR investigation,” Biochemistry, vol. 48, no. 27, pp. 6361–6368, 2009. View at: Publisher Site  Google Scholar
 J. Hu, T. Asbury, S. Achuthan et al., “Backbone structure of the amantadineblocked transmembrane domain M2 proton channel from influenza A virus,” Biophysical Journal, vol. 92, no. 12, pp. 4335–4343, 2007. View at: Publisher Site  Google Scholar
 F. A. Kovacs and T. A. Cross, “Transmembrane fourhelix bundle of influenza A M2 protein channel: structural implications from helix tilt and orientation,” Biophysical Journal, vol. 73, no. 5, pp. 2511–2517, 1997. View at: Publisher Site  Google Scholar
 F. A. Kovacs, J. K. Denny, Z. Song, J. R. Quine, and T. A. Cross, “Helix tilt of the M2 transmembrane peptide from influenza A virus: an intrinsic property,” Journal of Molecular Biology, vol. 295, no. 1, pp. 117–125, 2000. View at: Publisher Site  Google Scholar
 J. F. Wang, S. Kim, F. Kovacs, and T. A. Cross, “Structure of the transmembrane region of the M2 protein H^{+} channel,” Protein Science, vol. 10, no. 11, pp. 2241–2250, 2001. View at: Publisher Site  Google Scholar
 C. Tian, P. F. Gao, L. H. Pinto, R. A. Lamb, and T. A. Cross, “Initial structural and dynamic characterization of the M2 protein transmembrane and amphipathic helices in lipid bilayers,” Protein Science, vol. 12, no. 11, pp. 2597–2605, 2003. View at: Publisher Site  Google Scholar
 S. D. Cady, C. Goodman, C. D. Tatko, W. F. DeGrado, and M. Hong, “Determining the orientation of uniaxially rotating membrane proteins using unoriented samples: a ^{2}H, ^{13}C, and ^{15}N solidstate NMR investigation of the dynamics and orientation of a transmembrane helical bundle,” Journal of the American Chemical Society, vol. 129, no. 17, pp. 5719–5729, 2007. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Raiker Witter et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.