Biomedical subjects

Marco De Vivo

Publications and source records attributed to Marco De Vivo.

9 recordsLinked to original sources

PubMed2007

Proton shuttles and phosphatase activity in soluble epoxide hydrolase.

Recently, a novel metal Mg2+-dependent phosphatase activity has been discovered in the N-terminal domain of the soluble epoxide hydrolase (sEH), opening a new branch of fatty acid metabolism and providing an additional site for drug targeting. Importantly, the sEH N-terminal fold belongs to the haloacid dehalogenase (HAD) superfamily, which comprises a vast majority of phosphotransferases. Herein, we present the results of a computational study of the sEH phosphatase activity, which includes classical molecular dynamics (MD) simulations and mixed quantum mechanical/molecular mechanics (QM/MM) calculations. On the basis of experimental results, a two-step mechanism has been proposed and herein investigated: (1) phosphoenzyme intermediate formation and (2) phosphoenzyme intermediate hydrolysis. Building on our earlier work, we now provide a detailed description of the reaction mechanism for the whole catalytic cycle along with its free energy profile. The present computations suggest metaphosphate-like transition states for these phosphoryl transfers. They also reveal that the enzyme promotes water deprotonation and facilitates shuttling of protons via a metal-ligand connecting water bridge (WB). These WB-mediated proton shuttles are crucial for the activation of the solvent nucleophile and for the stabilization of the leaving group. Moreover, due to the conservation of structural features in the N-terminal catalytic site of sEH and other members of the HAD superfamily, we suggest a generalization of our findings to these other metal-dependent phosphatases.

Catalysis↗

PubMed2006

Modeling the charge distribution at metal sites in proteins for molecular dynamics simulations.

Almost half of the proteome of living organisms is constituted of metalloproteins. Unfortunately, the ability of the current generation of molecular dynamics pairwise-additive forcefields to properly describe metal pockets is severely lacking due to the intrinsic difficulty of handling polarization and charge transfer contributions. In order to improve the description of metalloproteins, a simple reparameterization strategy is proposed herein that does not involve artificial constraints. Specifically, a non-bonded quantum mechanical-based model is used to capture the mean polarization and charge transfer contributions to the interatomic forces within the metal site. The present approach is demonstrated to provide enough accuracy to maintain the integrity of the metal pocket for a variety of metalloproteins during extended (multi-nanosecond) molecular dynamics simulations. The method enables the sampling of small conformational changes and the relaxation of local frustrations in NMR structures.

Binding Sites↗

PubMed2006

Role of phosphorylated Thr160 for the activation of the CDK2/Cyclin A complex.

The enzymatic activity of the CDK2/Cyclin A complex increases upon the specific phosphorylation of Thr160@CDK2. In the present study, we have performed a comparative molecular dynamics (MD) study of models of the complex CDK2/Cyclin A/Substrate, which differ for the presence or absence of the phosphate group bound to Thr160. The models are based on two X-ray structures available for CDK2/CyclinA and pCDK2/CyclinA/Substrate complexes. In this way, we analyze the influence of the phosphorylated Thr160 (pThr160) on both the flexibility of CDK2 activation loop (AL) and substrate binding in CDK2. Our calculations point to a decreased flexibility of the AL in the phosphorylated model, in fairly good agreement with experimental data, and to a key role of pThr160 for substrate recognition and stability. Multiple alignments of the CDKs sequences point to the very high conservation of the AL sequence among the CDKs, thus extending our results to all CDKs.

Binding Sites↗

PubMed2006

Metadynamics as a tool for exploring free energy landscapes of chemical reactions.

The metadynamics or hills method is a relatively new molecular dynamics technique aimed to enhance the sampling of separated regions in phase space and map out the underlying free energy landscape as a function of a small number of order parameters or collective variables. The high efficiency allows for the application of metadynamics in combination with first principles dynamics methods, in particular with Car-Parrinello molecular dynamics, to study processes in which changes in the electronic structure play a dominant role, such as chemical reactions. The option to choose several independent collective variables is important to tackle complex and concerted transformations that lack an obvious a priori choice for a single reaction coordinate. In this Account, we discuss the role of metadynamics in the search of transition states, local minima, reaction paths, free energy profiles, and reaction coordinates among a growing list of alternative methods.

Journal Article↗

PubMed2005

Computational study of phosphatase activity in soluble epoxide hydrolase: high efficiency through a water bridge mediated proton shuttle.

Recently, a new branch of fatty acid metabolism has been opened by the novel phosphatase activity found in the N-terminal domain of the, hence bifunctional, soluble epoxide hydrolase (sEH). Importantly, this finding has also provided a new site for drug targeting in sEH's activity regulation. Classical MD and hybrid Car-Parrinello QM/MM calculations have been performed to investigate the reaction mechanism of the phosphoenzyme intermediate formation in the first step of the catalysis. The results support a concerted multi-event reaction mechanism: (1) a dissociative in-line nucleophilic substitution for the phosphoryl transfer reaction; (2) a double proton transfer involved in the formation of a good leaving group in the transition state. The presence of a water bridge in the substrate/enzyme complex allowed an efficient proton shuttle, showing its key role in speeding up the catalysis. The calculated free energy of the favored catalytic pathway is approximately 19 kcal/mol, in excellent agreement with experimental data.

Animals↗

PubMed2005

A theoretical DFT investigation of the lysozyme mechanism: computational evidence for a covalent intermediate pathway.

A theoretical DFT(B3LYP) investigation of the catalytic cycle of lysozyme has provided further evidence for a mechanism involving a glycosil-enzyme covalent intermediate, in agreement with recent experimental data. This type of intermediate has been located along two different pathways. Along the favored path the retention of the anomeric configuration of the peptidoglycan NAM unit involved in the reaction, is the result of two subsequent inversions at the C(1) carbon. The other path involves the opening of the pyranose ring and a nucleophilic attack on the prochiral carbonyl group of the open aldehyde, restoring the original anomeric configuration. No evidence has been found for a pathway characterized by the formation of an oxocarbenium ion (stabilized by resonance and electrostatic interactions), as suggested in the most popular mechanistic schemes.

Catalysis↗

PubMed2004

A computational study of the binding of propidium to the peripheral anionic site of human acetylcholinesterase.

Combined docking and molecular dynamics (MD) simulations were carried out in order to investigate the binding mode of propidium at the human acetylcholinesterase (HuAChE) peripheral site. Two different docking protocols followed by cluster analyses were performed, allowing the identification of five high-populated and low-energy configuration families. To dynamically explore the behavior of the ligand at the peripheral HuAChE binding site, six complexes (five low-energy and one high-energy) were submitted to 2.5 ns of MD simulations. The representative propidium/HuAChE binding modes were chosen on the basis of both the docking energy score and the dynamic stability of the complexes throughout the MD simulations. The most stable poses of propidium at HuAChE PAS were similar to those experimentally determined with the murine enzyme. We therefore suggest that the present modeling protocol might be used in the dynamic investigation of the interactions of a small-molecule inhibitor with a surface-like binding site of a protein. Finally, because of the biological relevance of the target studied here, the present results can be of interest for the rational design of molecules potentially useful in the treatment of the Alzheimer's disease.

Acetylcholinesterase↗

PubMed2004

Synthesis, biological activity, and three-dimensional quantitative structure-activity relationship model for a series of benzo[c]quinolizin-3-ones, nonsteroidal inhibitors of human steroid 5alpha-reductase 1.

New 5alpha-reductase 1 (5alphaR-1) inhibitors were designed to complete a consistent set of analogues suitable for a 3D QSAR study. These compounds were synthesized by a modification of the aza-Robinson annulation, further functionalized by Pd-catalyzed cross-coupling processes, and were tested with human 5alphaR-1 expressed in Chinese hamster ovary 1827 cells. It turned out that the potency of the resulting inhibitors was strongly dependent on the type of substitution at the 8 position, with the IC(50) values ranging from 8.1 to 1050 nM. The construction of this homogeneous set of molecules allowed a 3D QSAR study. In particular, comparative molecular field analysis (CoMFA) was used to correlate the potency of the inhibitors with their physicochemical features. Highly accurate evaluations of the atomic point charges were carried out by means of quantum chemical calculations at the DFT/B3LYP level of theory followed by the RESP fitting procedure. It turned out that increasing the reliability of electrostatic parameters greatly affected the statistical results of the QSAR analysis. The 3D QSAR model proposed could be very useful in the further development of 5alphaR-1 inhibitors, which are suitable candidates to be evaluated as drugs in the treatment of 5alphaR-1 related diseases such as acne and alopecia in men and hirsutism in women.

3-Oxo-5-alpha-Steroid 4-Dehydrogenase↗

PubMed2003

Density functional study of the enzymatic reaction catalyzed by a cyclin-dependent kinase.

Density functional theory (DFT) calculations were carried out to study the molecular mechanism of the phosphoryl transfer reaction catalyzed by cyclin-dependent kinases (CDKs). The DFT study presented here shows that CDKs catalyze the phosphoryl transfer reaction from ATP to the serine substrate through a single step mechanism with a SN2-like transition state.

Catalysis↗