Subdiffusive-Brownian crossover in membrane proteins: a generalized Langevin equation-based approach.

In this work, we propose a generalized Langevin equation-based model to describe the lateral diffusion of a protein in a lipid bilayer. The memory kernel is represented in terms of a viscous (instantaneous) and an elastic (noninstantaneous) component modeled through a Dirac δ function and a three-parameter Mittag-Leffler type function, respectively. By imposing a specific relationship between the parameters of the three-parameter Mittag-Leffler function, the different dynamical regimes-namely ballistic, subdiffusive, and Brownian, as well as the crossover from one regime to another-are retrieved. Within this approach, the transition time from the ballistic to the subdiffusive regime and the spectrum of relaxation times underlying the transition from the subdiffusive to the Brownian regime are given. The reliability of the model is tested by comparing the mean-square displacement derived in the framework of this model and the mean-square displacement of a protein diffusing in a membrane calculated through molecular dynamics simulations.

Hydrodynamics of immiscible binary fluids with viscosity contrast: a multiparticle collision dynamics approach.

We present a multiparticle collision dynamics (MPC) implementation of layered immiscible fluids A and B of different shear viscosities separated by planar interfaces. The simulated flow profile for imposed steady shear motion and the time-dependent shear stress functions are in excellent agreement with our continuum hydrodynamics results for the composite fluid. The wave-vector dependent transverse velocity auto-correlation functions (TVAF) in the bulk-fluid regions of the layers decay exponentially, and agree with those of single-phase isotropic MPC fluids. In addition, we determine the hydrodynamic mobilities of an embedded colloidal sphere moving steadily parallel or transverse to a fluid-fluid interface, as functions of the distance from the interface. The obtained mobilities are in good agreement with hydrodynamic force multipoles calculations, for a no-slip sphere moving under creeping flow conditions near a clean, ideally flat interface. The proposed MPC fluid-layer model can be straightforwardly implemented, and it is computationally very efficient. Yet, owing to the spatial discretization inherent to the MPC method, the model can not reproduce all hydrodynamic features of an ideally flat interface between immiscible fluids.

Open-Boundary Molecular Mechanics/Coarse-Grained Framework for Simulations of Low-Resolution G-Protein-Coupled Receptor-Ligand Complexes.

G-protein-coupled receptors (GPCRs) constitute as much as 30% of the overall proteins targeted by FDA-approved drugs. However, paucity of structural experimental information and low sequence identity between members of the family impair the reliability of traditional docking approaches and atomistic molecular dynamics simulations for in silico pharmacological applications. We present here a dual-resolution approach tailored for such low-resolution models. It couples a hybrid molecular mechanics/coarse-grained (MM/CG) scheme, previously developed by us for GPCR-ligand complexes, with a Hamiltonian-based adaptive resolution scheme (H-AdResS) for the solvent. This dual-resolution approach removes potentially inaccurate atomistic details from the model while building a rigorous statistical ensemble-the grand canonical one-in the high-resolution region. We validate the method on a well-studied GPCR-ligand complex, for which the 3D structure is known, against atomistic simulations. This implementation paves the way for future accurate in silico studies of low-resolution ligand/GPCRs models.

Molecular Dynamics Simulations of the [2Fe-2S] Cluster-Binding Domain of NEET Proteins Reveal Key Molecular Determinants That Induce Their Cluster Transfer/Release.

The NEET proteins are a novel family of iron-sulfur proteins characterized by an unusual three cysteine and one histidine coordinated [2Fe-2S] cluster. Aberrant cluster release, facilitated by the breakage of the Fe-N bond, is implicated in a variety of human diseases, including cancer. Here, the molecular dynamics in the multi-microsecond timescale, along with quantum chemical calculations, on two representative members of the family (the human NAF-1 and mitoNEET proteins), show that the loss of the cluster is associated with a dramatic decrease in secondary and tertiary structure. In addition, the calculations provide a mechanism for cluster release and clarify, for the first time, crucial differences existing between the two proteins, which are reflected in the experimentally observed difference in the pH-dependent cluster reactivity. The reliability of our conclusions is established by an extensive comparison with the NMR data of the solution proteins, in part measured in this work.

Open Boundary Simulations of Proteins and Their Hydration Shells by Hamiltonian Adaptive Resolution Scheme.

The recently proposed Hamiltonian adaptive resolution scheme (H-AdResS) allows the performance of molecular simulations in an open boundary framework. It allows changing, on the fly, the resolution of specific subsets of molecules (usually the solvent), which are free to diffuse between the atomistic region and the coarse-grained reservoir. So far, the method has been successfully applied to pure liquids. Coupling the H-AdResS methodology to hybrid models of proteins, such as the molecular mechanics/coarse-grained (MM/CG) scheme, is a promising approach for rigorous calculations of ligand binding free energies in low-resolution protein models. Toward this goal, here we apply for the first time H-AdResS to two atomistic proteins in dual-resolution solvent, proving its ability to reproduce structural and dynamic properties of both the proteins and the solvent, as obtained from atomistic simulations.

Computational metallomics of the anticancer drug cisplatin.

Cisplatin, cis-diamminedichlorido-platinum(II), is an important therapeutic tool in the struggle against different tumors, yet it is plagued with the emergence of resistance mechanisms after repeated administrations. This hampers greatly its efficacy. Overcoming resistance problems requires first and foremost an integrated and systematic understanding of the structural determinants and molecular recognition processes involving the drug and its cellular targets. Here we review a strategy that we have followed for the last few years, based on the combination of modern tools from computational chemistry with experimental biophysical methods. Using hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) simulations, validated by spectroscopic experiments (including NMR, and CD), we have worked out for the first time at atomic level the structural determinants in solution of platinated cellular substrates. These include the copper homeostasis proteins Ctr1, Atox1, and ATP7A. All of these proteins have been suggested to influence the pre-target resistance mechanisms. Furthermore, coupling hybrid QM/MM simulations with classical Molecular Dynamics (MD) and free energy calculations, based on force field parameters refined by the so-called "Force Matching" procedure, we have characterized the structural modifications and the free energy landscape associated with the recognition between platinated DNA and the protein HMGB1, belonging to the chromosomal high-mobility group proteins HMGB that inhibit the repair of platinated DNA. This may alleviate issues relative to on-target resistance process. The elucidation of the mechanisms by which tumors are sensitive or refractory to cisplatin may lead to the discovery of prognostic biomarkers. The approach reviewed here could be straightforwardly extended to other metal-based drugs.

Adaptation of Extremophilic Proteins with Temperature and Pressure: Evidence from Initiation Factor 6.

In this work, we study dynamical properties of an extremophilic protein, Initiation Factor 6 (IF6), produced by the archeabacterium Methanocaldococcus jannascii, which thrives close to deep-sea hydrothermal vents where temperatures reach 80 °C and the pressure is up to 750 bar. Molecular dynamics simulations (MD) and quasi-elastic neutron scattering (QENS) measurements give new insights into the dynamical properties of this protein with respect to its eukaryotic and mesophilic homologue. Results obtained by MD are supported by QENS data and are interpreted within the framework of a fractional Brownian dynamics model for the characterization of protein relaxation dynamics. IF6 from M. jannaschii at high temperature and pressure shares similar flexibility with its eukaryotic homologue from S. cerevisieae under ambient conditions. This work shows for the first time, to our knowledge, that the very common pattern of corresponding states for thermophilic protein adaptation can be extended to thermo-barophilic proteins. A detailed analysis of dynamic properties and of local structural fluctuations reveals a complex pattern for "corresponding" structural flexibilities. In particular, in the case of IF6, the latter seems to be strongly related to the entropic contribution given by an additional, C-terminal, 20 amino-acid tail which is evolutionary conserved in all mesophilic IF6s.

Structural predictions of neurobiologically relevant G-protein coupled receptors and intrinsically disordered proteins.

G protein coupled receptors (GPCRs) and intrinsic disordered proteins (IDPs) are key players for neuronal function and dysfunction. Unfortunately, their structural characterization is lacking in most cases. From one hand, no experimental structure has been determined for the two largest GPCRs subfamilies, both key proteins in neuronal pathways. These are the odorant (450 members out of 900 human GPCRs) and the bitter taste receptors (25 members) subfamilies. On the other hand, also IDPs structural characterization is highly non-trivial. They exist as dynamic, highly flexible structural ensembles that undergo conformational conversions on a wide range of timescales, spanning from picoseconds to milliseconds. Computational methods may be of great help to characterize these neuronal proteins. Here we review recent progress from our lab and other groups to develop and apply in silico methods for structural predictions of these highly relevant, fascinating and challenging systems.

Structural biology of cisplatin complexes with cellular targets: the adduct with human copper chaperone atox1 in aqueous solution.

Cisplatin is one of the most used anticancer drugs. Its cellular influx and delivery to target DNA may involve the copper chaperone Atox1 protein. Although the mode of binding is established by NMR spectroscopy measurements in solution-the Pt atom binds to Cys12 and Cys15 while retaining the two ammine groups-the structural determinants of the adduct are not known. Here a structural model by hybrid Car-Parrinello density functional theory-based QM/MM simulations is provided. The platinated site minimally modifies the fold of the protein. The calculated NMR and CD spectral properties are fully consistent with the experimental data. Our in silico/in vitro approach provides, together with previous studies, an unprecedented view into the structural biology of cisplatin-protein adducts.

Platination of the copper transporter ATP7A involved in anticancer drug resistance.

The clinical efficacy of the widely used anticancer drug cisplatin is severely limited by the emergence of resistance. This is related to the drug binding to proteins such as the copper influx transporter Ctr1, the copper chaperone Atox1, and the copper pumps ATP7A and ATP7B. While the binding modes of cisplatin to the first two proteins are known, the structural determinants of platinated ATP7A/ATP7B are lacking. Here we investigate the interaction of cisplatin with the first soluble domain of ATP7A. First, we establish by ESI-MS and (1)H, (13)C, and (15)N NMR that, in solution, the adduct is a monomer in which the sulfur atoms of residues Cys19 and Cys22 are cis-coordinated to the [Pt(NH3)2](2+) moiety. Then, we carry out hybrid Car-Parrinello QM/MM simulations and computational spectroscopy calculations on a model adduct based on the NMR structure of the apo protein and featuring the experimentally determined binding mode of the metal ion. These calculations show quantitative agreement with CD spectra and (1)H, (13)C, and (15)N NMR chemical shifts, thus providing a quantitative molecular view of the 3D binding mode of cisplatin to ATP7A. Importantly, the same comparison rules out a variety of alternative models with different coordination modes, that we explored to test the robustness of the computational approach. Using this combined in silico-in vitro approach we provide here for the first time a quantitative 3D atomic view of the platinum binding to the first soluble domain of ATP7A.

From NMR relaxation to fractional Brownian dynamics in proteins: results from a virtual experiment.

In a recent simulation study [J. Chem. Phys. 2010, 133, 145101], it has been shown that the time correlation functions probed by nuclear magnetic resonance (NMR) relaxation spectroscopy of proteins are well described by a fractional Brownian dynamics model, which accounts for the wide spectrum of relaxation rates characterizing their internal dynamics. Here, we perform numerical experiments to explore the possibility of using this model directly in the analysis of experimental NMR relaxation data. Starting from a molecular dynamics simulation of the 266 residue protein 6PGL in explicit water, we construct virtual (15)N R(1), R(2), and NOE relaxation rates at two different magnetic fields, including artificial noise, and test how far the parameters obtained from a fit of the model to the virtual experimental data coincide with those obtained from an analysis of the MD time correlation functions that have been used to construct these data. We show that in most cases, close agreement is found. Acceptance or rejection of parameter values obtained from relaxation rates are discussed on a physical basis, therefore avoiding overfitting.

Conformational transition of DNA bound to Hfq probed by infrared spectroscopy.

Hfq is a bacterial protein involved in RNA metabolism. Besides this, Hfq's role in DNA restructuring has also been suggested. Since this mechanism remains unclear, we examined the DNA conformation upon Hfq binding by combining vibrational spectroscopy and neutron scattering. Our analysis reveals that Hfq, which preferentially interacts with deoxyadenosine rich sequences, induces partial opening of dA-dT sequences accompanied by sugar repuckering of the dA strand and hence results in a heteronomous A/B duplex. Sugar repuckering is probably correlated with a global dehydration of the complex. By taking into account Hfq's preferential binding to A-tracts, which are commonly found in promoters, potential biological implications of Hfq binding to DNA are discussed.

Fractional protein dynamics seen by nuclear magnetic resonance spectroscopy: Relating molecular dynamics simulation and experiment.

We propose a fractional Brownian dynamics model for time correlation functions characterizing the internal dynamics of proteins probed by NMR relaxation spectroscopy. The time correlation functions are represented by a broad distribution of exponential functions which are characterized by two parameters. We show that the model describes well the restricted rotational motion of N-H vectors in the amide groups of lysozyme obtained from molecular dynamics simulation and that reliable predictions of experimental relaxation rates can be obtained on that basis.

Protein dynamics from a NMR perspective: networks of coupled rotators and fractional Brownian dynamics.

Nuclear magnetic resonance (NMR) has proven to be the most valuable tool for investigating internal dynamics of proteins. In this perspective, the interpretation of NMR relaxation data eventually relies on a model of the motions. In this article, we propose to compare two radically different approaches that aim at describing internal dynamics in proteins. It is shown that the correlation functions predicted by a network of coupled rotators can be interpreted in terms of a heuristic approach based on fractional Brownian dynamics for each of the vectors in the network. Our results are interpreted in terms of the probability distributions of relaxation modes in both processes, the median of which turns out to be the relevant quantity for the comparison of both models.

Influence of pressure on the slow and fast fractional relaxation dynamics in lysozyme: a simulation study.

The article reports on a molecular dynamics simulation study of the influence of moderate, nondenaturing pressure on the slow and fast internal relaxation dynamics of lysozyme. The model parameters of the fractional Ornstein-Uhlenbeck process are used to quantify the changes. We find that the nonexponential character for diffusive motions on time scales above 10 ps is enhanced and that the diffusion processes are slowed down. The diffusive motions on the subpicosecond time scale appear, in contrast, accelerated, whereas the nonexponential character is not altered by pressure. We attribute these findings to the different natures of slow and fast relaxation processes, which are characterized by structural rearrangements and collisions, respectively. The analyses are facilitated by the use of spatially resolved relaxation rate spectra.

Estimating the influence of finite instrumental resolution on elastic neutron scattering intensities from proteins.

Recent experimental and simulation studies show that the fractional Ornstein-Uhlenbeck process describes well the single particle motions in internal protein dynamics. Here the authors use this model to estimate the influence of finite instrumental resolution on elastic neutron scattering intensities from hydrated protein powders. They give, in particular, an estimation of the attenuation factor for the observed atomic position fluctuations, assuming a Gaussian and a triangular resolution function.

Diffusive dynamics of water in tert-butyl alcohol/water mixtures.

Quasielastic neutron scattering has been used to investigate the dynamical behavior of H(2)O in water/tert-butyl alcohol solutions. The measurements were made at fixed temperature (293 K) as a function of tert-butyl alcohol molar fraction, x, in the range 0-0.042. The data have been compared to those of pure water in the temperature range 269-293 K. The effect of tert-butyl alcohol addition on water dynamics is equivalent to that obtained by lowering the temperature of pure water by an amount proportional to the alcohol concentration. The temperature dependence of the diffusivity parameters in pure water and their concentration dependence in tert-butyl alcohol/water solutions can be rescaled to a common curve attributing to each solution a concentration-dependent "structural temperature" lower than the actual thermodynamic one. These results can be understood in terms of Stillinger's picture of water structuring and of other more recent theoretical pictures that emphasize the influence of the geometrical properties of hydrogen bond networks on water mobility.