Bioscience Biotechnology Research Communications

An International  Peer Reviewed Refereed Open Access Journal

P-ISSN: 0974-6455 E-ISSN: 2321-4007

Bioscience Biotechnology Research Communications

An Open Access International Journal

N. A. Moghaddam1, Saharnaz Ahmadi2 and Reza Rasoolzadeh*3

1Department of Biology, Central Tehran Branch, Islamic Azad University, Tehran, Iran

2Department of Biology, Science and Research Branch, Islamic Azad University, Tehran, Iran

3Young Researchers and Elites club, Science and Research Branch, Islamic Azad University, Tehran, Iran

Article Publishing History

Received: 11/08/2016

Accepted After Revision: 20/09/2016

ABSTRACT:

The importance of ionic channels is due to the passage of ions across the cell membrane which is based on electrochemical gradients. The structure of ionic channels often includes one or several central cores which makes up the pore. The direct electron transfer between the enzyme and unmodified electrode is usually prohibited due to shielding of the redox active sites by the protein shells. Monte Carlo simulation have been used to investigate protein folding pathways with some success. Monte Carlo was originally developed for calculating equilibrium properties of physical systems .In calculations we optimized the geometry and defined Potential Energy of the nanotube structure by performing molecular mechanics calculation using MM+ force field, if too large a time step is used in Monte Carlo simulation, it is possible to have a basic instability in the equations that result in a molecule blowing apart, we need small time steps to preserve integration accuracy, however in the Monte Carlo time step 50 femtoseconds (0.05ps) was appropriate. next step we calculated the Vibrational modes of the tube by applying the semi-empirical molecular orbital method. In this paper, we have studies the stability of CNT-Amino acids clusters using by semi-empirical method and investigation of vibrational frequencies and electrical properties. In the more the potential energy increases the more the conductivity of nanochannels decreases and we chose the least energy among nanotube and amino acid complexes. Also the more energy we use, the more conductivity we will have; therefore, we choose the complex which conducts the most current.

KEYWORDS:

Amino Acids- Cnt- Membranes Proteins- Monte Carlo- Semi Empirical

Download this article as:

Copy the following to cite this article:

Moghaddam N. A, Ahmadi S, Rasoolzadeh R. Amino Acid Binding to Nanotube: Simulation of Membrane Protein Channels by Computational Methods. Biosc.Biotech.Res.Comm. 2016;9(3).


Copy the following to cite this URL:

Moghaddam N. A, Ahmadi S, Rasoolzadeh R. Amino Acid Binding to Nanotube: Simulation of Membrane Protein Channels by Computational Methods. Biosc.Biotech.Res.Comm. 2016;9(3). Available from: https://bit.ly/2Nl2xDR


Introduction

Protein helices which make up the pore have consisted of four distinct subunits or one subunit which includes repetitive parts. Any disorder in protein-made channels causes paroxysm attacks. For instance, we can mention neuromuscular diseases as one type of these illnesses. These diseases are called disorders of ionic canals.

The function of channels is to allow selectivity and specificity for a Variety of molecular species transport across the cell membrane (Hornig et al. 2013, Rui et al. 2012). These channels, included: a) ligand-gated channels, b) voltage-gated channels, c) Second messenger gated channels, d) mechanosensitive channels, e) Gap junctions: porins not gated (Mollaamin et al. 2010, Lacroix et al. 2014).

The first step in understanding the physical mechanism of potassium transport through this protein nanopore is the determination of the water molecular distribution along the axial length of the pore Ion channels which are membrane proteins that mediation flux between the outside of the cell through a small, water-filled hole in the membrane-a pore. Ion- Selective pores were originally proposed to explain separate Components of Na+, K+ and leak currents in the classic experiments of Hodgkin and Huxley (Chone 2002 , Lee. AG. 2004). Potassium channels are the most diverse group of the ion channel family (Hornig et al. 2013). The recent determination of the crystallographic structure a bacterial K+ channel from Streptomyces lividans (KcsA) (Syeda et al. 2012) has provided the molecular basic for understanding the physical mechanisms controlling ionic selectivity, permeation, and transport through Various types of K+ channels. (Mantegazza and Catterall 2012, Rui et al. 2012).

In all cases, the functional K+ channel is tetramer (Mollaamin et al. 2010), typically of four identical Subunits folded around a central Port (Syeda et al. 2012). Voltage – gated potassium (Kv) channels are members of the voltage – gated ion channel superfamily (Hornig et al. 2013, Syeda et al. 2012), which is important for initation and propagation of action potentials in excitable cells. They are composed of four identical or homologous Subunits, each containing six transmembrane segments: S1-S6. Segments: S1-S4 form the voltage- sensing domain (VSD), and segments S5 and S6 Connected by the P loop, which is involved in ion selectivity, Comprise the pore- forming domain (PD) S4 has four gating – charge- carrying arginines (R1-R4) spaced at intervals of three amino acid residues, which are highly conserved and are thought to play a key role in coupling changes in membrane Voltage to opening and closing of the pore (Mantegazza and Catterall 2012, Mollaamin et al. 2010). In the Kv channels 13 electronic charges across the membrane electrical field per channel between the closed and open states (Lacroix et al. 2014, Vladimir et al. 2006). Arginine residues interacting with lipid phosphate groups play an important role in stabilizing the voltage- Sensor domain of the KvAP channel within a bilayer. Simulations of the bacterial potassium channel kcsA reveal specific interactions of phosphatidylglycerol with an acidic lipid-binding site the interface between adjacent protein monomers.

Molecular and langevin dynamics simulation as well as Monte Carlo simulation have been used to investigate protein folding pathways with some success. The metropolis Monte Carlo was originally developed for calculating equilibrium properties of physical systems (Sonoda et al. 2011, Sansom et al.2005, Jafari-Dehkordi et al. 2015).

The metropolis algorithm performs a sample of the configuration space of system starting from a random conformation and repeating a large number of steps. Molecular dynamics simulation is one the most promising approaches for solving the protein folding problem .in this method we observe the time behavior of atoms of the system. In MD simulation, new positions of atoms are calculated by numerical integration of newton’s equation of motion (Mollaamin et al. 2011, Cooke and Schiller 2008).

The repentance of the existence of buckminsterfullerene CNT (Kroto et al. 1985), theoretical speculation about carbon clusters (Mollaamin et al. 2011) over 36 years was finally verified. Since then, this beautiful molecule has attracted ever more attention of theoretical and experimental scientists. Some chemists began to focus their research on the chemistry of this molecule, but real fullerene chemistry began only after 1990. described a method for preparing macroscopic quantities of CNT (Oh et al. 2011).

Direct electron transfers between the electrode and the redox enzyme is very important for fundamental studies and construction of biosensors (Liang et al. 2010, Albareda-Sirvent and Hart 2001). However, the direct electron transfer between the enzyme and unmodified electrode is usually prohibited due to shielding of the redox active sites by the protein shells (Wang et al.2008, Mollaamin et al. 2010).

Therefore, several studies have been made to enhance the electron transfer. Mediators are widely used to access the redox center of an enzyme and then to act as the charge carriers. Mediators can minimize the effects of interferences, lower the operating potential of the electrodes, and improve the linear response range and sensitivity of the sensor (Zhang et al. 2005). Use of carbon nanotubes (CNTs) as mediators has attracted increasing attention in recent years. Comparing with traditional carbon electrodes, CNTs show unique properties, such as good conductivity, high chemical stability, and catalytic activities towards many electrochemical reactions (Mollaamin et al. 2010, Wang et al. 2008, Liang et al. 2010, Manso et al. 2007, Jurkschat et al. 2006). More importantly, it is possible to bring the nanotubes close to the redox centers of the proteins (Gooding et al. 2003, Liu et al. 2003).

Material And Methods

Many studies have shown that the carbon nanotubes possess remarkable mechanical and physical properties leading to many potential applications such as fluid transport, fluid storage at nanoscale, and nanodevices for drug delivery Since controlled experiments at the nanometer scale are very difficult, the simulation techniques have been widely and successfully used to investigate the mechanical property, wave propagation and resonant frequency (Natsukia et al. 2008).The vibration of molecules is best described using a quantum mechanical approach. A harmonic oscillator does not exactly describe molecular vibrations. Bond stretching is better described by a Morse potential and conformational changes have sine-wave-type behavior. However, the harmonic oscillator description is very useful as an approximate treatment for low vibrational quantum numbers (Youky et al. 2009).

A harmonic oscillator approximation is most widely used for computing molecular vibrational frequencies because more accurate methods require very large amounts of CPU time. Frequencies computed with the Hartree-Fock approximation and a quantum harmonic oscillator approximation tends to be 10% too high due to the harmonic oscillator approximation and lack of electron correlation (Fernandez et al. 2006). The high-frequency oscillations encountered using flexible water models are of the order of 3500 cm-1 which are somewhat larger than the CNT vibrational modes of 1500 cm-1 (Walther et al. 2001).

Hence, for this case study, the use of the flexible water and solvent model. Vibrational frequencies from semi-empirical calculations tend to be qualitative in (which) they reproduce the general trend mentioned in the Results here. However, the actual values are erratic. Some values will be close, whereas others are often too high. However, PM3 is generally more accurate than AM1.

Since periodic boundary conditions cannot be adopted, first principles calculations of finite-length SWCNTs are only affordable to relatively small systems with C atom number less than 300 within our present computational ability (Ma and Guo (2008).

The molecular mechanics method using the MM+ force field, and the Austin Model 1 (AM1) (Dewar et al. 1985) and Parameterized Model number 3 (PM3) (Stewart 1989) semi-empirical method within the Restricted Hartree–Fock (RHF) formalism are sufficient to study carbon systems (Erkoc 2004). In 1989, Stewart improved the techniques of parameterization and published PM3, which gave lower average errors than AM1, are sufficient to study carbon systems, mainly for the enthalpies of formation (Stewart 1989).

Table 1b: Calculated Properties of C60 and binding to amino acids.

(A) Electric Potential, (B) Electrostatic Properties, (C) Total atomic charges

CNT-His-Pro CNT-Lys-Arg
Atomic number A B C Atomic number A B C
N1 -18.421 10.12 -0.502 N1 -18.340 6.28 -0.593
C2 -14.757 8.93 -0.054 C2 -14.753 7.43 -0.065
C3 -14.713 7.67 0.321 C3 -14.769 7.96 0.178
O4 -22.380 7.92 -0.414 O4 -22.375 7.54 -0.398
C5 -14.800 9.29 -0.253 C5 -14.751 8.66 -0.255
C6 -14.798 10.82 -0.276 C6 -14.732 9.86 -0.282
C7 -14.769 11.22 -0.140 C7 -14.652 11.10 -0.120
N8 -18.322 6.66 -0.608 N8 -18.228 12.30 -0.649
C9 -14.690 7.78 -0.063 C9 -14.496 13.51 0.863
C10 -14.702 9.04 0.182 N10 -18.215 13.80 -0.728
O11 -22.316 9.14 -0.392 N11 -18.218 14.54 -0.741
C12 -14.694 8.25 -0.288 C17 -14.771 6.44 -0.159
C13 14.619 9.55 0.276 C70 -14.761 5.38 0.052
N14 -18.192 10.31 -0.679 N72 -18.513 9.18 -0.514
C15 -14.633 10.36 0.149 C73 -14.797 8.59 -0.123
C16 -14.564 11.42 0.411 C74 -14.728 7.48 0.371
N17 -18.187 11.46 -0.679 O75 -22.350 7.22 -0.447
C19 -14.698 5.48 0.186 C76 -14.795 9.74 -0.236
C73 -14.774 6.36 -0.217 C77 -14.778 9.46 -0.234
C78 -14.718 10.81 -0.271
C79 -14.637 10.84 -0.207
N80 -18.155 12.24 -0.680

Table 1a: Calculated Properties of CNT and binding to amino acids. (A) Electric Potential, (B) Electrostatic Properties, (C) Total atomic charges
CNT (B)

CNT-Ala-Gly

CNT-Val-Ala
Atomic number A B C Atomic number A B C Atomic number A B C
C1 -14.562 5.91621 0.049 N1 -18.293 6.77 -0.446 N1 -18.28 6.52 -0.25
C2 -14.562 5.91621 0.050 C2 -14.653 7.79 -0.107 C2 -14.65 7.76 0.02
C5 -14.559 5.91623 0.053 C3 -14.621 9.10 0.337 C3 -14.64 8.55 0.33
C6 -14.560 5.91621 0.051 O4 -22.238 9.32 -0.285 O4 -22.25 8.26 -0.33
C9 -14.56 5.91624 0.052 C5 -14.714 8.26 -0.413 C5 -14.71 8.78 -0.12
C10 -14.560 5.91625 0.052 C14 -14.744 6.07 -0.207 C6 -14.73 9.20 -0.42
C51 -14.567 5.91618 0.068 C56 -14.657 5.75 0.027 C7 -14.74 10.07 -0.40
C52 -14.568 5.91619 0.061 N66 -18.404 9.75 -0.406 C13 -14.69 5.52 -0.24
C55 -14.592 5.91619 0.040 C67 -14.735 8.57 -0.240 C66 -14.73 6.28 -0.10
C56 -14.588 5.91617 0.045 C68 -14.690 7.36 0.358 N68 -18.39 8.64 -0.63
C59 -14.604 5.91634 0.022 O69 -22.367 7.65 -0.429 C69 -14.72 7.89 -0.11
C60 -14.607 5.91634 0.025 C70 -14.69 7.51 0.31
O71 -22.35 7.89 -0.40
C72 -14.77 8.95 -0.37
CNT- Ala-Leu CNT-Ser-Thr CNT-Gln-Asn
Atomic

Number

A B C Atomic number A B C Atomic number A B C
N1 -18.301 7.33 -0.393 N1 -18.30 7.14 -0.50 N1 -18.264 6.30 -0.469
C2 -14.673 7.85 -0.188 C2 -14.67 8.09 -0.10 C2 -14.664 7.48 -0.047
C3 -14.651 9.29 0.322 C3 -14.64 8.44 0.36 C3 -14.666 7.53 0.339
O4 -22.260 9.85 -0.310 O4 -22.23 7.99 -0.25 O4 -22.257 6.68 -0.386
C5 -14.765 7.20 -0.399 C5 -14.67 9.46 -0.01 C5 -14.728 8.71 -0.236
C6 -14.683 7.85 0.135 O6 -22.31 10.44 -0.59 C6 -14.733 9.94 -0.351
C60 -14.776 5.48 -0.172 C8 -14.71 5.85 0.30 C7 -14.733 11.19 0.513
N66 -18.382 8.94 -0.368 C61 -14.69 5.97 0.03 O8 -22.361 11.23 -0.453
C67 -14.714 7.86 -0.116 N67 -18.42 8.89 -0.46 N9 -18.329 12.34 -0.718
C68 -14.729 6.50 0.3 62 N68 -14.70 8.05 -0.12 C15 -14.776 6.32 0.020
O69 -22.364 6.59 -0.452 C69 -14.68 6.94 0.04 C69 -14.649 5.53 0.101
C70 -14.743 8.23 -0.231 O70 -22.30 6.94 -0.34 N70 -18.412 9.61 -0.423
C71 -14.735 7.38 -0.133 C71 -14.64 7.63 -0.03 C71 -14.726 8.39 -0.136
C72 -14.771 7.57 -0.389 O72 -22.17 6.91 -0.47 C72 -14.714 7.20 0.175
C73 -14.756 8.17 -0.395 C73 -14.74 8.989 -0.38 O73 -22.346 7.53 -0.320
C74 -14.728 8.61 -0.344
C75 -14.619 7.57 0.561
O76 -22.335 6.55 -0.470
N77 -18.309 8.02 -0.711

Purpose of C60 is nanotube that includes 60 carbons. All calculations presented here were performed with semi-empirical Molecular mechanics (MM+) (Table 2).

Table 2: Semi empirical Calculations for CNT and conjunction to amino acids.
  C60 Nanotube C60-Ala-Gly C60-Ala-Leu C60–Val-Ala C60-Gln-Asn C60–His-Pro C60-Lys-Arg C60-Phe-Tyr C60-Ser-Thr C60-Trp-Pro
Total Energy (kcal/mol) -207257 -254702 -273475 -188038 -298292 -291066 739779.4 789947 573814.9 831669.4
Binding Energy (kcal/mol) -37784.3 -45141.8 -50207.3 32104.26 -52071.7 -52672.9 993483.9 1046278 804308.2 1083965
Electronic Energy (kcal/mol) -2069240 -2722240 -3106366 -2914485 -3373457 -3256055 -2520644 -2590608 -2239658 -2639509
Rotation Frequency 352.31 361.96 358.71 3773.92 352.64 369.47 251.51 1816.71 -3825.72 1864.92
Rotation Intensity 0.064 1.285 0.347 8.368 0.315 1.969 1.6 2.669 2.892 5.153
∆E -2283.99 -2739.14 -2895.54 -2857.42 -3154.71 -3082.21 -3235.42 -3315.58 -2967.73 -3226.15
Dipole moment 4.1475 11.759 3.2852 8.2113 14.002 32.3063 33.2251 6.3174 6.7668 8.1059

In the first step of the calculations we optimized the geometry and defined Potential Energy of the nanotube structure by performing molecular mechanics calculation using MM+ force field, if too large a time step is used in Monte Carlo simulation, it is possible to have a basic instability in the equations that result in a molecule blowing apart, we need small time steps to preserve integration accuracy, however in the Monte Carlo time step 50 femtoseconds (0.05ps) was appropriate (Safi Najafabadi et al. 2015). In the next step we calculated the Vibrational modes of the tube by applying the semi-empirical molecular orbital method by the Hyperchem-7 package program (Hyperchem 7.0 2007).

Results And Discussion

The resulting method was denoted, and in a sense, it is the best set of parameters (or at least a good local minimum) for the given set of experimental data. The optimization process, however, still requires some human intervention in selecting the experimental data and assigning appropriate weight factors to each set of data. As a reference Table1 is the result of semi-empirical computation using both method AM1& PM3.At the first glance in Table1, it can be observed by increasing dielectrics, normal modes will move to upper normal modes ratio and Vibrational frequencies resulted from semi-empirical Calculations tend to be qualitative. Therefore, by increasing dielectric the higher frequency will be gained in which semi-empirical methods will have the same operating procedure (Mollaamin and Monajjemi 2015).

The PM3 and AM1 methods are also more popular than other semi-empirical methods due to the availability of algorithms for including salvation effects in these calculations. There are also some known strengths and limitations of PM3. Overall heats of formation are more accurate than with AM1. Hypervalent molecules are also predicted more accurately. On average, PM3 predicts energies more accurately than AM1.The heats of formation are more accurate than AM1 or PM3 depending on the nature of the system and information desired, they will often give the most accurate obtainable results for organic molecules with semi-empirical methods. On average, PM3 predicts energies and geometries better than AM1 (Christensen et al. 2016).

There is energy of interaction in between solvent and solute. Therefore, the solute properties dependent on energy, such as geometry, total energy and vibrational frequencies depend on the solvent. The presence of a solvent, particularly a polar solvent, can also stabilize charge separation within the molecule. This not only changes the energy, but also results in a shift in the electron density and associated properties. In reality, these are the result of the quantum mechanical interaction between solvent and solute, which must be averaged over all possible attitudes of solvent molecules consistent with the principles of statistical mechanics. The results of the CNT-amino acids cluster simulation can be used to analyze the energetic aspects which are associated with the process of introducing a CNT fullerene from the gas phase into different amino acids (Table 2 and Fig 1).

Figure 1: Plotting of Calculated Properties of CNT and binding to amino acids. Figure 1: Plotting of Calculated Properties of CNT and binding to amino acids.

The net result clearly indicates that the process of introducing a CNT to different amino acids energetically is remarkable. Whereas a molecular mechanics potential used to characterize the response of a SWCNT which allows long range interactions between atoms, results in Molecular Mechanics force field which indicates that potential energy is maximum and also the potential energy will increase.

Conclusion

To reconstruct membranous proteins, we used simulated nanotubes in order to transfer ions across the membrane and transfer of ions was done successfully. In this study the more the potential energy increases the more the conductivity of nanochannels decreases and we chose the least energy among nanotube and amino acid complexes. And also the more energy we use, the more conductivity we will have; therefore, we choose the complex which conducts the most current. was observed most changes of Rotation Intensity and Rotation Frequency in the CNT-ALA-VAL and CNT-SER-THR complexes. Also we found maximum of Dipole moment (Debye) in the CNT-LYS-ARG and CNT-HIS-PRO complexes. Minimum of ∆E for the CNT-PHE-TYR obtained. This way we can simulate the channels which have hereditary defects and are not efficiently and can observe the fundamental cure of the diseases.

Acknowledgement

Support by the Young Researchers and Elites club, Science and Research Branch, Islamic Azad University, Tehran, Iran. is greatly acknowledged and we appreciate the Dr Fahimeh Sadat Vajedi in Kashan University.

References

Albareda-Sirvent, M. and A. L. Hart (2002). Preliminary estimates of lactic and malic acid in wine using electrodes printed from inks containing sol–gel precursors. Sensors and Actuators B: Chemical Vol. 87 No.1: Pages 73-81.

Chone, S. (2002). Potassim channel structures. Nature Reviews Neuroscience Vol. 3:115-121.

Christensen, A. S., T. Kubar, Q. Cui and M. Elstner (2016). Semiempirical Quantum Mechanical Methods for Noncovalent Interactions for Chemical and Biochemical Applications. Chemical Review Vol. 116 No. 9: Pages 5301–5337.

Cooke, B. and S. C. Schiller (2008). Statistical Prediction and Molecular Dynamics Simulation. Biophysical Journal Vol. 95 No. 10: Pages 4497–4511.

Dewar, M. J. S., E. G. Zoebisch, E. F. Healy, J. J. P. Stewart (1985). Development and use of quantum mechanical molecular models. 76. AM1: a new general purpose quantum mechanical molecular model. Journal of The American. Chemical Society Vol. 107 No.13: Pages 3902-3909.

Fernandez, I. R., H. Fangohr and A. Bhaskar (2006). Normal modes of carbon nanotubes: similarities and differences with their continuum counterpart. Journal of Physics: Conference Series Vol.26 No.1: Pages 131–134.

Gooding, J. J., R. Wibowo, J. Q. Liu, W. R. Yang, D. Losic, S. Orbons, F. J. Mearns, J. G. Shapter, D. B. Hibbert (2003). Protein electrochemistry using aligned carbon nanotube arrays. Journal of the American Chemical Society Vol. 125 No. 30: Pages 9006-9007.

Hornig, S., I. Ohmert, D. Trauner, Ch. Ader, M. Baldus and O. Pongs (2013). Tetraphenylporphyrin derivative specifically blocks members of the voltage-gated potassium channel subfamily Kv1. Channels Vol. 7 No 6: Pages 473–482.

Hyperchem 7.0 (2007). Quantum Biochemistry, Hype cube, Inc, Gainesville, FL. U.S.A. Vol. 103: Pages 1233-1237.

Hyperchem 8.0 (2007). Quantum Biochemistry, Hype cube, Inc, Gainesville, FL. U.S.A Vol. 103: Pages 1233-1237.

Jafari-Dehkordi, Sahar., Z. Aghili, S. Ahmadi, S. Jabbari, I. Rezazadeh, R. Hasani and R. Rasoolzadeh (2015). Molecular Mechanics Investigation of Diffrent Temperature Effects on Bacillus licheniformis á-amylase: A Computational Study. Journal of pure and applied microbiology Vol. 9 No. 1: Pages 607-611.

Jurkschat, K., S. J. Wilkins, C. J. Salter, H. C. Leventis, G. G. Wildgooses, L. Jiang, T. G. Jones, A. Crossley and R. G. Compton (2006). Multiwalled carbon nanotubes with molybdenum dioxide nanoplugs-New chemical nanoarchitectures by electrochemical modification. Small Vol. 2 No.1: Pages 95-98.

Kolinski, A. and J. Skolnick (2004). Reduced models of proteins and their applications. Polymers Vol. 45: Pages 511-524.

Kroto, H. W., J. R. Heath, S. C. O’Brien, R. F. Curl and R. E. Smalley (1985). C60-buckminsterfullerene. Nature Vol. 318 No. 36: Pages 162-163.

Lacroix, J. J., H. Clark HydeF.V. Campos and F. Bezanilla (2014). Moving gating charges through the gating pore in a Kv channel voltage sensor. PNAS Vol. 111, No. 19: Pages 1950-1959.

Lee. AG. (2004). How lipids affect the activities of integral membrane proteins. Biochimica et Biophysica Acta Vol.1666 No.1-2: Pages 62-87.

Liang. R. , M. Deng , S. Cui , H. Chen and J. Qiu (2010). Direct electrochemistry and electrocatalysis of myoglobin immobilized on zirconia/multi-walled carbon nanotube nanocomposite. Materials Research Bulletin Vol. 45: Pages 1855-1860.

Liu, J., A. Chou, W. Rahmat, M. N. Paddon-Row and J. J. Gooding (2003). Achieving Direct Electrical Connection to Glucose Oxidase using Aligned Single Walled Carbon Nanotube Arrays. Electroanalysis Vol. 17 No. 1: Pages 38-46.

Ma, Sh. and W. Guo (2008). Size-dependent polarizabilities of finite-length single-walled carbon nanotubes, Physics Letters A Vol. 372 No.27-28: Pages 4835–4838.

Manso, J., M. L. Mean, P. Ya´n˜ez-Seden˜o and J. Pingarro´n (2007). Electrochemical biosensors based on colloidal gold–carbon nanotubes composite electrodes, Journal of Electroanalytical Chemistry Vol. 603 No.1: Pages 1-7.

Mantegazza, M. and W. A. Catterall (2012). Voltage-Gated Na+ Channels: Structure, Function, and Pathophysiology. Jasper’s Basic Mechanisms of the Epilepsies.

Mollaamin, F., I. Layali, A. R. Ilkhani and M. Monajjemi (2010). Nanomolecular simulation of the voltage–gated potassium channel protein by gyration radius study. African Journal of Microbiology Research Vol. 4 No. 24: Pages 2795-2803.

Mollaamin, F., T. Nejadsattari and I. Layali (2011). Energy study at different solvents for potassium Channel Protein by Monte Carlo,Molecular and Langevin Dynamics Simulations . Journal of Physical and Theoretical Chemistry Vol. 8 No. 1: Pages 23-32.

Mollaamin, F., Sh. Momeni, M. Movahedi and M. Monajjemi (2011). Transportation of single wall carbon nanotube (SWCNT)through the cell membrane. African Journal of Microbiology Research. Vol. 5 No. 24: Pages 4175-4181.

Mollaamin, F., K. Shahani poor, T. Nejadsattari and M. Monajjemi (2010). Bio-modeling of active site in oxidized azurin using by computational methods. African Journal Microbiology Research Vol. 4 No.20Pages 2098-2108.

Mollaamin, F. and M. Monajjemi (2015). Harmonic Linear Combination and Normal Mode Analysis of Semiconductor Nanotubes Vibrations. Journal of Computational and Theoretical Nanoscience Vol. 12 No.6: Pages 1030–1039.

Natsukia, T., Q. Q. Ni and E. Morinobu (2008). Analysis of the vibration characteristics of double-walled carbon nanotubes. Carbon. Vol.46 No. 12: Pages 1570-1573.

Oh, W. C., Z. D. Meng and W. B. Ko (2011). Preparation of Fullerene Function Materials under Ultrasonic Irradiation: Review. Journal of Photocatalysis Science Vol. 2 No. 1: Pages 39-46.

Rui, H., M. Rivera and W. Im (2012). Protein Dynamics and Ion Traffic in Bacterioferritin. Biochemistry Vol. 51 No. 49: Pages 9900–9910.

Safi Najafabadi, A., S. Ahmadi, M. Mohammadian Fardin, R. Rasoolzadeh and F.S. Vajedi (2015). Molecular Simulation of GABA(A) Receptor to Study of Effects on Nervous Stimulants Inhibitory & Blood Pressure; A Nano Molecular Modeling of GABARAP. Biosciences Biotechnology Research Asia. Vol. 12 No.1: Pages 419-424.

Sakir Erkoc, S. (2004). Does tubular structure of carbon form only from graphine sheet? Physica E: Low-dimensional Systems and Nanostructures Vol. 25 No.1: Pages 69–77.

Sansom, M. S., P. J. Bond, S. S. Deol, A. Grottesi, S. Haider and Z. A Sands (2005). Molecular simulations and lipid-protein interactions: potassium channels and other membrane proteins. Biochemical Society Transactions Vol. 33 No. 5: Pages 916-920.

Sonoda, Y., S. Newstead, N.J. Hu, Y. Alguel, E. Nji, K. Beis, S. Yashiro, C. Lee, J. Leung, A. D. Cameron, B. Byrne, S. Iwata and D. Drew (2011). Benchmarking membrane protein detergent stability for improving throughput of high-resolution. X-ray structures. Structure Vol. 19 No. 1: Pages 17-25.

Stewart, J. J. P. (1989). Optimization of parameters for semiempirical methods I. Method. Journal of Computational Chemistry Vol. 10 No.2: Pages 209-220.

Syeda, R., J. S. Santosb, M. Montal and H. Bayley (2012). Tetrameric assembly of KvLm K+ channels with defined numbers of voltage sensors. PNAS Vol. 109 No. 42: Pages 16917–16922.

Vladimir Y. Y., D. Baker and W. A. Catterall (2006). Voltage sensor conformations in the open and closed states in ROSETTA structural models of K+ channels Vol. 103 No. 19: Pages 7292-7297. Walther, J. H., R. Jaffe, T. Halicioglu and P. Koumoutsakos (2001). Carbon Nanotubes in Water: Structural Characteristics and Energetics. The Journal of. Physical Chemistry B Vol. 105 No.41: Pages 9980-9987.

Wang, Z., M. Li, P. Su, Y. Zhang, Y. Shen, D. Han, A. Ivaska and L. Niu (2008). Direct electron transfers of horseradish peroxidase and its electrocatalysis based on carbon nanotube/thionine/gold. Electrochemistry Communications Vol.10 No.2: Pages 306-310.

Youky, Ono., N. Kusuno, K. Kusakabe and N. Suzuki (2009). Conduction Properties of a Distorted Buckled Carbon Nanotube in the Vibrational Normal Mode. Japanese Journal of Applied Physics Vol. 48 No.4: Pages 045001-045007.

Zhang, L., R. Yuan, Y. Chai and X. Li (2007). Investigation of the electrochemical and electrocatalytic behavior of positively charged gold nanoparticle and l-cysteine film on an Auelectrode. Analytica Chimica Acta Vol. 596 No.1: Pages 99-105.