Simple biophysics underpins collective conformations of the intrinsically disordered proteins of the Nuclear Pore Complex.

Nuclear Pore Complexes (NPCs) are key cellular transporter that control nucleocytoplasmic transport in eukaryotic cells, but its transport mechanism is still not understood. The centerpiece of NPC transport is the assembly of intrinsically disordered polypeptides, known as FG nucleoporins, lining its passageway. Their conformations and collective dynamics during transport are difficult to assess in vivo. In vitro investigations provide partially conflicting results, lending support to different models of transport, which invoke various conformational transitions of the FG nucleoporins induced by the cargo-carrying transport proteins. We show that the spatial organization of FG nucleoporin assemblies with the transport proteins can be understood within a first principles biophysical model with a minimal number of key physical variables, such as the average protein interaction strengths and spatial densities. These results address some of the outstanding controversies and suggest how molecularly divergent NPCs in different species can perform essentially the same function.

Calculation of iron transport through human H-chain ferritin.

Influx of ferrous ions from the cytoplasm through 3-fold pores in the shell of ferritin protein is computed using a 3-dimensional Poisson-Nernst-Planck electrodiffusion model, with inputs such as the pore structure and the diffusivity profile of permeant Fe(2+) ions extracted from all-atom molecular dynamics (MD) simulations. These calculations successfully reproduce experimental estimates of the transit time of Fe(2+) through the ferritin coat, which is on the millisecond time scale and hence much too long to be directly simulated via all-atom MD. This is also much longer than the typical time scale for ion transit in standard membrane spanning ion channels whose pores bear structural similarity to that of the 3-fold ferritin pore. The slow time scale for Fe(2+) transport through ferritin pores is traced to two features that distinguish the ferritin pore system from standard ion channels, namely, (i) very low concentration of cytoplasmic Fe(2+) under physiological conditions and (ii) very small internal diffusion coefficients for ions inside the ferritin pore resulting from factors that include the divalent nature of Fe(2+) and two rings of negatively charged amino acids surrounding a narrow geometric obstruction within the ferritin pore interior.

Metal binding sites of human H-chain ferritin and iron transport mechanism to the ferroxidase sites: a molecular dynamics simulation study.

We study via all atom classical molecular dynamics (MD) simulation the process of uptake of ferrous ions (Fe(2+)) into the human ferritin protein and the catalytic ferroxidase sites via pores ("channels") in the interior of the protein. We observe that the three-fold hydrophilic channels serve as the main entrance pathway for the Fe(2+) ions. The binding sites along the ion pathway are investigated. Two strong binding sites, at the Asp131 and Glu134 residues and two weak binding sites, at the His118 and Cys130 are observed inside the three-fold channel. We also identify an explicit pathway for an ion exiting the channel into the central core of the protein as it moves to the ferroxidase site. The diffusion of an Fe(2+) ion from the inner opening of the channel to a ferroxidase site located in the interior region of the protein coat is assisted by Thr135, His136 and Tyr137. The Fe(2+) ion binds preferentially to site A of the ferroxidase site.

Energetics and ion permeation characteristics in a glutamate-gated chloride (GluCl) receptor channel.

An invertebrate glutamate-gated chloride channel (GluCl) has recently been crystallized in an open-pore state. This channel is homologous to the human Cys-loop receptor family of pentameric ligand-gated ion channels, including anion-selective GlyR and GABAR and cation-selective nAChR and 5HT(3). We implemented molecular dynamics (MD) in conjunction with an elastic network model to perturb the X-ray structure of GluCl and investigated the open channel stability and its ion permeation characteristics. Our study suggests that TM2 helical tilting may close GluCl near the hydrophobic constriction L254 (L9'), similar to its cation-selective homologues. Ion permeation characteristics were determined by Brownian dynamics simulations using a hybrid MD/continuum electrostatics approach to evaluate the free energy profiles for ion transport. Near the selectivity filter region (P243 or P-2'), the free energy barrier for Na(+) transport is over 4 k(B)T higher than that for Cl(-), indicating anion selectivity of the channel. Furthermore, three layers of positivity charged rings in the extracellular domain also contribute to charge selectivity and facilitate Cl(-) permeability over Na(+). Collectively, the charge selectivity of GluCl may be determined by overall electrostatic and ion dehydration effects, perhaps not deriving from a single region of the channel (the selectivity filter region near the intracellular entrance).

Calcium inhibits paracellular sodium conductance through claudin-2 by competitive binding.

Claudins form paracellular pores at the tight junction in epithelial cells. Profound depletion of extracellular calcium is well known to cause loosening of the tight junction with loss of transepithelial resistance. However, moderate variations in calcium concentrations within the physiological range can also regulate transepithelial permeability. To investigate the underlying molecular mechanisms, we studied the effects of calcium on the permeability of claudin-2, expressed in an inducible MDCK I cell line. We found that in the physiological range, calcium acts as a reversible inhibitor of the total conductance and Na(+) permeability of claudin-2, without causing changes in tight junction structure. The effect of calcium is enhanced at low Na(+) concentrations, consistent with a competitive effect. Furthermore, mutation of an intrapore negatively charged binding site, Asp-65, to asparagine partially abrogated the inhibitory effect of calcium. This suggests that calcium competes with Na(+) for binding to Asp-65. Other polyvalent cations had similar effects, including La(3+), which caused severe and irreversible inhibition of conductance. Brownian dynamics simulations demonstrated that such inhibition can be explained if Asp-65 has a relatively high charge density, thus favoring binding of Ca(2+) over that of Na(+), reducing Ca(2+) permeation by inhibiting its dissociation from this site, and decreasing Na(+) conductance through repulsive electrostatic interaction with Ca(2+). These findings may explain why hypercalcemia inhibits Na(+) reabsorption in the proximal tubule of the kidney.

In silico models for the human alpha4beta2 nicotinic acetylcholine receptor.

The neuronal alpha4beta2 nicotinic acetylcholine receptor (nAChR) is one of the most widely expressed nAChR subtypes in the brain. Its subunits have high sequence identity (54 and 46% for alpha4 and beta2, respectively) with alpha and beta subunits in Torpedo nAChR. Using the known structure of the Torpedo nAChR as a template, the closed-channel structure of the alpha4beta2 nAChR was constructed through homology modeling. Normal-mode analysis was performed on this closed structure and the resulting lowest frequency mode was applied to it for a "twist-to-open" motion, which increased the minimum pore radius from 2.7 to 3.4 A and generated an open-channel model. Nicotine could bind to the predicted agonist binding sites in the open-channel model but not in the closed one. Both models were subsequently equilibrated in a ternary lipid mixture via extensive molecular dynamics (MD) simulations. Over the course of 11 ns MD simulations, the open channel remained open with filled water, but the closed channel showed a much lower water density at its hydrophobic gate comprised of residues alpha4-V259 and alpha4-L263 and their homologous residues in the beta2 subunits. Brownian dynamics simulations of Na+ permeation through the open channel demonstrated a current-voltage relationship that was consistent with experimental data on the conducting state of alpha4beta2 nAChR. Besides establishment of the well-equilibrated closed- and open-channel alpha4beta2 structural models, the MD simulations on these models provided valuable insights into critical factors that potentially modulate channel gating. Rotation and tilting of TM2 helices led to changes in orientations of pore-lining residue side chains. Without concerted movement, the reorientation of one or two hydrophobic side chains could be enough for channel opening. The closed- and open-channel structures exhibited distinct patterns of electrostatic interactions at the interface of extracellular and transmembrane domains that might regulate the signal propagation of agonist binding to channel opening. A potential prominent role of the beta2 subunit in channel gating was also elucidated in the study.

A rigid-body Newtonian propagation scheme based on instantaneous decomposition into rotation and translation blocks.

The rotation and translation block (RTB) method of Durand et al. [Biopolymers 34, 759 (1994)] and Tama et al. [Proteins 41, 1 (2000)] provides an appealing way to calculate low-frequency normal modes of large biomolecules by restricting the space of motions to exclude internal motions of preselected rigid fragments within the molecule. These fragments are modeled essentially as rigid bodies and the need to calculate high-frequency relative motions of the atoms that form them is obviated in a natural way. Here we extend the RTB approach into a method for computing the classical (Newtonian) dynamics of a biomolecule, or any large molecule, with effective rigid-body constraints applied to a prechosen set of internal molecular fragments. This method, to be termed RTB dynamics, is easy to implement, conserves the total energy of the system, does not require the construction of the matrix of second spatial derivatives of the potential-energy function (Hessian matrix), and can be used to compute the classical dynamics of a system moving in an arbitrary anharmonic force field. An elementary numerical application to signal propagation in the small membrane-bound polypeptide gramicidin-A is presented for illustration purposes.