Differential dynamics of the serotonin(1A) receptor in membrane bilayers of varying cholesterol content revealed by all atom molecular dynamics simulation

The serotonin(1A) receptor belongs to the superfamily of G protein-coupled receptors (GPCRs) and is a potential drug target in neuropsychiatric disorders. The receptor has been shown to require membrane cholesterol for its organization, dynamics and function. Although recent work suggests a close interaction of cholesterol with the receptor, the structural integrity of the serotonin(1A) receptor in the presence of cholesterol has not been explored. In this work, we have carried out all atom molecular dynamics simulations, totaling to 3s, to analyze the effect of cholesterol on the structure and dynamics of the serotonin(1A) receptor. Our results show that the presence of physiologically relevant concentration of membrane cholesterol alters conformational dynamics of the serotonin(1A) receptor and, on an average lowers conformational fluctuations. Our results show that, in general, transmembrane helix VII is most affected by the absence of membrane cholesterol. These results are in overall agreement with experimental data showing enhancement of GPCR stability in the presence of membrane cholesterol. Our results constitute a molecular level understanding of GPCR-cholesterol interaction, and represent an important step in our overall understanding of GPCR function in health and disease.

Folding of Protein L with Implications for Collapse in the Denatured State Ensemble

A fundamental question in protein folding is whether the coil to globule collapse transition occurs during the initial stages of folding (burst phase) or simultaneously with the protein folding transition. Single molecule fluorescence resonance energy transfer (FRET) and small-angle X-ray scattering (SAXS) experiments disagree on whether Protein L collapse transition occurs during the burst phase of folding. We study Protein L folding using a coarse-grained model and molecular dynamics simulations. The collapse transition in Protein L is found to be concomitant with the folding transition. In the burst phase of folding, we find that FRET experiments overestimate radius of gyration, R-g, of the protein due to the application of Gaussian polymer chain end-to-end distribution to extract R-g from the FRET efficiency. FRET experiments estimate approximate to 6 angstrom decrease in R-g when the actual decrease is approximate to 3 angstrom on guanidinium chloride denaturant dilution from 7.5 to 1 M, thereby suggesting pronounced compaction in the protein dimensions in the burst phase. The approximate to 3 angstrom decrease is close to the statistical uncertainties of the R-g data measured from SAXS experiments, which suggest no compaction, leading to a disagreement with the FRET experiments. The transition-state ensemble (TSE) structures in Protein L folding are globular and extensive in agreement with the Psi-analysis experiments. The results support the hypothesis that the TSE of single domain proteins depends on protein topology and is not stabilized by local interactions alone.

Insights into the Structural Dynamics of Nucleocytoplasmic Transport of tRNA by Exportin-t

Exportin-t (Xpot) transports mature 5'- and 3'-end processed tRNA from the nucleus to the cytoplasm by associating with a small G-protein Ran (RAs-related nuclear protein), in the nucleus. The release of tRNA in cytoplasm involves RanGTP hydrolysis. Despite the availability of crystal structures of nuclear and cytosolic forms of Xpot, the molecular details regarding the sequential events leading to tRNA release and subsequent conformational changes occurring in Xpot remain unknown. We have performed a combination of classical all-atom and accelerated molecular dynamics simulations on a set of complexes involving Xpot to study a range of features including conformational flexibility of free and cargo-bound Xpot and functionally critical contacts between Xpot and its cargo. The systems investigated include free Xpot and its different complexes, bound either to Ran (GTP/GDP) or tRNA or both. This approach provided a statistically reliable estimate of structural dynamics of Xpot after cargo release. The mechanistic basis for Xpot opening after cargo release has been explained in terms of dynamic structural hinges, about which neighboring region could be displaced to facilitate the nuclear to cytosolic state transition. Post-RanGTP hydrolysis, a cascade of events including local conformational change in RanGTP and loss of critical contacts at Xpot/tRNA interface suggest factors responsible for eventual release of tRNA. The level of flexibility in different Xpot complexes varied depending on the arrangement of individual HEAT repeats. Current study provides one of the most comprehensive and robust analysis carried out on this protein using molecular dynamics schemes.

Self-organized criticality in dynamics without branching

We demonstrate the phenomenon of self-organized criticality (SOC) in a simple random walk model described by a random walk of a myopic ant, i.e., a walker who can see only nearest neighbors. The ant acts on the underlying lattice aiming at uniform digging, i.e., reduction of the height profile of the surface but is unaffected by the underlying lattice. In one, two, and three dimensions we have explored this model and have obtained power laws in the time intervals between consecutive events of "digging." Being a simple random walk, the power laws in space translate to power laws in time. We also study the finite size scaling of asymptotic scale invariant process as well as dynamic scaling in this system. This model differs qualitatively from the cascade models of SOC.

Shock propagation in gas dynamics: explicit form of higher order compatibility conditions

With the use of tensor analysis and the method of singular surfaces, an infinite system of equations can be derived to study the propagation of curved shocks of arbitrary strength in gas dynamics. The first three of these have been explicitly given here. This system is further reduced to one involving scalars only. The choice of dependent variables in the infinite system is quite important, it leads to coefficients free from singularities for all values of the shock strength.

Experiments on a plume with off-source heating: Implications for cloud fluid dynamics

We report here on a series of laboratory experiments on plumes, undertaken with the object of simulating the effect of the heat release that occurs in clouds on condensation of water vapor. The experimental technique used for this purpose relies on ohmic heating generated in an electrically conducting plume fluid subjected to a suitable alternating voltage across specified axial stations in the plume flow [Bhat et al., 1989]. The present series of experiments achieves a value of the Richardson number that is toward the lower end of the range that characteristics cumulus clouds. It is found that the buoyancy enhancement due to heating disrupts the eddy structures in the flow and reduces the dilution owing to entrainment of ambient fluid that would otherwise have occurred in the central region of the plume. Heating also reduces the spread rate of the plume, but as it accelerates the flow as well, the overall specific mass flux in the plume does not show a very significant change at the heat input employed in the experiment. However, there is some indication that the entrainment rate (proportional to the streamwise derivative of the mass flux) is slightly higher immediately after heat injection and slightly lower farther downstream. The measurements support a previous proposal for a cloud scenario [Bhat and Narasimha, 1996] and demonstrate how fresh insights into certain aspects of the fluid dynamics of clouds may be derived from the experimental techniques employed here.

Analysis on conservation of disulphide bonds and their structural features in homologous protein domain families

Background: Disulphide bridges are well known to play key roles in stability, folding and functions of proteins. Introduction or deletion of disulphides by site-directed mutagenesis have produced varying effects on stability and folding depending upon the protein and location of disulphide in the 3-D structure. Given the lack of complete understanding it is worthwhile to learn from an analysis of extent of conservation of disulphides in homologous proteins. We have also addressed the question of what structural interactions replaces a disulphide in a homologue in another homologue. Results: Using a dataset involving 34,752 pairwise comparisons of homologous protein domains corresponding to 300 protein domain families of known 3-D structures, we provide a comprehensive analysis of extent of conservation of disulphide bridges and their structural features. We report that only 54% of all the disulphide bonds compared between the homologous pairs are conserved, even if, a small fraction of the non-conserved disulphides do include cytoplasmic proteins. Also, only about one fourth of the distinct disulphides are conserved in all the members in protein families. We note that while conservation of disulphide is common in many families, disulphide bond mutations are quite prevalent. Interestingly, we note that there is no clear relationship between sequence identity between two homologous proteins and disulphide bond conservation. Our analysis on structural features at the sites where cysteines forming disulphide in one homologue are replaced by non-Cys residues show that the elimination of a disulphide in a homologue need not always result in stabilizing interactions between equivalent residues. Conclusion: We observe that in the homologous proteins, disulphide bonds are conserved only to a modest extent. Very interestingly, we note that extent of conservation of disulphide in homologous proteins is unrelated to the overall sequence identity between homologues. The non-conserved disulphides are often associated with variable structural features that were recruited to be associated with differentiation or specialisation of protein function.

Identification of actin as an estradiol 17-beta-stimulated protein in the human placenta

Addition of estradiol 17-beta to first trimester human placental minces resulted in an increased synthesis of a protein of apparent molecular weight 45 kDa. The specific involvement of estrogen in the stimulation of this protein was established by demonstrating a reduction in the level of this protein by the addition of CCS 16949 A, an inhibitor of aromatase, a key enzyme in the biosynthesis of estradiol 17-beta and ICI 182,780, an estrogen receptor antagonist. The protein was purified to homogeneity and N-terminal sequencing of two of the internal peptides obtained by enzymatic digestion of the protein, as well as the absence of a free N-terminal indicated that it could be actin. This was confirmed by Western blotting using commercially available actin antiserum. The role of estradiol 17-beta in the stimulation of actin synthesis in human placenta was also established by monitoring the quantitative inhibition of DNase I by actin.

Cell Dynamics of Model Proteins

We design rapidly folding sequences by assigning the strongest couplings to the contacts present in a target native state in a two dimensional model of heteropolymers. The pathways to folding and their dependence on the temperature are illustrated via a mapping of the dynamics into motion within the space of the maximally compact cells.

The effect of pH on the unfolding pathway and stability of ribosome-inactivating protein abrin-II

The effect of pH on the unfolding pathway acid the stability of the toxic protein abrin-II have been studied by increasing denaturant concentrations of guanidine hydrochloride and by monitoring the change in 8,1-anilino naphthalene sulfonic acid (ANS) fluorescence upon binding to the hydrophobic sites of the protein. Intrinsic protein fluorescence, far and near UV-circular dichroism (CD) spectroscopy and ANS binding studies reveal that the unfolding of abrin-II occurs through two intermediates at pH 7.2 and one intermediate at pH 4.5. At pH 7.2, the two subunits A and B of abrin-II unfold sequentially. The native protein is more stable at pH 4.5 than at pH 7.2. However, the stability of the abrin-II A-subunit is not affected by a change in pH. These observations may assist in an understanding of the physiologically relevant transmembrane translocation of the toxin.

Thermodynamic characterization of the conformational stability of the homodimeric protein, pea lectin

The conformational stability of the homodimeric pea lectin was determined by both isothermal urea-induced and thermal denaturation in the absence and presence of urea. The denaturation profiles were analyzed to obtain the thermodynamic parameters associated with the unfolding of the protein. The data not only conform to the simple A(2) double left right arrow 2U model of unfolding but also are well described by the linear extrapolation model for the nature of denaturant-protein interactions. In addition, both the conformational stability (Delta G(s)) and the Delta C-p for the protein unfolding is quite high, at about 18.79 kcal/ mol and 5.32 kcal/(mol K), respectively, which may be a reflection of the relatively larger size of the dimeric molecule (M-r 49 000) and, perhaps, a consequent larger buried hydrophobic core in the folded protein. The simple two-state (A(2) double left right arrow 2U) nature of the unfolding process, with the absence of any monomeric intermediate, suggests that the quaternary interactions alone may contribute significantly to the conformational stability of the oligomer-a point that may be general to many oligomeric proteins.

Self-Absorption and Focusing Effects in Carrier Relaxation Dynamics under Depletion Conditions

The effect of a one-dimensional field (1) on the self-absorption characteristics and (2) when we have a finite numerical aperture for the objective lens that focuses the laser beam on the solid are considered here. Self-absorption, in particular its manifestation as an inner filter for the emitted signal, has been observed in luminescence experiments. Models for this effect exist and have been analyzed, but only in the absence of space charge. Using our previous results on minority carrier relaxation in the presence of a field, we obtain expressions incorporating inner filter effects. Focusing of a light beam on the sample, by an objective lens, results in a three-dimensional source and consequently a three-dimensional continuity equation to be solved for the minority carrier concentration. Assuming a one-dimensional electric field and employing Fourier-Bessel transforms, we recast the problem of carrier relaxation and solve the same via an identity that relates it to solutions obtained in the absence of focusing effects. The inner filter effect as well as focusing introduces new time scales in the problem of carrier relaxation. The interplay between the electric field and the parameters which characterize these effects and the consequent modulation of the intensity and time scales of carrier decay signals are analyzed and discussed.

Vortex dynamics at rf frequencies in Bi2Sr2CaCu2O8 single crystals

We have studied the low magnetic field high temperature region of the H-T phase diagram of Bi2Sr2CaCu2O8 single crystals using the technique of non-resonant rf response at a frequency of 20 MHz. With H(rf)parallel to a, H parallel to c, the isothermal magnetic field scans below T-c show that the frequency f(H) of the tank circuit decreases continuously with increase in H before saturating at H similar to H-D(T). Such a decrease in f(H) reflects increasing rf penetration into the weakly screened region between CuO bilayers. The saturation of f(H) at its lowest value for H similar to H-D(T) indicates complete rf penetration land hence the disappearance of field dependence) due to the vanishing of the screening rf currents I-rf(c) in those regions or equivalently when the phase coherence between adjacent superconducting layers vanishes. Therefore H,(T) represents the decoupling of the adjacent superconducting bilayers, and hence also a 3D to 2D decoupling transition of the vortex structure. Simultaneous monitoring of the field dependent rf power dissipation P(H) shows a maximum in dP/dH at H-D(T). The observed H-D(T) line in many crystals is in excellent agreement with the (l/t-1) behavior proposed for decoupling.

Dynamics of a dilute sheared inelastic fluid. I. Hydrodynamic modes and velocity correlation functions

The hydrodynamic modes and the velocity autocorrelation functions for a dilute sheared inelastic fluid are analyzed using an expansion in the parameter epsilon=(1-e)(1/2), where e is the coefficient of restitution. It is shown that the hydrodynamic modes for a sheared inelastic fluid are very different from those for an elastic fluid in the long-wave limit, since energy is not a conserved variable when the wavelength of perturbations is larger than the ``conduction length.'' In an inelastic fluid under shear, there are three coupled modes, the mass and the momenta in the plane of shear, which have a decay rate proportional to k(2/3) in the limit k -> 0, if the wave vector has a component along the flow direction. When the wave vector is aligned along the gradient-vorticity plane, we find that the scaling of the growth rate is similar to that for an elastic fluid. The Fourier transforms of the velocity autocorrelation functions are calculated for a steady shear flow correct to leading order in an expansion in epsilon. The time dependence of the autocorrelation function in the long-time limit is obtained by estimating the integral of the Fourier transform over wave number space. It is found that the autocorrelation functions for the velocity in the flow and gradient directions decay proportional to t(-5/2) in two dimensions and t(-15/4) in three dimensions. In the vorticity direction, the decay of the autocorrelation function is proportional to t(-3) in two dimensions and t(-7/2) in three dimensions.

Dynamics of a dilute sheared inelastic fluid. II. The effect of correlations

The effect of correlations on the viscosity of a dilute sheared inelastic fluid is analyzed using the ring-kinetic equation for the two-particle correlation function. The leading-order contribution to the stress in an expansion in epsilon=(1-e)(1/2) is calculated, and it is shown that the leading-order viscosity is identical to that obtained from the Green-Kubo formula, provided the stress autocorrelation function in a sheared steady state is used in the Green-Kubo formula. A systemmatic extension of this to higher orders is also formulated, and the higher-order contributions to the stress from the ring-kinetic equation are determined in terms of the terms in the Chapman-Enskog solution for the Boltzmann equation. The series is resummed analytically to obtain a renormalized stress equation. The most dominant contributions to the two-particle correlation function are products of the eigenvectors of the conserved hydrodynamic modes of the two correlated particles. In Part I, it was shown that the long-time tails of the velocity autocorrelation function are not present in a sheared fluid. Using those results, we show that correlations do not cause a divergence in the transport coefficients; the viscosity is not divergent in two dimensions, and the Burnett coefficients are not divergent in three dimensions. The equations for three-particle and higher correlations are analyzed diagrammatically. It is found that the contributions due to the three-particle and higher correlation functions to the renormalized viscosity are smaller than those due to the two-particle distribution function in the limit epsilon -> 0. This implies that the most dominant correlation effects are due to the two-particle correlations.