Noniterative Doubles Corrections to the Random Phase and Higher Random Phase Approximations:Singlet and Triplet Excitation Energies

The second-order noniterative doubles-corrected random phase approximation (RPA) method has been extended to triplet excitation energies and the doubles-corrected higher RPA method as well as a shifted version for calculating singlet and triplet excitation energies are presented here for the first time. A benchmark set consisting of 20 molecules with a total of 117 singlet and 71 triplet excited states has been used to test the performance of the new methods by comparison with previous results obtained with the second-order polarization propagator approximation (SOPPA) and the third order approximate coupled cluster singles, doubles and triples model CC3. In general, the second-order doubles corrections to RPA and HRPA significantly reduce both the mean deviation as well as the standard deviation of the errors compared to the CC3 results. The accuracy of the new methods approaches the accuracy of the SOPPA method while using only 10–60% of the calculation time.</p

Simple models of protein folding and of non--conventional drug design

While all the information required for the folding of a protein is contained in its amino acid sequence, one has not yet learned how to extract this information to predict the three--dimensional, biologically active, native conformation of a protein whose sequence is known. Using insight obtained from simple model simulations of the folding of proteins, in particular of the fact that this phenomenon is essentially controlled by conserved (native) contacts among (few) strongly interacting ("hot"), as a rule hydrophobic, amino acids, which also stabilize local elementary structures (LES, hidden, incipient secondary structures like $\alpha$--helices and $\beta$--sheets) formed early in the folding process and leading to the postcritical folding nucleus (i.e., the minimum set of native contacts which bring the system pass beyond the highest free--energy barrier found in the whole folding process) it is possible to work out a succesful strategy for reading the native structure of designed proteins from the knowledge of only their amino acid sequence and of the contact energies among the amino acids. Because LES have undergone millions of years of evolution to selectively dock to their complementary structures, small peptides made out of the same amino acids as the LES are expected to selectively attach to the newly expressed (unfolded) protein and inhibit its folding, or to the native (fluctuating) native conformation and denaturate it. These peptides, or their mimetic molecules, can thus be used as effective non--conventional drugs to those already existing (and directed at neutralizing the active site of enzymes), displaying the advantage of not suffering from the uprise of resistance

Theoretical estimates of the anapole magnetizabilities of C4H4X2 cyclic molecules for X=O, S, Se, and Te

Calculations have been carried out for C4H4X2 cyclic molecules, with X=O, S, Se, and Te, characterized by the presence of magnetic-field induced toroidal electron currents and associated orbital anapole moments. The orbital anapole induced by a static nonuniform magnetic field B, with uniform curl C =∇× B, is rationalized via a second-rank anapole magnetizability tensor aαβ , defined as minus the second derivative of the second-order interaction energy with respect to the components Cα and Bβ. The average anapole magnetizability a equals −χ, the pseudoscalar obtained by spatial averaging of the dipole-quadrupole magnetizability χα,βγ . It has different sign for D and L enantiomeric systems and can therefore be used for chiral discrimination. Therefore, in an isotropic chiral medium, a homogeneous magnetic field induces an electronic anapole Aα, having the same magnitude, but opposite sign, for two enantiomorphs.Fil: Pagola, Gabriel Ignacio. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; ArgentinaFil: Ferraro, Marta Beatriz. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; ArgentinaFil: Provasi, Patricio Federico. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Nordeste. Instituto de Modelado e Innovación Tecnológica. Universidad Nacional del Nordeste. Facultad de Ciencias Exactas Naturales y Agrimensura. Instituto de Modelado e Innovación Tecnologica; ArgentinaFil: Pelloni, Stefano. Universidad de Modena y Reggio Emilia. Departamento de Química; ItaliaFil: Lazzeretti, Paolo. Universidad de Modena y Reggio Emilia. Departamento de Química; Italia; Itali

Ligand-Induced Modulation of the Free-Energy Landscape of G Protein-Coupled Receptors Explored by Adaptive Biasing Techniques

Extensive experimental information supports the formation of ligand-specific conformations of G protein-coupled receptors (GPCRs) as a possible molecular basis for their functional selectivity for signaling pathways. Taking advantage of the recently published inactive and active crystal structures of GPCRs, we have implemented an all-atom computational strategy that combines different adaptive biasing techniques to identify ligand-specific conformations along pre-determined activation pathways. Using the prototypic GPCR β2-adrenergic receptor as a suitable test case for validation, we show that ligands with different efficacies (either inverse agonists, neutral antagonists, or agonists) modulate the free-energy landscape of the receptor by shifting the conformational equilibrium towards active or inactive conformations depending on their elicited physiological response. Notably, we provide for the first time a quantitative description of the thermodynamics of the receptor in an explicit atomistic environment, which accounts for the receptor basal activity and the stabilization of different active-like states by differently potent agonists. Structural inspection of these metastable states reveals unique conformations of the receptor that may have been difficult to retrieve experimentally

Biapenem Inactivation by B2 Metallo β-Lactamases: Energy Landscape of the Post-Hydrolysis Reactions

<div><h3>Background</h3><p>The first line of defense by bacteria against <em>β</em>-lactam antibiotics is the expression of β-lactamases, which cleave the amide bond of the β-lactam ring. In the reaction of biapenem inactivation by B2 metallo β-lactamases (MβLs), after the β-lactam ring is opened, the carboxyl group generated by the hydrolytic process and the hydroxyethyl group (common to all carbapenems) rotate around the C5–C6 bond, assuming a new position that allows a proton transfer from the hydroxyethyl group to C2, and a nucleophilic attack on C3 by the oxygen atom of the same side-chain. This process leads to the formation of a bicyclic compound, as originally observed in the X-ray structure of the metallo β-lactamase CphA in complex with product.</p> <h3>Methodology/Principal Findings</h3><p>QM/MM and metadynamics simulations of the post-hydrolysis steps in solution and in the enzyme reveal that while the rotation of the hydroxyethyl group can occur in solution or in the enzyme active site, formation of the bicyclic compound occurs primarily in solution, after which the final product binds back to the enzyme. The calculations also suggest that the rotation and cyclization steps can occur at a rate comparable to that observed experimentally for the enzymatic inactivation of biapenem only if the hydrolysis reaction leaves the N4 nitrogen of the β-lactam ring unprotonated.</p> <h3>Conclusions/Significance</h3><p>The calculations support the existence of a common mechanism (in which ionized N4 is the leaving group) for carbapenems hydrolysis in all MβLs, and suggest a possible revision of mechanisms for B2 MβLs in which the cleavage of the β-lactam ring is associated with or immediately followed by protonation of N4. The study also indicates that the bicyclic derivative of biapenem has significant affinity for B2 MβLs, and that it may be possible to obtain clinically effective inhibitors of these enzymes by modification of this lead compound.</p> </div

The Integrin Receptor in Biologically Relevant Bilayers: Insights from Molecular Dynamics Simulations

Integrins are heterodimeric (αβ) cell surface receptors that are potential therapeutic targets for a number of diseases. Despite the existence of structural data for all parts of integrins, the structure of the complete integrin receptor is still not available. We have used available structural data to construct a model of the complete integrin receptor in complex with talin F2–F3 domain. It has been shown that the interactions of integrins with their lipid environment are crucial for their function but details of the integrin/lipid interactions remain elusive. In this study an integrin/talin complex was inserted in biologically relevant bilayers that resemble the cell plasma membrane containing zwitterionic and charged phospholipids, cholesterol and sphingolipids to study the dynamics of the integrin receptor and its effect on bilayer structure and dynamics. The results of this study demonstrate the dynamic nature of the integrin receptor and suggest that the presence of the integrin receptor alters the lipid organization between the two leaflets of the bilayer. In particular, our results suggest elevated density of cholesterol and of phosphatidylserine lipids around the integrin/talin complex and a slowing down of lipids in an annulus of ~30 Å around the protein due to interactions between the lipids and the integrin/talin F2–F3 complex. This may in part regulate the interactions of integrins with other related proteins or integrin clustering thus facilitating signal transduction across cell membranes

Ligand-Dependent Conformations and Dynamics of the Serotonin 5-HT2A Receptor Determine Its Activation and Membrane-Driven Oligomerization Properties

From computational simulations of a serotonin 2A receptor (5-HT2AR) model complexed with pharmacologically and structurally diverse ligands we identify different conformational states and dynamics adopted by the receptor bound to the full agonist 5-HT, the partial agonist LSD, and the inverse agonist Ketanserin. The results from the unbiased all-atom molecular dynamics (MD) simulations show that the three ligands affect differently the known GPCR activation elements including the toggle switch at W6.48, the changes in the ionic lock between E6.30 and R3.50 of the DRY motif in TM3, and the dynamics of the NPxxY motif in TM7. The computational results uncover a sequence of steps connecting these experimentally-identified elements of GPCR activation. The differences among the properties of the receptor molecule interacting with the ligands correlate with their distinct pharmacological properties. Combining these results with quantitative analysis of membrane deformation obtained with our new method (Mondal et al, Biophysical Journal 2011), we show that distinct conformational rearrangements produced by the three ligands also elicit different responses in the surrounding membrane. The differential reorganization of the receptor environment is reflected in (i)-the involvement of cholesterol in the activation of the 5-HT2AR, and (ii)-different extents and patterns of membrane deformations. These findings are discussed in the context of their likely functional consequences and a predicted mechanism of ligand-specific GPCR oligomerization

Is autoimmunity the Achilles' heel of cancer immunotherapy?

The emergence of immuno-oncology as the first broadly successful strategy for metastatic cancer will require clinicians to integrate this new pillar of medicine with chemotherapy, radiation, and targeted small-molecule compounds. Of equal importance is gaining an understanding of the limitations and toxicities of immunotherapy. Immunotherapy was initially perceived to be a relatively less toxic approach to cancer treatment than other available therapies-and surely it is, when compared to those. However, as the use of immunotherapy becomes more common, especially as first- and second-line treatments, immunotoxicity and autoimmunity are emerging as the Achilles' heel of immunotherapy. In this Perspective, we discuss evidence that the occurrence of immunotoxicity bodes well for the patient, and describe mechanisms that might be related to the induction of autoimmunity. We then explore approaches to limit immunotoxicity, and discuss the future directions of research and reporting that are needed to diminish it

Analysis of Isotope Effects in NMR One-bond Indirect Nuclear Spin-spin Coupling Constants in Terms of Localized Molecular Orbitals

We recently showed, by analyzing contributions from localized molecular orbitals, that the anomalous deuterium isotope effect in the one-bond indirect nuclear spin–spin coupling constant of methane, also called the unexpected differential sensitivity, can be explained by the transfer of s-orbital character from the stretched bond to the other unchanged bonds [ChemPhysChem, 2008, 9, 1259]. We now extend this analysis of isotope effects to the molecules BH4−, NH4+, SiH4, H2O and NH3 in order to test our conclusions on a wider rage of XH4 compounds and to investigate whether the lone-pair orbitals are really responsible for the absence of a similar effect in water and ammonia as proposed earlier [J. Chem. Phys., 2000, 113, 3121].Fil: Provasi, Patricio Federico. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Nordeste. Instituto de Modelado e Innovación Tecnológica. Universidad Nacional del Nordeste. Facultad de Ciencias Exactas Naturales y Agrimensura. Instituto de Modelado e Innovación Tecnologica; ArgentinaFil: Sauer, Stephan P. A.. Universidad de Copenhague. Departamento de Química; Dinamarc