Prediction of Peptide Binding to Major Histocompatibility II Receptors with Molecular Mechanics and Semi-Empirical Quantum Mechanics Methods

Methods for prediction of the binding of peptides to major histocompatibility complex (MHC) II receptors are examined, using literature values of IC50 as a benchmark. Two sets of IC50 data for closely structurally related peptides based on hen egg lysozyme (HEL) and myelin basic protein (MBP) are reported first. This shows that methods based on both molecular mechanics and semi-empirical quantum mechanics can predict binding with good-to-reasonable accuracy, as long as a suitable method for estimation of solvation effects is included. A more diverse set of 22 peptides bound to HLA-DR1 provides a tougher test of such methods, especially since no crystal structure is available for these peptide-MHC complexes. We therefore use sequence based methods such as SYFPEITHI and SVMHC to generate possible binding poses, using a consensus approach to determine the most likely anchor residues, which are then mapped onto the crystal structure of an unrelated peptide bound to the same receptor. This analysis shows that the MM/GBVI method performs particularly well, as does the AMBER94 forcefield with Born solvation model. Indeed, MM/GBVI can be used as an alternative to sequence based methods in generating binding poses, leading to still better accuracy

Giardia: a pathogen or commensal for children in high-prevalence settings?

Giardia is a common intestinal parasite worldwide, and infection can be associated with clear and sometimes persistent symptomatology. However, in children in high prevalence settings, it is not associated with or is perhaps even protective against acute diarrhea, and the association with long-term outcomes has been difficult to discern

Differences in both prevalence and titre of specific immunoglobulin E among children with asthma in affluent and poor communities within a large town in Ghana.

Background Reports from several African countries have noted an increasing prevalence of asthma in areas of extensive urbanization. Objective To investigate the relevance of allergen-specific sensitization and body mass index (BMI) to asthma/wheezing and exercise-induced bronchospasm (EIB) among children from affluent and poorer communities within a large town in Ghana. Methods Children with physician-diagnosed asthma and/or current wheezing aged 9-16 years (n=99; cases) from three schools with differing socio-economic backgrounds [urban affluent (UA), urban poor (UP) or suburban/rural (SR)] were recruited from a cross-sectional study (n=1848) in Kumasi, Ghana, and matched according to age, sex and area of residence with non-asthmatic/non-wheezy controls. We assayed sera for IgE antibodies to mite, cat, dog, cockroach, Ascaris and galactose-α-1,3-galactose. Results Children from the UA school had the lowest total serum IgE. However, cases from the UA school had a higher prevalence and mean titre of sIgE to mite (71.4%, 21.2IU/mL) when compared with controls (14.3%, 0.8IU/mL) or cases from UP (30%, 0.8IU/mL) and SR community (47.8%, 1.6IU/mL). While similar findings were observed with EIB in the whole population, among cases there was no difference in IgE antibody prevalence or titre between children with or without EIB. BMI was higher among UA children with and without asthma; in UP and SR communities, children with EIB (n=14) had a significantly higher BMI compared with children with asthma/wheezing without EIB (n=38) (18.2 vs. 16.4, respectively, P<0.01). Conclusions and Clinical Relevance In the relatively affluent school, asthma/wheezing and EIB were associated with high titre IgE antibodies to mite, decreased total IgE, and increased BMI. This contrasted with children in the urban poor school and suggests that changes relevant to a Western model of childhood asthma can occur within a short geographical distance within a large city in Africa. © 2011 Blackwell Publishing Ltd

Quantum chemical molecular dynamics and metadynamics simulation of aluminium binding to amyloid-β and related peptides

We report semi-empirical tight-binding simulations of the interaction between Al(III) and biologically relevant peptides. The GFN2-XTB method is shown to accurately reproduce previously reported and density functional theory (DFT)-calculated geometries of model systems. Molecular dynamics simulations based on this method are able to sample peptide flexibility over timescales of up to nanoseconds, but these timescales are insufficient to explore potential changes in metal–peptide binding modes. To achieve this, metadynamics simulations using root mean square deviation as a collective variable were employed. With suitably chosen biasing potentials, these are able to efficiently explore diverse coordination modes, for instance, through Glu and/or Asp residues in a model peptide. Using these methods, we find that Al(III) binding to the N-terminal sequence of amyloid-β is highly fluxional, with all acidic sidechains and several backbone oxygens participating in coordination. We also show that such simulations could provide a means to predict a priori possible binding modes as a precursor to longer, atomistic simulations

Stereoelectronic control of photophysics:red and yellow axial and equatorial anomers of a rhenium-quinoline complex

A novel quinoline-substituted pyrimidine ligand forms two different coloured complexes upon reaction with Re(CO)(5)Br. These compounds display distinct photophysical properties that are dictated by their stereochemistry

Molecular dynamics simulations of copper binding to N-terminus mutants of amyloid-β

We report results of molecular dynamic (MD) simulations on N-terminus mutants of the copper-bound, amyloid-β (Aβ) peptide. Eight structures of Aβ were modelled, including seven mutant peptides in addition to the unaltered wild-type (WT). Trajectories analysed for each individual system were all approximately 1.4 μs in length, yielding a total of over 11 μs in total. The impact of these mutations are marked and varied compared to the wild-type peptide, including effects on secondary structure, stability and conformational changes. Each system showed differing levels of stability with some showing consistent, compact conformations whereas others displayed more flexible structures. Contrasts between comparable mutations at similar sites, such as A2T/A2V and D7H/D7N, show the location as well as the type of mutation have effects on protein structure observed in Ramachandran plots. We also report notable changes in peptide structure at residues remote to the site of substitution showing these mutations influence the entirety of Aβ. Salt-bridge profiles show this most clearly: addition or removal of charged residues affecting all salt-bridge interactions present in WT, even those remote from the site of mutation. Effects on secondary structure differ between mutations, most notably a change in incidence of β-strand, which has been linked to enhanced aggregational properties for the peptide. GFN2-xTB semi-empirical calculations show clear differences in binding energies of the copper-centre for each system

Structural variability of 4f and 5f thiocyanate complexes and dissociation of uranium(III)–thiocyanate bonds with increased ionicity

A series of complexes [Et4N][Ln(NCS)4(H2O)4] (Ln = Pr, Tb, Dy, Ho, Yb) have been structurally characterized, all showing the same structure, namely a distorted square antiprismatic coordination geometry, and the Ln–O and Ln–N bond lengths following the expected lanthanide contraction. When the counterion is Cs+, a different structural motif is observed and the eight-coordinate complex Cs5[Nd(NCS)8] isolated. The thorium compounds [Me4N]4[Th(NCS)7(NO3)] and [Me4N]4[Th(NCS)6(NO3)2] have been characterized, and high coordination numbers are also observed. Finally, attempts to synthesize a U(III) thiocyanate compound has been unsuccessful; from the reaction mixture, a heterocycle formed by condensation of five MeCN solvent molecules, possibly promoted by U(III), was isolated and structurally characterized. To rationalize the inability to isolate U(III) thiocyanate compounds, thin-layer cyclic voltammetry and IR spectroelectrochemistry have been utilized to explore the cathodic behavior of [Et4N]4[U(NCS)8] and [Et4N][U(NCS)5(bipy)2] along with a related uranyl compound [Et4N]3[UO2(NCS)5]. In all examples, the reduction triggers a rapid dissociation of [NCS]− ions and decomposition. Interestingly, the oxidation chemistry of [Et4N]3[UO2(NCS)5] in the presence of bipy gives the U(IV) compound [Et4N]4[U(NCS)8], an unusual example of a ligand-based oxidation triggering a metal-based reduction. The experimental results have been augmented by a computational investigation, concluding that the U(III)–NCS bond is more ionic than the U(IV)–NCS bond

Exploring the impact of mutation and post-translational modification on α-Synuclein: Insights from molecular dynamics simulations with and without copper

We report molecular dynamics simulations of two modifications to α-Synuclein, namely A53T mutation and phosphorylation at Ser129, which have been observed in Parkinson's disease patients. Both modifications are close to known metal binding sites, so as well as each modified peptide we also study Cu(II) bound to N-terminal and C-terminal residues. We show that A53T is predicted to cause increased β-sheet content of the peptide, with a persistent β-hairpin between residues 35–55 particularly notable. Phosphorylation has less effect on secondary structure but is predicted to significantly increase the size of the peptide, especially when bound to Cu(II), which is ascribed to reduced interaction of C-terminal sequence with central non-amyloid component. In addition, estimate of binding free energy to Cu(II) indicates A53T has little effect on metal-ion affinity, whereas phosphorylation markedly enhances the strength of binding. We suggest that the predicted changes in spatial extent and secondary structure of α-Synuclein may have implications for aggregation into Lewy bodies

How Cu(II) binding affects structure and dynamics of α-synuclein revealed by molecular dynamics simulations

We report accelerated molecular dynamics simulations of α-Synuclein and its complex with two Cu(II) ions bound to experimentally determined binding sites. Adding two Cu(II) ions, one bound to the N-terminal region and one to the C-terminus, decreases size and flexibility of the peptide while introducing significant new contacts within and between N-terminus and non-Aβ component (NAC). Cu(II) ions also alter the pattern of secondary structure within the peptide, inducing more and longer-lasting elements of secondary structure such as β-strands and hairpins. Free energy surfaces, obtained from reweighting the accelerated molecular dynamics boost potential, further demonstrate the restriction on size and flexibility that results from binding of copper ions

Computational investigation of copper-mediated conformational changes in α-synuclein dimer †

We report molecular dynamics simulation of dimers of α-synuclein, the peptide closely associated with onset of Parkinson's disease, both as metal-free dimer and with inter-chain bridging provided by Cu(ii) ions. Our investigation reveals that the presence of copper-induced inter-chain bridging not only stabilizes α-synuclein dimers, but also leads to enhanced β-sheet formation at critical regions within the N-terminal and NAC regions of the protein. These contacts are larger and longer-lived in the presence of copper, and as a result each peptide chain is more extended and less flexible than in the metal-free dimer. The persistence of these inter-peptide contacts underscores their significance in stabilising the dimers, potentially influencing the aggregation pathway. Moreover, the increased flexibility in the two termini, as well as the absence of persistent contacts in the metal-free dimer, correlates with the presence of amorphous aggregates. This phenomenon is known to mitigate fibrillation, while their absence in the metal-bound dimer suggests an increased propensity to form fibrils in the presence of copper ions