Explicit/implicit partitioning and a new explicit form of the generalized alpha method

Most finite element packages use the Newmark algorithm for time integration of structural dynamics. Various algorithms have been proposed to better optimize the high frequency dissipation of this algorithm. Hulbert and Chung proposed both implicit and explicit forms of the generalized alpha method. The algorithms optimize high frequency dissipation effectively, and despite recent work on algorithms that possess momentum conserving/energy dissipative properties in a non-linear context, the generalized alpha method remains an efficient way to solve many problems, especially with adaptive timestep control. However, the implicit and explicit algorithms use incompatible parameter sets and cannot be used together in a spatial partition, whereas this can be done for the Newmark algorithm, as Hughes and Liu demonstrated, and for the HHT-alpha algorithm developed from it. The present paper shows that the explicit generalized alpha method can be rewritten so that it becomes compatible with the implicit form. All four algorithmic parameters can be matched between the explicit and implicit forms. An element interface between implicit and explicit partitions can then be used, analogous to that devised by Hughes and Liu to extend the Newmark method. The stability of the explicit/implicit algorithm is examined in a linear context and found to exceed that of the explicit partition. The element partition is significantly less dissipative of intermediate frequencies than one using the HHT-alpha method. The explicit algorithm can also be rewritten so that the discrete equation of motion evaluates forces from displacements and velocities found at the predicted mid-point of a cycle. Copyright (C) 2003 John Wiley Sons, Ltd.

A partial velocity approach to subcycling structural dynamics

Subcycling, or the use of different timesteps at different nodes, can be an effective way of improving the computational efficiency of explicit transient dynamic structural solutions. The method that has been most widely adopted uses a nodal partition. extending the central difference method, in which small timestep updates are performed interpolating on the displacement at neighbouring large timestep nodes. This approach leads to narrow bands of unstable timesteps or statistical stability. It also can be in error due to lack of momentum conservation on the timestep interface. The author has previously proposed energy conserving algorithms that avoid the first problem of statistical stability. However, these sacrifice accuracy to achieve stability. An approach to conserve momentum on an element interface by adding partial velocities is considered here. Applied to extend the central difference method. this approach is simple. and has accuracy advantages. The method can be programmed by summing impulses of internal forces, evaluated using local element timesteps, in order to predict a velocity change at a node. However, it is still only statistically stable, so an adaptive timestep size is needed to monitor accuracy and to be adjusted if necessary. By replacing the central difference method with the explicit generalized alpha method. it is possible to gain stability by dissipating the high frequency response that leads to stability problems. However. coding the algorithm is less elegant, as the response depends on previous partial accelerations. Extension to implicit integration, is shown to be impractical due to the neglect of remote effects of internal forces acting across a timestep interface. (C) 2002 Elsevier Science B.V. All rights reserved.

Biochemical characterization of two mutants of human pyruvate dehydrogenase, F205L and T231A of the E1 alpha subunit

Mutations in the E1alpha subunit of the pyruvate dehydrogenase multienzyme complex may result in congenital lactic acidosis, but little is known about the consequences of these mutations at the enzymatic level. Here we characterize two mutants (F205L and T231A) of human pyruvate dehydrogenase in vitro, using the enzyme expressed in Escherichia coli. Wild-type and mutant proteins were purified successfully and their kinetic parameters were measured. F205L shows impaired binding of the thiamin diphosphate cofactor, which may explain why patients carrying this mutation respond to high-dose vitamin B-1 therapy. T231A has very low activity and a greatly elevated K-m for pyruvate, and this combination of effects would be expected to result in severe lactic acidosis. The results lead to a better understanding of the consequences of these mutations on the functional and structural properties of the enzyme, which may lead to improved therapies for patients carrying these mutations.

Structure of the HERG K+ channel S5P extracellular linker - Role of an amphipathic alpha-helix in C-type inactivation

The HERG K+ channel has very unusual kinetic behavior that includes slow activation but rapid inactivation. These features are critical for normal cardiac repolarization as well as in preventing lethal ventricular arrhythmias. Mutagenesis studies have shown that the extracellular peptide linker joining the fifth transmembrane domain to the pore helix is critical for rapid inactivation of the HERG K+ channel. This peptide linker is also considerably longer in HERG K+ channels, 40 amino acids, than in most other voltage-gated K+ channels. In this study we show that a synthetic 42-residue peptide corresponding to this linker region of the HERG K+ channel does not have defined structural elements in aqueous solution; however, it displays two well defined helical regions when in the presence of SDS micelles. The helices correspond to Trp(585)-Ile(593) and Gly(604)-Tyr(611) of the channel. The Trp(585)-Ile(593) helix has distinct hydrophilic and hydrophobic surfaces. The Gly(604)-Tyr(611) helix corresponds to an N-terminal extension of the pore helix. Electrophysiological studies of HERG currents following application of exogenous S5P peptides show that the amphipathic helix in the S5P linker interacts with the pore region of the channel in a voltage-dependent manner.

Insights into the structural basis for zinc inhibition of the glycine receptor

Histidines 107 and 109 in the glycine receptor ( GlyR) alpha(1) subunit have previously been identified as determinants of the inhibitory zinc-binding site. Based on modeling of the GlyR alpha(1) subunit extracellular domain by homology to the acetylcholine-binding protein crystal structure, we hypothesized that inhibitory zinc is bound within the vestibule lumen at subunit interfaces, where it is ligated by His(107) from one subunit and His(109) from an adjacent subunit. This was tested by co-expressing alpha(1) subunits containing the H107A mutation with alpha(1) subunits containing the H109A mutation. Although sensitivity to zinc inhibition is markedly reduced when either mutation is individually incorporated into all five subunits, the GlyRs formed by the co-expression of H107A mutant subunits with H109A mutant subunits exhibited an inhibitory zinc sensitivity similar to that of the wild type alpha(1) homomeric GlyR. This constitutes strong evidence that inhibitory zinc is coordinated at the interface between adjacent alpha(1) subunits. No evidence was found for beta subunit involvement in the coordination of inhibitory zinc, indicating that a maximum of two zinc-binding sites per alpha(1)beta receptor is sufficient for maximal zinc inhibition. Our data also show that two zinc-binding sites are sufficient for significant inhibition of alpha(1) homomers. The binding of zinc at the interface between adjacent alpha(1) subunits could restrict intersubunit movements, providing a feasible mechanism for the inhibition of channel activation by zinc.

Biochemical and structural characterization of β-lactamases tem-180 and tem-201

With the constant development of new antibiotics, selective pressure is a force to reckon when investigating antibiotic resistance. Although advantageous for medical treatments, it leads to increasing resistance. It is essential to use more potent and toxic antibiotics. Enzymes capable of hydrolyzing antibiotics are among the most common ways of resistance and TEM variants have been detected in several resistant isolates. Due to the rapid evolution of these variants, complex phenotypes have emerged and the need to understand their biological activity becomes crucial. To investigate the biochemical properties of TEM-180 and TEM-201 several computational methodologies have been used, allowing the comprehension of their structure and catalytic activity, which translates into their biological phenotype. In this work we intent to characterize the interface between these proteins and the several antibiotics used as ligands. We performed explicit solvent molecular dynamics (MD) simulations of these complexes and studied a variety of structural and energetic features. The interfacial residues show a distinct behavior when in complex with different antibiotics. Nevertheless, it was possible to identify some common Hot Spots among several complexes – Lys73, Tyr105 and Glu166. The structural changes that occur during the Molecular Dynamic (MD) simulation lead to the conclusion that these variants have an inherent capacity of adapting to the various antibiotics. This capability might be the reason why they can hydrolyze antibiotics that have not been described until now to be degraded by TEM variants. The results obtained with computational and experimental methodologies for the complex with Imipenem have shown that in order to this type of enzymes be able to acylate the antibiotics, they need to be capable to protect the ligand from water molecules.

Structural characterization and immunogenicity in wild-type and immune tolerant mice of degraded recombinant human interferon Alpha2b

Purpose: This study was conducted to study the influence of protein structure on the immunogenicity in wild-type and immune tolerant mice of well-characterized degradation products of recombinant human interferon alpha2b (rhIFNα2b). Methods: RhIFNα2b was degraded by metal-catalyzed oxidation (M), cross-linking with glutaraldehyde (G), oxidation with hydrogen peroxide (H), and incubation in a boiling water bath (B). The products were characterized with UV absorption, circular dichroism and fluorescence spectroscopy, gel permeation chromatography, reverse-phase high-pressure liquid chromatography, sodium dodecyl sulfate polyacrylamide gel electrophoresis, Western blotting, and mass spectrometry. The immunogenicity of the products was evaluated in wild-type mice and in transgenic mice immune tolerant for hIFNα2. Serum antibodies were detected by enzyme-linked immunosorbent assay or surface plasmon resonance. Results: M-rhIFNα2b contained covalently aggregated rhIFNα2b with three methionines partly oxidized to methionine sulfoxides. G-rhIFNα2b contained covalent aggregates and did not show changes in secondary structure. H-rhIFNα2b was only chemically changed with four partly oxidized methionines. B-rhIFNα2b was largely unfolded and heavily aggregated. Nontreated (N) rhIFNα2b was immunogenic in the wild-type mice but not in the transgenic mice, showing that the latter were immune tolerant for rhIFNα2b. The anti-rhIFNα2b antibody levels in the wild-type mice depended on the degradation product: M-rhIFNα2b > H-rhIFNα2b ∼ N-rhIFNα2b ≫ B-rhIFNα2b; G-rhIFNα2b did not induce anti-rhIFNα2b antibodies. In the transgenic mice, only M-rhIFNα2b could break the immune tolerance. Conclusions: RhIFNα2b immunogenicity is related to its structural integrity. Moreover, the immunogenicity of aggregated rhIFNα2b depends on the structure and orientation of the constituent protein molecules and/or on the aggregate size.

Structural investigation of the Bacillus subtilis morphogenic factor RodZ

A thesis to obtain a Master degree in Structural and Functional Biochemistry

The fourth extracellular loop of the alpha subunit of Na,K-ATPase. Functional evidence for close proximity with the second extracellular loop.

Na,K-ATPase is the main active transport system that maintains the large gradients of Na(+) and K(+) across the plasma membrane of animal cells. The crystal structure of a K(+)-occluding conformation of this protein has been recently published, but the movements of its different domains allowing for the cation pumping mechanism are not yet known. The structure of many more conformations is known for the related calcium ATPase SERCA, but the reliability of homology modeling is poor for several domains with low sequence identity, in particular the extracellular loops. To better define the structure of the large fourth extracellular loop between the seventh and eighth transmembrane segments of the alpha subunit, we have studied the formation of a disulfide bond between pairs of cysteine residues introduced by site-directed mutagenesis in the second and the fourth extracellular loop. We found a specific pair of cysteine positions (Y308C and D884C) for which extracellular treatment with an oxidizing agent inhibited the Na,K pump function, which could be rapidly restored by a reducing agent. The formation of the disulfide bond occurred preferentially under the E2-P conformation of Na,K-ATPase, in the absence of extracellular cations. Using recently published crystal structure and a distance constraint reproducing the existence of disulfide bond, we performed an extensive conformational space search using simulated annealing and showed that the Tyr(308) and Asp(884) residues can be in close proximity, and simultaneously, the SYGQ motif of the fourth extracellular loop, known to interact with the extracellular domain of the beta subunit, can be exposed to the exterior of the protein and can easily interact with the beta subunit.

Ab initio modeling and molecular dynamics simulation of the alpha 1b-adrenergic receptor activation.

This work describes the ab initio procedure employed to build an activation model for the alpha 1b-adrenergic receptor (alpha 1b-AR). The first version of the model was progressively modified and complicated by means of a many-step iterative procedure characterized by the employment of experimental validations of the model in each upgrading step. A combined simulated (molecular dynamics) and experimental mutagenesis approach was used to determine the structural and dynamic features characterizing the inactive and active states of alpha 1b-AR. The latest version of the model has been successfully challenged with respect to its ability to interpret and predict the functional properties of a large number of mutants. The iterative approach employed to describe alpha 1b-AR activation in terms of molecular structure and dynamics allows further complications of the model to allow prediction and interpretation of an ever-increasing number of experimental data.

Vaniprevir with pegylated interferon alpha-2a and ribavirin in treatment-naïve patients with chronic hepatitis C: a randomized phase II study.

Vaniprevir (MK-7009) is a macrocyclic hepatitis C virus (HCV) nonstructural protein 3/4A protease inhibitor. The aim of the present phase II study was to examine virologic response rates with vaniprevir in combination with pegylated interferon alpha-2a (Peg-IFN-α-2a) plus ribavirin (RBV). In this double-blind, placebo-controlled, dose-ranging study, treatment-naïve patients with HCV genotype 1 infection (n = 94) were randomized to receive open-label Peg-IFN-α-2a (180 μg/week) and RBV (1,000-1,200 mg/day) in combination with blinded placebo or vaniprevir (300 mg twice-daily [BID], 600 mg BID, 600 mg once-daily [QD], or 800 mg QD) for 28 days, then open-label Peg-IFN-α-2a and RBV for an additional 44 weeks. The primary efficacy endpoint was rapid viral response (RVR), defined as undetectable plasma HCV RNA at week 4. Across all doses, vaniprevir was associated with a rapid two-phase decline in viral load, with HCV RNA levels approximately 3 log(10) IU/mL lower in vaniprevir-treated patients, compared to placebo recipients. Rates of RVR were significantly higher in each of the vaniprevir dose groups, compared to the control regimen (68.8%-83.3% versus 5.6%; P < 0.001 for all comparisons). There were numerically higher, but not statistically significant, early and sustained virologic response rates with vaniprevir, as compared to placebo. Resistance profile was predictable, with variants at R155 and D168 detected in a small number of patients. No relationship between interleukin-28B genotype and treatment outcomes was demonstrated in this study. The incidence of adverse events was generally comparable between vaniprevir and placebo recipients; however, vomiting appeared to be more common at higher vaniprevir doses. CONCLUSION: Vaniprevir is a potent HCV protease inhibitor with a predictable resistance profile and favorable safety profile that is suitable for QD or BID administration.

Qualités psychométriques de la version française de la TAS-20 et prévalence de l'alexithymie chez 264 adolescents tout-venant [The 20-item Toronto alexithymia scale: structural validity, internal consistency and prevalence of alexithymia in a Swiss adolescent sample]

L'objectif de cette étude est d'examiner la structure factorielle et la consistance interne de la TAS-20 sur un échantillon d'adolescents (n = 264), ainsi que de décrire la distribution des caractéristiques alexithymiques dans cet échantillon. La structure à trois facteurs de la TAS-20 a été confirmée par notre analyse factorielle confirmatoire. La consistance interne, mesurée à l'aide d'alpha de Cronbach, est acceptable pour le premier facteur (difficulté à identifier les sentiments (DIF)), bonne pour le second (difficulté à verbaliser les sentiments (DDF)), mais en revanche, faible pour le troisième facteur (pensées orientées vers l'extérieur (EOT)). Les résultats d'une Anova mettent en évidence une tendance linéaire indiquant que plus l'âge augmente plus le niveau d'alexithymie (score total TAS-20), la difficulté à identifier les sentiments et les pensées orientées vers l'extérieur diminuent. En ce qui concerne la prévalence de l'alexithymie, on remarque en effet que 38,5 % des adolescents de moins de 16 ans sont considérés comme alexithymiques, contre 30,1 % des 16-17 ans et 22 % des plus de 17 ans. Notre étude indique donc que la TAS-20 est un instrument adéquat pour évaluer l'alexithymie à l'adolescence, tout en suggérant quelques précautions étant donné l'aspect développemental de cette période.

Syntaxin1a variants lacking an N-peptide or bearing the LE mutation bind to Munc18a in a closed conformation.

In neurons, soluble N-ethylmaleimide-sensitive factor attachment receptor (SNARE) proteins drive the fusion of synaptic vesicles to the plasma membrane through the formation of a four-helix SNARE complex. Members of the Sec1/Munc18 protein family regulate membrane fusion through interactions with the syntaxin family of SNARE proteins. The neuronal protein Munc18a interacts with a closed conformation of the SNARE protein syntaxin1a (Syx1a) and with an assembled SNARE complex containing Syx1a in an open conformation. The N-peptide of Syx1a (amino acids 1-24) has been implicated in the transition of Munc18a-bound Syx1a to Munc18a-bound SNARE complex, but the underlying mechanism is not understood. Here we report the X-ray crystal structures of Munc18a bound to Syx1a with and without its native N-peptide (Syx1aΔN), along with small-angle X-ray scattering (SAXS) data for Munc18a bound to Syx1a, Syx1aΔN, and Syx1a L165A/E166A (LE), a mutation thought to render Syx1a in a constitutively open conformation. We show that all three complexes adopt the same global structure, in which Munc18a binds a closed conformation of Syx1a. We also identify a possible structural connection between the Syx1a N-peptide and SNARE domain that might be important for the transition of closed-to-open Syx1a in SNARE complex assembly. Although the role of the N-peptide in Munc18a-mediated SNARE complex assembly remains unclear, our results demonstrate that the N-peptide and LE mutation have no effect on the global conformation of the Munc18a-Syx1a complex.

Sequence conservation in Plasmodium falciparum alpha-helical coiled coil domains proposed for vaccine development.

BACKGROUND: The availability of the P. falciparum genome has led to novel ways to identify potential vaccine candidates. A new approach for antigen discovery based on the bioinformatic selection of heptad repeat motifs corresponding to alpha-helical coiled coil structures yielded promising results. To elucidate the question about the relationship between the coiled coil motifs and their sequence conservation, we have assessed the extent of polymorphism in putative alpha-helical coiled coil domains in culture strains, in natural populations and in the single nucleotide polymorphism data available at PlasmoDB. METHODOLOGY/PRINCIPAL FINDINGS: 14 alpha-helical coiled coil domains were selected based on preclinical experimental evaluation. They were tested by PCR amplification and sequencing of different P. falciparum culture strains and field isolates. We found that only 3 out of 14 alpha-helical coiled coils showed point mutations and/or length polymorphisms. Based on promising immunological results 5 of these peptides were selected for further analysis. Direct sequencing of field samples from Papua New Guinea and Tanzania showed that 3 out of these 5 peptides were completely conserved. An in silico analysis of polymorphism was performed for all 166 putative alpha-helical coiled coil domains originally identified in the P. falciparum genome. We found that 82% (137/166) of these peptides were conserved, and for one peptide only the detected SNPs decreased substantially the probability score for alpha-helical coiled coil formation. More SNPs were found in arrays of almost perfect tandem repeats. In summary, the coiled coil structure prediction was rarely modified by SNPs. The analysis revealed a number of peptides with strictly conserved alpha-helical coiled coil motifs. CONCLUSION/SIGNIFICANCE: We conclude that the selection of alpha-helical coiled coil structural motifs is a valuable approach to identify potential vaccine targets showing a high degree of conservation.

The fourth transmembrane segment of the Na,K-ATPase alpha subunit: a systematic mutagenesis study.

The Na,K-ATPase is a major ion-motive ATPase of the P-type family responsible for many aspects of cellular homeostasis. To determine the structure of the pathway for cations across the transmembrane portion of the Na,K-ATPase, we mutated 24 residues of the fourth transmembrane segment into cysteine and studied their function and accessibility by exposure to the sulfhydryl reagent 2-aminoethyl-methanethiosulfonate. Accessibility was also examined after treatment with palytoxin, which transforms the Na,K-pump into a cation channel. Of the 24 tested cysteine mutants, seven had no or a much reduced transport function. In particular cysteine mutants of the highly conserved "PEG" motif had a strongly reduced activity. However, most of the non-functional mutants could still be transformed by palytoxin as well as all of the functional mutants. Accessibility, determined as a 2-aminoethyl-methanethiosulfonate-induced reduction of the transport activity or as inhibition of the membrane conductance after palytoxin treatment, was observed for the following positions: Phe(323), Ile(322), Gly(326), Ala(330), Pro(333), Glu(334), and Gly(335). In accordance with a structural model of the Na,K-ATPase obtained by homology modeling with the two published structures of sarcoplasmic and endoplasmic reticulum calcium ATPase (Protein Data Bank codes 1EUL and 1IWO), the results suggest the presence of a cation pathway along the side of the fourth transmembrane segment that faces the space between transmembrane segments 5 and 6. The phenylalanine residue in position 323 has a critical position at the outer mouth of the cation pathway. The residues thought to form the cation binding site II ((333)PEGL) are also part of the accessible wall of the cation pathway opened by palytoxin through the Na,K-pump.