Effect of knee joint angle on motor unit synchronization

Activity of the vasti has been argued to vary through knee range of movement due to changes in passive support of the patellofemoral joint and the relative contribution of these muscles to knee extension. Efficient function of the knee is dependent on optimal control of the patellofemoral joint, largely through coordinated activity of the medial and lateral quadriceps. Motor unit synchronization may provide a mechanism to coordinate the activity of vastus medialis (VMO) and vastus lateralis (VL), and may be more critical in positions of reduced passive support for the patellofemoral joint (i.e., full extension). Therefore, the aim of this study was to determine whether the degree of motor unit synchronization between the vasti muscles is dependent on joint angle. Electromyographic (EMG) recordings of single motor unit action potentials (MUAPs) were made from VMO and multiunit recordings from VL during isometric contractions of the quadriceps at 0 degrees, 30 degrees, and 60 degrees of knee flexion. The degree of synchronization between motor unit firing was evaluated by identification of peaks in the rectified EMG averages of VL, triggered from MUA-Ps in VMO. The proportion of cases in which there was a significant peak in the triggered averages was calculated. There was no significant difference in the degree of synchronization between the vasti at different knee angles (p = 0.57). These data suggest that this basic coordinative mechanism between the vasti muscles is controlled consistently throughout knee range of motion, and is not augmented at specific angles where the requirement for dynamic control of stability is increased. (D 2006 Orthopaedic Research Society. Published by Wiley Periodicals, Inc.

Spectroscopic characterization of copper(II) binding to the immunosuppressive drug mycophenolic acid

Mycophenolic acid (MPA) is a drug that has found widespread use as an immunosuppressive agent which limits rejection of transplanted organs. Optimal use of this drug is hampered by gastrointestinal side effects which can range in severity. One mechanism by which MPA causes gastropathy may involve a direct interaction between the drug and gastric phospholipids. To combat this interaction we have investigated the potential of MPA to coordinate Cu(II), a metal which has been used to inhibit gastropathy associated with use of the NSAID indomethacin. Using a range of spectroscopic techniques we show that Cu(II) is coordinated to two MPA molecules via carboxylates and, at low pH, water ligands. The copper complex formed is stable in solution as assessed by mass spectrometry and H-1 NMR diffusion experiments. Competition studies with glycine and albumin indicate that the copper-MPA complex will release Cu(II) to amino acids and proteins thereby allowing free MPA to be transported to its site of action. Transfer to serum albumin proceeds via a Cu(MPA)(albumin) ternary complex. These results raise the possibility that copper complexes of MPA may be useful in a therapeutic situation.

Automatic initialisation of 3D deformable models for cartilage segmentation

Deformable models are a highly accurate and flexible approach to segmenting structures in medical images. The primary drawback of deformable models is that they are sensitive to initialisation, with accurate and robust results often requiring initialisation close to the true object in the image. Automatically obtaining a good initialisation is problematic for many structures in the body. The cartilages of the knee are a thin elastic material that cover the ends of the bone, absorbing shock and allowing smooth movement. The degeneration of these cartilages characterize the progression of osteoarthritis. The state of the art in the segmentation of the cartilage are 2D semi-automated algorithms. These algorithms require significant time and supervison by a clinical expert, so the development of an automatic segmentation algorithm for the cartilages is an important clinical goal. In this paper we present an approach towards this goal that allows us to automatically providing a good initialisation for deformable models of the patella cartilage, by utilising the strong spatial relationship of the cartilage to the underlying bone.

3D statistical shape models to embed spatial relationship information

This paper presents the creation of 3D statistical shape models of the knee bones and their use to embed information into a segmentation system for MRIs of the knee. We propose utilising the strong spatial relationship between the cartilages and the bones in the knee by embedding this information into the created models. This information can then be used to automate the initialisation of segmentation algorithms for the cartilages. The approach used to automatically generate the 3D statistical shape models of the bones is based on the point distribution model optimisation framework of Davies. Our implementation of this scheme uses a parameterized surface extraction algorithm, which is used as the basis for the optimisation scheme that automatically creates the 3D statistical shape models. The current approach is illustrated by generating 3D statistical shape models of the patella, tibia and femoral bones from a segmented database of the knee. The use of these models to embed spatial relationship information to aid in the automation of segmentation algorithms for the cartilages is then illustrated.

MHC class I bound to an immunodominant Theileria parva epitope demonstrates unconventional presentation to T cell receptors

T cell receptor (TCR) recognition of peptide-MHC class I (pMHC) complexes is a crucial event in the adaptive immune response to pathogens. Peptide epitopes often display a strong dominance hierarchy, resulting in focusing of the response on a limited number of the most dominant epitopes. Such T cell responses may be additionally restricted by particular MHC alleles in preference to others. We have studied this poorly understood phenomenon using Theileria parva, a protozoan parasite that causes an often fatal lymphoproliferative disease in cattle. Despite its antigenic complexity, CD8+ T cell responses induced by infection with the parasite show profound immunodominance, as exemplified by the Tp1(214-224) epitope presented by the common and functionally important MHC class I allele N*01301. We present a high-resolution crystal structure of this pMHC complex, demonstrating that the peptide is presented in a distinctive raised conformation. Functional studies using CD8+ T cell clones show that this impacts significantly on TCR recognition. The unconventional structure is generated by a hydrophobic ridge within the MHC peptide binding groove, found in a set of cattle MHC alleles. Extremely rare in all other species, this feature is seen in a small group of mouse MHC class I molecules. The data generated in this analysis contribute to our understanding of the structural basis for T cell-dependent immune responses, providing insight into what determines a highly immunogenic p-MHC complex, and hence can be of value in prediction of antigenic epitopes and vaccine design.

Flexibility of short-strand RNA in aqueous solution as revealed by molecular dynamics simulation:are A′-RNA and A-RNA distinct conformational structures?

We use molecular dynamics simulations to compare the conformational structure and dynamics of a 21-base pair RNA sequence initially constructed according to the canonical A-RNA and A'-RNA forms in the presence of counterions and explicit water. Our study aims to add a dynamical perspective to the solid-state structural information that has been derived from X-ray data for these two characteristic forms of RNA. Analysis of the three main structural descriptors commonly used to differentiate between the two forms of RNA namely major groove width, inclination and the number of base pairs in a helical twist over a 30 ns simulation period reveals a flexible structure in aqueous solution with fluctuations in the values of these structural parameters encompassing the range between the two crystal forms and more. This provides evidence to suggest that the identification of distinct A-RNA and A'-RNA structures, while relevant in the crystalline form, may not be generally relevant in the context of RNA in the aqueous phase. The apparent structural flexibility observed in our simulations is likely to bear ramifications for the interactions of RNA with biological molecules (e.g. proteins) and non-biological molecules (e.g. non-viral gene delivery vectors). © CSIRO 2009.

The design of novel antifolates against pneumocystis carinii

The pneumonia caused by Pneumocystis carinii is ultimately responsible for the death of many acquired immunodeficiency syndrome (AIDS) patients. Large doses of trimethoprim and pyrimethamine in combination with a sulphonamide and/or pentamidine suppress the infection but produce serious side-effects and seldom prevent recurrence after treatment withdrawal. However, the partial success of the aforementioned antifolates, and also trimetrexate used alone, does suggest dihydrofolate reductase (DHFR) as a target for the development of antipneumocystis agents. From the DHFR inhibitory activities of 3'-substituted pyrimethamine analogues it was suggested that the 3'-(3'',3''-dimethyltriazen-1''-yl) substituent may be responsible for the greater activity for the P.carinii over the mammalian enzyme. Crystallographic and molecular modeling studies revealed considerable geometrical and electronic differences between the triazene and the chemically related formamidine functions that may account for the differences in DHFR inhibitory profiles. Structural and electronic parameters calculated for a series of 3'-(3'',3''-disubstitutedtriazen-1''-yl) pyrimethamine analogues did not correlate with the DHFR inhibitory activities. However, the in vitro screening against P.carinii DHFR revealed that the 3''-hydroxyethyl-3''-benzyl analogue was the most active and selective. Models of the active sites of human and P.carinii DHFRs were constructed using DHFR sequence and structural homology data which had identified key residues involved in substrate and cofactor binding. Low energy conformations of the 3'',3''-dimethyl and 3''-hydroxyethyl-3''-benzyle analogues, determined from nuclear magnetic resonance studies and theoretical calculations, were docked by superimposing the diaminopyrimidine fragment onto a previously docked pyrimethamine analogue. Enzyme kinetic data supported the 3''-hydroxyethyl-3''-benzyl moiety being located in the NADPH binding groove. The 3''-benzyl substituent was able to locate to within 3 AA of a valine residue in the active site of P.carinii DHFR thereby producing a hydrophobic contact. The equivalent residue in human DHFR is threonine, more hydrophilic and less likely to be involved in such a contact. This difference may account for the greater inhibitory activity this analogue has for P.carinii DHFR and provide a basis for future drug design. From an in vivo model of PCP in immunosuppressed rats it was established that the 3"-hydroxyethyl-3"-benzyl analogue was able to reduce the.P.carinii burden more effectively with increasing doses, without causmg any visible signs of toxicity. However, equivalent doses were not as effective as pentamidine, a current treatment of choice for Pneumocystis carinii pneumonia.

Peptide length significantly influences in vitro affinity for MHC class II molecules

Class II Major Histocompatibility Complex (MHC) molecules have an open-ended binding groove which can accommodate peptides of varying lengths. Several studies have demonstrated that peptide flanking residues (PFRs) which lie outside the core binding groove can influence peptide binding and T cell recognition. By using data from the AntiJen database we were able to characterise systematically the influence of PFRs on peptide affinity for MHC class II molecules.

Hydrotalcite intercalated siRNA:computational characterization of the interlayer environment

Using molecular dynamics (MD) simulations, we explore the structural and dynamical properties of siRNA within the intercalated environment of a Mg:Al 2:1 Layered Double Hydroxide (LDH) nanoparticle. An ab initio force field (Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies: COMPASS) is used for the MD simulations of the hybrid organic-inorganic systems. The structure, arrangement, mobility, close contacts and hydrogen bonds associated with the intercalated RNA are examined and contrasted with those of the isolated RNA. Computed powder X-ray diffraction patterns are also compared with related LDH-DNA experiments. As a method of probing whether the intercalated environment approximates the crystalline or rather the aqueous state, we explore the stability of the principle parameters (e.g., the major groove width) that differentiate both A- and A'- crystalline forms of siRNA and contrast this with recent findings for the same siRNA simulated in water. We find the crystalline forms remain structurally distinct when intercalated, whereas this is not the case in water. Implications for the stability of hybrid LDH-RNA systems are discussed.

Design, synthesis and biological evaluation of potential antibacterial oligonucleotides

Purine and pyrimidine triplex-forming oligonucleotides (TFOs), as potential antibacterial agents, were designed to bind by Hoogsteen and reverse Hoogsteen hydrogen bonds in a sequence specific manner in the major groove of genomic DNA at specific polypurine sites within the gyrA gene of E. coli and S. pneumoniae. Sequences were prepared by automated synthesis, with purification and characterisation determined by high performance liquid chromatograpy, capillary electrophoresis and mass spectrometry. Triplex stability was assessed using melting curves where the binding of the third strand to the duplex target, was assessed over a temperature range of 0-80°C, and at pH 6.4 and 7.2. The most successful of the unmodified TFOs (6) showed a Tm value of 26 °C at both pH values with binding via reverse Hoogsteen bonds. Binding to genomic DNA was also demonstrated by spectrofluorimetry, using fluorescein-labelled TFOs, from which dissociation constants were determined. Modifications in the form of 5mC, 5' acridine attachment, phosphorothioation, 2'-0-methylation and phosphoramidation, were made in order to. increase Tm values. Phosphoramidate modification was the most with increased Tm values of 42°C. However, the final purity of these sequences was poor due to their difficult syntheses. FACS (fluorescent activated cell sorting) analysis was used to determine the potential uptake of a fluorescently labelled analogue of 6 via passive, coJd shock mediated, and anionic liposome aided, uptake. This was established at 20°C and 37°C. At both temperatures anionic lipid-mediated uptake produced unrivalled fluorescence, equivalent to 20 and 43% at 20 and 37°C respectively. Antibacterial activity of each oligonucleotide was assessed by viable count anaJysis relying on passive uptake, cold shocking techniques, chlorpromazine-mediated uptake, and, cationic and anionic lipid-aided uptake. All oligonucleotides were assessed for their ability to enhance uptake, which is a major barrier to the effectiveness of these agents. Compound 6 under cold shocking conditions produced the greatest consistent decline in colony forming units per ml. Results for this compound were sometimes variable indicating inconsistent uptake by this particular assay method.

Extracellular loops 1 and 3 and their associated transmembrane regions of the calcitonin receptor-like receptor are needed for CGRP receptor function

The first and third extracellular loops (ECL) of G protein-coupled receptors (GPCRs) have been implicated in ligand binding and receptor function. This study describes the results of an alanine/leucine scan of ECLs 1 and 3 and loop-associated transmembrane (TM) domains of the secretin-like GPCR calcitonin receptor-like receptor which associates with receptor activity modifying protein 1 to form the CGRP receptor. Leu195Ala, Val198Ala and Ala199Leu at the top of TM2 all reduced aCGRP-mediated cAMP production and internalization; Leu195Ala and Ala199Leu also reduced aCGRP binding. These residues form a hydrophobic cluster within an area defined as the "minor groove" of rhodopsin-like GPCRs. Within ECL1, Ala203Leu and Ala206Leu influenced the ability of aCGRP to stimulate adenylate cyclase. In TM3, His219Ala, Leu220Ala and Leu222Ala have influences on aCGRP binding and cAMP production; they are likely to indirectly influence the binding site for aCGRP as well as having an involvement in signal transduction. On the exofacial surfaces of TMs 6 and 7, a number of residues were identified that reduced cell surface receptor expression, most noticeably Leu351Ala and Glu357Ala in TM6. The residues may contribute to the RAMP1 binding interface. Ile360Ala impaired aCGRP-mediated cAMP production. Ile360 is predicted to be located close to ECL2 and may facilitate receptor activation. Identification of several crucial functional loci gives further insight into the activation mechanism of this complex receptor system and may aid rational drug design.

Identification of the HLA-DM/HLA-DR interface

Human leukocyte antigen (HLA)-DM is a critical participant in antigen presentation that catalyzes the dissociation of the Class II-associated Invariant chain-derived Peptide (CLIP) from the major histocompatibility complex (MHC) Class II molecules. There is competition amongst peptides for access to an MHC Class II groove and it has been hypothesised that DM functions as a 'peptide editor' that catalyzes the replacement of one peptide for another within the groove. It is established that the DM catalyst interacts directly with the MHC Class II but the precise location of the interface is unknown. Here, we combine previously described mutational data with molecular docking and energy minimisation simulations to identify a putative interaction site of >4000A2 which agrees with known point mutational data for both the DR and DM molecule. The docked structure is validated by comparison with experimental data and previously determined properties of protein-protein interfaces. A possible dissociation mechanism is suggested by the presence of an acidic cluster near the N terminus of the bound peptide.

Modulation of interleukin-6 and matrix metalloproteinase 2 expression in human fibroblast-like synoviocytes by functional ionotropic glutamate receptors

Objective. Patients with rheumatoid arthritis (RA) have increased concentrations of the amino acid glutamate in synovial fluid. This study was undertaken to determine whether glutamate receptors are expressed in the synovial joint, and to determine whether activation of glutamate receptors on human synoviocytes contributes to RA disease pathology. Methods. Glutamate receptor expression was examined in tissue samples from rat knee joints and in human fibroblast-like synoviocytes (FLS). FLS from 5 RA patients and 1 normal control were used to determine whether a range of glutamate receptor antagonists influenced expression of the proinflammatory cytokine interleukin-6 (IL-6), enzymes involved in matrix degradation and cytokine processing (matrix metalloproteinase 2 [MMP-2] and MMP-9), and the inhibitors of these enzymes (tissue inhibitor of metalloproteinases 1 [TIMP-1] and TIMP-2). IL-6 concentrations were determined by enzyme-linked immunosorbent assay, MMP activity was measured by gelatin zymography, and TIMP activity was determined by reverse zymography. Fluorescence imaging of intracellular calcium concentrations in live RA FLS stimulated with specific antagonists was used to reveal functional activation of glutamate receptors that modulated IL-6 or MMP-2. Results. Ionotropic and metabotropic glutamate receptor subunit mRNA were expressed in the patella, fat pad, and meniscus of the rat knee and in human articular cartilage. Inhibition of N-methyl-D-aspartate (NMDA) receptors in RA FLS increased proMMP-2 release, whereas non-NMDA ionotropic glutamate receptor antagonists reduced IL-6 production by these cells. Stimulation with glutamate, NMDA, or kainate (KA) increased intracellular calcium concentrations in RA FLS, demonstrating functional activation of specific ionotropic glutamate receptors. Conclusion. Our findings indicate that activation of NMDA and KA glutamate receptors on human synoviocytes may contribute to joint destruction by increasing IL-6 expression. © 2007, American College of Rheumatology.

Predicting class II MHC-Peptide binding:a kernel based approach using similarity scores

Background - Modelling the interaction between potentially antigenic peptides and Major Histocompatibility Complex (MHC) molecules is a key step in identifying potential T-cell epitopes. For Class II MHC alleles, the binding groove is open at both ends, causing ambiguity in the positional alignment between the groove and peptide, as well as creating uncertainty as to what parts of the peptide interact with the MHC. Moreover, the antigenic peptides have variable lengths, making naive modelling methods difficult to apply. This paper introduces a kernel method that can handle variable length peptides effectively by quantifying similarities between peptide sequences and integrating these into the kernel. Results - The kernel approach presented here shows increased prediction accuracy with a significantly higher number of true positives and negatives on multiple MHC class II alleles, when testing data sets from MHCPEP [1], MCHBN [2], and MHCBench [3]. Evaluation by cross validation, when segregating binders and non-binders, produced an average of 0.824 AROC for the MHCBench data sets (up from 0.756), and an average of 0.96 AROC for multiple alleles of the MHCPEP database. Conclusion - The method improves performance over existing state-of-the-art methods of MHC class II peptide binding predictions by using a custom, knowledge-based representation of peptides. Similarity scores, in contrast to a fixed-length, pocket-specific representation of amino acids, provide a flexible and powerful way of modelling MHC binding, and can easily be applied to other dynamic sequence problems.

Statistical deconvolution of enthalpic energetic contributions to MHC-peptide binding affinity

Background - MHC Class I molecules present antigenic peptides to cytotoxic T cells, which forms an integral part of the adaptive immune response. Peptides are bound within a groove formed by the MHC heavy chain. Previous approaches to MHC Class I-peptide binding prediction have largely concentrated on the peptide anchor residues located at the P2 and C-terminus positions. Results - A large dataset comprising MHC-peptide structural complexes was created by re-modelling pre-determined x-ray crystallographic structures. Static energetic analysis, following energy minimisation, was performed on the dataset in order to characterise interactions between bound peptides and the MHC Class I molecule, partitioning the interactions within the groove into van der Waals, electrostatic and total non-bonded energy contributions. Conclusion - The QSAR techniques of Genetic Function Approximation (GFA) and Genetic Partial Least Squares (G/PLS) algorithms were used to identify key interactions between the two molecules by comparing the calculated energy values with experimentally-determined BL50 data. Although the peptide termini binding interactions help ensure the stability of the MHC Class I-peptide complex, the central region of the peptide is also important in defining the specificity of the interaction. As thermodynamic studies indicate that peptide association and dissociation may be driven entropically, it may be necessary to incorporate entropic contributions into future calculations.