Exploring potential energy surfaces to reach saddle points above convex regions.

Saddle points on high-dimensional potential energy surfaces (PES) play a determining role in the activated dynamics of molecules and materials. Building on approaches dating back more than 50 years, many open-ended transition-state search methods have been developed to follow the direction of negative curvature from a local minimum to an adjacent first-order saddle point. Despite the mathematical justification, these methods can display a high failure rate: using small deformation steps, up to 80% of the explorations can end up in a convex region of the PES, where all directions of negative curvature vanish, while if the deformation is aggressive, a similar fraction of attempts lead to saddle points that are not directly connected to the initial minimum. In high-dimension PES, these reproducible failures were thought to only increase the overall computational cost, without having any effect on the methods' capacity to find all saddle points surrounding a minimum. Using activation-relaxation technique nouveau (ARTn), we characterize the nature of the PES around minima, considerably expanding on previous knowledge. We show that convex regions can lie on activation pathways and that not exploring beyond them can introduce significant bias in the saddle-point search. We introduce an efficient approach for traversing the convex regions, almost eliminating exploration failures, while multiplying by almost 10 the number of identified unique and connected saddle points as compared to the standard ARTn, thus underlining the importance of correctly handling convex regions for completeness of saddle point explorations.

Structural and thermodynamical properties of early human amylin oligomers using replica exchange molecular dynamics: mutation effect of three key residues F15, H18 and F23.

Human islet amyloid polypeptide (hIAPP) is a 37-residue polypeptide, considered to be the main component of the pancreatic islet amyloid associated with type 2 diabetes and is one of the most amyloidogenic polypeptides known. Although the structure of hIAPP fibrils has already been obtained, structures of early oligomers and the mechanism of ß-sheet formation remain poorly understood. Herein, we characterize the atomic structure and the thermodynamics of the 14-37 residue fragment of hIAPP wild-type and mutated dimers and trimers. More precisely, three key residues, F15, H18 and F23, thought to affect the aggregation process, are the focus of this numerical study using replica exchange molecular dynamics coupled with the OPEP coarse-grained protein force field. Our simulations show that the oligomerization process takes place through the formation of anti-parallel ß-sheets most probably between C-terminal regions. Two main characteristics are associated with the onset of the fibrillation process. First, the sequence matching between the central (20-29) and C-terminal (30-37) regions, at intra and inter-molecular levels, helps in stabilizing the secondary structure and facilitates intermolecular interactions. Second, hydrophobic residues I26 and L27 are likely to promote the capture of further oligomeric structures and thus facilitate fibril elongation. Histidine mutation should have a more pronounced effect in the N-terminal region while phenylalanine mutations do not seem to prevent amyloid formation since central/C-terminal interactions are conserved through other contacts.

Strain-driven diffusion process during silicon oxidation investigated by coupling density functional theory and activation relaxation technique.

The reaction of oxygen molecules on an oxidized silicon model-substrate is investigated using an efficient potential energy hypersurface exploration that provides a rich picture of the associated energy landscape, energy barriers, and insertion mechanisms. Oxygen molecules are brought in, one by one, onto an oxidized silicon substrate, and accurate pathways for sublayer oxidation are identified through the coupling of density functional theory to the activation relaxation technique nouveau, an open-ended unbiased reaction pathway searching method, allowing full exploration of potential energy surface. We show that strain energy increases with O coverage, driving the kinetics of diffusion at the Si/SiO2 interface in the interfacial layer and deeper into the bulk: at low coverage, interface reconstruction dominates while at high coverage, oxygen diffusion at the interface or even deeper into the bottom layers is favored. A changing trend in energetics is observed that favors atomic diffusions to occur at high coverage while they appear to be unlikely at low coverage. Upon increasing coverage, strain is accumulated at the interface, allowing the oxygen atom to diffuse as the strain becomes large enough. The observed atomic diffusion at the interface releases the accumulated strain, which is consistent with a layer-by-layer oxidation growth.

Investigating the kinetics of the formation of a C Cottrell atmosphere around a screw dislocation in bcc iron: a mixed-lattice atomistic kinetic Monte-Carlo analysis.

We present a mixed-lattice atomistic kinetic Monte-Carlo algorithm (MLKMC) that integrates a rigid-lattice AKMC approach with the kinetic activation-relaxation technique (k-ART), an off-lattice/self-learning AKMC. This approach opens the door to study large and complex systems adapting the cost of identification and evaluation of transition states to the local environment. To demonstrate its capacity, MLKMC is applied to the problem of the formation of a C Cottrell atmosphere decorating a screw dislocation in α-Fe. For this system, transitions that occur near the dislocation core are searched by k-ART, while transitions occurring far from the dislocation are computed before the simulation starts using the rigid-lattice AKMC. This combination of the precision of k-ART and the speed of the rigid-lattice makes it possible to follow the onset of the C Cottrell atmosphere and to identify interesting mechanisms associated with its formation.

Interaction between interstitial carbon atoms and a ½ 〈1 1 1〉 self-interstitial atoms loop in an iron matrix: a combined DFT, off lattice KMC and MD study.

A static and kinetic study of the interaction between a 19 ½ 〈1 1 1〉 self-interstitial atoms loop and C atoms in body-centred cubic iron is presented in this work. An empirical potential matching the density functional theory calculations is used to study the static properties of the system. The usual kinetic Monte-Carlo (KMC) on-lattice restriction is not valid when the material is highly distorted, especially in the presence of a dislocation loop. Therefore, the dynamics of the system are investigated using both molecular dynamics simulations and k-ART, a self-learning/off-lattice atomic kinetic monte-carlo. The presented work is thus a full study of the C-loop and the C2-loop systems. A good agreement is observed between the statics and the kinetics (e.g. the discovery of a zone of stability of the C atom around the Fe cluster where the C can almost freely move), even though the kinetics show some unexpected behaviours of the studied systems. The pinning time of the loop induced by the C atoms is also estimated.

ARTIST: an activated method in internal coordinate space for sampling protein energy landscapes.

We present the first applications of an activated method in internal coordinate space for sampling all-atom protein conformations, the activation-relaxation technique for internal coordinate space trajectories (ARTIST). This method differs from all previous internal coordinate-based studies aimed at folding or refining protein structures in that conformational changes result from identifying and crossing well-defined saddle points connecting energy minima. Our simulations of four model proteins containing between 4 and 47 amino acids indicate that this method is efficient for exploring conformational space in both sparsely and densely packed environments, and offers new perspectives for applications ranging from computer-aided drug design to supramolecular assembly.

Increased basal expression of hepatic Cyp1a1 and Cyp1a2 genes in inbred mice selected for susceptibility to acetaminophen-induced hepatotoxicity.

Susceptibility to acetaminophen-induced hepatotoxicity was found to vary widely in an outbred colony of Swiss Webster mice. Some acetaminophen-treated male mice showed a significant elevation in serum levels of the hepatic enzyme alanine aminotransferase at a normally non-hepatotoxic oral dose. A selective breeding program over 17 generations produced inbred mice which were either susceptible or nonsusceptible to the hepatotoxic effects of acetaminophen. Liver microsomes from the susceptible group showed a statistically significant increase in the ability to metabolize acetaminophen to a reactive intermediate which covalently binds N-acetylcysteine. Microsomal cytochrome P450 activities associated with CYP1A2 (acetanilide 4-hydroxylation and methoxyresorufin O-demethylase) were significantly increased in the susceptible group. Ethoxyresorufin O-deethylase activity, associated with both CYP1A1 and CYP1A2, was also significantly elevated in this group. Further examination of both CYP1A isoforms revealed that hepatic CYP1A1 and CYP1A2 mRNA and protein levels were significantly elevated in animals from the susceptible group. In vivo caffeine 3-demethylation, which is associated with CYP1A2 activity, co-segregated with acetaminophen susceptibility and showed a significant positive correlation (r = 0.626, p < 0.005) with CYP1A2 mRNA expression in animals from both the susceptible and nonsusceptible groups. The co-segregation of elevated basal Cyp1a1 and CYP1a2 gene expression levels in animals selected for susceptibility to acetaminophen-induced hepatotoxicity suggested a common heritable basis for regulation of basal expression of both of these CYP1A isoforms. This was supported by the correlated expression of both CYP1A mRNAs within individual mice (r = 0.644, p < 0.02).

Expression of matrix metalloproteinases and their inhibitors in human brain tumors.

Sixty human brain tumors, classified according to the New World Health Organization (WHO) classification including, grade I schwannomas, meningiomas and pilocytic astrocytomas, grade II astrocytomas, grade III anaplastic astrocytomas, grade IV glioblastomas, grade III anaplastic oligodendrogliomas and grade IV glioblastomas and lung and melanoma metastases were analyzed for the expression of three matrix metalloproteinases (MMPs), two tissue inhibitors of MMPs (TIMPs) and for MMP activity. Some correlation was found between MMP expression and the degree of malignancy. Western blotting analysis revealed a more uniform pattern of distribution of MMP-2 (gelatinase A) than of MMP-9 (gelatinase B) and MMP-12 (metalloelastase) among tumors. MMP-9 levels were found to be significantly higher in grade III anaplastic astrocytomas and anaplastic oligodendrogliomas than those in grade I schwannomas and meningiomas. Anaplastic astrocytomas and Grade IV glioblastomas expressed significantly higher levels MMP-12 than grade I meningiomas. All sixty tumors showed a similar pattern of activity in zymography, proMMP-9 being the major species detected. Interestingly, TIMP-1 and TIMP-2 expression levels were especially low in tumors of grade II and grade III but significantly higher in tumors of grade I, particularly in schwannomas. Taken together, these data suggest that: 1) a balance between MMPs and TIMPs has an important role to play in human brain tumors; 2) TIMP expression may be valuable markers for tumor malignancy.

Perioperative variability of binding of lidocaine, quinidine, and propranolol after cardiac operations.

This study examined the effect of changes in plasma concentrations of total protein, albumin, alpha 1-acid glycoprotein, and free fatty acids occurring after heart operations on the protein binding of chemically basic drugs. Plasma protein and free fatty acid concentrations were measured simultaneously with in vitro determinations of the protein binding of lidocaine, quinidine, and propranolol: immediately before operation, immediately on weaning from cardiopulmonary bypass, on arrival in the recovery room, and 12, 24, 72, and 120 hours postoperatively. Initial decreases in the concentrations of all proteins were followed by a rise in alpha 1-acid glycoprotein to 254% of baseline at 72 to 120 hours. The free fractions of drug were initially increased to 168% of baseline for lidocaine, 206% for quinidine, and 200% for propranolol and fell progressively with time, reaching sustained troughs of 65% for lidocaine, 50% for quinidine, and 57% for propranolol at 72 to 120 hours. Regression analysis indicated a major influence of changing alpha 1-acid glycoprotein concentrations on free fractions of all three drugs, with a smaller effect of albumin that reached statistical significance only for lidocaine. There were no significant perioperative changes in plasma concentrations of free fatty acids when the in vitro effects of heparin were controlled. In conclusion, sequential changes in plasma protein concentrations after cardiac operations predictably alter the protein binding of lidocaine, quinidine, and propranolol and should be considered when interpreting total plasma drug concentrations.

Expression of matrix metalloproteinases and their inhibitors in human brain tumors.

Sixty human brain tumors, including grade I meningiomas, schwannomas, and pilocytic astrocytomas, grade II astrocytomas, grade III anaplastic astrocytomas and oligodendrogliomas, and grade IV glioblastomas and lung and melanoma metastases were analyzed for expression of four matrix metalloproteinases (MMPs), two tissue inhibitors of MMPs (TIMPs), and MMP activity. No marked correlation was found between MMP expression and the degree of malignancy. Western blotting analysis revealed a more uniform pattern of distribution of MMP-2 (gelatinase A) than of MMP-9 (gelatinase B) and MMP-12 (metalloelastase) among tumors. All 60 tumors showed a similar pattern of activity in zymography, MMP-2 being the major species detected. Interestingly, TIMP-1 and TIMP-2 expression levels were low in tumors of grade III but significantly higher in tumors of grade I, particularly schwannomas. Altogether, these data suggest that: (1) the balance between MMP-2 and TIMP-2 is important in human brain tumors; and (2) TIMP expression may be a valuable marker for tumor malignancy.

Matrix metalloproteinases and their inhibitors in human pituitary tumors.

OBJECTIVE: To determine the expression of matrix metalloproteinases (MMP)-1, -2, and -3 and the tissue inhibitors of metalloproteinases (TIMP)-1, -2, and -3 in 12 tissue samples from normal pituitary glands and in 28 human pituitary tumors ranging from Grade 0 to Grade IV, and to establish a correlation between the level of expression of MMPs and TIMPs and the tumor grade. METHODS: The expression of MMPs and TIMPs was determined by Western blotting. MMP activity was detected by gelatin zymography. RESULTS: MMPs were expressed in the majority of tumors, and their levels of expression were unrelated to tumor grade or to their invasive phenotype. Some correlation was observed between MMP activity detected by zymography and tumor grade. TIMP-2 and TIMP-3 were poorly expressed in high-grade tumors and strongly expressed in normal pituitary glands and in the majority of low-grade tumors. CONCLUSION: No correlation could be established between the invasive potential of tumors and MMP-1, -2, and -3 expression levels. Some correlation was observed between MMP activity detected by zymography and tumor grade. A good inverse correlation was observed between TIMP-2 and TIMP-3 expression levels and tumor grade. These data suggest that monitoring the expression of TIMP-2 and TIMP-3 or gelatinolytic activity could be of prognostic value.

Sampling activated mechanisms in proteins with the activation-relaxation technique.

The activated dynamics of proteins occur on time scales of milliseconds and longer. Standard all-atom molecular dynamics simulations are limited to much shorter times, of the order of tens of nanoseconds. Therefore, many activated mechanisms that are crucial for long-time dynamics will not be observed in such molecular dynamics simulation; different methods are required. Here, we describe in detail the activation-relaxation technique (ART) that generates directly activated mechanisms. The method is defined in the configurational energy landscape and defines moves in a two step fashion: (a) a configuration is first brought from a local minimum to a nearby first-order saddle point (the activation); and (b) the configuration is relaxed to a new metastable state (the relaxation). The method has already been applied to a wide range of problems in condensed matter, including metallic glasses, amorphous semiconductors and silica glass. We review the algorithm in detail, discuss some previously published results and present simulations of activated mechanisms for a two-helix bundle protein using an all-atom energy function.

Matrix proteinase inhibition by AE-941, a multifunctional antiangiogenic compound.

BACKGROUND: Matrix metalloproteinases (MMPs) play an important role in tissue remodeling under normal physiological and pathological conditions and are thus attractive targets for both diagnostic and therapeutic purposes. Here, we examined the effect of AE-941, an orally bioavailable standardized extract made of cartilage that shows significant antiangiogenic and antimetastatic properties in vivo, on the activity of various members of the MMP family. MATERIALS AND METHODS: The effect of AE-941 on the activity of MMPs was assessed by fluorimetric assays and by substrate gel zymography. RESULTS: AE-941 markedly inhibits the gelatinolytic activity of MMP-2 and to a lesser extent those of MMP-1, MMP-7, MMP-9 and MMP-13. AE-941 also inhibited the elastinolytic activities of MMP-2 and MMP-9 as well as MMP-12 (metalloelastase), porcine pancreatic elastase (PPE), and human leukocyte elastase (HLE). Western blot analysis revealed the presence within AE-941 of immunoreactive TIMP-like proteins, suggesting that these proteins may be at least partly responsible for the observed MMP inhibition. CONCLUSIONS: Taken together, these results demonstrate that AE-941 contains TIMP-like proteins that could be responsible for the specific inhibition of MMPs. Given the recent studies suggesting the presence within this compound of specific inhibitor(s) of endothelial cell proliferation, AE-941 appears as a pleotropic agent able to interfere with several biochemical steps leading to angiogenesis and to other physiopathological conditions. Since AE-941 is currently under Phase III clinical investigations, these findings are also of considerable importance for our understanding of its anticancer properties.

Quantitation of acetazolamide in plasma by high-performance liquid chromatography.

A method for estimating acetazolamide concentrations in human plasma is described. Buffered plasma (pH 4.8) containing chlorothiazide as an internal standard is extracted twice with ethyl acetate. The extract is evaporated, redissolved, and chromatographed on silica gel with hexane-chloroform-methanol-acetic acid (65:25:10:0.25) as the mobile phase. The extraction efficiencies were greater than 90%, the coefficients of variation at 1 and 30 micrograms/ml of plasma were 3.5 and 2.0%, respectively, and the calibration curves were linear and had an intercept of essentially zero. The suitability of the method for pharmacokinetic studies was verified in a normal volunteer dosed with 250-mg (solution) and 500-mg (sustained-release tablet) acetazolamide formulations.

Chromatographic and spectroscopic characterization of sulphur-bound dimetridazole and ranidazole derivatives.

The 5-nitroimidazoles, dimetridazole and ronidazole, two important veterinary drugs, were reacted under reductive conditions with the sulfhydryl-containing substrates cysteine and glutathione to yield 5-amino-4-S-substituted imidazoles. After purification by reversed-phase liquid chromatography (RP-LC), the four adducts were characterized by RP-LC with photodiode array detection using conditions where their parent drugs were not eluted from the column. Structural identification was conducted by spectroscopic techniques, mainly 1-dimensional and 2-dimensional NMR. While the dimetridazole adducts were found to be monosubstituted at the C-4 position, the two ronidazole products contained two units of the sulfhydryl substrate, located at the C-4 and C-6 positions.