Interactive 3D Complement Research Articles

Interaction of Piscidin-1 with zwitterionic versus anionic membranes: a comparative molecular dynamics study

Plasma membrane of each micro-organism has a unique set of lipid composition as a consequence of the environmental adaptation or a response to exposure to antimicrobial peptides (AMPs) as antibiotic agents. Understanding the relationship between lipid composition and action of antimicrobial peptides or considering how different lipid bilayers respond to AMPs may help us design more effective peptide drugs in the future. In this contribution, we intend to elucidate how two currently used membrane models, namely palmitoyl-oleoyl-phosphtidylglycerol (POPG) and 1-palmitoyl-oleoyl-glycero-phosphocholine (POPC), respond to antimicrobial peptide Piscidin-1 (Pis-1).The computed density profile of the peptide as it moves from the bulk solvent toward the membrane core suggests that Pis-1 penetrates into the POPG bilayer less than the POPC membrane. Furthermore, we showed that the two model membranes used in this study have different behavior in the presence of Pis-1. Hence, we suggest that membrane composition could be an important factor in determining lytic ability of peptide drugs to kill a unique bacterial species.An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:37

Han ethnicity-specific type 2 diabetic treatment from traditional Chinese medicine?

Insulin-degrading enzyme (IDE) gene is one of the type 2 diabetes mellitus susceptibility genes specific to the Han Chinese population. IDE, a zinc-metalloendopeptidase, is a potential target for controlling insulin degradation. Potential lead compounds for IDE inhibition were identified from traditional Chinese medicine (TCM) through virtual screening and evaluation of their pharmacokinetic properties of absorption, distribution, metabolism, excretion, and toxicity. Molecular dynamics (MD) simulation was performed to validate the stability of complexes from docking simulation. The top three TCM compounds, dihydrocaffeic acid, isopraeroside IV, and scopolin, formed stable H–bond interactions with key residue Asn139, and were linked to active pocket residues His108, His112, and Glu189 through zinc. Torsion angle trajectories also indicated some stable interactions for each ligand with IDE. Molecular level analysis revealed that the TCM candidates might affect IDE through competitive binding to the active site and steric hindrance. Structural feature analysis reveals that high amounts of hydroxyl groups and carboxylic moieties contribute to anchor the ligand within the complex. Hence, we suggest the top three TCM compounds as potential inhibitor leads against IDE protein to control insulin degradation for type 2 diabetes mellitus.An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:29

Investigation on the site-selective binding of bovine serum albumin by erlotinib hydrochloride

The purpose of this study was to investigate the site-selective binding of erlotinib hydrochloride (ET), a targeted anticancer drug, to bovine serum albumin (BSA) through 1H NMR, spectroscopic, thermodynamic, and molecular modeling methods. The fluorescence quenching of BSA by ET was a result of the formation of BSA–ET complex with high binding affinity. The site marker competition study combined with isothermal titration calorimetry experiment revealed that ET binds to site II of BSA mainly through hydrogen bond and van der Waals force. Molecular docking was further applied to define the specific binding site of ET to BSA. The conformation of BSA was changed in the presence of ET, revealed by synchronous fluorescence, circular dichroism, and three-dimensional fluorescence spectroscopy results. Further, NMR analysis of the complex revealed that the binding capacity contributed by the aromatic protons in the binding site of BSA might be greater than the aliphatic protons. An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:26

Conformational dynamics of full-length inducible human Hsp70 derived from microsecond molecular dynamics simulations in explicit solvent

Human 70 kDa heat shock protein (hHsp70) is an ATP-dependent chaperone and is currently an important target for developing new drugs in cancer therapy. Knowledge of the conformations of hHsp70 is central to understand the interactions between its nucleotide-binding domain (NBD) and substrate-binding domain (SBD) and is a prerequisite to design inhibitors. The conformations of ADP-bound (or nucleotide-free) hHsp70 and ATP-bound hHsp70 was investigated by using unbiased all-atom molecular dynamics (MD) simulations of homology models of hHsp70 in explicit solvent on a timescale of .5 and 2.7 μs, respectively. The conformational heterogeneity of hHsp70 was analyzed by computing effective free-energy landscapes (FELs) and distance distribution between selected pair of residues. These theoretical data were compared with those extracted from single-molecule Förster resonance energy transfer (FRET) experiments and to small-angle X-ray scattering (SAXS) data of Hsp70 homologs. The distance between a pair of residues in FRET is systematically larger than the distance computed in MD which is interpreted as an effect of the size and of the dynamics of the fluorescent probes. The origin of the conformational heterogeneity of hHsp70 in the ATP-bound state is due to different binding modes of the helix B of the SBD onto the NBD. In the ADP-bound (or nucleotide-free) state, it arises from the different closed conformations of the SBD and from the different positions of the SBD relative to the NBD. In each nucleotide-binding state, Hsp70 is better represented by an ensemble of conformations on a μs timescale corresponding to different local minima of the FEL. An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:30

Investigation of silent information regulator 1 (Sirt1) agonists from Traditional Chinese Medicine

Silent information regulator 1 (Sirt1), a class III nicotinamide adenine dinucleotide dependent histone deacetylases, is important in cardioprotection, neuroprotection, metabolic disease, calorie restriction, and diseases associated with aging. Traditional Chinese Medicine (TCM) compounds from TCM Database@Taiwan (http://tcm.cmu.edu.tw/) were employed for screening potent Sirt1 agonists, and molecular dynamics (MD) simulation was implemented to simulate ligand optimum docking poses and protein structure under dynamic conditions. TCM compounds such as (S)-tryptophan-betaxanthin, 5-O-feruloylquinic acid, and RosA exhibited good binding affinity across different computational methods, and their drug-like potential were validated by MD simulation. Docking poses indicate that the carboxylic group of the three candidates generated H-bonds with residues in the protein chain from Ser441 to Lys444 and formed H-bond, π–cation interactions, or hydrophobic contacts with Phe297 and key active residue, His363. During MD, stable π–cation interactions with residues Phe273 or Arg274 were formed by (S)-tryptophan-betaxanthin and RosA. All candidates were anchored to His363 by stable π- or H-bonds. Hence, we propose (S)-tryptophan-betaxanthin, 5-O-feruloylquinic acid, and RosA as potential lead compounds that can be further tested in drug development process for diseases associated with agingAn animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:28

Remarkable disparity in mechanical response among the extracellular domains of type I and II cadherins

Cadherins, a large family of calcium-dependent adhesion molecules, are critical for intercellular adhesion. While crystallographic structures for several cadherins show clear structural similarities, their relevant adhesive strengths vary and their mechanisms of adhesion between types I and II cadherin subfamilies are still unclear. Here, stretching of cadherins was explored experimentally by atomic force microscopy and computationally by steered molecular dynamics (SMD) simulations, where partial unfolding of the E-cadherin ectodomains was observed. The SMD simulations on strand-swapping cadherin dimers displayed similarity in binding strength, suggesting contributions of other mechanisms to explain the strength differences of cell adhesion in vivo. Systematic simulations on the unfolding of the extracellular domains of type I and II cadherins revealed diverse pathways. However, at the earliest stage, a remarkable similarity in unfolding was observed for the various type I cadherins that was distinct from that for type II cadherins. This likely correlates positively with their distinct adhesive properties, suggesting that the initial forced deformation in type I cadherins may be involved in cadherin-mediated adhesion. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:25

Insight into TPMT∗23 mutation mis-folding using molecular dynamics simulation and protein structure analysis

Thiopurine S-methyltransferase (TPMT) is an important enzyme that metabolizes thiopurine drugs. This enzyme exhibits a large number of interindividual polymorphism. TPMT∗23 polymorphism has been reported in a few cases in the world in co-dominance with TPMT∗3A. The phenotype has been reported to affect enzyme activity in vivo and in vitro. Its underlying structural basis is not clarified yet. In our study, the wild type (WT) protein structure was analyzed and the amino acids bordering water channels in thiopurine sites were identified. Molecular dynamics of both the WT and TPMT∗23 mutation was carried out. In addition, the effects of this mutation, especially on the thiopurine site which is closed with a pincer like mechanism, were investigated. We focused on explaining how a locally occurred A167G substitution propagated through hydrogen bonds alteration to induce structural modification which affects both thiopurine and S-adenosylmethionine receptors. Finally, a genetic prediction of mutation functional consequences has been conducted confirming altered activity. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:20

Molecular dynamics simulation to investigate the impact of disulfide bond formation on conformational stability of chicken cystatin I66Q mutant

Chicken cystatin (cC) mutant I66Q is located in the hydrophobic core of the protein and increases the propensity for amyloid formation. Here, we demonstrate that under physiological conditions, the replacement of Ile with the Gln in the I66Q mutant increases the susceptibility for the disulfide bond Cys71–Cys81 to be reduced when compared to the wild type (WT) cC. Molecular dynamics (MD) simulations under conditions favoring cC amyloid fibril formation are in agreement with the experimental results. MD simulations were also performed to investigate the impact of disrupting the Cys71–Cys81 disulfide bond on the conformational stability of cC at the atomic level, and highlighted major disruption to the cC appendant structure. Domain swapping and extensive unfolding has been proposed as one of the possible mechanisms initiating amyloid fibril formation by cystatin. Our in silico studies suggest that disulfide bond formation between residues Cys95 and Cys115 is necessary to maintain conformational stability of the I66Q mutant following breakage of the Cys71–Cys81 disulfide bridge. Subsequent breakage of disulfide bond Cys95–Cys115 resulted in large structural destabilization of the I66Q mutant, which increased the α–β interface distance and expanded the hydrophobic core. These experimental and computational studies provide molecular-level insight into the relationship between disulfide bond formation and progressive unfolding of amyloidogenic cC mutant I66Q. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:23

Identification of structural motifs in the E2 glycoprotein of Chikungunya involved in virus–host interaction

Chikungunya fever is one of the reemerging vector-borne diseases. It has become a major global health problem especially in the developing countries. There are no vaccines or specific antiviral drugs available to date. This study reports small molecule inhibitors of envelope glycoprotein 2 (E2 glycoprotein) which are predicted based on Chikungunya virus–host interactions. E2 glycoprotein of Chikungunya virus interacts at 216 residue of the host receptor protein which plays a vital role in initiating infection. Understanding the structural aspects of E2 glycoprotein is crucial to develop specific inhibitors to prevent the virus binding from host receptors. In silico method was adopted to predict the sequence motifs of envelope protein, as the method like yeast two hybrid system is laborious, time consuming, and costly. The E2 glycoprotein structure of the Indian isolate was modeled using two templates (2XFC and 3JOC) and then validated. The class III PDZ domain binding motif was found to be identified at 213–216 amino acids. The corresponding peptide structures which recognize the PDZ domain binding motif were identified by the literature search and were used for generating five point pharmacophore model (ADDDR) containing acceptor, donor and aromatic ring features. Databases such as Asinex, TosLab and Maybridge were searched for the matches for the predicted pharmacophore model. Two compounds were identified as lead molecules as their glide score is > 5 kcal/mol. Since the pharmacophore model is developed based on Chikungunya virus–host interaction, it can be used for designing promising antiviral lead compounds for the treatment of Chikungunya fever.An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:21

Insights into the drug resistance induced by the BaDHPS mutations: molecular dynamic simulations and MM/GBSA studies

Dihydropteroate synthase (DHPS) is essential for the folic acid biosynthetic pathway in prokaryotes; the mutation forms for DHPS are found to be relative to the urgent drug resistance problems. In our study, the Bacillus anthracis DHPS (BaDHPS) was selected for molecular dynamics and binding free energy studies to investigate the biochemistry behaviors of the wild-type and mutation form BaDHPS proteins (D184N and K220Q). It is found that the conformational change of the ligand dihydropteroate sulfathiazole binding site in mutation D184N and K220Q systems is mainly attributed from the Loop 1, Loop 2, and Loop 7 regions, and the binding free energy of these mutation systems is lower than that of the wild-type system. Additionally, some important hydrogen bonds of the mutation systems are disrupted during the simulations. But the shortening of the distance between residue Thr67 and the ligand would cause significant change of the binding pose in the K220Q system. These studies of DHPS family will be helpful for further drug resistance investigations. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:24

Molecular dynamics simulations of the thermal stability of tteRBP and ecRBP

Molecular dynamics simulations were performed for investigating the thermal stability of the extremely thermophilic Thermoanaerobacter tengcongensis ribose binding protein (tteRBP) and the mesophilic homologous Escherichia coli ribose binding protein (ecRBP). The simulations for the two proteins were carried out under the room temperature (300 K) and the optimal activity temperature (tteRBP 375 K and ecRBP 329 K), respectively. The comparative analyses of the trajectories show that the two proteins have stable overall structures at the two temperatures; further analyses indicate that they both have strong side-chain interactions and different backbone flexibilities at the different temperatures. The tteRBP 375 K and ecRBP 329 K have stronger internal motion and higher flexibility than tteRBP 300 K and ecRBP 300 K, respectively, it is noted that the flexibility of tteRBP is much higher than that of ecRBP at the two temperatures. Therefore, tteRBP 375 K can adapt to high temperature due to its higher flexibility of backbone. Combining with the researches by Cuneo et al., it is concluded that the side-chain interactions and flexibility of backbone are both the key factors to maintain thermal stability of the two proteins. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:22

Crystal structure of the CN-hydrolase SA0302 from the pathogenic bacterium Staphylococcus aureus belonging to the Nit and NitFhit Branch of the nitrilase superfamily

The nitrilases include a variety of enzymes with functional specificities of nitrilase, amidase, and hydrolase reactions. The crystal structure of the uncharacterized protein SA0302 from the pathogenic microorganism Staphylococcus aureus is solved at 1.7 Å resolution. The protein contains 261 amino acids and presents a four-layer αββα sandwich with a chain topology similar to that of a few known CN-hydrolase folds. In the crystal, the proteins are arranged as dimers whose monomers are related by a pseudo twofold rotation symmetry axis. Analysis of the sequences and structures of CN-hydrolases with known 3D structures shows that SA0302 definitely is a member of Branch 10 (Nit and NitFhit) of the nitrilase superfamily. Enzyme activities and substrate specificities of members of this branch are not yet characterized, in contrast to those of the members of Branches 1–9. Although the sequence identities between Branch 10 members are rather low, less than 30%, five conserved regions are common in this subfamily. Three of them contain functionally important catalytic residues, and the two other newly characterized ones are associated with crucial intramolecular and intermolecular interactions. Sequence homology of the area near the active site shows clearly that the catalytic triad of SA0302 is Glu41-Lys110-Cys146. We suggest also that the active site includes a fourth residue, the closely located Glu119. Despite an extensive similarity with other Nit-family structural folds, SA0302 displays an important difference. Protein loop 111–122, which follows the catalytic Lys110, is reduced to half the number of amino acids found in other Nit-family members. This leaves the active site fully accessible to solvent and substrates. We have identified conservative sequence motifs around the three core catalytic residues, which are inherent solely to Branch 10 of the nitrilase superfamily. On the basis of these new sequence fingerprints, 10 previously uncharacterized proteins also could be assigned to this hydrolase subfamily. An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:19

Docking and molecular dynamics studies of peptide inhibitors of ornithine decarboxylase: a rate-limiting enzyme for the metabolism of Fusarium solani

Fusarium solani causes a wide variety of diseases in plants. Polyamine biosynthesis is responsible for the growth and pathogenicity of the fungus. The initial step of this pathway involves the decarboxylation of ornithine to putrescine, and is catalyzed by the enzyme ornithine decarboxylase (ODC). Inhibiting this process may be a promising approach for the management of fungal disease in various crops. Therefore, there is a need to develop inhibitors of ODC that have higher binding capacity than ornithine. Fifteen peptides were designed and modeled based on physicochemical properties of residues in the active site of ODC. The peptide GLIWGNGPF showed the highest dock score. It is assumed that the de novo design of peptides could be a potential approach to inhibit polyamine biosynthesis. Molecular dynamics studies make an important contribution to understanding the effect of the binding of peptides and the stability of an ODC-peptide complex system. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:8.

Protein flexibility and conformational states of Leishmania antigen eIF-4A: identification of a novel plausible protein adjuvant using comparative genomics and molecular modeling

Recent homology modeling studies have identified specific residues (epitope) of the Leishmania RNA helicase protein (LmeIF) that stimulates production of IL-12 cytokine. However, question remains concerning how LmeIF’s N-terminal moiety initiates adjuvant effects. Extensive molecular modeling combining the normal mode analysis (NMA) and molecular dynamics simulations, in the present study, has demonstrated that the LmeIF structure may exist in two different forms corresponding to the extended and collapsed (closed) states of the entire structure. The computational results showed that the two domains of the LmeIF structure tend to undergo large fluctuations in a concerted fashion and have strong effect on the solvent accessible surface of the epitope situated on the N-terminal structure. The conformational freedom of the C-terminal domains may explain why the entire LmeIF protein is not as active as the N-terminal moiety. Thereafter, a comparative genome analysis with subsequent homology modeling and molecular electrostatic potential (MEP) techniques allowed us to predict a novel and plausible RNA helicase (LI-helicase) from the Listeria source with adjuvant property as observed for the Leishmania eIF-4A protein. The structural folding and MEP maps revealed similar topologies of the epitope of both LmeIF and LI-helicase proteins and striking identity in the local disposition of the charged groups. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:7

Evidence-based docking of the urease activation complex

1Ureases require accessory proteins for their activation and proper function. In Klebsiella aerogenes, UreD, UreF, UreG, and UreE are sequentially complexed to UreABC as required for its activation. Until now, only low-resolution structures are available for this activation complex. To circumvent such limitation, our work intends to provide an atomic-level model for the (UreABC–UreDFG)3 complex from K. aerogenes, by employing comparative modeling associated to sequential macromolecular dockings, validated through small-angle X-ray scattering profiles and comparison with results from cross-linking, mutagenesis, and pull-down experiments. Additionally, normal mode analyses of the obtained complex supported the characterization of the elevated flexibility of both UreD–UreF dimer and (UreABC–UreDFG)3 oligomer, explaining the previously observed diffuse binding of UreD to the apoenzyme. The model shown here is the first atomic-level depiction of this complex, a required step for the unraveling of the urease activation process. 1Both authors share senior authorship. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:6

Influence of divalent magnesium ion on DNA: molecular dynamics simulation studies

A large amount of experimental evidence is available on the effect of magnesium ions on the structure and stability of DNA double helix. Less is known, however, on how these ions affect the stability and dynamics of the molecule. The static time average pictures from X-ray structures or the quantum chemical energy minimized structures lack understanding of the dynamic DNA–ion interaction. The present work addresses these questions by molecular dynamics simulation studies on two DNA duplexes and their interaction with magnesium ions. Results show typical B-DNA character with occasional excursions to deviated states. We detected expected stability of the duplexes in terms of backbone conformations and base pair parameter by the CHARMM-27 force field. Ion environment analysis shows that Mg2+ retains the coordination sphere throughout the simulation with a preference for major groove over minor. An extensive analysis of the influence of the Mg2+ ion shows no evidence of the popular predictions of groove width narrowing by dipositive metal ion. The major groove atoms show higher occupancy and residence time compared to minor groove for magnesium, where no such distinction is found for the charge neutralizing Na+ ions. The determining factor of Mg2+ ion’s choice in DNA binding site evolves as the steric hindrance faced by the bulky hexahydrated cation where wider major groove gets the preference. We have shown that in case of binding of Mg2+ to DNA non electrostatic contributions play a major role. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:5

Dominant-negative effects in prion diseases: insights from molecular dynamics simulations on mouse prion protein chimeras

Mutations in the prion protein (PrP) can cause spontaneous prion diseases in humans (Hu) and animals. In transgenic mice, mutations can determine the susceptibility to the infection of different prion strains. Some of these mutations also show a dominant-negative effect, thus halting the replication process by which wild type mouse (Mo) PrP is converted into Mo scrapie. Using all-atom molecular dynamics (MD) simulations, here we studied the structure of HuPrP, MoPrP, 10 Hu/MoPrP chimeras, and 1 Mo/sheepPrP chimera in explicit solvent. Overall, ∼2 μs of MD were collected. Our findings suggest that the interactions between α1 helix and N-terminal of α3 helix are critical in prion propagation, whereas the β2–α2 loop conformation plays a role in the dominant-negative effect. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:4.

Elucidation by NMR solution of neurotensin in small unilamellar vesicle environment: molecular surveys for neurotensin receptor recognition

Neurotensin (NT) is a tridecapeptide, hormone in the periphery and neurotransmitter in the brain. We used high-resolution nuclear magnetic resonance (NMR) to resolve the three-dimensional structure of NT in a small unilamellar vesicle (SUV) environment. We demonstrate that if the dynamic of the association–dissociation processes of peptide to SUV binding is rapid enough, structural determination can be obtained by solution NMR experiments. Thus, according to the global dynamic of the system, SUVs seem to be an effective model to mimic biological membranes, especially since the lipid composition can be modified or sterols may be added to closely mimic the biological membranes studied. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:2

An insight to the dynamics of conserved water-mediated salt bridge interaction and interdomain recognition in hIMPDH isoforms

Inosine monophosphate dehydrogenase (IMPDH) is involved in de novo biosynthesis pathway of guanosine nucleotide. Type II isoform of this enzyme is selectively upregulated in lymphocytes and chronic myelogenous leukemia (CML) cells, and is an excellent target for antileukemic agent. The molecular dynamics simulation results (15 ns) of three unliganded 1B3O, 1JCN, and 1JR1 structures have clearly revealed that IN, IC (N- and C-terminal of catalytic domains) and C1, C2 (cystathionine-beta-synthase-1 and 2) domains of IMPDH enzyme have been stabilized by six conserved water (center) mediated salt bridge interactions. These conserved water molecules could be involved in interdomain or intradomain recognition, intradomain coupling, and charge transfer processes. The binding propensity of cystathionine-beta-synthase domain to catalytic domain (through conserved water-mediated salt bridges) has provided a new insight to the biochemistry of IMPDH. Stereospecific interaction of IN with C2 domain through conserved water molecule (K109–WII 1–D215/D216) is observed to be unique in the simulated structure of hIMPDH-II. The geometrical/structural consequences and topological feature around the WII 1 water center may be utilized for isoform specific inhibitor design for CML cancer. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:1

Gating and conduction of nano-channel forming proteins: a computational approach

Monitoring conformational changes in ion channels is essential to understand their gating mechanism. Here, we explore the structural dynamics of four outer membrane proteins with different structures and functions in the slowest nonzero modes of vibration. Normal mode analysis was performed on the modified elastic network model of channel in the membrane. According to our results, when membrane proteins were analyzed in the dominant mode, the composed pores, TolC and α-hemolysin showed large motions at the intramembrane β-barrel region while, in other porins, OmpA and OmpF, largest motions observed in the region of external flexible loops. A criterion based on equipartition theorem was used to measure the possible amplitude of vibration in channel forming proteins. The current approach complements theoretical and experimental techniques including HOLE, Molecular Dynamics (MD), and voltage clamp used to address the channel’s structure and dynamics and provides the means to conduct a theoretical simultaneous study of the structure and function of the channel. An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:3