The human Islet amyloid polypeptide (hIAPP or amylin) is a 37-residue peptide hormone that is normally cosecreted with insulin by the pancreatic β-cells. In patients with type 2 diabetes, hIAPP ...deposits as amyloid fibrils in the extracellular spaces of the pancreatic islets. Recent experimental studies show that the intramolecular disulfide bond between Cys2 and Cys7 plays a central role in the process of fibril formation. However, the effect of the disulfide bond on the intrinsic structural properties of monomeric hIAPP is yet to be determined. In this study, we characterize the atomic structure and the thermodynamics of full-length hIAPP in the presence and absence of a disulfide bond using extensive combined Hamiltonian and temperature replica exchange molecular dynamics simulations (HT-REMD) with a coarse grained protein force field. Our simulations show that HT-REMD is more efficient in sampling than temperature REMD. On the basis of a total simulation time of 28 μs, we find that, although native hIAPP (in the presence of a disulfide bond) essentially adopts a disordered conformation in solution, consistent with the signal measured by ultraviolet-circular dichroism (UV-CD) spectroscopy, it also transiently samples α-helical structure for residues 5−16. In comparison with the N-terminal region, the C-terminal region is highly disordered and populates a much lesser content of isolated β-strand conformation for residues 22−26 and 30−35. Moreover, the absence of the disulfide bond greatly decreases the extent of helix formed throughout residues 5−9 in favor of random coil and β-sheet structure. Implications of the stabilization of N-terminal helical structure by disulfide bond on the initialization of hIAPP amyloid formation are discussed.
In this paper, we apply the Matteoli−Mansoori empirical formula for the pair correlation function of simple fluids obeying the Lennard-Jones potential to calculate reduced self-diffusion coefficients ...on the basis of the modified free volume theory. The self-diffusion coefficient thus computed as functions of temperature and density is compared with the molecular dynamics simulation data and the self-diffusion coefficient obtained by the modified free volume theory implemented with the Monte Carlo simulation method for the pair correlation function. We show that the Matteoli−Mansoori empirical formula yields sufficiently accurate self-diffusion coefficients in the supercritical regime, provided that the minimum free volume activating diffusion is estimated with the classical turning point of binary collision at the mean relative kinetic energy 3k B T/2, where k B is the Boltzmann constant and T is the temperature. In the subcritical regime, the empirical formula yields qualitatively correct, but lower values for the self-diffusion coefficients compared with computer simulation values and those from the modified free volume theory implemented with the Monte Carlo simulations for the pair correlation function. However, with a slightly modified critical free volume, the results can be made quite acceptable.
In previous work on the density fluctuation theory of transport coefficients of liquids, it was necessary to use empirical self-diffusion coefficients to calculate the transport coefficients (e.g., ...shear viscosity of carbon dioxide). In this work, the necessity of empirical input of the self-diffusion coefficients in the calculation of shear viscosity is removed, and the theory is thus made a self-contained molecular theory of transport coefficients of liquids, albeit it contains an empirical parameter in the subcritical regime. The required self-diffusion coefficients of liquid carbon dioxide are calculated by using the modified free volume theory for which the generic van der Waals equation of state and Monte Carlo simulations are combined to accurately compute the mean free volume by means of statistical mechanics. They have been computed as a function of density along four different isotherms and isobars. A Lennard-Jones site−site interaction potential was used to model the molecular carbon dioxide interaction. The density and temperature dependence of the theoretical self-diffusion coefficients are shown to be in excellent agreement with experimental data when the minimum critical free volume is identified with the molecular volume. The self-diffusion coefficients thus computed are then used to compute the density and temperature dependence of the shear viscosity of liquid carbon dioxide by employing the density fluctuation theory formula for shear viscosity as reported in an earlier paper (J. Chem. Phys. 2000, 112, 7118). The theoretical shear viscosity is shown to be robust and yields excellent density and temperature dependence for carbon dioxide. The pair correlation function appearing in the theory has been computed by Monte Carlo simulations.
In recent years, much effort has focused on the early stages of aggregation and the formation of amyloid oligomers. Aggregation processes for these proteins are complex and their non-equilibrium ...nature makes any experimental study very difficult. Under these conditions, simulations provide a useful alternative for understanding the dynamics of the early stages of oligomerization. Here, we focus on the non-Aβ amyloid component (NAC) of the monomer, dimer, and trimer of α-synuclein, an important 35-residue sequence involved in the aggregation and fibrillation of this protein associated with Parkinson's disease. Using Hamiltonian and temperature replica exchange molecular dynamics simulations combined with the coarse grained Optimized Potential for Efficient peptide structure Prediction potential, we identify the role of the various regions and the secondary structures for the onset of oligomerization. For this sequence, we clearly observe the passage from α-helix to β-sheet, a characteristic transition of amyloid proteins. More precisely, we find that the NAC monomer is highly structured with two α-helical regions, between residues 2-13 and 19-25. As the dimer and trimer form, β-sheet structures between residues 2-14 and 26-34 appear and rapidly structure the system. The resulting conformations are much more structured than similar dimers and trimers of β-amyloid and amylin proteins and yet display a strong polymorphism at these early stages of aggregation. In addition to its inherent experimental interest, comparison with other sequences shows that NAC could be a very useful numerical model for understanding the onset of aggregation.
Lambert-Eaton myasthenic syndrome (LEMS) is an autoimmune-mediated neuromuscular disease thought to be caused by autoantibodies against P/Q-type voltage-gated calcium channels (VGCCs), which attack ...and reduce the number of VGCCs within transmitter release sites (active zones; AZs) at the neuromuscular junction (NMJ), resulting in neuromuscular weakness. However, patients with LEMS also have antibodies to other neuronal proteins, and about 15% of patients with LEMS are seronegative for antibodies against VGCCs. We hypothesized that a reduction in the number of P/Q-type VGCCs alone is not sufficient to explain LEMS effects on transmitter release. Here, we used a computational model to study a variety of LEMS-mediated effects on AZ organization and transmitter release constrained by electron microscopic, pharmacological, immunohistochemical, voltage imaging, and electrophysiological observations. We show that models of healthy AZs can be modified to predict the transmitter release and short-term facilitation characteristics of LEMS and that in addition to a decrease in the number of AZ VGCCs, disruption in the organization of AZ proteins, a reduction in AZ number, a reduction in the amount of synaptotagmin, and the compensatory expression of L-type channels outside the remaining AZs are important contributors to LEMS-mediated effects on transmitter release. Furthermore, our models predict that antibody-mediated removal of synaptotagmin in combination with disruption in AZ organization alone could mimic LEMS effects without the removal of VGCCs (a seronegative model). Overall, our results suggest that LEMS pathophysiology may be caused by a collection of pathological alterations to AZs at the NMJ, rather than by a simple loss of VGCCs.
We used a computational model of the active zone (AZ) in the mammalian neuromuscular junction to investigate Lambert-Eaton myasthenic syndrome (LEMS) pathophysiology. This model suggests that disruptions in presynaptic active zone organization and protein content (particularly synaptotagmin), beyond the simple removal of presynaptic calcium channels, play an important role in LEMS pathophysiology.
The action potential (AP) waveform controls the opening of voltage-gated calcium channels and contributes to the driving force for calcium ion flux that triggers neurotransmission at presynaptic ...nerve terminals. Although the frog neuromuscular junction (NMJ) has long been a model synapse for the study of neurotransmission, its presynaptic AP waveform has never been directly studied, and thus the AP waveform shape and propagation through this long presynaptic nerve terminal are unknown. Using a fast voltage-sensitive dye, we have imaged the AP waveform from the presynaptic terminal of male and female frog NMJs and shown that the AP is very brief in duration and actively propagated along the entire length of the terminal. Furthermore, based on measured AP waveforms at different regions along the length of the nerve terminal, we show that the terminal is divided into three distinct electrical regions: A beginning region immediately after the last node of Ranvier where the AP is broadest, a middle region with a relatively consistent AP duration, and an end region near the tip of nerve terminal branches where the AP is briefer. We hypothesize that these measured changes in the AP waveform along the length of the motor nerve terminal may explain the proximal-distal gradient in transmitter release previously reported at the frog NMJ.
The AP waveform plays an essential role in determining the behavior of neurotransmission at the presynaptic terminal. Although the frog NMJ is a model synapse for the study of synaptic transmission, there are many unknowns centered around the shape and propagation of its presynaptic AP waveform. Here, we demonstrate that the presynaptic terminal of the frog NMJ has a very brief AP waveform and that the motor nerve terminal contains three distinct electrical regions. We propose that the changes in the AP waveform as it propagates along the terminal can explain the proximal-distal gradient in transmitter release seen in electrophysiological studies.
Abstract only
Background
Claudin‐2 is a paracellular channel. Na and water flow through it driven by an osmotic or Na concentration gradient, suggesting that Na and water pass through a common pore ...and that their fluxes are frictionally coupled. We investigated the mechanism by molecular dynamics simulations.
Methods
The claudin pore was modeled either as a double cone or a cylindrical nanotube with narrowest diameter of 6.5 Å and 6 negatively charged sites in the middle of the pore (each ‐0.2e). A hydrostatic pressure gradient was applied via an additional constant force in the channel axis direction to a layer of water molecules parallel to the membrane interface. Simulations were performed using the CHARMM27 force field and the TIP3P water model, with full electrostatics was included.
Results
With an osmotic gradient, water molecules oriented within the pore in 2 water wires and exhibited directional flow with Pf= 4.7 x 10‐13 cm
3
/s and a water:Na flux ratio of 233. Increasing the magnitude of the negative charge increased permeability only when the charges were smeared evenly along the pore surface. Application of a Na concentration gradient, also drove water flow with water:Na flux of 14.6. Increasing the pore diameter to 7Å also increased water flow significantly by allowing 3 water wires to form.
Conclusions
Our findings confirm that water can permeate through claudin‐2 via a water wire mechanism, demonstrate that water permeability is dependent on the magnitude and distribution of intrapore charges and on the pore diameter, and provide a molecular explanation for the frictional coupling of Na and water fluxes within the pore.
The Amyloid-beta protein is related to Alzheimer's disease, and various experiments have shown that oligomers as small as the dimer are cytotoxic. Two alloforms are mainly produced: Aβ(1-40) and ...Aβ(1-42). They have very different oligomer distributions, and it was recently suggested, from experimental studies, that this variation may originate from structural differences in their dimer structures. Little structural information is available on the Aβ dimer, however, and to complement experimental observations, we simulated the folding of the wild-type Aβ(1-40) and Aβ(1-42) dimers as well as the mutated Aβ(1-40)(D23N) dimer using an accurate coarse-grained force field coupled to Hamiltonian-temperature replica exchange molecular dynamics. The D23N variant impedes the salt-bridge formation between D23 and K28 seen in the wild-type Aβ, leading to very different fibrillation properties and final amyloid fibrils. Our results show that the Aβ(1-42) dimer has a higher propensity than the Aβ(1-40) dimer to form β-strands at the central hydrophobic core (residues 17-21) and at the C-terminal (residues 30-42), which are two segments crucial to the oligomerization of Aβ. The free energy landscape of the Aβ(1-42) dimer is also broader and more complex than that of the Aβ(1-40) dimer. Interestingly, D23N also impacts the free energy landscape by increasing the population of configurations with higher β-strand propensities when compared against Aβ(40). In addition, while Aβ(1-40)(D23N) displays a higher β-strand propensity at the C-terminal, its solvent accessibility does not change with respect to the wild-type sequence. Overall, our results show the strong impact of the two amino acids Ile41-Ala42 and the salt-bridge D23-K28 on the folding of the Aβ dimer.