Crystal Structure of Staphylococcus aureus Metallopeptidase (Sapep) Reveals Large Domain Motions between the Manganese-bound and Apo-states

Proteases belonging to the M20 family are characterized by diverse substrate specificity and participate in several metabolic pathways. The Staphylococcus aureus metallopeptidase, Sapep, is a member of the aminoacylase-I/M20 protein family. This protein is a Mn2+-dependent dipeptidase. The crystal structure of this protein in the Mn2+-bound form and in the open, metal-free state suggests that large interdomain movements could potentially regulate the activity of this enzyme. We note that the extended inactive conformation is stabilized by a disulfide bond in the vicinity of the active site. Although these cysteines, Cys(155) and Cys(178), are not active site residues, the reduced form of this enzyme is substantially more active as a dipeptidase. These findings acquire further relevance given a recent observation that this enzyme is only active in methicillin-resistant S. aureus. The structural and biochemical features of this enzyme provide a template for the design of novel methicillin-resistant S. aureus-specific therapeutics.

Dynamic coupling between the LID and NMP domain motions in the catalytic conversion of ATP and AMP to ADP by adenylate kinase

The catalytic conversion of adenosine triphosphate (ATP) and adenosine monophosphate (AMP) to adenosine diphosphate (ADP) by adenylate kinase (ADK) involves large amplitude, ligand induced domain motions, involving the opening and the closing of ATP binding domain (LID) and AMP binding domain (NMP) domains, during the repeated catalytic cycle. We discover and analyze an interesting dynamical coupling between the motion of the two domains during the opening, using large scale atomistic molecular dynamics trajectory analysis, covariance analysis, and multidimensional free energy calculations with explicit water. Initially, the LID domain must open by a certain amount before the NMP domain can begin to open. Dynamical correlation map shows interesting cross-peak between LID and NMP domain which suggests the presence of correlated motion between them. This is also reflected in our calculated two-dimensional free energy surface contour diagram which has an interesting elliptic shape, revealing a strong correlation between the opening of the LID domain and that of the NMP domain. Our free energy surface of the LID domain motion is rugged due to interaction with water and the signature of ruggedness is evident in the observed root mean square deviation variation and its fluctuation time correlation functions. We develop a correlated dynamical disorder-type theoretical model to explain the observed dynamic coupling between the motion of the two domains in ADK. Our model correctly reproduces several features of the cross-correlation observed in simulations. (C) 2011 American Institute of Physics. doi:10.1063/1.3516588]

Understanding the role of domain-domain linkers in the spatial orientation of domains in multi-domain proteins

Inter-domain linkers (IDLs)' bridge flanking domains and support inter-domain communication in multi-domain proteins. Their sequence and conformational preferences enable them to carry out varied functions. They also provide sufficient flexibility to facilitate domain motions and, in conjunction with the interacting interfaces, they also regulate the inter-domain geometry (IDG). In spite of the basic intuitive understanding of the inter-domain orientations with respect to linker conformations and interfaces, we still do not entirely understand the precise relationship among the three. We show that IDG is evolutionarily well conserved and is constrained by the domain-domain interface interactions. The IDLs modulate the interactions by varying their lengths, conformations and local structure, thereby affecting the overall IDG. Results of our analysis provide guidelines in modelling of multi-domain proteins from the tertiary structures of constituent domain components.

Role of Water in the Enzymatic Catalysis: Study of ATP plus AMP -> 2ADP Conversion by Adenylate Kinase

The catalytic conversion ATP + AMP -> 2ADP by the enzyme adenylate kinase (ADK) involves the binding of one ATP. molecule to the LID domain and one AMP molecule to the NMP domain. The latter is followed by a. phosphate transfer and then the release of two ADP molecules. We have computed a novel two-dimensional configurational free energy surface (2DCFES), with one reaction coordinate each for the LID and the NMP domain motions, while considering explicit water interactions. Our computed 2DCFES clearly reveals the existence of a stable half-open half-closed (HOHC) intermediate stale of the enzyme. Cycling of the enzyme through the HOHC state reduces the conformational free energy barrier for. the reaction by about 20 kJ/mol. We find that the stability of the HOHC state (missed in all earlier studies with implicit solvent model) is largely because of the increase of specific interactions of the polar amino acid side chains with water, particularly with the arginine and the histidine residues. Free energy surface of the LID domain is rather rugged, which can conveniently slow down LID's conformational motion, thus facilitating a new substrate capture after the product release in the catalytic cycle.

Crystallographic Structure of Phosphofructokinase-2 from Escherichia coli in Complex with Two ATP Molecules. Implications for Substrate Inhibition

Phosphofructokinase-1 and -2 (Pfk-1 and Pfk-2, respectively) from Escherichia coli belong to different homologous superfamilies. However, in spite of the lack of a common ancestor, they share the ability to catalyze the same reaction and are inhibited by the substrate MgATP. Pfk-2, an ATP-dependent 6-phosphofructokinase member of the ribokinase-like superfamily, is a homodimer of 66 kDa subunits whose oligomerization state is necessary for catalysis and stability. The presence of MgATP favors the tetrameric form of the enzyme. In this work, we describe the structure of Pfk-2 in its inhibited tetrameric form, with each subunit bound to two ATP molecules and two Mg ions. The present structure indicates that substrate inhibition occurs due to the sequential binding of two MgATP molecules per subunit, the first at the usual site occupied by the nucleotide in homologous enzymes and the second at the allosteric site, making a number of direct and Mg-mediated interactions with the first. Two configurations are observed for the second MgATP, one of which involves interactions with Tyr23 from the adjacent subunit in the dimer and the other making an unusual non-Watson-Crick base pairing with the adenine in the substrate ATP. The oligomeric state observed in the crystal is tetrameric, and some of the structural elements involved in the binding of the Substrate and allosteric ATPs are also participating in the dimer-dimer interface. This structure also provides the grounds to compare analogous features of the nonhomologous phosphofructokinases from E. coli. (C) 2008 Elsevier Ltd. All rights reserved.

The Crystal Complex of Phosphofructokinase-2 of Escherichia coli with Fructose-6-phosphate KINETIC AND STRUCTURAL ANALYSIS OF THE ALLOSTERIC ATP INHIBITION

Substrate inhibition by ATP is a regulatory feature of the phosphofructokinases isoenzymes from Escherichia coli (Pfk-1 and Pfk-2). Under gluconeogenic conditions, the loss of this regulation in Pfk-2 causes substrate cycling of fructose-6-phosphate (fructose-6-P) and futile consumption of ATP delaying growth. In the present work, we have broached the mechanism of ATP-induced inhibition of Pfk-2 from both structural and kinetic perspectives. The crystal structure of Pfk-2 in complex with fructose-6-P is reported to a resolution of 2 angstrom. The comparison of this structure with the previously reported inhibited form of the enzyme suggests a negative interplay between fructose-6-P binding and allosteric binding of MgATP. Initial velocity experiments show a linear increase of the apparent K(0.5) for fructose-6-P and a decrease in the apparent k(cat) as a function of MgATP concentration. These effects occur simultaneously with the induction of a sigmoidal kinetic behavior (n(H) of approximately 2). Differences and resemblances in the patterns of fructose-6-P binding and the mechanism of inhibition are discussed for Pfk-1 and Pfk-2, as an example of evolutionary convergence, because these enzymes do not share a common ancestor.

Structural role of the active-site metal in the conformation of Trypanosoma brucei phosphoglycerate mutase

Phosphoglycerate mutases (PGAMs) participate in both the glycolytic and the gluconeogenic pathways in reversible isomerization of 3-phosphoglycerate and 2-phosphoglycerate. PGAMs are members of two distinct protein families: enzymes that are dependent on or independent of the 2,3-bisphosphoglycerate cofactor. We determined the X-ray structure of the monomeric Trypanosoma brucei independent PGAM (TbiPGAM) in its apoenzyme form, and confirmed this observation by small angle X-ray scattering data. Comparing the TbiPGAM structure with the Leishmania mexicana independent PGAM structure, previously reported with a phosphoglycerate molecule bound to the active site, revealed the domain movement resulting from active site occupation. The structure reported here shows the interaction between Asp319 and the metal bound to the active site, and its contribution to the domain movement. Substitution of the metal-binding residue Asp319 by Ala resulted in complete loss of independent PGAM activity, and showed for the first time its involvement in the enzymes function. As TbiPGAM is an attractive molecular target for drug development, the apoenzyme conformation described here provides opportunities for its use in structure-based drug design approaches.

Time-domain simulation of wave-induced ship motions by a Rankine panel method

In this thesis, a numerical program has been developed to simulate the wave-induced ship motions in the time domain. Wave-body interactions have been studied for various ships and floating bodies through forced motion and free motion simulations in a wide range of wave frequencies. A three-dimensional Rankine panel method is applied to solve the boundary value problem for the wave-body interactions. The velocity potentials and normal velocities on the boundaries are obtained in the time domain by solving the mixed boundary integral equations in relation to the source and dipole distributions. The hydrodynamic forces are calculated by the integration of the instantaneous hydrodynamic pressures over the body surface. The equations of ship motion are solved simultaneously with the boundary value problem for each time step. The wave elevation is computed by applying the linear free surface conditions. A numerical damping zone is adopted to absorb the outgoing waves in order to satisfy the radiation condition for the truncated free surface. A numerical filter is applied on the free surface for the smoothing of the wave elevation. Good convergence has been reached for both forced motion simulations and free motion simulations. The computed added-mass and damping coefficients, wave exciting forces, and motion responses for ships and floating bodies are in good agreement with the numerical results from other programs and experimental data.

Dynamic state estimation for identifying earthquake support motions in instrumented structures

The problem of identification of multi-component and (or) spatially varying earthquake support motions based on measured responses in instrumented structures is considered. The governing equations of motion are cast in the state space form and a time domain solution to the input identification problem is developed based on the Kalman and particle filtering methods. The method allows for noise in measured responses, imperfections in mathematical model for the structure, and possible nonlinear behavior of the structure. The unknown support motions are treated as hypothetical additional system states and a prior model for these motions are taken to be given in terms of white noise processes. For linear systems, the solution is developed within the Kalman filtering framework while, for nonlinear systems, the Monte Carlo simulation based particle filtering tools are employed. In the latter case, the question of controlling sampling variance based on the idea of Rao-Blackwellization is also explored. Illustrative examples include identification of multi-component and spatially varying support motions in linear/nonlinear structures.

Phase mapping of domain kinetics in lithium niobate by digital holographic interferometry

The phase mapping of domain kinetics under the uniform steady-state electric field is achieved and investigated in the LiNbO3 crystals by digital holographic interferometry. We obtained the sequences of reconstructed three-dimensional and two-dimensional wave-field phase distributions during the electric poling in the congruent and near stoichiometric LiNbO3 crystals. The phase mapping of individual domain nucleation and growth in the two crystals are obtained. It is found that both longitudinal and lateral domain growths are not linear during the electric poling. The phase mapping of domain wall motions in the two crystals is also obtained. Both the phase relaxation and the pinning-depinning mechanism are observed during the domain wall motion. The residual phase distribution is observed after the high-speed domain wall motion. The corresponding analyses and discussions are proposed to explain the phenomena.

Tracking magnetic bright point motions through the solar atmosphere

High-cadence, multiwavelength observations and simulations are employed for the analysis of solar photospheric magnetic bright points (MBPs) in the quiet Sun. The observations were obtained with the Rapid Oscillations in the Solar Atmosphere (ROSA) imager and the Interferometric Bidimensional Spectrometer at the Dunn Solar Telescope. Our analysis reveals that photospheric MBPs have an average transverse velocity of approximately 1 km s-1, whereas their chromospheric counterparts have a slightly higher average velocity of 1.4 km s-1. Additionally, chromospheric MBPs were found to be around 63 per cent larger than the equivalent photospheric MBPs. These velocity values were compared with the output of numerical simulations generated using the muram code. The simulated results were similar, but slightly elevated, when compared to the observed data. An average velocity of 1.3 km s-1 was found in the simulated G-band images and an average of 1.8 km s-1 seen in the velocity domain at a height of 500 km above the continuum formation layer. Delays in the change of velocities were also analysed. Average delays of ˜4 s between layers of the simulated data set were established and values of ˜29 s observed between G-band and Ca ii K ROSA observations. The delays in the simulations are likely to be the result of oblique granular shock waves, whereas those found in the observations are possibly the result of a semi-rigid flux tube.