The hinge-bending mode of a lysozyme-inhibitor complex

The hinge-bending mode of hen egg white lysozyme is studied by a constrained minimization technique. Results with and without a bound inhibitor, tri-N-acetyl-glucosamine, are obtained. The frequency of the mode with the inhibitor is found to be 4.3 cm^?1, in contrast to 3.0 cm^?1 for the free enzyme. Also, the hinge-bending angle with the lowest energy is shifted 10° towards a more closed cleft in the bound species. The main contribution to these differences arise from interactions with the residues lining the cleft and those on the back side of it. Structural details that account for the energetics are presented. The method of calculation is somewhat different from a previous study [J. A. McCammon, B. R. Gelin, M. Karplus & P. G. Wolynes, (1976) Nature 262 , 325–326] to reduce the likelihood of artifacts in the results.     

Domain motions in bacteriophage T4 lysozyme: A comparison between molecular dynamics and crystallographic data

A comparison of a series of extended molecular dynamics (MD) simulations of bacteriophage T4 lysozyme in solvent with X-ray data is presented. Essential dynamics analyses were used to derive collective fluctuations from both the simulated trajectories and a distribution of crystallographic conformations. In both cases the main collective fluctuations describe domain motions. The protein consists of an N- and C-terminal domain connected by a long helix. The analysis of the distribution of crystallographic conformations reveals that the N-terminal helix rotates together with either of these two domains. The main domain fluctuation describes a closure mode of the two domains in which the N-terminal helix rotates concertedly with the C-terminal domain, while the domain fluctuation with second largest amplitude corresponds to a twisting mode of the two domains, with the N-terminal helix rotating concertedly with the N-terminal domain. For the closure mode, the difference in hinge-bending angle between the most open and most closed X-ray structure along this mode is 49 degrees. In the MD simulation that shows the largest fluctuation along this mode, a rotation of 45 degrees was observed. Although the twisting mode has much less freedom than the closure mode in the distribution of crystallographic conformations, experimental results suggest that it might be functionally important. Interestingly, the twisting mode is sampled more extensively in all MD simulations than it is in the distribution of X-ray conformations. Proteins 31:116–127, 1998. © 1998 Wiley-Liss, Inc.     

The essential dynamics of thermolysin: Confirmation of the hinge-bending motion and comparison of simulations in vacuum and water

Comparisons of the crystal structures of thermolysin and the thermolysin-like protease produced by B. cereus have recently led to the hypothesis that neutral proteases undergo a hinge-bending motion. We have investigated this hypothesis by analyzing molecular dynamics simulations of thermolysin in vacuum and water, using the essential dynamics method. This method is able to extract large concerted atomic motions of biological importance from a molecular dynamics trajectory. The analysis of the thermolysin trajectories indeed revealed a large rigid body hinge-bending motion of the Nterminal and C-terminal domains, similar to the motion hypothesized from the crystal structure comparisons. In addition, it appeared that the essential dynamics properties derived from the vacuum simulation were similar to those derived from the solvent simulation. © 1995 Wiley-Liss, Inc.     

Site-Specific Structural Variations Accompanying Tubular Assembly of the HIV-1 Capsid Protein

The 231-residue capsid (CA) protein of human immunodeficiency virus type 1 (HIV-1) spontaneously self-assembles into tubes with a hexagonal lattice that is believed to mimic the surface lattice of conical capsid cores within intact virions. We report the results of solid-state nuclear magnetic resonance (NMR) measurements on HIV-1 CA tubes that provide new information regarding changes in molecular structure that accompany CA self-assembly, local dynamics within CA tubes, and possible mechanisms for the generation of lattice curvature. This information is contained in site-specific assignments of signals in two- and three-dimensional solid-state NMR spectra, conformation-dependent ¹⁵N and ¹³C NMR chemical shifts, detection of highly dynamic residues under solution NMR conditions, measurements of local variations in transverse spin relaxation rates of amide ¹H nuclei, and quantitative measurements of site-specific ¹⁵N–¹⁵N dipole–dipole couplings. Our data show that most of the CA sequence is conformationally ordered and relatively rigid in tubular assemblies and that structures of the N-terminal domain (NTD) and the C-terminal domain (CTD) observed in solution are largely retained. However, specific segments, including the N-terminal β-hairpin, the cyclophilin A binding loop, the inter-domain linker, segments involved in intermolecular NTD–CTD interactions, and the C-terminal tail, have substantial static or dynamical disorder in tubular assemblies. Other segments, including the 3₁₀-helical segment in CTD, undergo clear conformational changes. Structural variations associated with curvature of the CA lattice appear to be localized in the inter-domain linker and intermolecular NTD–CTD interface, while structural variations within NTD hexamers, around local 3-fold symmetry axes, and in CTD–CTD dimerization interfaces are less significant.     

Direct Measurement of Hydration-Related Dynamic Changes in Lysozyme using Inelastic Neutron Scattering Spectroscopy

Abstract

Inelastic neutron scattering spectroscopy is used to investigate dynamic changes in lysozyme powder at two different low D,0 hydrations (0.07g D²,O/g protein and 0,20 g D²,O/g protein). In the higher hydration sample, the inelastic scattering between 0.8 and 4.0 cm^?1 energy transfer is increased and the elastic scattering is decreased. The decreased elastic scattering suggests increased atomic amplitudes of motion and the increased 0.8 to 4.0 cm^?1 scattering suggests increased motions in this frequency range. Comparison with normal mode models of lysozyme dynamics shows that the inelastic difference occurs in the frequency region predicted for the lowest frequency, largest amplitude, global modes of the molccule[M. Levitt, C. Sanderand P. S. Stern, J. Mol. Biol. 181. 423 (1985). B.Brooks and M.Karplus.Prot. Natl Acad. Sci (U.S.A) 82. 4995 (1985), R.E. Bruccoleri, M. Karplus and J.A. McCammon, Biopolymers 25 1767 (1986)]. Our results are consistent with a model in which an increased number of low frequency global modes are present in the higher hydrated sample.     

E. coli 5'-nucleotidase undergoes a hinge-bending domain rotation resembling a ball-and-socket motion

Structures of nine independent conformers of E. coli 5'-nucleotidase (5'-NT) have been analyzed using four different crystal forms. These data show that the two-domain protein undergoes an unusual 96 degrees hinge-bending domain rotation. Structures of the open and closed forms with substrates and inhibitors reveal that the substrate moves by approximately 25 A with the large domain rotation into the catalytic site. The domain motions derived from a comparison of the nine conformations agree well with motions obtained from a normal mode analysis in that all independent domain rotations are around axes that are roughly located in the plane which includes the domain centers and the hinge. Two residues, Lys355 and Gly356, form the core of the hinge region and undergo a large change of the main-chain torsion angles. The hinge-bending movement observed for 5'-nucleotidase differs markedly from a classical hinge-bending closure motion which involves an opening of the substrate or ligand-binding cleft between two domains. In contrast, the movement observed in 5'-nucleotidase resembles that of a ball-and-socket joint. The smaller C-terminal domain rotates approximately around its center such that the residues at the domain interface move in a sliding motion along the interface. Few direct interdomain contacts and a layer of water molecules between the two domains facilitate the sliding motion.     

An Effective Enzyme Interacting with Poly (dT-dA)· Poly (dT-dA): A Dynamic Enhancer-Repressor Action

Abstract

The Green function technique is used to study the open hydrogen bond probability of poly(dT-dA)·poly(dT-dA) when an effective enzyme is attached to the helix. The DNA interstrand hydrogen bond mean motion and probability of fluctuating to an open state depends on the internal vibrational frequency of the enzyme. An enzyme with internal frequency of 80 cm ^?1 reduces hydrogen bond motion and the resulting probability of hydrogen bond fluctuational opening. An enzyme with internal frequency of 72 cm ^?1 increases hydrogen bond motion and the probability of hydrogen bond breaking.     

Phosphorylation triggers domain separation in the DNA binding response regulator NarL

DNA binding proteins of two-component signal transduction systems in microorganisms are activated by phosphorylation through an unknown mechanism. NarL is an example from the nitrate/nitrite signal transduction system of Escherichia coli. NarL consists of N- and C-terminal domains, the latter of which contains the DNA binding elements. To explore the mechanism of activation, single nitroxide side chains were introduced, one at a time, at nine different sites throughout the C-terminal domain to monitor the tertiary structure and the status of the surface in contact with the N-terminal domain. In addition, three pairs of doubly labeled proteins were prepared to monitor the interdomain distance using the magnetic dipolar interaction. The results of these site-directed spin-labeling studies reveal that phosphorylation at a distant site in the N-terminal domain triggers domain separation, likely by a hinge-bending motion. This in turn presents key elements of the C-terminal domain for docking to the DNA target in the configuration described in the recent crystal structure. The data also imply that a single conformation of unphosphorylated NarL exists in solution, and there is no detectable equilibrium between the closed and open conformations.     

Investigation of Domain Motions in Bacteriophage T4 Lysozyme

Abstract

Hinge-bending in T4 lysozyme has been inferred from single amino acid mutant crystalline allomorphs by Matthews and coworkers. This raises an important question: are the different conformers in the unit cell artifacts of crystal packing forces, or do they represent different solution state structures? The objective of this theoretical study is to determine whether domain motions and hinge-bending could be simulated in T4 lysozyme using molecular dynamics. An analysis of a 400 ps molecular dynamics simulation of the 164 amino acid enzyme T4 lysozyme is presented. Molecular dynamics calculations were computed using the Discover software package (Biosym Technologies). All hydrogen atoms were modeled explicitly with the inclusion of all 152 crystallographic waters at a temperature of 300 K. The native T4 lysozyme molecular dynamics simulation demonstrated hinge-bending in the protein. Relative domain motions between the N-terminal and C-terminal domains were evident. The enzyme hinge bending sites resulted from small changes in backbone atom conformations over several residues rather than rotation about a single bound. Two hinge loci were found in the simulation. One locus comprises residues 8–14 near the C-terminal of the A helix; the other site, residues 77–83 near the C-terminal of the C helix. Comparison of several snapshot structures from the dynamics trajectory clearly illustrates domain motions between the two lysozyme lobes. Time correlated atomic motions in the protein were analyzed using a dynamical cross-correlation map. We found a high degree of correlated atomic motions in each of the domains and, to a lesser extent, anticorrelated motions between the two domains. We also found that the hairpin loop in the N-terminal lobe (residues 19–24) acted as a mobile ‘flap’ and exhibited highly correlated dynamic motions across the cleft of the active site, especially with residue 142.     

Protein Dynamics of the Sensor Protein HemAT as Probed by Time-Resolved Step-Scan FTIR Spectroscopy

The heme-based aerotactic transducer (HemAT) is an oxygen-sensor protein consisting of a sensor and a signaling domain in the N- and C-terminal regions, respectively. Time-resolved step-scan FTIR spectroscopy was employed to characterize protein intermediate states obtained by photolysis of the carbon monoxide complexes of sensor-domain, full-length HemAT, and the Y70F (B-helix), L92A (E-helix), T95A (E-helix), and Y133F (G-helix) HemAT mutants. We assign the spectral components to discrete substructures, which originate from a helical structure that is solvated (1638 cm^?1) and a native helix that is protected from solvation by interhelix tertiary interactions (1654 cm^?1). The full-length protein is characterized by an additional amide I absorbance at 1661 cm^?1, which is attributed to disordered structure suggesting that further protein conformational changes occur in the presence of the signaling domain in the full-length protein. The kinetics monitored within the amide I absorbance of the polypeptide backbone in the sensor domain exhibit two distinct relaxation phases (t₁ = 24 and t₂ = 694 μs), whereas that of the full-length protein exhibits monophasic behavior for all substructures in a time range of t = 1253–2090 μs. These observations can be instrumental in monitoring helix motion and the role of specific mutants in controlling the dynamics in the communication pathway from the sensor to the signaling domain. The kinetics observed for the amide I relaxation for the full-length protein indicate that the discrete substructures within full-length HemAT, unlike those of the sensor domain, relax independently.     

Quaternary structures and low frequency molecular vibrations of haems of deoxy and oxyhaemoglobin studied by resonance Raman scattering

The influence of quaternary structure on the low frequency molecular vibrations of the haem within deoxyhaemoglobin (deoxy Hb) and Oxyhaemoglobin (oxy Hb) was studied by resonance Raman scattering. The FeO₂ stretching frequency was essentially identical between the high affinity (R) state (Hb A) and low affinity (T) state (Hb Kansas and Hb M Milwaukee with inositol hexaphosphate). However in deoxy Hb, only one of the polarized lines showed an appreciable frequency shift upon switch of quaternary structure, i.e. 215 to 218 cm^?1 for the T state (Hb A, des-His(146β) Hb, and des-Arg(141α) Hb (pH 6.5)) and 220 to 221 cm^?1 for the R state (des-Arg(141α) Hb (pH 9.0), des-His(146β)-Arg(141α) Hb and NES des-Arg(141α) Hb). Based on the observed ⁵⁴Fe isotopic frequency shift of the corresponding Raman lines of deoxy Hb A (214 → 217 cm^?1), of deoxy NES des-Arg Hb (220 → 223 cm^?1), of the protoporphyrinato-Fe(II)-(2-methylimidazole) complex in the ferrous high spin state (207 → 211 cm^?1) and of deoxymyoglobin (220 → 222 cm^?1) (Kitagawa et al., 1979), and on substitution of perdeuterated for protonated 2-methylimidazole in the deoxygenated picket fence complex (T_pivPP)Fe²⁺ (2-MeIm) (209 → 206 cm^?1), and on the results of normal co-ordinates calculation carried out previously, we proposed that the 216 cm^?1 line of deoxy Hb is associated primarily with the FeN_ε(HisF8) stretching mode and accordingly that the FeN_ε(HisF8) bond is stretched in the T state due to a strain exerted by globin.     

Raman spectroscopic study on low-frequency collective modes in self-associates of guanosine monophosphates

Low-frequency Raman spectra of the self-associates of guanosine monophosphates (GMPs) Na₂ · 5′GMP, K₂ · 5′GMP, Na₂ · 3′GMP, and K₂ · 3′GMP, and polyribonucleic acid K · poly(rG), were obtained. In acidic gels and dried fibers, GMP molecules are known to form helical stacks of hydrogen-bonded tetramers. Some low-frequency collective modes specific to the helically stacked structures were observed. We examined the dependence of these modes on counterions and water content. The lowest frequency mode at ca. 20 cm^?1 is sensitive to the water content of the sample and is clearly visible in solid-state samples, so it works as a marker band of the environmental condition of the helices. The intensity and the peak frequency of this mode in solid-state samples depend on the helical structure and counterions. The broad peaks in the vicinity of 100 cm^?1 are influenced by cations and are independent of water content.     

The N-terminal Domain of the Drosophila Mitochondrial Replicative DNA Helicase Contains an Iron-Sulfur Cluster and Binds DNA

The metazoan mitochondrial DNA helicase is an integral part of the minimal mitochondrial replisome. It exhibits strong sequence homology with the bacteriophage T7 gene 4 protein primase-helicase (T7 gp4). Both proteins contain distinct N- and C-terminal domains separated by a flexible linker. The C-terminal domain catalyzes its characteristic DNA-dependent NTPase activity, and can unwind duplex DNA substrates independently of the N-terminal domain. Whereas the N-terminal domain in T7 gp4 contains a DNA primase activity, this function is lost in metazoan mtDNA helicase. Thus, although the functions of the C-terminal domain and the linker are partially understood, the role of the N-terminal region in the metazoan replicative mtDNA helicase remains elusive. Here, we show that the N-terminal domain of Drosophila melanogaster mtDNA helicase coordinates iron in a 2Fe-2S cluster that enhances protein stability in vitro. The N-terminal domain binds the cluster through conserved cysteine residues (Cys⁶⁸, Cys⁷¹, Cys¹⁰², and Cys¹⁰⁵) that are responsible for coordinating zinc in T7 gp4. Moreover, we show that the N-terminal domain binds both single- and double-stranded DNA oligomers, with an apparent K_d of ∼120 nm. These findings suggest a possible role for the N-terminal domain of metazoan mtDNA helicase in recruiting and binding DNA at the replication fork.     

Cloning and characterization of mitochondrial methionyl-tRNA synthetase from a pathogenic fungi Candida albicans

A genomic sequence encoding mitochondrial methionyl-tRNA synthetase (MetRS) was determined from a pathogenic fungi Candida albicans. The gene is distinct from that encoding the cytoplasmic MetRS. The encoded protein consists of 577 amino acids (aa) and contains the class I defining sequences in the N-terminal domain and the conserved anticodon-binding amino acid, Trp, in the C-terminal domain. This protein showed the highest similarity with the mitochondrial MetRSs of Saccharomyces cerevisiae and Shizosaccharomyces pombe. The mitochondrial MetRSs of these fungi were distinguished from their cytoplasmic forms. The protein lacks the zinc binding motif in the N-terminal domain and the C-terminal dimerization appendix that are present in MetRSs of several other species. Escherichia coli tRNA^Met was a substrate for the encoded protein as determined by genetic complementation and in vitro aminoacylation reaction. This cross-species aminoacylation activity suggests the conservation of interaction mode between tRNA^Met and MetRS.     

Examination of the secondary structure of proteins by deconvolved FTIR spectra

Fourier transform ir (FTIR) spectra of 21 globular proteins have been obtained at 2 cm^?1 resolution from 1600 to 1700 cm^?1 in deuterium oxide solution. Fourier self-deconvolution was applied to all spectra, revealing that the amide I band of each protein except casein consists of six to nine components. The components are observed at 11 well-defined frequencies, although all proteins do not exhibit components at every characteristic frequency. The root mean square (RMS) deviation of 124 individual values from the 11 average characteristic frequencies is 1.9 cm^?1. The observed components are assigned to helical segments, extended beta-segments, unordered segments, and turns. Segments with similar structures do not necessarily exhibit band components with identical frequencies. For instance, the lower frequency beta-structure band can vary within a range of approximately 15 cm^?1. The relative areas of the individual components of the deconvolved spectra were determined by a Gauss–Newton, iterative curve-fitting procedure that assumed Gaussian band envelopes for the deconvolved components. The measured areas were used to estimate the percentage of helix and beta-structure for each of 21 globular proteins. The results are in good general agreement with values derived from x-ray data by Levitt and Greer. The RMS deviation between 22 values (alpha- and beta-content of 11 beta-rich proteins measured by both techniques) is 2.5 percentage points; the maximum absolute deviation is 4 percentage points.     

VCD spectroscopic properties of the β‐hairpin forming miniprotein CLN025 in various solvents

Electronic and vibrational circular dichroism are often used to determine the secondary structure of proteins, because each secondary structure has a unique spectrum. Little is known about the vibrational circular dichroic spectroscopic features of the β‐hairpin. In this study, the VCD spectral features of a decapeptide, YYDPETGTWY (CLN025), which forms a stable β‐hairpin that is stabilized by intramolecular weakly polar interactions and hydrogen bonds were determined. Molecular dynamics simulations and ECD spectropolarimetry were used to confirm that CLN025 adopts a β‐hairpin in water, TFE, MeOH, and DMSO and to examine differences in the secondary structure, hydrogen bonds, and weakly polar interactions. CLN025 was synthesized by microwave‐assisted solid phase peptide synthesis with N^α‐Fmoc protected amino acids. The VCD spectra displayed a (?,+,?) pattern with bands at 1640 to 1656 cm^?1, 1667 to 1687 cm^?1, and 1679 to 1686 cm^?1 formed by the overlap of a lower frequency negative couplet and a higher frequency positive couplet. A maximum IR absorbance was observed at 1647 to 1663 cm^?1 with component bands at 1630 cm^?1, 1646 cm^?1, 1658 cm^?1, and 1675 to 1680 cm^?1 that are indicative of the β‐sheet, random meander, either random meander or loop and turn, respectively. These results are similar to the results of others, who examined the VCD spectra of β‐hairpins formed by ^DPro‐Xxx turns and indicated that observed pattern is typical of β‐hairpins. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 442–450, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Vibrational Analysis of Phosphorothioate DNA: II. The POS Group in the Model Compound Dimethyl Phosphorothioate [(CH3O)2(POS)]-

Abstract

The results of Raman and Infrared (IR) spectroscopic investigations on the vibrational modes of dimethyl phosphorothioate (DMPS) anion, [(CH₃O)₂(POS)]^?, are reported. Ab initio calculations of the vibrational modes, the IR and Raman spectra and the interatomic force constants of DMPS were performed. A normal mode calculation was performed and the results were used to calculate the potential energy distribution for the vibrational modes. This analysis shows that in DMPS the P-S stretching mode has a frequency of about 630 cm^?1 and an angle bending mode involving the sulfur atom has a frequency of about 440 cm^?1. The proposed vibrational mode assignments will serve as marker bands in the conformational studies of phosphorothioate oligonucleotides which play a central role in the novel antisense therapeutic paradigm.     

Hinge-bending deformation of enzyme observed by microwave dielectric measurement

A dielectric relaxation peak due to an intramolecular motion in the active site of trypsin was first observed in aqueous solution below the freezing temperature of bulk water by a time domain reflectometry method. If trypsin inhibitor is added to the solution, it vanishes. It is suggested that the motion observed is a hinge-bending deformation giving rise to the enzymatic activity of trypsin, which is prohibited by linkage of the trypsin inhibitor. Relaxation time of the motion is 3 × 10² ns at − 10°C. © 1996 John Wiley & Sons, Inc.     

Hierarchical extensibility in the PEVK domain of skeletal-muscle titin

Titin is the main determinant of passive muscle force. Physiological extension of titin derives largely from its PEVK (Pro-Glu-Val-Lys) domain, which has a different length in different muscle types. Here we characterized the elasticity of the full-length, human soleus PEVK domain by mechanically manipulating its contiguous, recombinant subdomain segments: an N-terminal (PEVKI), a middle (PEVKII), and a C-terminal (PEVKIII) one third. Measurement of the apparent persistence lengths revealed a hierarchical arrangement according to local flexibility: the N-terminal PEVKI is the most rigid and the C-terminal PEVKIII is the most flexible segment within the domain. Immunoelectron microscopy supported the hierarchical extensibility within the PEVK domain. The effective persistence lengths decreased as a function of ionic strength, as predicted by the Odijk-Skolnick-Fixman model of polyelectrolyte chains. The ionic strength dependence of persistence length was similar in all segments, indicating that the residual differences in the elasticity of the segments derive from nonelectrostatic mechanisms.