Analysis of Hydrogen Bonding and Stability of Protein Secondary Structures in Molecular Dynamics Simulation

Molecular Dynamics (MD) simulations provide an atomic level account of the molecular motions and have proven to be immensely useful in the investigation of the dynamical structure of proteins. Once an MD trajectory is obtained, specific interactions at the molecular level can be directly studied by setting up appropriate combinations of distance and angle monitors. However, if a study of the dynamical behavior of secondary structures in proteins becomes important, this approach can become unwieldy. We present herein a method to study the dynamical stability of secondary structures in proteins, based on a relatively simple analysis of backbone hydrogen bonds. The method was developed for studying the thermal unfolding of beta-lactamases, but can be extended to other systems and adapted to study relevant properties.

Edge-sharing bioctahedral diruthenium(III) complexes containing methoxy and carboxylate bridges:x-ray structures, redox behavior, and core stability

Edge-sharing bioctahedral (ESBO) complexes [Ru-2(OMe)(O2CC6H4-p-X)3(1-MeIm)(4)](ClO4)2 (X = OMe (1a), Me (1b)) and [Ru-2(O2CC6H4-P-X)(4)(1-MeIm)(4)](ClO4)(2) (X = OMe (2a), Me (2b)) are prepared by reacting Ru2Cl(O(2)CR)(4) with 1-methylimidazole (1-MeIm) in methanol followed by treatment with NaClO4. Complex 2a and the PF6- salt (1a') of 1a have been structurally characterized. Crystal data for 1a.1.5MeCN. 0.5Et(2)O: triclinic, P (1) over bar, a = 13.125(2) Angstrom, b = 15.529(3) Angstrom, c 17.314(5) Angstrom, a; 67.03(2)degrees, beta 68.05(2)degrees, gamma = 81.38(1)degrees, V 3014(1) Angstrom(3), Z = 2. Crystal data for 2a: triclinic, P (1) over bar, a 8.950(1) Angstrom, b = 12.089(3) Angstrom, c = 13.735(3) Angstrom, alpha 81.09(2)degrees, beta = 72.27(1)degrees, gamma = 83.15(2)degrees, V = 1394(1) Angstrom(3), Z = 1. The complexes consist of a diruthenium(III) unit held by two monoatomic and two three-atom bridging ligands. The 1-MeIm ligands are at the terminal sites of the [Ru-2(mu-L)(eta(1):mu-O(2)CR)(eta(1):eta(1):mu-O(2)CR)(2)](2+) core having a Ru-Ru single bond (L = OMe or eta(1)-O(2)CR). The Ru-Ru distance and the Ru-O-Ru angle in the core of 1a' and 2a are 2.49 Angstrom and similar to 76 degrees. The complexes undergo one-electron oxidation and reduction processes in MeCN-0.1 M TBAP to form mixed-valence diruthenium species with Ru-Ru bonds of orders 1.5 and 0.5, respectively.

Gradient solid electrolytes for thermodynamic measurements: System Na2CO3-Na2SO4

The thermodynamic properties of Na2CO3-Na2SO4 solid solution with hexagonal structure have been measured in the temperature range of 873 to 1073 K, using a composite-gradient solid electrolyte. The cell used can be represented as The composite-gradient solid electrolyte consisted of pure Na2CO3 at one extremity and the solid solution under study at the other, with variation in composition across the electrolyte. A CO2 + O2 + Ar gas mixture was used to fix the chemical potential of sodium at each electrode. The Nernstian response of the cell to changes in partial pressures of CO2 and O2 at the electrodes has been demonstrated. The activity of Na2CO3 in the solid solution was measured by two techniques. In the first method, the electromotive force (emf) of the cell was measured with the same CO2 + O2 + Ar mixture at both electrodes. The resultant emf is directly related to the activity of Na2CO3 at the solid solution electrode. By the second approach, the activity was calculated from the difference in compositions Of CO2 + O2 + Ar mixtures at the two electrodes required to produce a null emf. Both methods gave identical results. The second method is more suitable for gradient solid electrolytes that exhibit significant electronic conduction. The activity of Na2CO3 exhibits positive deviation from Raoult's law. The excess Gibbs' energy of mixing of the solid solution can be represented using a subregular solution model such as the following: DELTAG(E) = X(1 - X)[6500(+/-200)X + 3320(+/-80)(1 - X)J mol-1 where X is the mole fraction of Na2CO3. By combining this information with the phase diagram, mixing properties of the liquid phase are obtained.

High-temperature stability of Ca2ZrSi4O12

The standard Gibbs free energy of formation of orthorhombic Ca2ZrSi4O12 from component oxides ZrO2 (monoclinic), CaO (rock salt), and SiO2 (quartz) has been determined in the temperature range 973 to 1273 K using a solid-state cell incorporating single-crystal CaF2 as the electrolyte: Delta G(f) degrees = -219930 + 11.77T (+/- 1500) J.mol(-1) This is the only quantitative information now available on the stability of Ca2ZrSi4O12.

Interactions of L-serine at the active site of serine hydroxymethyltransferases: induction of thermal stability

Serine hydroxymethyltransferase (SHMT), EC 2.1.2.1, exhibits broad substrate and reaction specificity. In addition to cleaving many 3-hydroxyamino acids to glycine and an aldehyde, the enzyme also catalyzed the decarboxylation, transamination and racemization of several substrate analogues of amino acids. To elucidate the mechanism of interaction of substrates, especially L-serine with the enzyme, a comparative study of interaction of L-serine with the enzyme from sheep liver and Escherichia coli, was carried out. The heat stability of both the enzymes was enhanced in the presence of serine, although to different extents. Thermal denaturation monitored by spectral changes indicated an alteration in the apparent T, of sheep liver and E. coli SHMTs from 55 +/- 1 degrees C to 72 +/- 3 degrees C at 40 mM serine and from 67 +/- 1 degrees C to 72 +/- 1 degrees C at 20 mM serine, respectively. Using stopped flow spectrophotometry k values of (49 +/- 5)(.)10(-3) s(-1) and (69 +/- 7).10(-3) s(-1) for sheep liver and E. coli enzymes were determined at 50 mM serine. The binding of serine monitored by intrinsic fluorescence and sedimentation velocity measurements indicated that there was no generalized change in the structure of both proteins. However, visible CD measurements indicated a change in the asymmetric environment of pyridoxal 5'-phosphate at the active site upon binding of serine to both the enzymes. The formation of an external aldimine was accompanied by a change in the secondary structure of the enzymes monitored by far UV-CD spectra. Titration microcalorimetric studies in the presence of serine (8 mM) also demonstrated a single class of binding and the conformational changes accompanying the binding of serine to the enzyme resulted in a more compact structure leading to increased thermal stability of the enzyme.

Phase relations in the system Cu─Ho─O and stability of Cu2Ho2O5

The phase relations in the system Cu-Ho-O have been determined at 1300 K using X-ray diffraction, optical microscopy, and electron microprobe analysis of samples equilibrated in evacuated quartz ampules and in pure oxygen. Only one ternary compound, Cu2Ho2O5, was found to be stable. The Gibbs free energy of formation of this compound has been measured using the solid-state cell Pt,Cu2O + Cu2Ho2O5 + Ho2O3/(Y2O3)ZrO2/CuO + Cu2O,Pt in the temperature range of 973 to 1350 K. For the formation of Cu2Ho2O5 from its binary component oxides, 2CuO(s) + Ho2O3(S) --> Cu2Ho2O5(s) DELTAG-degrees = 11190 - 13.8T(+/- 120) J-mol-1 Since the formation is endothermic, CU2Ho2O5 becomes thermodynamically unstable with respect to CuO and Ho2O3 below 810 K. When the oxygen partial pressure over Cu2Ho2O5 is lowered, it decomposes according to the reaction 2Cu2Ho2O5(s) --> 2Ho2O3(s) + 2Cu2O(S) + O2(g) for which the equilibrium oxygen potential is given by DELTAmu(O2) = - 238510 + 160.2T(+/- 450) J.mol-1 The decomposition temperature at an oxygen partial pressure of 1.52 x 10(4) Pa was measured using a combined DTA-TGA apparatus. Based on these results, an oxygen potential diagram for the system Cu-Ho-O at 1300 K is presented.

Stability of copper tolerance in Thiobacillus ferrooxidans

A strain of Thiobacillus ferrooxidans MAL-4-1 was adapted to grow at higher concentrations of copper by repeated subculturing in the presence of increasing levels of added cupric ions in 9K medium. The strains adapted to copper were found to be more efficient in bioleaching of copper from concentrates. When copper tolerant strains were back cultured repeatedly in 9K medium without cupric ions, the initially developed metal tolerance was observed to be lost. This indicates that the copper tolerance developed is stress-dependent and not a permanent trait of the adapted strain.

Structure and stability of human telomeric sequence

Telomeric DNA of a variety of vertebrates including humans contains the tandem repeat d(TTAGGG)(n). We have investigated the structural properties of the human telomeric repeat oligonucleotide models d(T(2)AG(3))(4), d(G(3)T(2)A)(3)G(3), and d(G(3)T(2)AG(3)) using CD, gel electrophoresis, and chemical probing techniques. The sequences d(G(3)T(2)A)(3)G(3) and d(T(2)AG(3))(4) assume an antiparallel G quartet structure by intramolecular folding, while the sequence d(G(3)T(2)AG(3)) also adopts an antiparallel G quartet structure but by dimerization of hairpins. In all the above cases, adenines are in the loop. The TTA loops are oriented at the same end of the G tetrad stem in the case of hairpin dimer. Further, the oligonucleotide D(G(3)T(2)AG(3)) forms a higher order structure by the association of two hairpin dimers via stacking of G tetrad planes. Here we show that N-7 of adenine in the hairpin dimer is Hoogsteen hydrogen-bonded. The partial reactivity of loop adenines with DEPC in d(T(2)AG(3))(4) suggests that the intramolecular G quartet structure is highly polymorphic and structures with different loop orientations and topologies are formed in solution. Intra- and interloop hydrogen bonding schemes for the TTA loops are proposed to account for the observed diethyl pyrocarbonate reactivities of adenines. Sodium-induced G quartet structures differ from their potassium-induced counterparts not only in stability but also in loop conformation and interactions. Thus, the overall structure and stability of telomeric sequences are modulated by the cation present, loop sequence, and the number of G tracts, which might be important for the telomere function.

Applicability of Butler's equation in interpreting the thermodynamic behavior of surfaces and adsorption in Fe-S-O melts

Expressions for various second-order derivatives of surface tension with respect to composition at infinite dilution in terms of the interaction parameters of the surface and those of the bulk phases of dilute ternary melts have been presented. A method of deducing the parameters, which consists of repeated differentiation of Butler's equations with subsequent application of the appropriate boundary conditions, has been developed. The present investigation calculates the surface tension and adsorption functions of the Fe-S-O melts at 1873 and 1923 K using the modified form of Butler's equations and the derived values for the surface interaction parameters of the system. The calculated values are found to be in good agreement with those of the experimental data of the system. The present analysis indicates that the energetics of the surface phase are considerably different from those of the bulk phase. The present research investigates a critical compositional range beyond which the surface tension increases with temperature. The observed increase in adsorption of sulfur with consequent desorption of oxygen as a function of temperature above the critical compositional range has been ascribed to the increase of activity ratios of oxygen to sulfur in the surface relative to those in the bulk phase of the system.

A thermodynamic criterion for selection of gas compositions for diamond deposition

Isoactivity lines for carbon with respect to diamond as the standard state have been calculated in the ternary system C-H-O at 1223 K to identify the diamond deposition domain. The gas composition is calculated by suppressing the formation of all condensed forms of carbon using the SOLGASMIX free-energy minimization program. Thirty six gas species were included in the calculation. From the gas composition, isoactivity lines are computed using recent data on the Gibbs energy of diamond. Except for activities less than 0.1, the isoactivity lines are almost linear on the C-H-O ternary diagram. Gas compositions which generate activity of diamond ranging from 1 to 100 at 1223 K fall inside a narrow wedge originating from the point representing CO. This wedge is very similar to the revised lens-shaped diamond growth domain identified by Bachman et al., using inputs from experiment. The small difference between the calculated and observed domains may be attributed to variation in the supersaturation required for diamond deposition with gas composition. The diamond solubility in the gas phase along the isoactivity line for a(di)=100 and P=6.7 kPa exhibits a minimum at 1280 K, which is close to the optimum temperature found experimentally. At higher supersaturations, non-diamond forms of carbon, including amorphous varieties, are expected. The results suggest that thermodynamic calculations can be useful for locating diamond growth domains in more complex CVD systems containing halogens, for which very little experimental data is available.

Integral treatment for the representation of thermodynamic properties in multicomponent systems using interaction parameters

The partial thermodynamic functions of the solvent component of a ternary system have been deduced in terms of the interaction parameters by integration of several series which emerge from the Maclaurin infinite series based on the integral property of the system and subjected to appropriate boundary conditions. The series integration shows that the resulting partial functions are suitable for interpreting the thermodynamic properties of the system and are independent of compositional paths. In the present analysis, the higher order terms of these series are found to make insignificant contributions.

Flap-lag-torsion stability in hover and forward flight with three-dimensional wake

The contributions of full-wake dynamics in trim analysis are demonstrated for finding the control inputs and periodic responses simultaneously, as well as in Floquet eigenanalysis for finding the damping levels. The equations of flap bending, lag bending, and torsion are coupled with a three-dimensional, finite state wake, and low-frequency (<1/rev) to high frequency (>1/rev) multiblade modes are considered. Full blade-wake dynamics is used in trim analysis and Floquet eigenanalysis. A uniform cantilever blade in trimmed flight is investigated over a range of thrust levels, advance ratios, number of blades, and blade torsional frequencies. The investigation includes the convergence characteristics of control inputs, periodic responses, and damping levels with respect to the number of spatial azimuthal harmonics and radial shape functions in the wake representation. It also includes correlation with the measured lag damping of a three-bladed untrimmed rotor. The parametric study shows the dominant influence of wake dynamics on control inputs, periodic responses, and damping levels, and wake theory generally improves the correlation.

Phase stability and dielectric properties of ceramics in the Pb(Zn1/3Nb2/3)O3-Pb(Ni1/3Nb2/3)O3-PbTio3 relaxor ferroelectric system

The perovskite structure in Pb(Zn1/3Nb2/3)O3 can be stabilized by the addition of Pb(Ni1/3Nb2/3)O3 and PbTiO3.Pb(Ni1/3Nb2/3)O3 assists in lowering the sintering temperature and shifting the Curie temperature of ceramics while PbTiO3 helps to optimize the dielectric properties. The phase stability and dielectric properties of several compositions in the Pb(Zn1/3Nb2/3)O3-Pb(Ni1/3Nb2/3)O3-PbTiO3 ternary relaxor ferroelectric system were investigated for possible capacitor applications. The effect of calcining and sintering temperature on the stability of perovskite phase in PZN rich compositions was studied extensively as a function of composition. The boundary line separating perovskite and mixed phases was determined for compositions near PZN. Several compositions can be sintered below 1050°C. The dielectric properties of compositions near the mixed phase boundary showed strong dependence on the percentage of pyrochlore phase. Compositions with a dielectric constant of 12.500 at room temperature have been identified which meet Z5T and Y5U specifications for dielectric constant and tan δ.

Interaction of sheep liver apo-serine hydroxymethyltransferase with pyridoxal-5'-phosphate: A physicochemical, kinetic, and thermodynamic study

Sheep liver serine hydroxymethyltransferase (EC 2.1.2.1) is a homotetramer of M(r) 213,000 requiring pyridoxal-5'-phosphate (PLP) as cofactor, Removal of PLP from the holoenzyme converted the enzyme to the apo form which, in addition to being inactive, was devoid of the characteristic absorption spectrum. Upon the addition of PLP to the apoenzyme, complete activity was restored and the visible absorption spectrum with a maximum at 425 nm was regained. The interaction of PLP with the apoenzyme revealed two phases of reaction with pseudo-first-order rate constants of 20 +/- 5 s(-1) and 12.2 +/- 2.0 x 10(-3) s(-1), respectively. However, addition of PLP to the apoenzyme did not cause gross conformational changes as evidenced by circular dichroic and fluorescence spectroscopy. Although conformationally apoenzyme and holoenzyme were indistinguishable, they had distinct apparent melting temperatures of 51 +/- 2 and 58 +/- 2 degrees C, respectively, and the reconstituted holoenzyme was thermally as stable as the native holoenzyme. These results suggested that there was no apparent difference in the secondary structure of holoenzyme, apoenzyme, and reconstituted holoenzyme, However, sedimentation analysis of the apoenzyme revealed the presence of two peaks of S-20,S-w values of 8.7 +/- 0.5 and 5.7 +/- 0.3 S, respectively. A similar pattern was observed when the apoenzyme was chromatographed on a calibrated Sephadex G-150 column. The first peak corresponded to the tetrameric form (M(r) 200,000 +/- 15,000) while the second peak had a M(r) of 130,000 +/- 10,000. Reconstitution experiments revealed that only the tetrameric form of the apoenzyme could be converted into an active holoenzyme while the dimeric form could not be reconstituted into an active enzyme. These results demonstrate that PLP plays an important role in maintaining the structural integrity of the enzyme by preventing the dissociation of the enzyme into subunits, in addition to its function in catalysis. (C) 1996 Academic Press, Inc.

Parallel computing concepts and methods for floquet analysis of helicopter trim and stability

Floquet analysis is widely used for small-order systems (say, order M < 100) to find trim results of control inputs and periodic responses, and stability results of damping levels and frequencies, Presently, however, it is practical neither for design applications nor for comprehensive analysis models that lead to large systems (M > 100); the run time on a sequential computer is simply prohibitive, Accordingly, a massively parallel Floquet analysis is developed with emphasis on large systems, and it is implemented on two SIMD or single-instruction, multiple-data computers with 4096 and 8192 processors, The focus of this development is a parallel shooting method with damped Newton iteration to generate trim results; the Floquet transition matrix (FTM) comes out as a byproduct, The eigenvalues and eigenvectors of the FTM are computed by a parallel QR method, and thereby stability results are generated, For illustration, flap and flap-lag stability of isolated rotors are treated by the parallel analysis and by a corresponding sequential analysis with the conventional shooting and QR methods; linear quasisteady airfoil aerodynamics and a finite-state three-dimensional wake model are used, Computational reliability is quantified by the condition numbers of the Jacobian matrices in Newton iteration, the condition numbers of the eigenvalues and the residual errors of the eigenpairs, and reliability figures are comparable in both the parallel and sequential analyses, Compared to the sequential analysis, the parallel analysis reduces the run time of large systems dramatically, and the reduction increases with increasing system order; this finding offers considerable promise for design and comprehensive-analysis applications.