3 Nsec molecular dynamics simulation of the protein ubiquitin and comparison with X-ray crystal and solution NMR structures.

Mainly due to computational limitations, past protein molecular dynamics simulations have rarely been extended to 300 psec; we are not aware of any published results beyond 350 psec. The present work compares a 3000 psec simulation of the protein ubiquitin with the available x-ray crystallographic and solution NMR structures. Aside from experimental structure availability, ubiquitin was studied because of its relatively small size (76 amino acids) and lack of disulfide bridges. An implicit solvent model was used except for explicit treatment of waters of crystallization. We found that the simulated average structure retains most of the character of the starting x-ray crystal structure. In two highly surface accessible regions, the simulation was not in agreement with the x-ray structure. In addition, there are six backbone-backbone hydrogen bonds that are in conflict between the solution NMR and x-ray crystallographic structures; two are bonds that the NMR does not locate, and four are ones that the two methods disagree upon the donor. Concerning these six backbone-backbone hydrogen bonds, the present simulation agrees with the solution NMR structure in five out-of-the six cases, in that if a hydrogen bond is present in the x-ray structure and not in the NMR structure, the bond breaks within 700 psec. Of the two hydrogen bonds that are found in the NMR structure and not in the x-ray structure, one forms at 1400 psec and the other forms rarely. The present results suggest that relatively long molecular dynamics simulations, that use protein x-ray crystal coordinates for the starting structure and a computationally efficient solvent representation, may be used to gain an understanding of conformational and dynamic differences between the solid-crystal and dilute-solution states.     

Interspecific variation of body shape and sexual dimorphism in three coexisting species of the genus <Emphasis Type="Italic">Petrotilapia</Emphasis> (Teleostei: Cichlidae) from Lake Malawi

Differences in color patterns have been the most used feature in describing cichlid species belonging to genus Petrotilapia from Lake Malawi. In this study, we quantified morphological variation in body shape within and among three coexisting Petrotilapia species using landmark-based geometric morphometric methods. Statistic analyses revealed significant body shape differences among species but not between sexes. Post hoc multiple comparisons based on Mahalanobis distances revealed that P. nigra was significantly different from P. genalutea and Petrotilapia sp., whereas the latter two were not significantly different. The splines generated showed that the most pronounced variation was in the head region, in which P. nigra had a relatively longer and deeper head than the other two. The most clear-cut distinction was in gape length; P. genalutea had the longest gape, followed by Petrotilapia sp., whereas P. nigra had the shortest gape. Body depth was shallower in P. nigra than the others. When comparing sexes by their centroid size, ANOVA revealed that males were bigger than females. Therefore, we conclude that color is not the only feature that can distinguish these congeners. We discuss the observed sexual dimorphism in terms of sexual selection and relate morphological variation among species to feeding behavior, which may help explain their coexistence in nature.     

Modeling Study on the Cleavage Step of the Self-Splicing Reaction in Group I Introns

Abstract

A three-dimensional model of the Tetrahymena thermophila group I intron is used to further explore the catalytic mechanism of the transphosphorylation reaction of the cleavage step. Based on the coordinates of the catalytic core model proposed by Michel and Westhof (Michel, F., Westhof, E.J. Mol. Biol. 216, 585–610 (1990)), we first converted their ligation step model into a model of the cleavage step by the substitution of several bases and the removal of helix P9. Next an attempt to place a trigonal bipyramidal transition state model in the active site revealed that this modified model for the cleavage step could not accommodate the transition state due to insufficient space. A lowering of PI helix relative to surrounding helices provided the additional space required. Simultaneously, it provided a better starting geometry to model the molecular contacts proposed by Pyle et al. (Pyle, A M., Murphy, F. L., Cech, T. R. Nature 358, 123–128. (1992)), based on mutational studies involving the J8/7 segment Two hydrated Mg²⁺ complexes were placed in the active site of the ribozyme model, using the crystal structure of the functionally similar Klenow fragment (Beese, L.S., Steitz, T.A. EMBO J. 10, 25–33 (1991)) as a guide. Thepresence of two metal ions in the active site of the intron differs from previous models, which incorporate one metal ion in the catalytic site to fulfill the postulated roles of Mg²⁺ in catalysis. The reaction profile is simulated based on a trigonal bipyramidal transition state, and the role of the hydrated Mg²⁺ complexes in catalysis is further explored using molecular orbital calculations.     

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.     

GluK1 antagonists from 6-(tetrazolyl)phenyl decahydroisoquinoline derivatives: In vitro profile and in vivo analgesic efficacy

We have explored the decahydroisoquinoline scaffold, bearing a phenyl tetrazole, as GluK1 antagonists with potential as oral analgesics. We have established the optimal linker atom between decahydroisoquinoline and phenyl rings and demonstrated an improvement of both the affinity for the GluK1 receptor and the selectivity against the related GluA2 receptor with proper phenyl substitution. In this Letter, we also disclose in vivo data that led to the discovery of LY545694·HCl, a compound with oral efficacy in two persistent pain models.     

Structure-activity relationship of a series of diaminoalkyl substituted benzimidazole as neuropeptide Y Y1 receptor antagonists

A series of benzimidazoles (4) was synthesized and evaluated in vitro as potent and selective NPY Y1 receptor antagonists. Substitution of the piperidine nitrogen of 4 with appropriate R groups resulted in compounds with more than 80-fold higher affinity at the Y receptor compared to the parent compound 5 (R = H). The most potent benzimidazole in this series was 21 (Ki = 0.052 nM).     

Molecular dynamics simulations of a protein-protein dimer: particle-mesh Ewald electrostatic model yields far superior results to standard cutoff model.

In this article we present two 1000 ps molecular dynamics simulations on the rat micro-glutathione S-transferase dimeric enzyme in complex with the product 1-(S-glutathionyl)-2,4-dinitrobenzene, in a periodic box with explicit solvent molecules, and investigate the effect of long-range electrostatics models on the structure and dynamics of the dimer and its components. One simulation used the standard cutoff method (10A), whilst the other used the particle-mesh Ewald (PME) method. We monitored the root mean-square atomic deviation (RMSD) from the initial crystal structure to examine the convergence of both simulations, as well as several other structural parameters such as the distance between active sites, rigid body rotation between domains in subunits, radius of gyration, B-factors, number of hydrogen bonds and salt bridges and solvent-accessible surface area. For example, with the PME method, the dimer structure remains much closer to the initial crystallographic structure with an average RMSD of 1.3A +/- 0.1A and 1.0A +/- 0.1A for all heavy and backbone atoms, respectively, in the last 200 ps; the respective values for the cutoff simulation are 4.7A +/- 0.3A and 4.2A +/- 0.3A. The large deviations observed in the cutoff simulation severely affected the stability of the enzyme dimer and its complex with the bound product. This finding is contrary to that found in a similar study of the monomeric protein ubiquitin [Fox, T. & Kollman, P. A. Proteins Struct. Func. Genet. 25, 315-334 (1996)]. Unlike the earlier published work, the present study provides evidence that the standard cutoff method is not generally valid for the study of protein complexes, or their subunits.     

Controlling the regiospecificity and coupling of cytochrome P450cam: T185F mutant increases coupling and abolishes 3-hydroxynorcamphor product.

Cytochrome P450cam (P450CIA1) catalyzes the hydroxylation of camphor and several substrate analogues such as norcamphor and 1-methyl-norcamphor. Hydroxylation was found experimentally at the 3, 5, and 6 positions of norcamphor, but only at the 5 and 6 positions of 1-methyl-norcamphor. In the catalytic cycle, the hydroxylation of substrate is coupled to the consumption of NADH. For camphor, the degree of coupling is 100%, but for both norcamphor and 1-methyl-norcamphor, the efficiency is dramatically lowered to 12% and 50%, respectively. Based on an examination of the active site of P450cam, it appeared that mutating position 185 might dramatically alter the product specificity and coupling of hydroxylation of norcamphor by P450cam. Analysis of molecular dynamics trajectories of norcamphor bound to the T185F mutant of cytochrome P450cam predicted that hydroxylation at the 3 position should be abolished and that the coupling should be dramatically increased. This mutant was constructed and the product profile and coupling experimentally determined. The coupling was doubled, and hydroxylation at the 3 position was essentially abolished. Both of these results are in agreement with the prediction.     

Applicability of PM3 to transphosphorylation reaction path: Toward designing a minimal ribozyme

A growing body of evidence shows that RNA can catalyze many of the reactions necessary both for replication of genetic material and the possible transition into the modern protein-based world. However, contemporary ribozymes are too large to have self-assembled from a prebiotic oligonucleotide pool. Still, it is likely that the major features of the earliest ribozymes have been preserved as molecular fossils in the catalytic RNA of today. Therefore, the search for a minimal ribozyme has been aimed at finding the necessary structural features of a modern ribozyme (Beaudry and Joyce, 1990). Both a three-dimensional model and quantum chemical calculations are required to quantitatively determine the effects of structural features of the ribozyme on the reaction it catalyzes. Using this model, quantum chemical calculations must be performed to determine quantitatively the effects of structural features on catalysis. Previous studies of the reaction path have been conducted at theab initio level, but these methods are limited to small models due to enormous computational requirements. Semiempirical methods have been applied to large systems in the past; however, the accuracy of these methods depends largely on the system under investigation. In the preent study we assess the validity of the MNDO/PM3 method on a simple model of the ribozyme-catalyzed reaction, or hydrolysis of phosphoric acid. We find that the results are qualitatively similar toab initio results using large basis sets. Therefore, PM3 is suitable for studying the reaction path of the ribozyme-catalyzed reaction.     

Correlation of Tm and sequence of DNA duplexes with delta H computed by an improved empirical potential method

T_m values for 20 DNA duplexes with different repeating base sequences provided the data base for developing a rational and relatively simple methodology for computing apparent enthalpies for the helix → coil transitions of DNA helices, ΔH _calc. The computational variables and their range of acceptable values were selected on the basis of physically plausible arguments. Over 350,000 different combinations of the variables were tested for degree of fit. It was therby possible to find a combination giving a high degree of linear fit between T_m and ΔH _calc (correlation coefficient, 0.99), with T_m values deviating (on average) from the regression line by only ±2.17°C. Most of this uncertainty is attributed to experimental limitations, although computational approximations also contribute. With ΔH _calc for the melting of each of the unique complementary dinucleotide fragments computed by the method developed, it is possible to estimate T_m and (relative) ΔH _calc reliable for the melting of any particular DNA [with base pairs G(I)·C and A·T] given only its base sequence. The ΔH_calc values for the complementary dinucleotide fragments, together with statistical considerations, make it apparent that T_m of DNAs with repeating base sequence show only an approximate linear dependence on G·C content because A·T and G·G pairs do not contribute to helix stability independently of the base-pair sequence in which they occur. In fact, the nearestneighbor stacking interactions are so significant that certain complementary dinucleotide fragment sequences with 0,50, and 100% G·C content have the same stability.