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.

A new pentameric structure of rotavirus NSP4 revealed by molecular replacement

The region spanning residues 95-146 of the rotavirus nonstructural protein NSP4 from the asymptomatic human strain ST3 has been purified and crystallized and diffraction data have been collected to a resolution of 2.6 angstrom. Several attempts to solve the structure by the molecular-replacement method using the available tetrameric structures of this domain were unsuccessful despite a sequence identity of 73% to the already known structures. A more systematic approach with a dimer as the search model led to an unexpected pentameric structure using the program Phaser. The various steps involved in arriving at this molecular-replacement solution, which unravelled a case of subtle variation between different oligomeric states unknown at the time of solving the structure, are presented in this paper.

Insights into the Substrate Specificity of a Thioesterase Rv0098 of Mycobacterium Tuberculosis through X-ray Crystallographic and Molecular Dynamics Studies

The crystal structure of Rv0098, a long-chain fatty acyl-CoA thioesterase from Mycobacterium tuberculosis with bound dodecanoic acid at the active site provided insights into the mode of substrate binding but did not reveal the structural basis of substrate specificities of varying chain length. Molecular dynamics studies demonstrated that certain residues of the substrate binding tunnel are flexible and thus modulate the length of the tunnel. The flexibility of the loop at the base of the tunnel was also found to be important for determining the length of the tunnel for accommodating appropriate substrates. A combination of crystallographic and molecular dynamics studies thus explained the structural basis of accommodating long chain substrates by Rv0098 of M. tuberculosis.

Crowding, molecular volume and plasticity: An assessment involving crystallography, NMR and simulations

The discrepancy between the X-ray and NMR structures of Mycobacterium tuberculosis peptidyl-tRNA hydrolase in relation to the functionally important plasticity of the molecule led to molecular dynamics simulations. The X-ray and the NMR studies along with the simulations indicated an inverse correlation between crowding and molecular volume. A detailed comparison of proteins for which X-ray and the NMR structures appears to confirm this correlation. In consonance with the reported results of the investigations in cellular compartments and aqueous solution, the comparison indicates that the crowding results in compaction of the molecule as well as change in its shape, which could specifically involve regions of the molecule important in function. Crowding could thus influence the action of proteins through modulation of the functionally important plasticity of the molecule. Selvaraj M, Ahmad R, Varshney U and Vijayan M 2012 Crowding, molecular volume and plasticity: An assessment involving crystallography, NMR and simulations. J. Biosci. 37 953-963] DOI 10.1007/s12038-012-9276-5

Structural and molecular basis of interaction of HCV non-structural protein 5A with human casein kinase 1 alpha and PKR

Background: Interaction of non-structural protein 5A (NS5A) of Hepatitis C virus (HCV) with human kinases namely, casein kinase 1 alpha (ck1 alpha) and protein kinase R (PKR) have different functional implications such as regulation of viral replication and evasion of interferon induced immune response respectively. Understanding the structural and molecular basis of interactions of the viral protein with two different human kinases can be useful in developing strategies for treatment against HCV. Results: Serine 232 of NS5A is known to be phosphorylated by human ck1 alpha. A structural model of NS5A peptide containing phosphoacceptor residue Serine 232 bound to ck1 alpha has been generated using the known 3-D structures of kinase-peptide complexes. The substrate interacting residues in ck1 alpha has been identified from the model and these are found to be conserved well in the ck1 family. ck1 alpha - substrate peptide complex has also been used to understand the structural basis of association between ck1 alpha and its other viral stress induced substrate, tumour suppressor p53 transactivation domain which has a crystal structure available. Interaction of NS5A with another human kinase PKR is primarily genotype specific. NS5A from genotype 1b has been shown to interact and inhibit PKR whereas NS5A from genotype 2a/3a are unable to bind and inhibit PKR efficiently. This is one of the main reasons for the varied response to interferon therapy in HCV patients across different genotypes. Using PKR crystal structure, sequence alignment and evolutionary trace analysis some of the critical residues responsible for the interaction of NS5A 1b with PKR have been identified. Conclusions: The substrate interacting residues in ck1 alpha have been identified using the structural model of kinase substrate peptide. The PKR interacting NS5A 1b residues have also been predicted using PKR crystal structure, NS5A sequence analysis along with known experimental results. Functional significance and nature of interaction of interferon sensitivity determining region and variable region 3 of NS5A in different genotypes with PKR which was experimentally shown are also supported by the findings of evolutionary trace analysis. Designing inhibitors to prevent this interaction could enable the HCV genotype 1 infected patients respond well to interferon therapy.

Multiple molecular alterations in phosphodegrons 1-3 within AAV2 capsid demonstrates higher hepatic gene transfer efficiency

The success of AAV2 mediated hepatic gene transfer in human trials for diseases such as hemophilia has been hampered by a combination of low transduction efficiency and a robust immune response directed against these vectors. We have previously shown that AAV2 is targeted for destruction in the cytoplasm by the host-cellular kinase/ubiquitination/proteasomal degradation machinery and modification of the serine(S)/threonine(T) kinase and lysine(K) targets on AAV capsid is beneficial. Thus targeted single mutations of S/T>A(S489A, S498A, T251A) and K>R (K532R) improved the efficiency of gene transfer in vivo as compared to wild type (WT)-AAV2 vectors (∼6-14 fold). In the present study, we evaluated if combined alteration of the phosphodegrons (PD), which are the phosphorylation sites recognized as degradation signals by ubiquitin ligases, improves further the gene transfer efficiency. Thus, we generated four multiple mutant vectors (PD: 1+3, S489A+K532R, PD: 1+3, S489A+K532R together with T251 residue which did not lie in any of the phosphodegrons but had shown increased transduction efficiency compared to the WT-AAV2 vector (∼6 fold) and was also conserved in 9 out of 10 AAV serotypes (AAV 1 to 10), PD: 1+3, S489A+K532R+S498A and a fourth combination of PD: 3, K532R+T251. We then evaluated them in vitro and in vivo and compared their gene transfer efficiency with either the WT-AAV2 or the best single mutant S489A-AAV2 vector. The novel multiple mutations on the AAV2 capsid did not affect the overall vector packaging efficiency. All the multiple AAV2 mutants showed superior transduction efficiency in HeLa cells in vitro when compared to either the WT (62-72% Vs 21%) or the single mutant S489A (62-72% Vs 50%) AAV2 vectors as demonstrated by FACS analysis (Fig. 1A). On hepatic gene transfer with 5x10^10 vgs per animal in C57BL/6 mice, all the multiple mutants showed increased transgene expression compared to either the WT-AAV2 (∼15-23 fold) or the S489A single mutant vector (∼2-3 fold) (Fig.1B and C). These novel multiple mutant AAV2 vectors also showed higher vector copy number in murine hepatocytes 4 weeks post transduction, as compared to the WT-AAV2 (∼5-6 Vs 1.4 vector copies/diploid genome) and further higher when compared to the single mutant S489A(∼5-6 fold Vs 3.8 fold) (Fig.1D). Further ongoing studies will demonstrate the therapeutic benefit of one or more of the multiple mutants vectors in preclinical models of hemophilia.

Analysis of designed beta-hairpin peptides: molecular conformation and packing in crystals

The crystal structures of several designed peptide hairpins have been determined in order to establish features of molecular conformations and modes of aggregation in the crystals. Hairpin formation has been induced using a centrally positioned (D)Pro-Xxx segment (Xxx = (L)Pro, Aib, Ac(6)c, Ala; Aib = alpha-aminoisobutyric acid; Ac(6)c = 1-aminocyclohexane-1-carboxylic acid). Structures of the peptides Boc-Leu-Phe-Val-(D)Pro-(L)Pro-Leu-Phe-Val-OMe (1), Boc-Leu-Tyr-Val-(D)Pro-(L)Pro-Leu-Phe-Val-OMe (2, polymorphic forms labeled as 2a and 2b), Boc-Leu-Val-Val-(D)Pro-(L)Pro-Leu-Val-Val-OMe (3), Boc-Leu-Phe-Val-(D)Pro-Aib-Leu-Phe-Val-OMe (4, polymorphic forms labeled as 4a and 4b), Boc-Leu-Phe-Val-(D)Pro-Ac(6)c-Leu-Phe-Val-OMe (5) and Boc-Leu-Phe-Val-(D)Pro-Ala-Leu-Phe-Val-OMe (6) are described. All the octapeptides adopt type II' beta-turn nucleated hairpins, stabilized by three or four cross-strand intramolecular hydrogen bonds. The angle of twist between the two antiparallel strands lies in the range of -9.8 degrees to -26.7 degrees. A detailed analysis of packing motifs in peptide hairpin crystals is presented, revealing three broad modes of association: parallel packing, antiparallel packing and orthogonal packing. An attempt to correlate aggregation modes in solution with observed packing motifs in crystals has been made by indexing of crystal faces in the case of three of the peptide hairpins. The observed modes of hairpin aggregation may be of relevance in modeling multiple modes of association, which may provide insights into the structure of insoluble polypeptide aggregates.

Molecular flexibility of Mycobacterium tuberculosis ribosome recycling factor and its functional consequences: An exploration involving mutants

Internal mobility of the two domain molecule of ribosome recycling factor (RRF) is known to be important for its action. Mycobacterium tuberculosis RRF does not complement E. coli for its deficiency of RRF (in the presence of E. coli EF-G alone). Crystal structure had revealed higher rigidity of the M. tuberculosis RRF due to the presence of additional salt bridges between domains. Two inter-domain salt bridges and one between the linker region and the domain containing C-terminal residues were disrupted by appropriate mutations. Except for a C-terminal deletion mutant, all mutants showed RRF activity in E. coli when M. tuberculosis EF-G was also co-expressed. The crystal structures of the point mutants, that of the C-terminal deletion mutant and that of the protein grown in the presence of a detergent, were determined. The increased mobility resulting from the disruption of the salt bridge involving the hinge region allows the appropriate mutant to weakly complement E. coli for its deficiency of RRF even in the absence of simultaneous expression of the mycobacterial EF-G. The loss of activity of the C-terminal deletion mutant appears to be partly due to the rigidification of the molecule consequent to changes in the hinge region.