Spectral properties of conformational motion in proteins

The conformational dynamics of large fragments of protein structure within the framework of a generalized model of bounded diffusion (GMBD) have been considered. This model is a development of the known ordinary model of bounded diffusion (MBD). The latter model assumes that the driving force of motion, which is at the same time a source of friction, has negligibly small correlation time (white noise). In contrast to the MBD, the GMBD takes into account the finite character of friction correlation time, i.e. memory friction effects. Two different mechanisms of friction for the fragment are considered: (1) friction by rapid density fluctuations due to the harmonic vibrational motion of atoms and (2) friction by slow encounters with surrounding fragments of protein structure. The present theory shows that at high frequencies the power decay of the spectrum which takes place in the MBD is replaced by an exponential one. In the narrow intermediate frequency range the results of the GMBD coincide with those of the MBD. At low frequencies the results of the two models differ quantitatively, though their qualitative behaviour is similar.     

Human thrombospondin's (TSP-1) C-terminal domain opens to interact with the CD-47 receptor: a molecular modeling study

Thrombospondin-1 (TSP-1) interaction with the membranous receptor CD-47 involves the peptide RFYVVMWK (4N-1) located in its C-terminal domain. However, the available X-ray structure of TSP-1 describes this peptide as completely buried into a hydrophobic pocket, preventing any interaction. Where classical standard methods failed, an appropriate approach combining normal mode analysis and an adapted protocol of energy minimization identified the large amplitude motions responsible of the partial solvent exposure of 4N-1. In agreement, the obtained model of the open TSP-1 was further used for protein-protein docking experiments against a homology model generated for CD-47. Considering the multiple applications of the CD-47 receptor as a target, our results open new pharmacological perspectives for the design of TSP-1:CD-47 inhibitors and CD-47 antagonists. We also suggest a common opening mechanism for proteins sharing the same fold as TSP-1. This work also suggests the usefulness of our approach in other topics in which predictions of protein-protein interactions are of importance.     

Intrinsic ligand binding properties of the human and bovine α,-interferon receptors

The Type I interferon receptor (IFN-αR) interacts with all IFN-αs, IFN-β and IFN-ω, and seems to be a multisubunit receptor. To investigate the role of a cloned receptor subunit (IFN-αR1), we have examined the intrinsic ligand binding properties of the bovine and human IFN-αR1 polypeptides expressed in Xenopus laevis oocytes. Albeit with different efficiencies, Xenopus oocytes expressing either the human or bovine IFN-αR1 polypeptide exhibit significant binding and formation of crosslinked complexes with human IFN-αA and IFN-αB. Thus, the IFN-αR1 polypeptide most likely plays a direct role in ligand binding.     

The impact of adding trunk motion to the interpretation of the role of joint moments during normal walking

Biomechanical model assumptions affect the interpretation of the role of the muscle or joint moments to the segmental power estimated by induced acceleration analysis (IAA). We evaluated the effect of modeling the pelvis and trunk segments as two separate segments (8 SM) versus as a single segment (7 SM) on the segmental power, support of the body, knee and hip extension acceleration produced by the joint moments during the stance phase of normal walking. Significant differences were observed in the contribution of the stance hip abductor and extensor moments to support, ipsilateral knee and hip acceleration, and ipsilateral thigh and upper body power. The primary finding was that the role of the stance hip moment in generating ipsilateral thigh and upper body power differed based on degrees of freedom in the model. Secondarily, the magnitude of contributions also differed. For example, the hip abductor and extensor moments showed greater contribution to support, hip and knee acceleration in the 8 SM. IAA and segment power analysis are sensitive to the degrees of freedom between the pelvis and trunk. There is currently no gold standard by which to evaluate the accuracy of IAA predictions. However, modeling the pelvis and trunk as separate segments is closer to the anatomical architecture of the body. An 8 SM appears to be more appropriate for estimating the role of joint moments, particularly to motion of more proximal segments during normal walking.     

大鼠运动病敏感性性别差异与相关脑区AVP及其V1b受体表达的关系

目的：检测不同性别大鼠旋转刺激后脑内相关区域精氨酸加压素（AVP）含量及V1b受体表达的变化,探讨AVP及受体参与运动病的可能机制。方法：给予SD大鼠30 min绕水平轴的旋转刺激,然后采用放免法检测相关脑区AVP含量,并通过荧光免疫组化方法测定相应脑区V1b受体的表达情况。结果：①在雌性大鼠,旋转刺激组各脑区AVP含量无显著性改变;对于雄性大鼠,对照组各检测脑区AVP含量高于雌性,旋转刺激组小脑、延髓内AVP含量的变化无显著性意义,但前脑、间脑、脑桥内AVP含量较对照组明显降低（P〈0.05）。②雌性大鼠视上核AVP的V1b受体表达阳性神经元数量旋转刺激组显著低于对照组（P〈0.05）,而前庭核、最后区V1b受体表达阳性神经元数量明显多于对照组（P〈0.05）;在雄性大鼠,旋转刺激组视上核与前庭核V1b受体表达阳性神经元数量无显著性改变,而最后区V1b受体表达阳性神经元数量有所增加（P〈0.05）,但增加幅度没有雌性大鼠明显。结论：前脑、间脑、脑桥内AVP含量与前庭核和最后区V1b受体表达及对旋转刺激反应的差异可能与运动病敏感性性别差异有关,并且前庭核、最后区可能是AVP-V1受体拮抗剂抗运动病作用的靶点。     

Spectral broadening of the Soret band in myoglobin: an interpretation by the full spectrum of low-frequency modes from a normal modes analysis

In this work the temperature dependence of the Soret band line shape in carbon-monoxy myoglobin is re-analyzed by using both the full correlator approach in the time domain and the frequency domain approach. The new analyses exploit the full density of vibrational states of carbon-monoxy myoglobin available from normal modes analysis, and avoid the artificial division of the entire set of vibrational modes coupled to the Soret transition into "high-frequency" and "low-frequency" subsets; the frequency domain analysis, however, makes use of the so-called short-times approximation, while the time domain one avoids it. Time domain and frequency domain analyses give very similar results, thus supporting the applicability of the short-times approximation to the analysis of hemeprotein spectra; in particular, they clearly indicate the presence of spectral heterogeneity in the Soret band of carbon-monoxy myoglobin. The analyses also show that a temperature dependence of the Gaussian width parameter steeper than the hyperbolic cotangent law predicted by the Einstein harmonic oscillator and/or a temperature dependence of inhomogeneous broadening are not sufficient to obtain quantitative information on the magnitude of an-harmonic contributions to the iron-heme plane motion. However, the dependence of the previous two quantities may be used to obtain semiquantitative information on the overall coupling of the Soret transition to the low-frequency modes and therefore on the dynamic properties of the heme pocket in different states of the protein.     

Intrinsic properties and regulation of Pannexin 1 channel

Pannexin 1 (Panx1) channels are generally represented as non-selective, large-pore channels that release ATP. Emerging roles have been described for Panx1 in mediating purinergic signaling in the normal nervous, cardiovascular, and immune systems, where they may be activated by mechanical stress, ionotropic and metabotropic receptor signaling, and via proteolytic cleavage of the Panx1 C-terminus. Panx1 channels are widely expressed in various cell types, and it is now thought that targeting these channels therapeutically may be beneficial in a number of pathophysiological contexts, such as asthma, atherosclerosis, hypertension, and ischemic-induced seizures. Even as interest in Panx1 channels is burgeoning, some of their basic properties, mechanisms of modulation, and proposed functions remain controversial, with recent reports challenging some long-held views regarding Panx1 channels. In this brief review, we summarize some well-established features of Panx1 channels; we then address some current confounding issues surrounding Panx1 channels, especially with respect to intrinsic channel properties, in order to raise awareness of these unsettled issues for future research.     

Maspin, the molecular bridge between the plasminogen activator system and beta1 integrin that facilitates cell adhesion

Maspin is a non-inhibitory serine protease inhibitor (serpin) that influences many cellular functions including adhesion, migration, and invasion. The underlying molecular mechanisms that facilitate these actions are still being elucidated. In this study we determined the mechanism by which maspin mediates increased MCF10A cell adhesion. Utilizing competition peptides and mutation analyses, we discovered two unique regions (amino acid residues 190-202 and 260-275) involved in facilitating the increased adhesion function of maspin. In addition, we demonstrate that the urokinase-type plasminogen activator (uPA)/uPA receptor (uPAR) complex is required for the localization and adhesion function of maspin. Finally, we showed that maspin, uPAR, and β1 integrin co-immunoprecipitate, suggesting a novel maspin-uPA-uPAR-β1 integrin mega-complex that regulates mammary epithelial cell adhesion.     

Influence of the secondary structure on the pore forming properties of synthetic alamethicin analogs: NMR and molecular modelling studies

Synthetic alamethicin analogs, in which all Aib residues had been replaced by Leu (L2) then proline 14 replaced by an alanine (L5), were studied in SDS micelles using circular dichroism and NMR spectroscopy. Nuclear Overhauser effects were used as constraints for molecular modelling. The structures determined for both peptides in SDS micelles were compared with those previously obtained in methanol in order to establish a secondary structure/ionophore activity relationship. Our results indicated that a shortening of peptide helices could be responsible for the observed decrease in ion channel lifetimes. However, the length of helices may not by itself explain the drastic destabilization of channels when Pro14 of alamethicin is replaced by Ala in L5. Indeed analysis of the helical wheel of L5 reveals heterogeneity in the amphipathicity depending on the medium. Thus, loss of amphipathicity seems to underly the observed destabilization of channels. © 1998 European Peptide Society and John Wiley & Sons, Ltd.     

Structure and hydration of the amylopectin trisaccharide building blocks--Synthesis, NMR, and molecular dynamics

To gain insight into the molecular details and hydration of amylopectin, the five constituting trisaccharides have been chemically synthesized as their methyl alpha-glycosides. All five trisaccharides were subjected to 950 MHz NMR spectroscopy for complete assignment and nanosecond molecular dynamics trajectories were calculated to study the structure and dynamics of the trisaccharides in aqueous solution. Systematic analysis of the simulation data revealed several examples of bridging water molecules playing an important role in the stabilization of specific amylopectin conformations, which was also supported by the experimental NMR data such as interresidue NOE's and heteronuclear scalar couplings between nuclei from neighboring residues. Although alpha-maltotriose, alpha-iso-maltotriose, alpha-panose and alpha-isopanose are relatively well characterized structures, the study also includes one less characterized trisaccharide with the structure alphaGlcp(1-->4)alphaGlcp(1-->6)alphaGlcp. This trisaccharide, tentatively labelled alpha-forkose, is located at the branch point of amylopectin, forking the amylopectin into two strands that align into double-helical segments. The results show that the conformation of alpha-forkose takes a natural bend form which fits well into the structure of the double-helical segment of amylopectin. As the only trisaccharide in this study the structure of alpha-forkose is not significantly influenced by the hydration. In contrast, alpha-isopanose takes a restricted, but rather extended form due to an exceptionally strong localized water density. The two homo-linkage oligomers, alpha-maltotriose and alpha-iso-maltotriose, showed to be the most extended and the most flexible trimers, respectively, providing regular structure for crystalline domains and maximum linker flexibility for amorphous domains.     

Bias-free separation of internal and overall motion of biomolecules

Collective internal motions are known to be important for the function of biological macromolecules. It has been discussed in the past whether the application of superimposing algorithms to remove the overall motion from a structural ensemble introduces artificial correlations between distant atoms. Here we present a new method to eliminate residual rotation and translation from cartesian modes derived from a normal mode analysis or from a principal component analysis. Bias-free separation is based on the idea that the addition of modes of pure rotation/translation can compensate the residual overall motion. Removal of overall motion must reduce the "total amount of motion" (TAM) in the mode. Our algorithm allows to back-calculate revised covariance matrices. The approach was applied to two model systems that show residual overall motion, when analyzed using all atoms as reference for the superimposing algorithm. In both cases, our algorithm was capable of eliminating residual covariances caused by the overall motion, while retaining internal covariances even for very distant atoms. A structural ensemble obtained for a 13-ns molecular dynamics simulation of the protein Ribonuclease T1 showed a covariance matrix of the corrected modes with significantly sharper contours after applying the bias-free separation.     

In silico approaches to evaluate the molecular properties of organophosphate compounds to inhibit acetylcholinesterase activity in housefly

Organophosphate compounds (OPC) have become the primary choice as insecticides and are widely used across the world. Additionally, OPCs were also commonly used as a chemical warfare agent that triggers a great challenge to public safety. Exposure of OPCs to human causes immediate excitation of cholinergic neurotransmission through transient elevation of synaptic acetylcholine (ACh) levels and accumulations. Likewise, prolonged exposure of OPCs can affect the processes in immune response, carbohydrate metabolism, cardiovascular toxicity, and several others. Studies revealed that the toxicity of OPCs was provoked by inhibition of acetylcholinesterase (AChE). Therefore, combined in silico approaches – pharmacophore-based 3D-QSAR model; docking and Molecular Dynamics (MD) – were used to assess the precise and comprehensive effects of series of known OP-derived compounds together with its ?log LD₅₀ values. The selected five-featured pharmacophore model – AAHHR.61 – displayed the highest correlation (R² = .9166), cross-validated coefficient (Q² = .8221), F = 63.2, Pearson-R = .9615 with low RMSE = .2621 values obtained using five component PLS factors. Subsequently, the well-validated model was then used as a 3D query to search novel OPCs using a high-throughput virtual screening technique. Simultaneously, the docking studies predicted the binding pose of the most active OPC in the MdAChE binding pocket. Additionally, the stability of docking was verified using MD simulation. The results revealed that OP22 and predicted lead compounds bound tightly to S315 of MdAChE through potential hydrogen bond interaction over time. Overall, this study might provide valuable insight into binding mode of OPCs and hit compounds to inhibit AChE in housefly.     

Alpha-crystallin and ATP facilitate the in vitro renaturation of xylanase: enhancement of refolding by metal ions

Alpha-crystallin is a multimeric protein that functions as a molecular chaperone and shares extensive structural homology to small heat shock proteins. For the functional in vitro analysis of alpha-crystallin, the xylanase Xyl II from alkalophilic thermophilic Bacillus was used as a model system. The mechanism of chaperone action of alpha-crystallin is less investigated. Here we studied the refolding of Gdn HCl-denatured Xyl II in the presence and absence of alpha-crystallin to elucidate the molecular mechanism of chaperone-mediated in vitro folding. Our results, based on intrinsic tryptophan fluorescence and hydrophobic fluorophore 8-anilino-1-naphthalene sulfonate binding studies, suggest that alpha-crystallin formed a complex with a putative molten globule-like intermediate in the refolding pathway of Xyl II. The alpha-crystallin.Xyl II complex exhibited no functional activity. Addition of ATP to the complex initiated the renaturation of Xyl II with 30%-35% recovery of activity. The nonhydrolyzable analog 5'-adenylyl imidodiphosphate (AMP-PNP) was capable of reconstitution of active Xyl II to a lesser extent than ATP. Although the presence of Ca(2+) was not required for the in vitro refolding of Xyl II, the renaturation yield was enhanced in its presence. Experimental evidence indicated that the binding of ATP to the alpha-crystallin.Xyl II complex brought about conformational changes in alpha-crystallin facilitating the dissociation of xylanase molecules. This is the first report of the enhancement of alpha-crystallin chaperone functions by metal ions.     

Study of the motion of magnetotactic bacteria

Motion of flagellate bacteria is considered from the point of view of rigid body mechanics. As a general case we consider a flagellate coccus magnetotactic bacterium swimming in a fluid in the presence of an external magnetic field. The proposed model generalizes previous approaches to the problem and allows one to access parameters of the motion that can be measured experimentally. The results suggest that the strong helical pattern observed in typical trajectories of magnetotactic bacteria can be a biological advantage complementary to magnetic orientation. In the particular case of zero magnetic interaction the model describes the motion of a non-magnetotactic coccus bacterium swimming in a fluid. Theoretical calculations based on experimental results are compared with the experimental track obtained by dark field optical microscopy. Correspondence to: H. G. P. Lins de Barros     

Mechanisms for allosteric activation of protease DegS by ligand binding and oligomerization as revealed from molecular dynamics simulations

A local perturbation of a protein may lead to functional changes at some distal site, a phenomenon denoted as allostery. Here, we study the allosteric control of a protease using molecular dynamics simulations. The system considered is the bacterial protein DegS which includes a protease domain activated on ligand binding to an adjacent PDZ domain. Starting from crystallographic structures of DegS homo‐trimers, we perform simulations of the ligand‐free and ‐bound state of DegS at equilibrium. Considering a single protomer only, the trimeric state was mimicked by applying restraints on the residues in contact with other protomers in the DegS trimer. In addition, the bound state was also simulated without any restraints to mimic the monomer. Our results suggest that not only ligand release but also disassembly of a DegS trimer inhibits proteolytic activity. Considering various observables for structural changes, we infer allosteric pathways from the interface with other protomers to the active site. Moreover, we study how ligand release leads to (i) catalytically relevant changes involving residues 199–201 and (ii) a transition from a stretched to a bent conformation for residues 217–219 (which prohibits proper substrate binding). Finally, based on ligand‐induced C_α shifts we identify residues in contact with other protomers in the DegS trimer that likely transduce the perturbation from ligand release from a given protomer to adjacent protomers. These residues likely play a key role in the experimentally known effect of ligand release from a protomer on the proteolytic activity of the other protomers. Proteins 2016; 84:1690–1705. © 2016 Wiley Periodicals, Inc.     

Intrinsic protein-protein interaction-mediated and chaperonin-assisted sequential assembly of stable bardet-biedl syndrome protein complex, the BBSome

The pleiotropic features of obesity, retinal degeneration, polydactyly, kidney abnormalities, cognitive impairment, hypertension, and diabetes found in Bardet-Biedl syndrome (BBS) make this disorder an important model disorder for identifying molecular mechanisms involved in common human diseases. To date, 16 BBS genes have been reported, seven of which (BBS1, 2, 4, 5, 7, 8, and 9) code for proteins that form a complex known as the BBSome. The function of the BBSome involves ciliary membrane biogenesis. Three additional BBS genes (BBS6, BBS10, and BBS12) have homology to type II chaperonins and interact with CCT/TRiC proteins and BBS7 to form a complex termed the BBS-chaperonin complex. This complex is required for BBSome assembly. Little is known about the process and the regulation of BBSome formation. We utilized point mutations and null alleles of BBS proteins to disrupt assembly of the BBSome leading to the accumulation of BBSome assembly intermediates. By characterizing BBSome assembly intermediates, we show that the BBS-chaperonin complex plays a role in BBS7 stability. BBS7 interacts with BBS2 and becomes part of a BBS7-BBS2-BBS9 assembly intermediate referred to as the BBSome core complex because it forms the core of the BBSome. BBS1, BBS5, BBS8, and finally BBS4 are added to the BBSome core to form the complete BBSome.     

Plasmin-sensitive dibasic sequences in the third fibronectin-like domain of L1-cell adhesion molecule (CAM) facilitate homomultimerization and concomitant integrin recruitment

L1 is a multidomain transmembrane neural recognition molecule essential for neurohistogenesis. While moieties in the immunoglobulin-like domains of L1 have been implicated in both heterophilic and homophilic binding, the function of the fibronectin (FN)-like repeats remains largely unresolved. Here, we demonstrate that the third FN-like repeat of L1 (FN3) spontaneously homomultimerizes to form trimeric and higher order complexes. Remarkably, these complexes support direct RGD-independent interactions with several integrins, including alpha(v)beta(3) and alpha(5)beta(1). A pep- tide derived from the putative C-C' loop of FN3 (GSQRKHSKRHIHKDHV(852)) also forms trimeric complexes and supports alpha(v)beta(3) and alpha(5)beta(1) binding. Substitution of the dibasic RK(841) and KR(845) sequences within this peptide or the FN3 domain limited multimerization and abrogated integrin binding. Evidence is presented that the multimerization of, and integrin binding to, the FN3 domain is regulated both by conformational constraints imposed by other domains and by plasmin- mediated cleavage within the sequence RK( downward arrow)HSK( downward arrow)RH(846). The integrin alpha(9)beta(1), which also recognizes the FN3 domain, colocalizes with L1 in a manner restricted to sites of cell-cell contact. We propose that distal receptor ligation events at the cell-cell interface may induce a conformational change within the L1 ectodomain that culminates in receptor multimerization and integrin recruitment via interaction with the FN3 domain.     

A Lagrange-based generalised formulation for the equations of motion of simple walking models

Simple 2D models of walking often approximate the human body to multi-link dynamic systems, where body segments are represented by rigid links connected by frictionless hinge joints. Performing forward dynamics on the equations of motion (EOM) of these systems can be used to simulate their movement. However, deriving these equations can be time consuming. Using Lagrangian mechanics, a generalised formulation for the EOM of n-link open-loop chains is derived. This can be used for single support walking models. This has an advantage over Newton–Euler mechanics in that it is independent of coordinate system and prior knowledge of the ground reaction force (GRF) is not required. Alternative strategies, such as optimisation algorithms, can be used to estimate joint activation and simulate motion. The application of Lagrange multipliers, to enforce motion constraints, is used to adapt this general formulation for application to closed-loop chains. This can be used for double support walking models. Finally, inverse dynamics are used to calculate the GRF for these general n-link chains. The necessary constraint forces to maintain a closed-loop chain, calculated from the Lagrange multipliers, are one solution to the indeterminate problem of GRF distribution in double support models. An example of this method’s application is given, whereby an optimiser estimates the joint moments by tracking kinematic data.