PubMed HealthSearch

Biomedical subjects

M Orozco

Publications and source records attributed to M Orozco.

At least 19 recordsLinked to original sources

Molecular dynamics simulations of the unfolding of barnase in water and 8 M aqueous urea.

Molecular dynamics simulations of barnase have been conducted both in water and in 8 M urea solution for 500 ps at 25 degrees C and for 2000 ps at 85 degrees C. The final structure of the aqueous simulation at room temperature matches closely the structure obtained by NMR and the experimentally observed protections from isotopic exchange. The comparison of the structures generated by the aqueous simulation at 85 degrees C reveals a trajectory composed of groups of geometrically related structures separated by narrow regions of rapid change in structure. The first of these regions displays changes in backbone rmsd to the crystal structure and solvent-accessible area suggestive of a transition state, while the properties observed during the final 300 ps of the simulation are consistent with a stable intermediate. These assignments were confirmed by calculation of the "progress along the reaction coordinate" phi-values using an empirical equation based on a linear response method. The pathway of unfolding defined in this fashion agrees well with the experimental results of site-directed mutagenesis in terms of secondary structure content of the transition state and the intermediate and reproduces the relative stability of the different elements of secondary structure. The results of the simulations in urea suggest a mechanism at the molecular level for its well-known enhancement of the denaturation of proteins. The analysis of radial distribution functions shows that the first solvation shell of the protein is enriched in urea relative to the bulk solvent. The displacement of water molecules allows greater exposure of hydrophobic side chains, as witnessed particularly in the analysis of solvent-accessible surface areas at the higher temperature. Almost all urea molecules in the first shell form at least one hydrogen bond with the protein. They provide a more favorable environment for accommodation of the remaining water molecules, and they facilitate the separation of secondary structure elements by acting as a bridge between groups previously forming intraprotein hydrogen bonds.

Bacterial Proteins

Role of tautomerism of 2-azaadenine and 2-azahypoxanthine in substrate recognition by xanthine oxidase.

The tautomerism of 2-azaadenine and 2-hypoxanthine has been examined in the gas phase and in aqueous solution. The tautomerism in the gas phase has been studied by means of semiempirical and ab initio quantum-mechanical computations, as well as density-functional calculations. The influence of the aqueous solvent on the relative stability between tautomers has been estimated from self-consistent reaction field calculations performed with different high-level continuum models. The results provide a detailed picture of the tautomeric preference for these purine bases. The importance of tautomerism in the substrate recognition by xanthine oxidase is discussed. Finally, the rate of oxidation of 2-azaadenine and 2-hypoxanthine by xanthine oxidase is discussed in terms of the recognition model at the enzyme active site.

Adenine

Helical preferences of alanine, glycine, and aminoisobutyric homopeptides.

The stability between helical conformations of homopeptides of alanine, glycine, and aminoisobutyric acid has been studied by means of quantum-mechanical methods. The influence of peptide length on the relative stability between helical conformations has also been analyzed by means of systematic studies for peptides of size up to 11 residues. Finally, the influence of the solvent has been examined by using self-consistent reaction field methods. The results provide a detailed picture of the modulation exerted by these factors on the helical preferences of these peptides.

Alanine

Tautomerism of xanthine and alloxanthine: a model for substrate recognition by xanthine oxidase.

Tautomerism of neutral xanthine and alloxanthine has been examined both in the gas phase and in aqueous solution. The tautomeric preference in the gas phase has been studied by means of semiempirical and ab initio quantum-mechanical computations with inclusion of correlation effects at the Møller-Plesset level, and from density-functional calculations. The influence of solvent on the relative stability between tautomers has been estimated from self-consistent reaction field calculations performed with different models. The results provide a detailed picture of tautomerism for these biologically relevant purine bases. The functional implications in the recognition by xanthine oxidase are analyzed from inspection of the interaction patterns of the most stable tautomeric forms. A model for the recognition of these purine derivatives in the enzyme binding site is discussed.

Binding Sites

Molecular dynamics study of the binding of elsamicin A to DNA.

Molecular mechanics (MM) and dynamics (MD) were used to examine the conformational flexibility of elsamicin A, an antitumour antibiotic, in aqueous solution. Furthermore, MM and MD simulations were performed to obtain the first three-dimensional model of the elsamicin-A-DNA complex. Comparison of MD simulations of free and DNA-bound elsamicin A provided insight into the mechanism of DNA cleavage. Moreover, the drug-DNA complexes were analyzed to investigate specific drug-DNA interactions, as well as the role of the solvent water molecules in the drug-DNA binding. On the basis of the results derived from these simulations, several issues related to the sequence-selective binding of elsamicin A to DNA are discussed.

Aminoglycosides

DNA sequence-specific reading by echinomycin: role of hydrogen bonding and stacking interactions.

The binding of echinomycin to DNA hexamers of the form GpApXpZpTpC, where the central XpZ step can be CpG, TpA, GpC, or ApT, has been studied by molecular modeling and molecular mechanics techniques. Interaction energies have also been calculated for the complexation of echinomycin with sequences containing the preferred central CpG step and different flanking base pairs. Besides, two more sets of sequences incorporating either 2,6-diaminopurine (DAP) or hypoxanthine in place of adenine or guanine, respectively, have been examined. The aim of this work was to evaluate the relative importance of hydrogen-bonding and stacking interactions in the association of echinomycin with DNA and further rationalize the experimental evidence. The results of these calculations are in consonance with available data from footprinting experiments and appear to support our previous hypothesis that, in addition to the crucial intermolecular hydrogen bonds in the central region, the stacking interactions involving the quinoxaline-2-carboxamide chromophores of the drug and the DNA base pairs play an important role in modulating the binding specificity of this bisintercalating antitumor antibiotic. This is most clearly seen when sequences with similar minor-groove environments are compared (e.g. CpI vs TpA or CpG vs TpDAP). The dipole moment of N-methylquinoxaline-2-carboxamide has been measured (mu = 4.15 +/- 0.03 D) and compares very well with the calculated value (mu = 4.14 D). The fact that G:C, I:C, A:T, and DAP:T base pairs are shown to be endowed with distinct van der Waals and electrostatic stacking properties with respect to this heteroaromatic ring system could have important implications for the design of novel DNA mono- and bis-intercalating agents.

2-Aminopurine

Binding of echinomycin to d(GCGC)2 and d(CCGG)2: distinct stacking interactions dictate the sequence-dependent formation of Hoogsteen base pairs.

Molecular dynamics simulations have been used to explore the behavior of the complexes of echinomycin with the DNA tetramers d(GCGC)2 and d(CCGG)2 in which the terminal bases have been paired according to either a Hoogsteen or a Watson-Crick hydrogen bonding scheme. The energy of the four resulting complexes has been monitored along the dynamics trajectories and the interaction energy between echinomycin and DNA has been decomposed into contributions arising from the planar aromatic systems and the depsipeptide part of the antibiotic. Our calculations predict a large increase in overall stabilization upon protonation of the terminal cytosines and subsequent Hoogsteen pair formation in the complex of echinomycin with d(GCGC)2 but not with d(CCGG)2, in agreement with the experimental evidence [Gao and Patel, Quart. Rev. Biophys. 22, 93-138 (1989)]. The conformational preferences appear to arise mainly from differential stacking interactions in which the electrostatic component is shown to play a dominant role. Differences in hydrogen bonding patterns are also found among the complexes and these are compared in relation to available crystal structures. The binding of echinomycin to DNA appears as a complex process involving many interrelated variables.

Base Composition

Idiopathic chronic hypertrophic craniocervical pachymeningitis: case report.

A 55-year-old woman with a unique form of chronic hypertrophic pachymeningitis involving the posterior fossa and upper cervical spine is reported. Unlike other cases previously described, the clinical picture was dominated by signs of increased intracranial pressure, lower cranial nerve disorders, and a progressive cervical radiculomyelopathy. The diagnosis was made by means of a contrast-enhanced magnetic resonance imaging scan and confirmed by histological examination of the excised dura. Surgical treatment with removal of the hypertrophic dura provided temporary relief, although the natural history of the disease was not modified. Exhaustive bacteriological and histopathological studies failed to identify a specific cause for this diffuse hypertrophy of the cranial and cervical dura. The literature is reviewed, and other histologically documented cases are discussed.

Cervical Vertebrae

Mechanism for the rotamase activity of FK506 binding protein from molecular dynamics simulations.

Molecular dynamics (MD) and free energy perturbation (FEP) methods are used to study the binding and mechanism of isomerization of a tetrapeptide (AcAAPFNMe) by FK506 binding protein (FKBP). Detailed structures are predicted for the complexes of FKBP with the peptide in both ground-state and transition-state forms. The results support a mechanism of catalysis by distortion, where a large number of nonbonded interactions act together to stabilize preferentially the twisted transition state. The two most important groups for the catalysis are suggested to be Trp59 and Asp37, but several other groups are identified as directly or indirectly involved in the binding and catalysis. However, the structural results do not support the notion that the keto oxygen of the immunosuppressive agents FK506 and rapamycin mimics the oxygen for the twisted peptide bond in the FKBP-transition-state complex.

Amino Acid Sequence

Biosynthesis of lipoic acid: characterization of the lipoic acid auxotrophs Escherichia coli W1485-lip2 and JRG33-lip9.

The abilities of the Escherichia coli lipoic acid auxotrophs W1485-lip2 and JRG33-lip9 to grow on succinate medium in the presence of octanoate, 8-mercaptooctanoate, or 6-mercaptooctanoate have been determined. Both organisms are mutated in lipA. Neither organism can use octanoate or 6-mercaptooctanoate for production of lipoate, but the lip2 allele can use 8-mercaptooctanoate. Chromosomal DNA from the auxotrophs was amplified by PCR using primers derived from the DNA sequence of wild-type lipA and then sequenced. Both mutants contain single G/C to A/T mutations in lipA, resulting in conversion of Ser307 into Phe in W1485-lip2 and Glu195 into Lys in JRG33-lip9. These results support the hypothesis that lipA is involved in the sulfur insertion step(s) of lipoate biosynthesis and indicate that it is possible to selectively block formation of the C8-S bond through suitable mutation in lipA.

Amino Acid Sequence

Experimental and modelling studies on the DNA cleavage by elsamicin A.

The ability of elsamicin A, an antitumour antibiotic, to cleave DNA in the presence of ferrous iron and reducing agents, has been analysed using experimental and theoretical approaches. Experimentally, the antibiotic causes DNA breakage in the presence of ferrous ions and a reducing agent. The DNA-cleaving activity appears to be partially blocked by the action of superoxide dismutase and catalase. These results indicate that the elsamicin aglycone moiety (chartarin) can be involved in the production of free radicals. We have performed a broad theoretical study based in the quantum-mechanical framework, which allow us to determine the redox properties of elsamicin that lead to the generation of radical species. Our results clearly show that elsamicin acts as a true catalyst in the production of superoxide radicals. Moreover, it is suggested that the oxidation/reduction mechanism of the aglycone moiety of elsamicin (a lactone), leading to DNA breakage, is different from the mechanism followed by other well-known anti-cancer drugs, whose chromophore is a quinone.

Aminoglycosides

On the suitability of semiempirical calculations as sources of force field parameters.

The suitability of Dewar's Hamiltonians as a source of bonded force field parameters is explored from the comparison analysis between up to 270 semiempirically derived force field parameters and experimentally derived values reported in some of the most popular force fields. From the statistical analysis of the results, some general conclusions about the semiempirical parametrization are formulated.

Mathematical Computing

Giant aneurysm of the high internal carotid artery: surgical treatment.

We report the successful surgical treatment of a patient with a giant aneurysm of the high internal carotid artery by resection of the aneurysm and end-to-end anastomosis at the skull base, after resection of the styloid process and section of the stylohyoid and digastric muscles. No shunt was used and no major neurologic deficits resulted. Thirty months after operation the repaired artery is patent, with no evidence of stenosis or residual aneurysm and with good distal perfusion. Surgical treatment of high carotid artery aneurysms is possible with rather simple techniques; excision of the aneurysm followed by restoration of continuity of the carotid ends should be the treatment of choice regardless of the location of the aneurysm.

Aneurysm

A quantum chemical study of the enzymatic deamination of benzoadenine derivatives. A theoretical model of the interactions occurring between nucleosides and the active site of adenosine deaminase.

A theoretical study is presented, where, by using both ab initio and semi-empirical methodologies, the properties of benzoadenine derivatives as substrates of adenosine deaminase are discussed. The results suggest that lin-benzoadenine and lin-benzoadenosine can be recognized with an affinity similar to that of adenosine, but only if they are introduced about 0.12 nm deeper inside the active site of the enzyme than the natural substrate adenosine. This fact implies the existence of non-linear hydrogen bonds inside the active site of adenosine deaminase. Ab initio molecular electrostatic potential values suggest that these hydrogen bonds can exist, and have stability similar to that of linear hydrogen bonds. Finally, the great rate of deamination of lin-benzoadenine, comparable with that of adenosine despite the absence of the ribose, is explained in the context of the hypothesis that the protonation at the N1 atom is the rate-determining step of the whole deamination reaction.

Adenosine

Quantum chemical study of the electronic and conformational characteristics of adenosine and 8-substituted derivatives: functional implications in the mechanism of reaction of adenosine deaminase.

A quantum chemical study of 10 substrates of adenosine deaminase is performed. The conformational preference around the glycosidic bond of several 8-substituted derivatives of adenosine is studied using semiempirical modified neglect of diatomic overlap (MNDO) and Austin model 1 (AM1) methods. All the compounds studied show preference for the anti conformation; the syn - anti energetic differences calculated are small and in excellent agreement with experimental data. A relationship between the ab initio molecular electrostatic potential minimum energy of N3 and the syn - anti energetic difference is found. A highly significant relationship is also found between the ab initio net charge over the purine and pyrimidine rings and the logarithm of the maximum rate of deamination (log Vm) of the nucleosides by adenosine deaminase. In contrast, no significant relationship is found between the anti preference of 8-substituted derivatives of adenosine and their log Vm of deamination.

Adenosine

A practical procedure for the determination of electrostatic charges of large molecules.

A practical procedure for the precise determination of electrostatic charges, which are evaluated by fitting the rigorous quantum mechanical molecular electrostatic potential to a monopole-monopole expression, is presented. The proposal of this procedure arises from the study of the minimum requirements necessary to obtain reliable electrostatic charges. Such a study is focused on: (i) the dependence of the electrostatic charges on the set of points where the quantum mechanical and the monopole-monopole molecular electrostatic potentials are fitted; thus, both the influence of the number of points and their distribution in layers located out of the van der Waals radii of the atoms are examined, and (ii) the reliability of the use of fractional models for the evaluation of electrostatic charges of large molecules. Results point out that the optimum number of points is defined by a density of points ranging from 0.45 to 0.60 points/A2 when four layers (separated by 0.2 A) are considered. Nevertheless, the use of only two layers (separated by 0.4 A) for large molecules is recommended, thus enabling one to obtain reliable charges at a reduced computational cost. Moreover, results justify the use of fractional models for the determination of electrostatic charges of extremely large molecules, even when aromatic structures exist.

Adenosine

Molecular-mechanics modelling of drug-DNA structures; the effects of differing dielectric treatment on helix parameters and comparison with a fully solvated structural model.

This study analyses the influence that the nature of the dielectric constant has on the final structures obtained from in vacuo molecular mechanics calculations on a drug-DNA complex and compares these structures with the energy minimised complex including explicit solvent molecules. Minimisations have been performed on a proflavine-decanucleotide structure, where the drug was intercalated at the d(CpG) site of the d(GpApTpApCpGpApTpApC) decamer duplex, using two expressions for the dielectric constant: a distance-independent, epsilon ij = EPS, and a distance-dependent, epsilon ij = EPS*Rij, form and for values of EPS from 1 to 8. Significantly different structures are obtained for the distance-independent and the distance-dependent expressions of the dielectric constant. The use of a distance-independent dielectric constant leads to distorted structures, which are very sensitive to slight changes in the value of EPS. The use of a distance-dependent dielectric constant leads to less distorted and more stable structures. The effects on helical parameters are analysed in detail. The structures obtained for different values of EPS (within the distance-dependent formalism) seem to converge for values of EPS equal to 4 or greater. Based on these results a distance-dependent form of the dielectric with an EPS value of 4 is recommended in order to produce reliable refined nucleic acid structures by molecular mechanics. These conclusions have been supported by molecular-mechanics minimisation of the same structure with the inclusion of explicit water molecules and counter-ions.

Base Composition

TEFOOL/2: a program for theoretical drug design on microcomputers.

TEFOOL/2, a program written in BASIC, is presented in this paper. The purpose of TEFOOL/2 is to provide people interested in drug design with an easy-to-handle program where some of the most important techniques in QSAR are included. The program permits the selection of the training series, performs regression calculations and searches for optimum substituents. The latter is achieved by using either a Hansch's strategy or geometrical procedures. The program is interactive and can be implemented on an IBM-PC or compatible microcomputer. Although TEFOOL/2 has been developed for its application in drug design studies, its great flexibility makes it suitable for application to any experimental design or optimization process.

Algorithms