Laboratory of Theoretical Biochemistry

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

Abstract Since the dawn of the computer age, scientists have designed devices to represent molecular structures and developed tools to simulate their dynamic behavior in silico. To this day, these tools remain central to our understanding of biomolecular phenomena. In contrast to other fields such as fluid mechanics or meteorology, the observation of molecular motions at the atomic level remains a major experimental challenge. Continuous advances in computer graphics and numerical computation, combined with the emergence of human–computer interaction approaches, led to the methodology of so‐called “interactive molecular simulations”, characterized by two main features. First, the possibility to visualize a running simulation in interactive time, that is, compatible with human perception. Second, the possibility to manipulate the simulation interactively by imposing a force, changing a biophysical property, or editing runtime parameters on the fly. Such simulations are still little used in computational biology, where it is more common to run a series of offline simulations and then visualize and analyze the results. However, interactive molecular simulation tools promise to handle time‐consuming tasks such as the modeling of particularly complex biomolecular structures more efficiently or to support approaches such as Rational Drug Design with regard to pharmaceutical applications.This article is categorized under: Computer and Information Science > Computer Algorithms and Programming Computer and Information Science > Visualization Software > Simulation Methods

<div> <div> <div> <p>In spite of the biological importance of the binding of Zn2+, Ca2+, and Mg2+ to the carboxylate group, cation-acetate binding affinities and binding modes remain actively debated. Here, we report the first use of Raman multivariate curve resolution (Raman-MCR) vibrational spectroscopy to obtain self-consistent free and bound metal acetate spectra and one-to-one binding constants, without the need to invoke any a priori assumptions regarding the shapes of the corresponding vibrational bands. The experimental results, combined with classical molecular dynamics simulations with a force field effectively accounting for electronic polarization via charge scaling and ab initio simulations, indicate that the measured binding constants pertain to direct (as opposed to water separated) ion pairing. The resulting binding constants do not scale with cation size, as </p><div> <div> <div> <p>the binding constant to Zn2+ is significantly larger than that to either Mg2+ or Ca2+, although Zn2+ and Mg2+ have similar radii that are about 25% smaller than Ca2+. Remaining uncertainties in the metal acetate binding free energies are linked to fundamental ambiguities associated with identifying the range of structures pertaining to non-covalently bound species. </p> </div> </div> </div> </div> </div> </div>

Electronic circular dichroism is one of the most used spectroscopic techniques for peptide and protein structural characterization. However, while valuable experimental spectra exist for α-helix, β-sheet and random coil secondary structures, previous studies showed important discrepancies for β-turns, limiting their use as a reference for structural studies. In this paper, we simulated circular dichroism spectra for the best-characterized β-turns in peptides, namely types I, II, I' and II'. In particular, by combining classical molecular dynamics simulations and state-of-the-art quantum time-dependent density functional theory (with the polarizable embedding multiscale model) computations, two common electronic circular dichroism patterns were found for couples of β-turn types (namely, type I/type II' and type II/type I'), at first for a minimal di-peptide model (Ace-Ala-Ala-NHMe), but also for all sequences tested with non-aromatic residues in the central positions. On the other hand, as expected, aromatic substitution causes important perturbations to the previously found ECD patterns. Finally, by applying suitable approximations, these patterns were subsequently rationalized based on the exciton chirality rule. All these results provide useful predictions and pave the way for a possible experimental characterization of β-turns based on circular dichroism spectroscopy.

The double-stranded RNA-binding domain (dsRBD) is a broadly distributed domain among RNA-maturing enzymes. Although this domain recognizes dsRNA's structures via a conserved canonical structure adopting an α1-β1β2β3-α2 topology, several dsRBDs can accommodate discrete structural extensions expanding further their functional repertoire. How these structural elements engage cooperative communications with the canonical structure and how they contribute to the dsRBD's overall folding are poorly understood. Here, we addressed these issues using the dsRBD of human dihydrouridine synthase-2 (hDus2) (hDus2-dsRBD) as a model. This dsRBD harbors N- and C-terminal extensions, the former being directly involved in the recognition of tRNA substrate of hDus2. These extensions engage residues that form a long-range hydrophobic network (LHN) outside the RNA-binding interface. We show by coarse-grain Brownian dynamics that the Nt-extension and its residues F359 and Y364 rigidify the major folding nucleus of the canonical structure via an indirect effect. hDus2-dsRBD unfolds following a two-state cooperative model, whereas both F359A and Y364A mutants, designed to destabilize this LHN, unfold irreversibly. Structural and computational analyses show that these mutants are unstable due to an increase in the dynamics of the two extensions favoring solvent exposure of α2-helix and weakening the main folding nucleus rigidity. This LHN appears essential for maintaining a thermodynamic stability of the overall system and eventually a functional conformation for tRNA recognition. Altogether, our findings suggest that functional adaptability of extended dsRBDs is promoted by a cooperative hydrophobic coupling between the extensions acting as effectors and the folding nucleus of the canonical structure.

Although the magnetic and transport properties of molecular junction systems composed of metalled porphyrins or phthalocyanines have been broadly studied in recent years, to date no studies have been devoted to evaluate the aforementioned properties in junction systems featuring metalled expanded porphyrins as active elements. The present work reports a detailed theoretical study of the magnetic and electronic transport properties of the recently synthesized dinuclear Cu(ii)-naphthoisoamethyrin complex (PyCu2). This is the first work on performing these kinds of studies using a magnetically coupled metallic expanded porphyrin as a molecular kernel. DFT and wave function-based methods have been used to determine the nature of the magnetic interaction between the metallic centres, characterized by the exchange coupling constant J, showing that although this was found to be weakly antiferromagnetic, after an exhaustive analysis it turns out that the coupling has a ferromagnetic nature with a value of J = 14.2 cm-1. Once the magnetic ground state of PyCu2 was rigorously established, the spin resolved transport properties of the device composed of the expanded porphyrin attached to two gold nano-wires were studied by means of the combination of DFT and the nonequilibrium Green's function formalism, in order to explore PyCu2 prospects as a possible spintronic device.

The nanomechanical properties of elastomeric proteins determine the elasticity of a variety of tissues. A widespread natural tactic to regulate protein extensibility lies in the presence of covalent disulfide bonds, which significantly enhance protein stiffness. The prevalent in vivo strategy to form disulfide bonds requires the presence of dedicated enzymes. Here we propose an alternative chemical route to promote non-enzymatic oxidative protein folding via disulfide isomerization based on naturally occurring small molecules. Using single-molecule force-clamp spectroscopy, supported by DFT calculations and mass spectrometry measurements, we demonstrate that subtle changes in the chemical structure of a transient mixed-disulfide intermediate adduct between a protein cysteine and an attacking low molecular-weight thiol have a dramatic effect on the protein’s mechanical stability. This approach provides a general tool to rationalize the dynamics of S-thiolation and its role in modulating protein nanomechanics, offering molecular insights on how chemical reactivity regulates protein elasticity.

The two hallmarks of Alzheimer's disease (AD) are the presence of neurofibrillary tangles (NFT) made of aggregates of the hyperphosphorylated tau protein and of amyloid plaques composed of amyloid-β (Aβ) peptides, primarily Aβ1-40 and Aβ1-42. Targeting the production, aggregation, and toxicity of Aβ with small molecule drugs or antibodies is an active area of AD research due to the general acceptance of the amyloid cascade hypothesis, but thus far all drugs targeting Aβ have failed. From a review of the recent literature and our own experience based on in vitro, in silico, and in vivo studies, we present some reasons to explain this repetitive failure.

A Ni-doped graphene nanosheet is used as an efficient catalyst for the reduction of NO by CO molecules.