Investigating the size and dynamics of voltage-gated sodium channel fenestrations: A molecular dynamics study

Eukaryotic voltage-gated sodium channels (VGSCs) are essential for the initiation and propagation of action potentials in electrically excitable cells, and are important pharmaceutical targets for the treatment of neurological disorders such as epilepsy, cardiac arrhythmias, and chronic pain. Evidence suggests that small, hydrophobic, VGSC-blocking drugs can gain access to binding residues within the central cavity of these channels by passing through lateral, lipid-filled “fenestrations” which run between the exterior of the protein and its central pore. Here, we use molecular dynamics simulations to investigate how the size and shape of fenestrations change over time in several bacterial VGSC models and a homology model of Nav1.4. We show that over the course of the simulations, the size of the fenestrations is primarily influenced by rapid protein motions, such as amino acid side-chain rotation, and highlight that differences between fenestration bottleneck-contributing residues are the primary cause of variations in fenestration size between the 6 bacterial models. In the eukaryotic channel model, 2 fenestrations are wide, but 2 are narrow due to differences in the amino acid sequence in the 4 domains. Lipid molecules are found to influence the size of the fenestrations by protruding acyl chains into the fenestrations and displacing amino acid side-chains. Together, the results suggest that fenestrations provide viable pathways for small, flexible, hydrophobic drugs.     

Characterizing fenestration size in sodium channel subtypes and their accessibility to inhibitors

Voltage-gated sodium channels (Nav) underlie the electrical activity of nerve and muscle cells. Humans have nine different subtypes of these channels, which are the target of small-molecule inhibitors commonly used to treat a range of conditions. Structural studies have identified four lateral fenestrations within the Nav pore module that have been shown to influence Nav pore blocker access during resting-state inhibition. However, the structural differences among the nine subtypes are still unclear. In particular, the dimensions of the four individual fenestrations across the Nav subtypes and their differential accessibility to pore blockers is yet to be characterized. To address this, we applied classical molecular dynamics simulations to study the recently published structures of Nav1.1, Nav1.2, Nav1.4, Nav1.5, and Nav1.7. Although there is significant variability in the bottleneck sizes of the Nav fenestrations, the subtypes follow a common pattern, with wider DI-II and DIII-IV fenestrations, a more restricted DII-III fenestration, and the most restricted DI-IV fenestration. We further identify the key bottleneck residues in each fenestration and show that the motions of aromatic residue sidechains govern the bottleneck radii. Well-tempered metadynamics simulations of Nav1.4 and Nav1.5 in the presence of the pore blocker lidocaine also support the DI-II fenestration being the most likely access route for drugs. Our computational results provide a foundation for future in vitro experiments examining the route of drug access to sodium channels. Understanding the fenestrations and their accessibility to drugs is critical for future analyses of diseases mutations across different sodium channel subtypes, with the potential to inform pharmacological development of resting-state inhibitors and subtype-selective drug design.     

Influence of lipids on the hydrophobic barrier within the pore of the TWIK-1 K2P channel

Several recent ion channel structures have revealed large side portals, or ‘fenestrations’ at the interface between their transmembrane helices that potentially expose the ion conduction pathway to the lipid core of the bilayer. In a recent study we demonstrated that functional activity of the TWIK-1 K2P channel is influenced by the presence of hydrophobic residues deep within the inner pore. These residues are located near the fenestrations in the TWIK-1 structure and promote dewetting of the pore by forming a hydrophobic barrier to ion conduction. During our previous MD simulations, lipid tails were observed to enter these fenestrations. In this addendum to that study, we investigate lipid contribution to the dewetting process. Our results demonstrate that lipid tails from both the upper and lower leaflets can occupy the fenestrations and partially penetrate into the pore. The lipid tails do not sterically occlude the pore, but there is an inverse correlation between the presence of water within the hydrophobic barrier and the number of lipids tails within the lining of the pore. However, dewetting still occurs in the absence of lipids tails, and pore hydration appears to be determined primarily by those side-chains lining the narrowest part of the pore cavity.     

Influence of lipids on the hydrophobic barrier within the pore of the TWIK-1 K2P channel

Several recent ion channel structures have revealed large side portals, or ‘fenestrations’ at the interface between their transmembrane helices that potentially expose the ion conduction pathway to the lipid core of the bilayer. In a recent study we demonstrated that functional activity of the TWIK-1 K2P channel is influenced by the presence of hydrophobic residues deep within the inner pore. These residues are located near the fenestrations in the TWIK-1 structure and promote dewetting of the pore by forming a hydrophobic barrier to ion conduction. During our previous MD simulations, lipid tails were observed to enter these fenestrations. In this addendum to that study, we investigate lipid contribution to the dewetting process. Our results demonstrate that lipid tails from both the upper and lower leaflets can occupy the fenestrations and partially penetrate into the pore. The lipid tails do not sterically occlude the pore, but there is an inverse correlation between the presence of water within the hydrophobic barrier and the number of lipids tails within the lining of the pore. However, dewetting still occurs in the absence of lipids tails, and pore hydration appears to be determined primarily by those side-chains lining the narrowest part of the pore cavity.     

BATMAS30: amino acid substitution matrix for alignment of bacterial transporters

Aligned amino acid sequences of three functionally independent samples of transmembrane (TM) transport proteins have been analyzed. The concept of TM-kernel is proposed as the most probable transmembrane region of a sequence. The average amino acid composition of TM-kernels differs from the published amino acid composition of transmembrane segments. TM-kernels contain more alanines, glycines, and less polar, charged, and aromatic residues in contrast to non-TM-proteins. There are also differences between TM-kernels of bacterial and eukaryotic proteins. We have constructed amino acid substitution matrices for bacterial TM-kernels, named the BATMAS (BActerial Transmembrane MAtrix of Substitutions) series. In TM-kernels, polar and charged residues, as well as proline and tyrosine, are highly conserved, whereas there are more substitutions within the group of hydrophobic residues, in contrast to non-TM-proteins that have fewer, relatively more conserved, hydrophobic residues. These results demonstrate that alignment of transmembrane proteins should be based on at least two amino acid substitution matrices, one for loops (e.g., the BLOSUM series) and one for TM-segments (the BATMAS series), and the choice of the TM-matrix should be different for eukaryotic and bacterial proteins.     

RNA editing modulates the binding of drugs and highly unsaturated fatty acids to the open pore of Kv potassium channels

The time course of inactivation of voltage‐activated potassium (Kv) channels is an important determinant of the firing rate of neurons. In many Kv channels highly unsaturated lipids as arachidonic acid, docosahexaenoic acid and anandamide can induce fast inactivation. We found that these lipids interact with hydrophobic residues lining the inner cavity of the pore. We analysed the effects of these lipids on Kv1.1 current kinetics and their competition with intracellular tetraethylammonium and Kvβ subunits. Our data suggest that inactivation most likely represents occlusion of the permeation pathway, similar to drugs that produce ‘open‐channel block’. Open‐channel block by drugs and lipids was strongly reduced in Kv1.1 channels whose amino acid sequence was altered by RNA editing in the pore cavity, and in Kv1.x heteromeric channels containing edited Kv1.1 subunits. We show that differential editing of Kv1.1 channels in different regions of the brain can profoundly alter the pharmacology of Kv1.x channels. Our findings provide a mechanistic understanding of lipid‐induced inactivation and establish RNA editing as a mechanism to induce drug and lipid resistance in Kv channels.     

Solution structure of omega-grammotoxin SIA,a gating modifier of P/Q and N-type Ca(2+) channel

omega-Grammotoxin SIA (GrTx) is a 36 amino acid residue protein toxin from spider venom that inhibits P/Q and N-type voltage-gated Ca(2+) channels by modifying voltage-dependent gating. We determined the three-dimensional structure of GrTx using NMR spectroscopy. The toxin adopts an "inhibitor cystine knot" motif composed of two beta-strands (Leu19-Cys21 and Cys30-Trp32) and a beta-bulge (Trp6, Gly7-Cys30) with a +2x, -1 topology, which are connected by four chain reversals. Although GrTx was originally identified as an inhibitor of voltage-gated Ca(2+) channel, it also binds to K(+) channels with lower affinity. A similar cross-reaction was observed for Hanatoxin1 (HaTx), which binds to the voltage-sensing domains of K(+) and Ca(2+) channels with different affinities. A detailed comparison of the GrTx and HaTx structures identifies a conserved face containing a large hydrophobic patch surrounded by positively charged residues. The slight differences in the surface shape, which result from the orientation of the surface aromatic residues and/or the distribution of the charged residues, may explain the differences in the binding affinity of these gating modifiers with different voltage-gated ion channels.     

Homology Models of the Tetramerization Domain of Six Eukaryotic Voltage-gated Potassium Channels Kv1.1-Kv1.6

Abstract

The homology models of the tetramerization (T1) domain of six eukaryotic potassium channels, Kv1.1-Kv1.6, were constructed based on the crystal structure of the Shaker T1 domain. The results of amino acid sequence alignment indicate that the T1 domains of these K⁺ channels are highly conserved, with the similarities varying from 77% between Shaker and Kv1.6 to 93% between Kv1.2 and Kv1.3. The homology models reveal that the T1 domains of these Kv channels exhibit similar folds as those of Shaker K⁺ channel. These models also show that each T1 monomer consists of three distinct layers, with N-terminal layer 1 and C- terminal layer 3 facing the cytoplasm and the membrane, respectively. Layer 2 exhibits the highest structural conservation because it is located around the central hydrophobic core. For each Kv channel, four identical subunits assemble into the homotetramer architecture around a four-fold axis through the hydrogen bonds and salt bridges formed by 15 highly conserved polar residues. The narrowest opening of the pore is formed by the four conserved residues corresponding to R115 of the Shaker T1 domain. The homology models of these Kv T1 domains provide particularly attractive targets for further structure-based studies.     

Homology models of the tetramerization domain of six eukaryotic voltage-gated potassium channels Kv1.1-Kv1.6

The homology models of the tetramerization (T1) domain of six eukaryotic potassium channels, Kv1.1-Kv1.6, were constructed based on the crystal structure of the Shaker T1 domain. The results of amino acid sequence alignment indicate that the T1 domains of these K+ channels are highly conserved, with the similarities varying from 77% between Shaker and Kv1.6 to 93% between Kv1.2 and Kv1.3. The homology models reveal that the T1 domains of these Kv channels exhibit similar folds as those of Shaker K+ channel. These models also show that each T1 monomer consists of three distinct layers, with N-terminal layer 1 and C-terminal layer 3 facing the cytoplasm and the membrane, respectively. Layer 2 exhibits the highest structural conservation because it is located around the central hydrophobic core. For each Kv channel, four identical subunits assemble into the homotetramer architecture around a four-fold axis through the hydrogen bonds and salt bridges formed by 15 highly conserved polar residues. The narrowest opening of the pore is formed by the four conserved residues corresponding to R115 of the Shaker T1 domain. The homology models of these Kv T1 domains provide particularly attractive targets for further structure-based studies.     

Bulky aromatic amino acids increase the antibacterial activity of 15-residue bovine lactoferricin derivatives.

A model peptide, FKCRRWQWRMKKLGA, residues 17-31 of bovine lactoferricin, has been subjected to structure-antibacterial activity relationship studies. The two Trp residues are very important for antibacterial activity, and analogue studies have demonstrated the significance of the size, shape and aromatic character of the side chains. In the current study we have replaced Trp residues in the model peptide with bulky aromatic amino acids to elucidate further the importance of size and shape. The counterproductive Cys residue in position 3 was also replaced by these aromatic amino acids. The largest aromatic amino acids employed resulted in the most active peptides. The peptides containing these hydrophobic residues were generally more active against Staphylococcus aureus than against Escherichia coli, indicating that the bacterial specificity as well as the antibacterial efficiency can be altered by employing large hydrophobic aromatic amino acid residues.     

Ion channel gating: insights via molecular simulations

Ion channels are gated, i.e. they can switch conformation between a closed and an open state. Molecular dynamics simulations may be used to study the conformational dynamics of ion channels and of simple channel models. Simulations on model nanopores reveal that a narrow (<4 A) hydrophobic region can form a functionally closed gate in the channel and can be opened by either a small (approximately 1 A) increase in pore radius or an increase in polarity. Modelling and simulation studies confirm the importance of hydrophobic gating in K channels, and support a model in which hinge-bending of the pore-lining M2 (or S6 in Kv channels) helices underlies channel gating. Simulations of a simple outer membrane protein, OmpA, indicate that a gate may also be formed by interactions of charged side chains within a pore, as is also the case in ClC channels.     

Monte Carlo studies on equilibrium globular protein folding. II. Beta-barrel globular protein models

In the context of dynamic Monte Carlo simulations on a model protein confined to a tetrahedral lattice, the interplay of protein size and tertiary structure, and the requirements for an all-or-none transition to a unique native state, are investigated. Small model proteins having a primary sequence consisting of a central bend neutral region flanked by two tails having an alternating hydrophobic/hydrophilic pattern of residues are seen to undergo a continuous transition to a beta-hairpin collapsed state. On increasing the length of the tails, the beta-hairpin structural motif is found to be in equilibrium with a four-member beta-barrel. Further increase of the tail length results in the shift of the structural equilibrium to the four-member beta-barrel. The random coil to beta-barrel transition is of an all-or-none character, but while the central turn is always the desired native bend, the location of the turns involving the two external strands is variable. That is, beta-barrels having the external stands that are two residues out of register are also observed in the transition region. Introduction into the primary sequence of two additional regions that are at the very least neutral toward turn formation produces an all-or-none transition to the unique, native, four-member beta-barrel. Various factors that can augment the stability of the native conformation are explored. Overall, these folding simulations strongly indicate that the general rules of globular protein folding are rather robust--namely, one requires a general pattern of hydrophobic/hydrophilic residues that allow the protein to have a well-defined interior and exterior and the presence of regions in the amino acid sequence that at the very least are locally indifferent to turn formation. Since no site-specific interactions between hydrophobic and hydrophilic residues are required to produce a unique four-member beta-barrel, these simulations strongly suggest that site specificity is involved in structural fine-tuning.     

Molecular Analysis of the Putative Inactivation Particle in the Inactivation Gate of Brain Type IIA Na+ Channels

Fast Na⁺ channel inactivation is thought to involve binding of phenylalanine 1489 in the hydrophobic cluster IFM in L_III-IV of the rat brain type IIA Na⁺ channel. We have analyzed macroscopic and single channel currents from Na⁺ channels with mutations within and adjacent to hydrophobic clusters in L_III-IV. Substitution of F1489 by a series of amino acids disrupted inactivation to different extents. The degree of disruption was closely correlated with the hydrophilicity of the amino acid at position 1489. These mutations dramatically destabilized the inactivated state and also significantly slowed the entry into the inactivated state, consistent with the idea that F1489 forms a hydrophobic interaction with a putative receptor during the fast inactivation process. Substitution of a phe residue at position 1488 or 1490 in mutants lacking F1489 did not restore normal inactivation, indicating that precise location of F1489 is critical for its function. Mutations of T1491 disrupted inactivation substantially, with large effects on the stability of the inactivated state and smaller effects on the rate of entry into the inactivated state. Mutations of several other hydrophobic residues did not destabilize the inactivated state at depolarized potentials, indicating that the effects of mutations at F1489 and T1491 are specific. The double mutant YY1497/8QQ slowed macroscopic inactivation at all potentials and accelerated recovery from inactivation at negative membrane potentials. Some of these mutations in L_III-IV also affected the latency to first opening, indicating coupling between L_III-IV and channel activation. Our results show that the amino acid residues of the IFM hydrophobic cluster and the adjacent T1491 are unique in contributing to the stability of the inactivated state, consistent with the designation of these residues as components of the inactivation particle responsible for fast inactivation of Na⁺ channels.     

A set of homology models of pore loop domain of six eukaryotic voltage-gated potassium channels Kv1.1-Kv1.6

Homology models of the pore loop domain of six eukaryotic potassium channels, Kv1.1-Kv1.6, were generated based on the crystallographic structure of KcsA. The results of amino acid sequence alignment indicate that these Kv channels are composed of two structurally and functionally independent domains: the N-terminal 'voltage sensor' domain and the C-terminal 'pore loop' domain. The homology models reveal that the pore loop domains of these Kv channels exhibit similar folds to those of KcsA. The structural features and specific packing of aromatic residues around the selectivity filter of these Kv channels are nearly identical to those of KcsA, whereas most of the structural variations occur in the turret as well as in the inner and outer helices. The distribution of polar and nonpolar side chains on the surfaces of the KcsA and Kv channels reveals that they exhibit a segregation of side chains common to most integral membrane proteins. As the hydrogen bond between Glu71 and Asp80 in KcsA plays an important role in stabilizing the channel, the substituted Val residue in the Kv family corresponding to Glu71 of KcsA stabilizes the channel by making hydrophobic contact with Tyr residue from the signature sequence of the selectivity filter. The homology models of these Kv channels provide particularly attractive subjects for further structure-based studies.     

Block of Brain Sodium Channels by Peptide Mimetics of the Isoleucine,Phenylalanine, and Methionine (IFM) Motif from the Inactivation Gate

Inactivation of sodium channels is thought to be mediated by an inactivation gate formed by the intracellular loop connecting domains III and IV. A hydrophobic motif containing the amino acid sequence isoleucine, phenylalanine, and methionine (IFM) is required for the inactivation process. Peptides containing the IFM motif, when applied to the cytoplasmic side of these channels, produce two types of block: fast block, which resembles the inactivation process, and slow, use-dependent block stimulated by strong depolarizing pulses. Fast block by the peptide ac-KIFMK-NH₂, measured on sodium channels whose inactivation was slowed by the α-scorpion toxin from Leiurus quinquestriatus (LqTx), was reversed with a time constant of 0.9 ms upon repolarization. In contrast, control and LqTx-modified sodium channels were slower to recover from use-dependent block. For fast block, linear peptides of three to six amino acid residues containing the IFM motif and two positive charges were more effective than peptides with one positive charge, whereas uncharged IFM peptides were ineffective. Substitution of the IFM residues in the peptide ac-KIFMK-NH₂ with smaller, less hydrophobic residues prevented fast block. The positively charged tripeptide IFM-NH₂ did not cause appreciable fast block, but the divalent cation IFM-NH(CH₂)₂NH₂ was as effective as the pentapeptide ac-KIFMK-NH₂. The constrained peptide cyclic KIFMK containing two positive charges did not cause fast block. These results indicate that the position of the positive charges is unimportant, but flexibility or conformation of the IFM-containing peptide is important to allow fast block. Slow, use-dependent block was observed with IFM-containing peptides of three to six residues having one or two positive charges, but not with dipeptides or phenylalanine-amide. In contrast to its lack of fast block, cyclic KIFMK was an effective use-dependent blocker. Substitutions of amino acid residues in the tripeptide IFM-NH₂ showed that large hydrophobic residues are preferred in all three positions for slow, use-dependent block. However, substitution of the large hydrophobic residue diphenylalanine or the constrained residues phenylglycine or tetrahydroisoquinoline for phe decreased potency, suggesting that this phe residue must be able to enter a restricted hydrophobic pocket during the binding of IFM peptides. Together, the results on fast block and slow, use-dependent block indicate that IFM peptides form two distinct complexes of different stability and structural specificity with receptor site(s) on the sodium channel. It is proposed that fast block represents binding of these peptides to the inactivation gate receptor, while slow, use-dependent block represents deeper binding of the IFM peptides in the pore.     

Molecular simulations study of novel 1,4‐dihydropyridines derivatives with a high selectivity for Cav3.1 calcium channel

1,4‐Dihydropyridines (DHPs) have been developed to treat hypertension, angina, and nerve system disease. They are thought to mainly target the L‐type calcium channels, but low selectivity prompts them to block Cav1.2 and Cav3.1 channels simultaneously. Recently, some novel DHPs with different hydrophobic groups have been synthesized and among them M12 has a higher selectivity for Cav3.1. However, the structural information about Cav3.1‐DHPs complexes is not available in the experiment. Thus, we combined homology modeling, molecular docking, molecular dynamics simulations, and binding free energy calculations to quantitatively elucidate the inhibition mechanism of DHPs. The calculated results indicate that our model is in excellent agreement with experimental results. On the basis of conformational analysis, we identify the main interactions between DHPs and calcium channels and further elaborate on the different selectivity of ligands from the micro perspective. In conjunction with energy distribution, we propose that the binding sites of Cav3.1‐DHPs is characterized by several interspersed hydrophobic amino acid residues on the IIIS6 and IVS6 segments. We also speculate the favorable function groups on prospective DHPs. Besides, our model provides important information for further mutagenesis experiments.     

Structure,Dynamics, and Substrate Specificity of the OprO Porin from Pseudomonas aeruginosa

The outer membrane (OM) of Gram-negative bacteria functions as a selective permeability barrier between cell and environment. For nutrient acquisition, the OM contains a number of channels that mediate uptake of small molecules by diffusion. Many of these channels are specific, i.e., they prefer certain substrates over others. In electrophysiological experiments, the OM channels OprP and OprO from Pseudomonas aeruginosa show a specificity for phosphate and diphosphate, respectively. In this study we use x-ray crystallography, free-energy molecular dynamics (MD) simulations, and electrophysiology to uncover the atomic basis for the different substrate specificity of these highly similar channels. A structural analysis of OprP and OprO revealed two crucial differences in the central constriction region. In OprP there are two tyrosine residues, Y62 and Y114, whereas the corresponding residues in OprO are phenylalanine F62 and aspartate D114. To probe the importance of these two residues in generating the different substrate specificities, the double mutants were generated in silico and in vitro. Applied-field MD simulations and electrophysiological experiments demonstrated that the double mutations interchange the phosphate and diphosphate specificities of OprP and OprO. Our findings outline a possible strategy to rationally design channel specificity by modification of a small number of residues that may be applicable to other pores as well.     

Helix packing moments reveal diversity and conservation in membrane protein structure

Helical membrane proteins are more tightly packed and the packing interactions are more diverse than those found in helical soluble proteins. Based on a linear correlation between amino acid packing values and interhelical propensity, we propose the concept of a helix packing moment to predict the orientation of helices in helical membrane proteins and membrane protein complexes. We show that the helix packing moment correlates with the helix interfaces of helix dimers of single pass membrane proteins of known structure. Helix packing moments are also shown to help identify the packing interfaces in membrane proteins with multiple transmembrane helices, where a single helix can have multiple contact surfaces. Analyses are described on class A G protein-coupled receptors (GPCRs) with seven transmembrane helices. We show that the helix packing moments are conserved across the class A family of GPCRs and correspond to key structural contacts in rhodopsin. These contacts are distinct from the highly conserved signature motifs of GPCRs and have not previously been recognized. The specific amino acid types involved in these contacts, however, are not necessarily conserved between subfamilies of GPCRs, indicating that the same protein architecture can be supported by a diverse set of interactions. In GPCRs, as well as membrane channels and transporters, amino acid residues with small side-chains (Gly, Ala, Ser, Cys) allow tight helix packing by mediating strong van der Waals interactions between helices. Closely packed helices, in turn, facilitate interhelical hydrogen bonding of both weakly polar (Ser, Thr, Cys) and strongly polar (Asn, Gln, Glu, Asp, His, Arg, Lys) amino acid residues. We propose the use of the helix packing moment as a complementary tool to the helical hydrophobic moment in the analysis of transmembrane sequences.     

Probing structural features of Alzheimer's amyloid-β pores in bilayers using site-specific amino acid substitutions

A current hypothesis for the pathology of Alzheimer's disease (AD) proposes that amyloid-β (Aβ) peptides induce uncontrolled, neurotoxic ion flux across cellular membranes. The mechanism of ion flux is not fully understood because no experiment-based Aβ channel structures at atomic resolution are currently available (only a few polymorphic states have been predicted by computational models). Structural models and experimental evidence lend support to the view that the Aβ channel is an assembly of loosely associated mobile β-sheet subunits. Here, using planar lipid bilayers and molecular dynamics (MD) simulations, we show that amino acid substitutions can be used to infer which residues are essential for channel structure. We created two Aβ(1-42) peptides with point mutations: F19P and F20C. The substitution of Phe19 with Pro inhibited channel conductance. MD simulation suggests a collapsed pore of F19P channels at the lower bilayer leaflet. The kinks at the Pro residues in the pore-lining β-strands induce blockage of the solvated pore by the N-termini of the chains. The cysteine mutant is capable of forming channels, and the conductance behavior of F20C channels is similar to that of the wild type. Overall, the mutational analysis of the channel activity performed in this work tests the proposition that the channels consist of a β-sheet rich organization, with the charged/polar central strand containing the mutation sites lining the pore, and the C-terminal strands facing the hydrophobic lipid tails. A detailed understanding of channel formation and its structure should aid studies of drug design aiming to control unregulated Aβ-dependent ion fluxes.     

Amino acid preferences of small, naturally occurring polypeptides

An analysis of amino acid composition of small, naturally occurring peptides ranging in size from 3 to 50 residues has been carried out. The purpose of the study is to determine whether differential trends in amino acid usage exist for small peptides compared to larger polypeptides and proteins. Results indicate that Cys, Trp, and Phe are substantially more frequent in peptides compared to their abundance in proteins at large. Aliphatic hydrophobic residues, particularly Leu and Ile, are somewhat underrepresented, while the frequency of Glu is significantly reduced. The shorter peptides are also more frequently neutral and become increasingly charged as their size increases.