Origin of functional diversity among tetrameric voltage-gated channels.

The aim of the present work is to relate functional differences of voltage-gated K(+) (K(v)), hyperpolarization-activated cyclic nucleotide-gated (HCN), and cyclic nucleotide gated (CNG) channels to differences in their amino acid sequences. By means of combined bioinformatic sequence analyses and homology modelling, we suggest that: (1) CNG channels are less voltage-dependent than K(v) channels since the charge of their voltage sensor, the S4 helix, is lower than that of K(v) channels and because of the presence of a conserved proline in the S4-S5 linker, which is quite likely to uncouple S4 from S5 and S6. (2) In HCN channels, S4 features a higher net positive charge with respect to K(v) channels and an extensive network of hydrophobic residues, which is quite likely to provide a tight coupling among S4 and the neighboring helices.

Using metadynamics to understand the mechanism of calmodulin/target recognition at atomic detail.

The ability of calcium-bound calmodulin (CaM) to recognize most of its target peptides is caused by its binding to two hydrophobic residues ('anchors'). In most of the CaM complexes, the anchors pack against the hydrophobic pockets of the CaM domains and are surrounded by fully conserved Met side chains. Here, by using metadynamics simulations, we investigate quantitatively the energetics of the final step of this process using the M13 peptide, which has a high affinity and spans the sequence of the skeletal myosin light chain kinase, an important natural CaM target.

Convergent dynamics in the protease enzymatic superfamily.

Proteases regulate various aspects of the life cycle in all organisms by cleaving specific peptide bonds. Their action is so central for biochemical processes that at least 2% of any known genome encodes for proteolytic enzymes. Here we show that selected proteases pairs, despite differences in oligomeric state, catalytic residues, and fold, share a common structural organization of functionally relevant regions which are further shown to undergo similar concerted movements.

Molecular simulation of the binding of nerve growth factor peptide mimics to the receptor tyrosine kinase A.

Nerve growth factor (NGF) mimics play an important role for therapies that target the receptor tyrosine kinase A (trkA). The N-terminal fragment of the NGF (N-term@NGF) was previously demonstrated to be an important determinant for affinity and specificity in the binding to trkA. Here we use a variety of computational tools (contact surface analysis and free energy predictions) to identify residues playing a key role for the binding to the receptor.

Polarization effects and charge transfer in the KcsA potassium channel.

The electronic structure of the selectivity filter of KcsA K(+) channel is investigated by density functional theory (DFT/BLYP) and QM/MM methods. The quantum part includes the selectivity filter, which is polarized by the electrostatic field of the environment treated with the Amber force field. The details of the electronic structure were investigated using the maximally localized Wannier function centers of charge and Bader's atoms in molecules charge analysis.

Deuterium isotope effects in A:T and A:U base pairs: a computational NMR study.

Recent measurements of trans-hydrogen bond deuterium isotope effects (DIEs) on 13C chemical shifts in nucleic acids (Vakonakis, I.; LiWang, A. C.

cAMP Modulation of the cytoplasmic domain in the HCN2 channel investigated by molecular simulations.

The hyperpolarization-activated cyclic nucleotide-modulated (HCN) cation channels are opened by membrane hyperpolarization, while their activation is modulated by the binding of cyclic adenosine monophosphate (cAMP) in the cytoplasm. Here we investigate the molecular basis of cAMP channel modulation by performing molecular dynamics simulations of a segment comprising the C-linker and the cyclic nucleotide binding domain (CNBD) in the presence and absence of cAMP, based on the available crystal structure of HCN2 from mouse. In presence of cAMP, the protein undergoes an oscillation of the quaternary structure on the order of 10 ns, not observed in the apoprotein.

Structure and function of vanadium haloperoxidases.

A quantum mechanics/molecular mechanics study of the resting state of the vanadium dependent chloroperoxidase from fungi Curvularia inaequalis and of the early intermediates of the halide oxidation is reported. The investigation of different protonation states indicates that the enzyme likely consists of an anionic H2VO4- vanadate moiety where one hydroxo group is in axial position. The calculations suggest that the hydrogen peroxide binding may not involve an initial protonation of the vanadate cofactor.

Duocarmycins binding to DNA investigated by molecular simulation.

Duocarmycins are a potent class of antitumor agents, whose activity arises through their covalent binding to adenine nucleobases of DNA.(1-3) Here, we perform molecular dynamics (MD) and hybrid Car-Parinello QM/MM simulations to investigate aspects of duocarmycin binding to the d(pGpApCpTpApApTpTpGpApC) oligonucleotide. We focus on the derivatives (+)-duocarmycin SA (DSA) and (+)-duocarmycin SI (DSI), for which structural information of the covalent complex with the oligonucleotide is available, as well as on the related, but less reactive, NBOC-duocarmycin SA (NBOC-DSA), interacting with the same oligonucleotide.

Binding of novel azole-bridged dinuclear platinum(II) anticancer drugs to DNA: insights from hybrid QM/MM molecular dynamics simulations.

Dinuclear Pt-containing compounds might be used to overcome the intrinsic and acquired cell resistance of widely used anticancer drugs such as cisplatin. Recently, the complexes [[cis-Pt(NH3)2]2(mu-OH)(mu-pz)](NO3)2 (with pz = pyrazolate) (1), [[cis-Pt(NH3)2]2(mu-OH)(mu-1,2,3-ta-N(1),N(2))](NO3)2 (with ta = 1,2,3-triazolate) (2), and the binding of 1 to d(CpTpCpTpG*pG*pTpCpTpCp) have been characterized. Here we provide the structural and electronic properties of the free drugs, of the intermediates of binding to guanine bases, and of the products, by performing DFT calculations.

Binding of phosphinate and phosphonate inhibitors to aspartic proteases: a first-principles study.

Phosphinate and phosphonate derivatives are potent inhibitors of aspartic proteases (APs). The affinity for the enzyme might be caused by the presence of low barrier hydrogen bonds between the ligand and the catalytic Asp dyad in the cleavage site. We have used density functional theory calculations along with hybrid quantum mechanics/molecular mechanics Car-Parrinello molecular dynamics simulations to investigate the hydrogen-bonding pattern at the binding site of the complexes of human immunodeficiency virus type-1 AP and the eukaryotic endothiapepsin and penicillopepsin.

Changes in ultraviolet absorbance and hence in protective efficacy against lipid peroxidation of organic sunscreens after UVA irradiation.

Owing to the spectral distribution of solar UV, the UVA component of sunlight is now believed to be the main cause of photoaging and photocarcinogenesis and is much more effective than UVB in inducing peroxidative damage. Consequently, most skin care cosmetic products now include UVA filters in their formulations along with UVB filters. These modern sunscreens should provide and maintain their initial absorbance, hence protection, throughout the entire period of exposure to sunlight.

Coarse-grained model of proteins incorporating atomistic detail of the active site.

We present a novel approach to explore the conformational space of globular proteins near their native state. It combines the advantages of coarse-grained models with those of all-atoms simulations, required to treat molecular recognition processes. The comparison between calculated structural properties with those obtained with all-atoms molecular dynamics simulations establishes the accuracy of the model.

Insights on HIV-1 Tat:P/CAF bromodomain molecular recognition from in vivo experiments and molecular dynamics simulations.

Structural and functional studies indicate that, through its bromodomain, the cellular acetyltransferase P/CAF binds the acetylated Tat protein of human immunodeficiency virus type 1 (HIV-1) and promotes transcriptional activation of the integrated provirus. Based on the NMR structure of P/CAF complexed with an acetylated Tat peptide, here we use molecular dynamics simulations to construct a model describing the interaction between full length Tat and the P/CAF bromodomain. Our calculations show that the protein-protein interface involves hydrophobic interactions between the P/CAF ZA loop and the Tat core domain.

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.

Unwinding the helical linker of calcium-loaded calmodulin: a molecular dynamics study.

The fold of calmodulin (CaM) consists of two globular domains connected by a helical segment (the linker), whose conformational properties play a crucial role for the protein's molecular recognition processes. Here we investigate the structural properties of the linker by performing a 11.5 ns molecular dynamics (MD) simulation of calcium-loaded human CaM in aqueous solution.

A homology model of the pore region of HCN channels.

HCN channels are activated by membrane hyperpolarization and regulated by cyclic nucleotides, such as cyclic adenosine-mono-phosphate (cAMP). Here we present structural models of the pore region of these channels obtained by using homology modeling and validated against spatial constraints derived from electrophysiological experiments. For the construction of the models we make two major assumptions, justified by electrophysiological observations: i), in the closed state, the topology of the inner pore of HCN channels is similar to that of K(+) channels.

Molecular basis of the allosteric mechanism of cAMP in the regulatory PKA subunit.

The second messenger cyclic Adenosine MonoPosphate (cAMP) mediates many biological process by interacting with structurally conserved nucleotide binding domains (cNBD's). Here, we use molecular dynamics simulations on RIIbeta-PKA, one of the best characterized members of the cNBD family, in presence and absence of cAMP. The results of our calculations are fully consistent with the available experimental data and suggest that the key factor of the cAMP allosteric mechanism in cNBDS's is the increased flexibility of the protein upon ligand release along with a mechanical coupling between helical segments.

Structural basis of gating of CNG channels.

Cyclic nucleotide-gated (CNG) ion channels, underlying sensory transduction in vertebrate photoreceptors and olfactory sensory neurons, require cyclic nucleotides to open. Here, we present structural models of the tetrameric CNG channel pore from bovine rod in both open and closed states, as obtained by combining homology modeling-based techniques, experimentally derived spatial constraints and structural patterns present in the PDB database. Gating is initiated by an anticlockwise rotation of the N-terminal region of the C-linker, which is then, transmitted through the S6 transmembrane helices to the P-helix, and in turn from this to the pore lumen, which opens up from 2 to 5A thus allowing for ion permeation.

Evolutionarily conserved functional mechanics across pepsin-like and retroviral aspartic proteases.

The biological function of the aspartic protease from HIV-1 has recently been related to the conformational flexibility of its structural scaffold. Here, we use a multistep strategy to investigate whether the same mechanism affects the functionality in the pepsin-like fold. (i) We identify the set of conserved residues by using sequence-alignment techniques.

Comparative analysis of HIV-1 Tat variants.

HIV-1 Tat protein is a crucial element for viral replication; therefore, its inhibition might be exploited against the AIDS infection. To gain insights on the natural variability of this protein, we present a comparative investigation on the relationship between the primary sequences and the experimentally available three-dimensional structures from the HIV-1 Tat variants Z2, BRU, and MAL. Our computational tools include sequence conservation algorithms, structural analysis, electrostatic modeling, and molecular dynamics (MD) simulations.

A proficient enzyme: insights on the mechanism of orotidine monophosphate decarboxylase from computer simulations.

Decarboxylation of orotidine 5'-monophosphate (Omp) to uridine 5'-monophosphate by orotidine 5'-monophosphate decarboxylase (ODCase) is currently the object of vivid debate. Here, we clarify its enzymatic activity with long time scale classical molecular dynamics and hybrid ab initio Car-Parrinello/molecular mechanics simulations. The lack of structural (experimental) information on the ground state of ODCase/Omp complex is overcome by a careful construction of the model and the analysis of three different strains of the enzyme.

Water-assisted reaction mechanism of monozinc beta-lactamases.

Hybrid Car-Parrinello QM/MM calculations are used to investigate the reaction mechanism of hydrolysis of a common beta-lactam substrate (cefotaxime) by the monozinc beta-lactamase from Bacillus cereus (BcII). The calculations suggest a fundamental role for an active site water in the catalytic mechanism. This water molecule binds the zinc ion in the first step of the reaction, expanding the zinc coordination number and providing a proton donor adequately oriented for the second step.