A finely balanced order-disorder equilibrium sculpts the folding-binding landscape of an antibiotic sequestering protein.

TipA, a MerR family transcription factor from Streptomyces lividans, promotes antibiotic resistance by sequestering broad-spectrum thiopeptide-based antibiotics, thus counteracting their inhibitory effect on ribosomes. TipAS, a minimal binding motif which is expressed as an isoform of TipA, harbors a partially disordered N-terminal subdomain that folds upon binding multiple antibiotics. The extent and nature of the underlying molecular heterogeneity in TipAS that shapes its promiscuous folding-function landscape is an open question and is critical for understanding antibiotic-sequestration mechanisms. Here, combining equilibrium and time-resolved experiments, statistical modeling, and simulations, we show that the TipAS native ensemble exhibits a pre-equilibrium between binding-incompetent and binding-competent substates, with the fully folded state appearing only as an excited state under physiological conditions. The binding-competent state characterized by a partially structured N-terminal subdomain loses structure progressively in the physiological range of temperatures, swells on temperature increase, and displays slow conformational exchange across multiple conformations. Binding to the bactericidal antibiotic thiostrepton follows a combination of induced-fit and conformational-selection-like mechanisms, via partial binding and concomitant stabilization of the binding-competent substate. These ensemble features are evolutionarily conserved across orthologs from select bacteria that infect humans, underscoring the functional role of partial disorder in the native ensemble of antibiotic-sequestering proteins belonging to the MerR family.

Theoretical Insights into MutY Glycosylase DNA Repair Mechanism.

Maintaining the integrity of the genome is fundamental to living organisms. To this end, nature developed several mechanisms to find and promptly repair DNA lesions. Among them, base excision repair (BER) enzymes evolved to efficiently carry out this task. Notably, the mechanisms allowing these proteins to search for, detect, and fix DNA damage on a biologically relevant time scale still remain partially unclear. By taking MutY, a BER enzyme implied in the repair of the 8-oxoguanine-adenine mismatches, as a model system, we shed some light on the repair mechanism through a theoretical-computational approach. First, we estimated the effect of the oxidation state of the MutY iron-sulfur cluster on the protein-DNA binding. Then, the redox thermodynamics of both the protein cluster and DNA nucleobases are calculated. Finally, the charge migration kinetics along the double strand bound to the enzyme has been evaluated. The rationalization of our results indicates that the search for DNA lesions is essentially dictated by the redox chemistry of the species involved, i.e., the iron-sulfur redox cofactor and the DNA bound to the enzyme.

Unveiling the Cleavage Mechanism of an RNA Model Compound on the whole pH Scale: Computations Meet Experiments in the Determination of Reaction Rates.

The knowledge of the mechanism of reactions occurring in solution is a primary research line both in the context of theoretical-computational chemistry and in the field of organic and bio-organic chemistry. Given the importance of the hydrolysis of nucleic acids in life-related phenomena, here we present a combined experimental and computational study on the cleavage of an RNA model compound. This phosphodiester features a cleavage rate strictly dependent on the pH with three different dependence domains. Such experimental evidence, highlighted by an in-depth kinetic investigation, unequivocally suggests a change in the reaction mechanism along the pH scale. In order to interpret the data and to explain the experimental behavior, we have applied a theoretical-computational procedure, involving a hybrid quantum/classical approach, able to model chemical reactions in complex environments, i. e. in solution. This study turns out to quantitatively reproduce the experimental data with accuracy and, in addition, provides useful mechanistic insight into the transesterification process of the investigated compound. The study indicates that the cleavage can occur through an A N D N ${A_N D_N }$ , an A N + D N ${A_N + D_N }$ , and a D N A N ${D_N A_N }$ mechanism depending on the pH values.

A Theoretical-Computational Study of Phosphodiester Bond Cleavage Kinetics as a Function of the Temperature.

The hydrolysis of the phosphodiester bond is an important chemical reaction involved in several biological processes. Here, we study the cleavage of this bond by means of a theoretical-computational method in a model system, the dineopentyl phosphate. By such an approach, we reconstructed the kinetics and related thermodynamics of this chemical reaction along an isochore. In particular, we evaluated the kinetic constants of all the reaction steps within a wide range of temperatures, mostly corresponding to conditions where no experimental measures are available due to the extremely slow kinetics. Our results, in good agreement with the experimental data, show the robustness of our theoretical-computational methodology which can be easily extended to more complex systems.

Role of the membrane anchor in the regulation of Lck activity.

Theoretical work suggests that collective spatiotemporal behavior of integral membrane proteins should be modulated by boundary lipids sheathing their membrane anchors. Here, we show evidence for this prediction while investigating the mechanism for maintaining a steady amount of the active form of integral membrane protein Lck kinase (LckA) by Lck trans-autophosphorylation regulated by the phosphatase CD45. We used super-resolution microscopy, flow cytometry, and pharmacological and genetic perturbation to gain insight into the spatiotemporal context of this process. We found that LckA is generated exclusively at the plasma membrane, where CD45 maintains it in a ceaseless dynamic equilibrium with its unphosphorylated precursor. Steady LckA shows linear dependence, after an initial threshold, over a considerable range of Lck expression levels. This behavior fits a phenomenological model of trans-autophosphorylation that becomes more efficient with increasing LckA. We then challenged steady LckA formation by genetically swapping the Lck membrane anchor with structurally divergent ones, such as that of Src or the transmembrane domains of LAT, CD4, palmitoylation-defective CD4 and CD45 that were expected to drastically modify Lck boundary lipids. We observed small but significant changes in LckA generation, except for the CD45 transmembrane domain that drastically reduced LckA due to its excessive lateral proximity to CD45. Comprehensively, LckA formation and maintenance can be best explained by lipid bilayer critical density fluctuations rather than liquid-ordered phase-separated nanodomains, as previously thought, with "like/unlike" boundary lipids driving dynamical proximity and remoteness of Lck with itself and with CD45.

ERAP1 and ERAP2 Haplotypes Influence Suboptimal HLA-B*27:05-Restricted Anti-Viral CD8+ T Cell Responses Cross-Reactive to Self-Epitopes.

The human leukocyte antigen (HLA)-B*27 family of alleles is strongly associated with ankylosing spondylitis (AS), a chronic inflammatory disorder affecting the axial and peripheral joints, yet some HLA-B*27 variants not associated with AS have been shown. Since no major differences in the ligandome of associated compared to not-associated alleles have emerged, a plausible hypothesis is that the quantity rather than the quality of the presented epitopes makes the difference. In addition, the Endoplasmic Reticulum AminoPeptidases (ERAPs) 1 and 2, playing a crucial role in shaping the HLA class I epitopes, act as strong AS susceptibility factors, suggesting that an altered peptidome might be responsible for the activation of pathogenic CD8+ T cells. In this context, we have previously singled out a B*27:05-restricted CD8+ T cell response against pEBNA3A (RPPIFIRRL), an EBV peptide lacking the B*27 classic binding motif. Here, we show that a specific ERAP1/2 haplotype negatively correlates with such response in B*27:05 subjects. Moreover, we prove that the B*27:05 allele successfully presents peptides with the same suboptimal N-terminal RP motif, including the self-peptide, pDYNEIN (RPPIFGDFL). Overall, this study underscores the cooperation between the HLA-B*27 and ERAP1/2 allelic variants in defining CD8+ T cell reactivity to suboptimal viral and self-B*27 peptides and prompts further investigation of the B*27:05 peptidome composition.

Modelling Complex Bimolecular Reactions in a Condensed Phase: The Case of Phosphodiester Hydrolysis.

(1) Background: the theoretical modelling of reactions occurring in liquid phase is a research line of primary importance both in theoretical-computational chemistry and in the context of organic and biological chemistry. Here we present the modelling of the kinetics of the hydroxide-promoted hydrolysis of phosphoric diesters. (2) Method: the theoretical-computational procedure involves a hybrid quantum/classical approach based on the perturbed matrix method (PMM) in conjunction with molecular mechanics. (3) Results: the presented study reproduces the experimental data both in the rate constants and in the mechanistic aspects (C-O bond vs. O-P bond reactivity). The study suggests that the basic hydrolysis of phosphodiesters occurs through a concerted ANDN mechanism, with no formation of penta-coordinated species as reaction intermediates. (4) Conclusions: the presented approach, despite the approximations, is potentially applicable to a large number of bimolecular transformations in solution and therefore leads the way to a fast and general method to predict the rate constants and reactivities/selectivities in complex environments.

Effect of Salts on the Conformational Dynamics of the Cytochrome P450 OleP.

Cytochrome P450 OleP catalytic activity is strongly influenced by its structural dynamic conformational behavior. Here, we combine equilibrium-binding experiments with all-atom molecular dynamics simulations to clarify how different environments affect OleP conformational equilibrium between the open and the closed-catalytic competent-forms. Our data clearly show that at high-ionic strength conditions, the closed form is favored, and, very interestingly, different mechanisms, depending on the chemistry of the cations, can be used to rationalize such an effect.

Atomic-Level View of the Functional Transition in Vertebrate Hemoglobins: The Case of Antarctic Fish Hbs.

Tetrameric hemoglobins (Hbs) are prototypal systems for studies aimed at unveiling basic structure-function relationships as well as investigating the molecular/structural basis of adaptation of living organisms to extreme conditions. However, a chronological analysis of decade-long studies conducted on Hbs is illuminating on the difficulties associated with the attempts of gaining functional insights from static structures. Here, we applied molecular dynamics (MD) simulations to explore the functional transition from the T to the R state of the hemoglobin of the Antarctic fish Trematomus bernacchii (HbTb). Our study clearly demonstrates the ability of the MD technique to accurately describe the transition of HbTb from the T to R-like states, as shown by a number of global and local structural indicators. A comparative analysis of the structural states that HbTb assumes in the simulations with those detected in previous MD analyses conducted on HbA (human Hb) highlights interesting analogies (similarity of the transition pathway) and differences (distinct population of intermediate states). In particular, the ability of HbTb to significantly populate intermediate states along the functional pathway explains the observed propensity of this protein to assume these structures in the crystalline state. It also explains some functional data reported on the protein that indicate the occurrence of other functional states in addition to the canonical R and T ones. These findings are in line with the emerging idea that the classical two-state view underlying tetrameric Hb functionality is probably an oversimplification and that other structural states play important roles in these proteins. The ability of MD simulations to accurately describe the functional pathway in tetrameric Hbs suggests that this approach may be effectively applied to unravel the molecular and structural basis of Hbs exhibiting peculiar functional properties as a consequence of the environmental adaptation of the host organism.

Absorption behavior of doxorubicin hydrochloride in visible region in different environments: a combined experimental and computational study.

The experimental absorption measurements in the interval 350-600 nm (Vis), molecular dynamics simulations, quantum-mechanics calculations and an advanced molecular treatment of simulation data are here combined to provide a complete picture of the absorption behavior in the visible portion of the electromagnetic spectrum of the doxorubicin hydrochloride (DX) molecule in different environments. By such an approach, we have shown that it is possible to characterize the effect of the environment on the DX absorption behavior - including the vibronic contributions - as well as to interpret such differences in terms of molecular electronic excited states, which are found to be strongly influenced by the environment.

Complexation and organization of doxorubicin on polystyrene sulfonate chains: impacts on doxorubicin dimerization and quenching.

Anthracycline doxorubicin hydrochloride (DX) is a positively charged fluorescent drug, which in water self-associates into non-fluorescent antiparallel dimers upon increasing concentration and/or ionic strength. The positive charge of DX allows for complexation with negatively charged polymers and drug carriers. The fluorescence of DX following complexation with polyanion polystyrene sulfonate (PSS) is studied here. The fluorescence emission of DX decreases in the presence of PSS, being almost completely quenched when the ratio (R) of PSS monomers-to-DX molecules is larger than 10. Increasing R values over 30 results in a progressive recovery of fluorescence. The circular dichroism of PSS-DX complexes shows inverted characteristic bands of DX dimers suggesting the presence of parallel dimers at a concentration of DX below dimerization in water. Molecular dynamics studies corroborate a preferential orientation of DX into parallel dimers when interacting with PSS and show that DX molecules interact with a binding pocket of PSS monomers rather than with one single monomer. Increasing the ionic strength results in a recovery of fluorescence without an apparent release of DX from the PSS-DX complex as shown by DOSY NMR. PSS acts as a template for concentrating DX, triggering dimerisation and orienting DX molecules with their charged groups facing the negatively charged PSS monomers.

Theoretical Evaluation of Sulfur-Based Reactions as a Model for Biological Antioxidant Defense.

Sulfur-containing amino acids, Methionine (Met) and Cysteine (Cys), are very susceptible to Reactive Oxygen Species (ROS). Therefore, sulfur-based reactions regulate many biological processes, playing a key role in maintaining cellular redox homeostasis and modulating intracellular signaling cascades. In oxidative conditions, Met acts as a ROS scavenger, through Met sulfoxide formation, while thiol/disulfide interchange reactions take place between Cys residues as a response to many environmental stimuli. In this work, we apply a QM/MM theoretical-computational approach, which combines quantum-mechanical calculations with classical molecular dynamics simulations to estimate the free energy profile for the above-mentioned reactions in solution. The results obtained, in good agreement with experimental data, show the validity of our approach in modeling sulfur-based reactions, enabling us to study these mechanisms in more complex biological systems.

Modeling Charge Transfer Reactions by Hopping between Electronic Ground State Minima: Application to Hole Transfer between DNA Bases.

In this paper, we extend the previously described general model for charge transfer reactions, introducing specific changes to treat the hopping between energy minima of the electronic ground state (i.e., transitions between the corresponding vibrational ground states). We applied the theoretical-computational model to the charge transfer reactions in DNA molecules which still represent a challenge for a rational full understanding of their mechanism. Results show that the presented model can provide a valid, relatively simple, approach to quantitatively study such reactions shedding light on several important aspects of the reaction mechanism.

Theoretical Modeling of Redox Potentials of Biomolecules.

The estimation of the redox potentials of biologically relevant systems by means of theoretical-computational approaches still represents a challenge. In fact, the size of these systems typically does not allow a full quantum-mechanical treatment needed to describe electron loss/gain in such a complex environment, where the redox process takes place. Therefore, a number of different theoretical strategies have been developed so far to make the calculation of the redox free energy feasible with current computational resources. In this review, we provide a survey of such theoretical-computational approaches used in this context, highlighting their physical principles and discussing their advantages and limitations. Several examples of these approaches applied to the estimation of the redox potentials of both proteins and nucleic acids are described and critically discussed. Finally, general considerations on the most promising strategies are reported.

A Simplified Treatment for Efficiently Modeling the Spectral Signal of Vibronic Transitions: Application to Aqueous Indole.

In this paper, we introduce specific approximations to simplify the vibronic treatment in modeling absorption and emission spectra, allowing us to include a huge number of vibronic transitions in the calculations. Implementation of such a simplified vibronic treatment within our general approach for modelling vibronic spectra, based on molecular dynamics simulations and the perturbed matrix method, provided a quantitative reproduction of the absorption and emission spectra of aqueous indole with higher accuracy than the one obtained when using the existing vibronic treatment. Such results, showing the reliability of the approximations employed, indicate that the proposed method can be a very efficient and accurate tool for computational spectroscopy.

Effects of Environmental and Electric Perturbations on the pKa of Thioredoxin Cysteine 35: A Computational Study.

Here we present a theoretical-computational study dealing with the evaluation of the pKa of the Cysteine residues in Thioredoxin (TRX) and in its complex with the Thioredoxin-interacting protein (TXNIP). The free energy differences between the anionic and neutral form of the Cysteine 32 and 35 have been evaluated by means of the Perturbed Matrix Method with classical perturbations due to both the environment and an exogenous electric field as provided by Molecular Dynamics (MD) simulations. The evaluation of the free energies allowed us to show that the effect of the perturbing terms is to lower the pKa of Cysteine 32 and Cysteine 35 with respect to the free amino-acid. On the other hand, in the complex TRX-TXNIP, our data show an enhanced stabilization of the neutral reduced form of Cys 35. These results suggest that external electric stimuli higher than 0.02 V/nm can modulate the Cysteine pKa, which can be connected to the tight regulation of the TRX acting as an antioxidant agent.

Quaternary Structure Transitions of Human Hemoglobin: An Atomic-Level View of the Functional Intermediate States.

Human hemoglobin (HbA) is one of the prototypal systems used to investigate structure-function relationships in proteins. Indeed, HbA has been used to develop the basic concepts of protein allostery, although the atomic-level mechanism underlying the HbA functionality is still highly debated. This is due to the fact that most of the three-dimensional structural information collected over the decades refers to the endpoints of HbA functional transition with little data available for the intermediate states. Here, we report molecular dynamics (MD) simulations by focusing on the relevance of the intermediate states of the protein functional transition unraveled by the crystallographic studies carried out on vertebrate Hbs. Fully atomistic simulations of the HbA T-state indicate that the protein undergoes a spontaneous transition toward the R-state. The inspection of the trajectory structures indicates that the protein significantly populates the intermediate HL-(C) state previously unraveled by crystallography. In the structural transition, it also assumes the intermediate states crystallographically detected in Antarctic fish Hbs. This finding suggests that HbA and Antarctic fish Hbs, in addition to the endpoints of the transitions, also share a similar deoxygenation pathway despite a distace of hundreds of millions of years in the evolution scale. Finally, using the essential dynamic sampling methodology, we gained some insights into the reverse R to T transition that is not spontaneously observed in classic MD simulations.

Exploiting Reaction-Diffusion Conditions to Trigger Pathway Complexity in the Growth of a MOF.

Coordination polymers (CPs), including metal-organic frameworks (MOFs), are crystalline materials with promising applications in electronics, magnetism, catalysis, and gas storage/separation. However, the mechanisms and pathways underlying their formation remain largely undisclosed. Herein, we demonstrate that diffusion-controlled mixing of reagents at the very early stages of the crystallization process (i.e., within ≈40 ms), achieved by using continuous-flow microfluidic devices, can be used to enable novel crystallization pathways of a prototypical spin-crossover MOF towards its thermodynamic product. In particular, two distinct and unprecedented nucleation-growth pathways were experimentally observed when crystallization was triggered under microfluidic mixing. Full-atom molecular dynamics simulations also confirm the occurrence of these two distinct pathways during crystal growth. In sharp contrast, a crystallization by particle attachment was observed under bulk (turbulent) mixing. These unprecedented results provide a sound basis for understanding the growth of CPs and open up new avenues for the engineering of porous materials by using out-of-equilibrium conditions.

Quantitative Characterization of Binding Pockets and Binding Complementarity by Means of Zernike Descriptors.

In this work, we describe the application of the Zernike formalism to quantitatively characterize the binding pockets of two sets of biologically relevant systems. Such an approach, when applied to molecular dynamics trajectories, is able to pinpoint the subtle differences between very similar molecular regions and their impact on the local propensity to ligand binding, allowing us to quantify such differences. The statistical robustness of our procedure suggests that it is very suitable to describe protein binding sites and protein-ligand interactions within a rigorous and well-defined framework.

Characterization of the catalytic flexible loop in the dihydroorotase domain of the human multi-enzymatic protein CAD.

The dihydroorotase (DHOase) domain of the multifunctional protein carbamoyl-phosphate synthetase 2, aspartate transcarbamoylase, and dihydroorotase (CAD) catalyzes the third step in the de novo biosynthesis of pyrimidine nucleotides in animals. The crystal structure of the DHOase domain of human CAD (huDHOase) revealed that, despite evolutionary divergence, its active site components are highly conserved with those in bacterial DHOases, encoded as monofunctional enzymes. An important element for catalysis, conserved from Escherichia coli to humans, is a flexible loop that closes as a lid over the active site. Here, we combined mutagenic, structural, biochemical, and molecular dynamics analyses to characterize the function of the flexible loop in the activity of CAD's DHOase domain. A huDHOase chimera bearing the E. coli DHOase flexible loop was inactive, suggesting the presence of distinctive elements in the flexible loop of huDHOase that cannot be replaced by the bacterial sequence. We pinpointed Phe-1563, a residue absolutely conserved at the tip of the flexible loop in CAD's DHOase domain, as a critical element for the conformational equilibrium between the two catalytic states of the protein. Substitutions of Phe-1563 with Ala, Leu, or Thr prevented the closure of the flexible loop and inactivated the protein, whereas substitution with Tyr enhanced the interactions of the loop in the closed position and reduced fluctuations and the reaction rate. Our results confirm the importance of the flexible loop in CAD's DHOase domain and explain the key role of Phe-1563 in configuring the active site and in promoting substrate strain and catalysis.