A-Site Residues Move Independently from P-Site Residues in all-Atom Molecular Dynamics Simulations of the 70S Bacterial Ribosome

The ribosome is a large macromolecular machine, and correlated motion between residues is necessary for coordinating function across multiple protein and RNA chains. We ran two all-atom, explicit solvent molecular dynamics simulations of the bacterial ribosome and calculated correlated motion between residue pairs by using mutual information. Because of the short timescales of our simulation (ns), we expect that dynamics are largely local fluctuations around the crystal structure. We hypothesize that residues that show coupled dynamics are functionally related, even on longer timescales. We validate our model by showing that crystallographic B-factors correlate well with the entropy calculated as part of our mutual information calculations. We reveal that A-site residues move relatively independently from P-site residues, effectively insulating A-site functions from P-site functions during translation

Dimensionality of Carbon Nanomaterials Determines the Binding and Dynamics of Amyloidogenic Peptides: Multiscale Theoretical Simulations

Experimental studies have demonstrated that nanoparticles can affect the rate of protein self-assembly, possibly interfering with the development of protein misfolding diseases such as Alzheimer's, Parkinson's and prion disease caused by aggregation and fibril formation of amyloid-prone proteins. We employ classical molecular dynamics simulations and large-scale density functional theory calculations to investigate the effects of nanomaterials on the structure, dynamics and binding of an amyloidogenic peptide apoC-II(60-70). We show that the binding affinity of this peptide to carbonaceous nanomaterials such as C60, nanotubes and graphene decreases with increasing nanoparticle curvature. Strong binding is facilitated by the large contact area available for π-stacking between the aromatic residues of the peptide and the extended surfaces of graphene and the nanotube. The highly curved fullerene surface exhibits reduced efficiency for π-stacking but promotes increased peptide dynamics. We postulate that the increase in conformational dynamics of the amyloid peptide can be unfavorable for the formation of fibril competent structures. In contrast, extended fibril forming peptide conformations are promoted by the nanotube and graphene surfaces which can provide a template for fibril-growth

Model for the Peptide-Free Conformation of Class II MHC Proteins

Background: Major histocompatibility complex proteins are believed to undergo significant conformational changes concomitant with peptide binding, but structural characterization of these changes has remained elusive. Methodology/Principal Findings: Here we use molecular dynamics simulations and experimental probes of protein conformation to investigate the peptide-free state of class II MHC proteins. Upon computational removal of the bound peptide from HLA-DR1-peptide complex, the a50-59 region folded into the P1-P4 region of the peptide binding site, adopting the same conformation as a bound peptide. Strikingly, the structure of the hydrophobic P1 pocket is maintained by engagement of the side chain of Phe a54. In addition, conserved hydrogen bonds observed in crystal structures between the peptide backbone and numerous MHC side chains are maintained between the a51-55 region and the rest of the molecule. The model for the peptide-free conformation was evaluated using conformationally-sensitive antibody and superantigen probes predicted to show no change, moderate change, or dramatic changes in their interaction with peptide-free DR1 and peptide-loaded DR1. The binding observed for these probes is in agreement with the movements predicted by the model. Conclusion/Significance: This work presents a molecular model for peptide-free class II MHC proteins that can help to interpret the conformational changes known to occur within the protein during peptide binding and release, and ca

An All-Atom Model of the Chromatin Fiber Containing Linker Histones Reveals a Versatile Structure Tuned by the Nucleosomal Repeat Length

In the nucleus of eukaryotic cells, histone proteins organize the linear genome into a functional and hierarchical architecture. In this paper, we use the crystal structures of the nucleosome core particle, B-DNA and the globular domain of H5 linker histone to build the first all-atom model of compact chromatin fibers. In this 3D jigsaw puzzle, DNA bending is achieved by solving an inverse kinematics problem. Our model is based on recent electron microscopy measurements of reconstituted fiber dimensions. Strikingly, we find that the chromatin fiber containing linker histones is a polymorphic structure. We show that different fiber conformations are obtained by tuning the linker histone orientation at the nucleosomes entry/exit according to the nucleosomal repeat length. We propose that the observed in vivo quantization of nucleosomal repeat length could reflect nature's ability to use the DNA molecule's helical geometry in order to give chromatin versatile topological and mechanical properties

Effectiveness and cost-effectiveness of 'BeweegKuur', a combined lifestyle intervention in the Netherlands: Rationale, design and methods of a randomized controlled trial

<p>Abstract</p> <p>Background</p> <p>Improving the lifestyle of overweight and obese adults is of increasing interest in view of its role in several chronic diseases. Interventions aiming at overweight or weight-related chronic diseases suffer from high drop-out rates. It has been suggested that Motivational Interviewing and more frequent and more patient-specific coaching could decrease the drop-out rate. 'BeweegKuur' is a multidisciplinary lifestyle intervention which offers three programmes for overweight persons. The effectiveness and the cost-effectiveness of intensively guided programmes, such as the 'supervised exercise programme' of 'BeweegKuur', for patients with high weight-related health risk, remain to be assessed. Our randomized controlled trial compares the expenses and effects of the 'supervised exercise programme' with those of the less intensively supervised 'start-up exercise programme'.</p> <p>Methods/Design</p> <p>The one-year intervention period involves coaching by a lifestyle advisor, a physiotherapist and a dietician, coordinated by general practitioners (GPs). The participating GP practices have been allocated to the interventions, which differ only in terms of the amount of coaching offered by the physiotherapist. Whereas the 'start-up exercise programme' includes several consultations with physiotherapists to identify barriers hampering independent exercising, the 'supervised exercise programme' includes more sessions with a physiotherapist, involving exercise under supervision. The main goal is transfer to local exercise facilities. The main outcome of the study will be the participants' physical activity at the end of the one-year intervention period and after one year of follow-up. Secondary outcomes are dietary habits, health risk, physical fitness and functional capacity. The economic evaluation will consist of a cost-effectiveness analysis and a cost-utility analysis. The primary outcome measures for the economic evaluation will be the physical activity and the number of quality-adjusted life years. Costs will be assessed from a societal perspective with a time horizon of two years. Additionally, a process evaluation will be used to evaluate the performance of the intervention and the participants' evaluation of the intervention.</p> <p>Discussion</p> <p>This study is expected to provide information regarding the additional costs and effects of the 'supervised exercise programme' in adults with very high weight-related health risk.</p> <p>Trial registration number</p> <p>ISRCTN: <a href="http://www.controlled-trials.com/ISRCTN46574304">ISRCTN46574304</a></p

The role of interfacial lipids in stabilizing membrane protein oligomers

Oligomerization of membrane proteins in response to lipid binding has a critical role in many cell-signalling pathways1 but is often difficult to define2 or predict3. Here we report the development of a mass spectrometry platform to determine simultaneously the presence of interfacial lipids and oligomeric stability and to uncover how lipids act as key regulators of membrane-protein association. Evaluation of oligomeric strength for a dataset of 125 α-helical oligomeric membrane proteins reveals an absence of interfacial lipids in the mass spectra of 12 membrane proteins with high oligomeric stability. For the bacterial homologue of the eukaryotic biogenic transporters (LeuT4, one of the proteins with the lowest oligomeric stability), we found a precise cohort of lipids within the dimer interface. Delipidation, mutation of lipid-binding sites or expression in cardiolipin-deficient Escherichia coli abrogated dimer formation. Molecular dynamics simulation revealed that cardiolipin acts as a bidentate ligand, bridging across subunits. Subsequently, we show that for the Vibrio splendidus sugar transporter SemiSWEET5, another protein with low oligomeric stability, cardiolipin shifts the equilibrium from monomer to functional dimer. We hypothesized that lipids are essential for dimerization of the Na+/H+ antiporter NhaA from E. coli, which has the lowest oligomeric strength, but not for the substantially more stable homologous Thermus thermophilus protein NapA. We found that lipid binding is obligatory for dimerization of NhaA, whereas NapA has adapted to form an interface that is stable without lipids. Overall, by correlating interfacial strength with the presence of interfacial lipids, we provide a rationale for understanding the role of lipids in both transient and stable interactions within a range of α-helical membrane proteins, including G-protein-coupled receptors

Structural Insights from Binding Poses of CCR2 and CCR5 with Clinically Important Antagonists: A Combined In Silico Study

Chemokine receptors are G protein-coupled receptors that contain seven transmembrane domains. In particular, CCR2 and CCR5 and their ligands have been implicated in the pathophysiology of a number of diseases, including rheumatoid arthritis and multiple sclerosis. Based on their roles in disease, they have been attractive targets for the pharmaceutical industry, and furthermore, targeting both CCR2 and CCR5 can be a useful strategy. Owing to the importance of these receptors, information regarding the binding site is of prime importance. Structural studies have been hampered due to the lack of X-ray crystal structures, and templates with close homologs for comparative modeling. Most of the previous models were based on the bovine rhodopsin and β2-adrenergic receptor. In this study, based on a closer homolog with higher resolution (CXCR4, PDB code: 3ODU 2.5 Å), we constructed three-dimensional models. The main aim of this study was to provide relevant information on binding sites of these receptors. Molecular dynamics simulation was done to refine the homology models and PROCHECK results indicated that the models were reasonable. Here, binding poses were checked with some established inhibitors of high pharmaceutical importance against the modeled receptors. Analysis of interaction modes gave an integrated interpretation with detailed structural information. The binding poses confirmed that the acidic residues Glu291 (CCR2) and Glu283 (CCR5) are important, and we also found some additional residues. Comparisons of binding sites of CCR2/CCR5 were done sequentially and also by docking a potent dual antagonist. Our results can be a starting point for further structure-based drug design

Snake Cytotoxins Bind to Membranes via Interactions with Phosphatidylserine Head Groups of Lipids

The major representatives of Elapidae snake venom, cytotoxins (CTs), share similar three-fingered fold and exert diverse range of biological activities against various cell types. CT-induced cell death starts from the membrane recognition process, whose molecular details remain unclear. It is known, however, that the presence of anionic lipids in cell membranes is one of the important factors determining CT-membrane binding. In this work, we therefore investigated specific interactions between one of the most abundant of such lipids, phosphatidylserine (PS), and CT 4 of Naja kaouthia using a combined, experimental and modeling, approach. It was shown that incorporation of PS into zwitterionic liposomes greatly increased the membrane-damaging activity of CT 4 measured by the release of the liposome-entrapped calcein fluorescent dye. The CT-induced leakage rate depends on the PS concentration with a maximum at approximately 20% PS. Interestingly, the effects observed for PS were much more pronounced than those measured for another anionic lipid, sulfatide. To delineate the potential PS binding sites on CT 4 and estimate their relative affinities, a series of computer simulations was performed for the systems containing the head group of PS and different spatial models of CT 4 in aqueous solution and in an implicit membrane. This was done using an original hybrid computational protocol implementing docking, Monte Carlo and molecular dynamics simulations. As a result, at least three putative PS-binding sites with different affinities to PS molecule were delineated. Being located in different parts of the CT molecule, these anion-binding sites can potentially facilitate and modulate the multi-step process of the toxin insertion into lipid bilayers. This feature together with the diverse binding affinities of the sites to a wide variety of anionic targets on the membrane surface appears to be functionally meaningful and may adjust CT action against different types of cells

Bak Conformational Changes Induced by Ligand Binding: Insight into BH3 Domain Binding and Bak Homo-Oligomerization

Recently we reported that the BH3-only proteins Bim and Noxa bind tightly but transiently to the BH3-binding groove of Bak to initiate Bak homo-oligomerization. However, it is unclear how such tight binding can induce Bak homo-oligomerization. Here we report the ligand-induced Bak conformational changes observed in 3D models of Noxa·Bak and Bim·Bak refined by molecular dynamics simulations. In particular, upon binding to the BH3-binding groove, Bim and Noxa induce a large conformational change of the loop between helices 1 and 2 and in turn partially expose a remote groove between helices 1 and 6 in Bak. These observations, coupled with the reported experimental data, suggest formation of a pore-forming Bak octamer, in which the BH3-binding groove is at the interface on one side of each monomer and the groove between helices 1 and 6 is at the interface on the opposite side, initiated by ligand binding to the BH3-binding groove

Insights on Glucocorticoid Receptor Activity Modulation through the Binding of Rigid Steroids

Background: The glucocorticoid receptor (GR) is a transcription factor that regulates gene expression in a ligand-dependent fashion. This modular protein is one of the major pharmacological targets due to its involvement in both cause and treatment of many human diseases. Intense efforts have been made to get information about the molecular basis of GR activity. Methodology/Principal Findings: Here, the behavior of four GR-ligand complexes with different glucocorticoid and antiglucocorticoid properties were evaluated. The ability of GR-ligand complexes to oligomerize in vivo was analyzed by performing the novel Number and Brightness assay. Results showed that most of GR molecules form homodimers inside the nucleus upon ligand binding. Additionally, in vitro GR-DNA binding analyses suggest that ligand structure modulates GRDNA interaction dynamics rather than the receptor's ability to bind DNA. On the other hand, by coimmunoprecipitation studies we evaluated the in vivo interaction between the transcriptional intermediary factor 2 (TIF2) coactivator and different GR-ligand complexes. No correlation was found between GR intranuclear distribution, cofactor recruitment and the homodimerization process. Finally, Molecular determinants that support the observed experimental GR LBD-ligand/TIF2 interaction were found by Molecular Dynamics simulation. Conclusions/Significance: The data presented here sustain the idea that in vivo GR homodimerization inside the nucleus can be achieved in a DNA-independent fashion, without ruling out a dependent pathway as well. Moreover, since at least one GR-ligand complex is able to induce homodimer formation while preventing TIF2 coactivator interaction, results suggest that these two events might be independent from each other. Finally, 21-hydroxy-6,19-epoxyprogesterone arises as a selective glucocorticoid with potential pharmacological interest. Taking into account that GR homodimerization and cofactor recruitment are considered essential steps in the receptor activation pathway, results presented here contribute to understand how specific ligands influence GR behavior. © 2010 Presman et al.Fil:Presman, D.M. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales; Argentina.Fil:Alvarez, L.D. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales; Argentina.Fil:Levi, V. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales; Argentina.Fil:Martí, M.A. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales; Argentina.Fil:Veleiro, A.S. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales; Argentina.Fil:Burton, G. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales; Argentina.Fil:Pecci, A. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales; Argentina