Studies on molecular evolution and structural features of double-headed inhibitors of ?-amylase and trypsin in plants

Plant seeds usually have high concentrations of proteinase and amylase inhibitors. These inhibitors exhibit a wide range of specificity, stability and oligomeric structure. In this communication, we report analysis of sequences that show statistically significant similarity to the double-headed alpha-amylase/trypsin inhibitor of ragi (Eleusine coracana). Our aim is to understand their evolutionary and structural features. The 14 sequences of this family that are available in the SWISSPROT database form three evolutionarily distinct branches. The branches relate to enzyme specificities and also probably to the oligomeric state of the proteins and not to the botanical class of the plant from which the enzymes are derived. This suggests that the enzyme specificities of the inhibitors evolved before the divergence of commercially cultivated cereals. The inhibitor sequences have three regions that display periodicity in hydrophobicity. It is likely that this feature reflects extended secondary structure in these segments. One of the most variable regions of the polypeptide corresponds to a loop, which is most probably exposed in the native structure of the inhibitors and is responsible for the inhibitory property.

Basis set dependence of the molecular electrostatic potential topography. A case study of substituted benzenes

The basis set dependence of the topographical structure of the molecular electrostatic potential (MESP), as well as the effect of substituents on the MESP distribution, has been investigated with substituted benzenes as test cases. The molecules are studied at HF-SCF 3�21G and 6�31G** levels, with a further MESP topographical investigation at the 3�21G, double-zeta, 6�31G*, 6�31G**, double-zeta polarized and triple-zeta polarized levels. The MESP critical points for a 3�21G optimized/6�31G** basis are similar to the corresponding 6�31G** optimized/6�31G** ones. More generally, the qualitative features of the MESP topography computed at the polarized level are independent of the level at which optimization is carried out. For a proper representation of oxygen lone pairs, however, optimization using a polarized basis set is required. The nature of the substituent drastically changes the MESP distribution over the phenyl ring. The values and positions of MESP minima indicate the most active site for electrophilic attack. This point is strengthened by a study of disubstituted benzenes.

Chain folding and A:T pairing in human telomeric DNA: a model-building and molecular dynamics study

The various types of chain folding and possible intraloop as well as interloop base pairing in human telomeric DNA containing d(TTAG(3)) repeats have been investigated by model-building, molecular mechanics, and molecular dynamics techniques. Model-building and molecular mechanics studies indicate that it is possible to build a variety of energetically favorable folded-back structures with the two TTA loops on same side and the 5' end thymines in the two loops forming TATA tetrads involving a number of different intraloop as well as interloop A:T pairing schemes. In these folded-back structures, although both intraloop and interloop Watson-Crick pairing is feasible, no structure is possible with interloop Hoogsteen pairing. MD studies of representative structures indicate that the guanine-tetraplex stem is very rigid and, while the loop regions are relatively much more flexible, most of the hydrogen bonds remain intact throughout the 350-ps in vacuo simulation. The various possible TTA loop structures, although they are energetically similar, have characteristic inter proton distances, which could give rise to unique cross-peaks in two-dimensional nuclear Overhauser effect spectroscopy (NOESY) experiments. These folded-back structures with A:T pairings in the loop region help in rationalizing the data from chemical probing and other biochemical studies on human telomeric DNA.

Molecular dynamics studies on nucleoside 2',3'-cyclic phosphates

2',3'-cyclic nucleotides are intermediates and substrates of Ribonuclease (RNase)-catalysed reactions. The characterization of the equilibrium conformation as well as the flexibility inherent in these molecules helps in understanding the enzymatic action of RNases. The present study explores parameters like phase angle, glycosydic torsion angle and hydrogen bond to find possible interrelationship between them through Molecular Dynamics (MD) simulations on 3'-GMP, 3'-UMP, A>p, G>p, U>p, C>p, GpA>p and UpA>p. Interesting results of the effect of cyclisation and other constraints such as hydrogen bond between certain groups on the equilibrium ribose conformation have emerged from this study.

Influence of glycosidic linkage on the nature of carbohydrate binding in beta-prism I fold lectins: An X-ray and molecular dynamics investigation on banana lectin-carbohydrate complexes

The three crystal structures reported here provide details of the interactions of mannose and the mannosyl-alpha-1,3-mannose component of a pentamannose with banana lectin and evidence for the binding of glucosyl-alpha-1,2-glucose to the lectin. The known structures involving the lectin include a complex with glucosyl-beta-1,3-glucose. Modeling studies on the three disaccharide complexes with the reducing end and the nonreducing end at the primary binding site are also provided here. The results of the Xray and modeling studies show that the disaccharides with an alpha-1,3 linkage prefer to have the nonreducing end at the primary binding site, whereas the reducing end is preferred at the site when the linkage is beta-1,3 in mannose/glucose-specific beta-prism I fold lectins. In the corresponding galactose-specific lectins, however, alpha-1,3-linked disaccharides cannot bind the lectin with the nonreducing end at the primary binding site on account of steric clashes with an aromatic residue that occurs only when the lectin is galactose-specific. Molecular dynamics simulations based on the known structures involving banana lectin enrich the information on lectin-carbohydrate interactions obtained from crystal structures. They demonstrate that conformational selection as well as induced fit operate when carbohydrates bind to banana lectin.

The stability of transmembrane helices: A molecular dynamics study on the isolated helices of bacteriorhodopsin

Bacteriorhodopsin (bR) continues to be a proven testing ground for the study of integral membrane proteins (IMPs). It is important to study the stability of the individual helices of bR, as they are postulated to exist as independently stable transmembmne helices (TMHs) and also for their utility as templates for modeling other IMPs with the postulated seven-helix bundle topology. Toward this purpose, the seven helices of bR have been studied by molecular dynamics simulation in this study. The suitability of using the backbone-dependent rotamer library of side-chain conformations arrived at from the data base of globular protein structures in the case TMHs has been tested by another set of ? helix simulations with the side-chain orientations taken from this library. The influence of the residue's net charge oil the helix stability was examined by simulating the helices III, IV, and VI (from both of the above sets of helices) with zero net charge on the side chains. The results of these 20 simulations demonstrate in general the stability of the isolated helices of bR in conformity with the two-stage hypothesis of IMP folding. However, the helices I, II, V, and VII are more stable than the other three helices. The helical nature of certain regions of III, IV, and VI are influenced by factors such as the net charge and orientation of several residues. It is seen that the residues Arg, Lys, Asp, and Glu (charged residues), and Ser, Thr, Gly, and Pro, play a crucial role in the stability of the helices of bR. The backbone-dependent rotamer library for the side chains is found to be suitable for the study of TMHs in IMP. (C) 1996 John Wiley & Sons, Inc.

Molecular-dynamics studies of systems of confined dumbbell molecules

We present the results of molecular-dynamics simulations of systems of dumbbell molecules confined by parallel molecular walls. We have carried out systematic studies of three cases: freezing, steady flows, and stick-slip friction. We find that the molecular orientational degrees of freedom cause the surface layers to deviate from a planar configuration. Nevertheless, steady flows, in a channel as narrow as 15 molecular sizes, display continuum behavior. A range of mechanisms in the dynamics of the freezing of a confined fluid is found, as a function of the wall-fluid interactions and the bond length of the dumbbell molecules. The simple order-disorder transition associated with stick-slip motion in the presence of a layer of monoatomic lubricant molecules is supplanted by more complex behavior due to rotational degrees of freedom of the diatomic molecules.

Quantum theoretical study of molecular mechanisms of mutation and cancer - A review

Several endogenous and exogenous chemical species, particularly the so-called reactive oxygen species (ROS) and reactive nitrogen oxide species (RNOS), attack deoxyribonucleic acid (DNA) in biological systems producing DNA lesions which hamper normal cell functioning and cause various diseases including mutation and cancer. The guanine (G) base of DNA among all the bases is most susceptible and certain modified guanines get involved in mispairing with other bases during DNA replication. The biological system repairs the abnormal base pairs, but those that are still left cause mutation and cancer. Anti-oxidants present in biological systems can scavenge the ROS and RNOS. Thus three types of molecular events occur in biological media: (i) DNA damage, (ii) DNA repair, and (iii) prevention of DNA damage by scavenging ROS and RNOS. Quantum mechanical methods may be used to unravel molecular mechanisms of such phenomena. Some recent quantum theoretical results obtained on these problems are reviewed here.

Effect of molecular association and reactivity on substitution in chromium(III)-salprn complexes

A new chromium(III)-Schiff base complex, [Cr(5-chlorosalprn)(H2O)(2)]ClO4, where salprn=N,N'-propylenebis(salicylideneimine) has been prepared and characterized by electrospray ionization mass spectrometric (ESIMS) analysis and other spectroscopic techniques. Single crystal X-ray data reveal that the complex assumes a trans-diaquo structure, [Cr(C17H18Cl2N2O4)]ClO4.H2O. The effect of phenyl ring substituents on the rate of formation of [O=Cr-V Schiff base](+) has been investigated. The bimolecular rate constant for the formation of O=Cr-V species by the [Cr(Schiff base)(H2O)(2)]ClO4, where the Schiff base=salprn, (1) and 5-chlorosalprn, (2) with PhOI was compared. In the case of (2) the rate was found to be faster by an order of magnitude at pH=4 compared to (1). The introduction of a chloro-substituent on the phenyl ring not only influences the rate of redox reactivity but also the pKa values of aquo ligands of the complexes, indicating the difference in the electronic environment around the metal ion in both (1) and (2).

A molecular dynamics study based post facto free energy analysis of the binding of bovine angiogenin with UMP and CMP ligands

Angiogenin is a protein belonging to the superfamily of RNase A. The RNase activity of this protein is essential for its angiogenic activity. Although members of the RNase A family carry out RNase activity, they differ markedly in their strength and specificity. In this paper, we address the problem of higher specificity of angiogenin towards cytosine against uracil in the first base binding position. We have carried out extensive nano-second level molecular dynamics(MD) computer simulations on the native bovine angiogenin and on the CMP and UMP complexes of this protein in aqueous medium with explicit molecular solvent. The structures thus generated were subjected to a rigorous free energy component analysis to arrive at a plausible molecular thermodynamic explanation for the substrate specificity of angiogenin.

Molecular insights into the mechanism of phenotypic tolerance to rifampicin conferred on mycobacterial RNA polymerase by MsRbpA

The protein MsRbpA from Mycobacterium smegmatis rescues RNA polymerase (RNAP) from the inhibitory effect of rifampicin (Rif). We have reported previously that MsRbpA interacts with the beta-subunit of RNAP and that the effect of MsRbpA on Rif-resistant (Rif(R)) RNAP is minimal. Here we attempted to gain molecular insights into the mechanism of action of this protein with respect to its role in rescuing RNAP from Rif-mediated transcription inhibition. Our experimental approach comprised multiple-round transcription assays, fluorescence spectroscopy, MS and surface plasmon resonance in order to meet the above objective. Based on our molecular studies we propose here that Rif is released from its binding site in the RNAP-Rif complex in the presence of MsRbpA. Biophysical studies reveal that the location of MsRbpA on RNAP is at the junction of the beta- and beta'-subunits, close to the Rif-binding site and the (i + 1) site on RNAP.

Tryptophan-water interaction in Monellin:Hydration patterns from molecular dynamics simulations

Femtosecond spectroscopy carried out earlier on Monellin and some other systems has given insights into the hydration dynamics of the proteins. In the present work, molecular dynamics simulations have been performed on Monellin to study the hydration dynamics. A method has been described to follow up the molecular events of the protein–water interactions in detail. The time constants of the survival correlation function match well with the reported experimental values. This validates the procedure, adapted here for Monellin, to investigate the hydration dynamics in general.

Quaternary association in beta-prism I fold plant lectins: Insights from X-ray crystallography, modelling and molecular dynamics

Dimeric banana lectin and calsepa, tetrameric artocarpin and octameric heltuba are mannose-specific beta-prism I fold lectins of nearly the same tertiary structure. MD simulations on individual subunits and the oligomers provide insights into the changes in the structure brought about in the protomers on oligomerization, including swapping of the N-terminal stretch in one instance. The regions that undergo changes also tend to exhibit dynamic flexibility during MD simulations. The internal symmetries of individual oligomers are substantially retained during the calculations. Energy minimization and simulations were also carried out on models using all possible oligomers by employing the four different protomers. The unique dimerization pattern observed in calsepa could be traced to unique substitutions in a peptide stretch involved in dimerization. The impossibility of a specific mode of oligomerization involving a particular protomer is often expressed in terms of unacceptable steric contacts or dissociation of the oligomer during simulations. The calculations also led to a rationale for the observation of a heltuba tetramer in solution although the lectin exists as an octamer in the crystal, in addition to providing insights into relations among evolution, oligomerization and ligand binding.