Solution Structure and Conformational Flexibility of a Polyketide Synthase Module.

Polyketide synthases (PKSs) are versatile C-C bond-forming enzymes that are broadly distributed in bacteria and fungi. The polyketide compound family includes many clinically useful drugs such as the antibiotic erythromycin, the antineoplastic epothilone, and the cholesterol-lowering lovastatin. Harnessing PKSs for custom compound synthesis remains an open challenge, largely because of the lack of knowledge about key structural properties. Particularly, the domains-well characterized on their own-are poorly understood in their arrangement, conformational dynamics, and interplay in the intricate quaternary structure of modular PKSs. Here, we characterize module 2 from the 6-deoxyerythronolide B synthase by small-angle X-ray scattering and cross-linking mass spectrometry with coarse-grained structural modeling. The results of this hybrid approach shed light on the solution structure of a cis-AT type PKS module as well as its inherent conformational dynamics. Supported by a directed evolution approach, we also find that acyl carrier protein (ACP)-mediated substrate shuttling appears to be steered by a nonspecific electrostatic interaction network.

Analysis and engineering of substrate shuttling by the acyl carrier protein (ACP) in fatty acid synthases (FASs).

Perturbations of domain-domain interactions impact the function of type I fatty acid synthases. We identify interface point mutations that modulate fatty acid chain lengths, and explain their effect in changes of domain-domain binding energetics. Engineering of similar interfaces in related megasynthases may be exploited for custom synthesis of natural products.

Empirical Conversion of pK a Values between Different Solvents and Interpretation of the Parameters: Application to Water, Acetonitrile, Dimethyl Sulfoxide, and Methanol.

An empirical conversion method (ECM) that transforms pK a values of arbitrary organic compounds from one solvent to the other is introduced. We demonstrate the method's usefulness and performance on pK a conversions involving water and organic solvents acetonitrile (MeCN), dimethyl sulfoxide (Me2SO), and methanol (MeOH). We focus on the pK a conversion from the known reference value in water to the other three organic solvents, although such a conversion can also be performed between any pair of the considered solvents. The ECM works with an additive parameter that is specific to a solvent and a molecular family (essentially characterized by a functional group that is titrated). We formally show that the method can be formulated with a single additive parameter, and that the extra multiplicative parameter used in other works is not required. The values of the additive parameter are determined from known pK a data, and their interpretation is provided on the basis of physicochemical concepts. The data set of known pK a values is augmented with pK a values computed with the recently introduced electrostatic transform method, whose validity is demonstrated. For a validation of our method, we consider pK a conversions for two data sets of titratable compounds. The first data set involves 81 relatively small molecules belonging to 19 different molecular families, with the pK a data available in all four considered solvents. The second data set involves 76 titratable molecules from 5 additional molecular families. These molecules are typically larger, and their experimental pK a values are available only in Me2SO and water. The validation tests show that the agreement between the experimental pK a data and the ECM predictions is generally good, with absolute errors often on the order of 0.5 pH units. The presence of a few outliers is rationalized, and observed trends with respect to molecular families are discussed.

Merging Structural Information from X-ray Crystallography, Quantum Chemistry, and EXAFS Spectra: The Oxygen-Evolving Complex in PSII.

Structural data of the oxygen-evolving complex (OEC) in photosystem II (PSII) determined by X-ray crystallography, quantum chemistry (QC), and extended X-ray absorption fine structure (EXAFS) analyses are presently inconsistent. Therefore, a detailed study of what information can be gained about the OEC through a comparison of QC and crystallographic structure information combined with the information from range-extended EXAFS spectra was undertaken. An analysis for determining the precision of the atomic coordinates of the OEC by QC is carried out. OEC model structures based on crystallographic data that are obtained by QC from different research groups are compared with one another and with structures obtained by high-resolution crystallography. The theory of EXAFS spectra is summarized, and the application of EXAFS spectra to the experimental determination of the structure of the OEC is detailed. We discriminate three types of parameters entering the formula for the EXAFS spectrum: (1) model-independent, predefined, and fixed; (2) model-dependent that can be computed or adjusted; and (3) model-dependent that must be adjusted. The information content of EXAFS spectra is estimated and is related to the precision of atomic coordinates and resolution power to discriminate different atom-pair distances of the OEC. It is demonstrated how a precise adjustment of atomic coordinates can yield a nearly perfect representation of the experimental OEC EXAFS spectrum, but at the expense of overfitting and losing the knowledge of the initial OEC model structure. Introducing a novel type of penalty function, it is shown that moderate adjustment of atomic coordinates to the EXAFS spectrum limited by constraints avoids overfitting and can be used to validate different OEC model structures. This technique is used to identify the OEC model structures whose computed OEC EXAFS spectra agree best with the measured spectrum. In this way, the most likely S-state and protonation pattern of the OEC for the most recent high-resolution crystal structure of PSII are determined. We find that the X-ray free-electron laser (XFEL) structure is indeed not significantly affected by exposure to XFEL pulses and thus results in a radiation-damage-free model of the OEC.

Computing pKa Values in Different Solvents by Electrostatic Transformation.

We introduce a method that requires only moderate computational effort to compute pKa values of small molecules in different solvents with an average accuracy of better than 0.7 pH units. With a known pKa value in one solvent, the electrostatic transform method computes the pKa value in any other solvent if the proton solvation energy is known in both considered solvents. To apply the electrostatic transform method to a molecule, the electrostatic solvation energies of the protonated and deprotonated molecular species are computed in the two considered solvents using a dielectric continuum to describe the solvent. This is demonstrated for 30 molecules belonging to 10 different molecular families by considering 77 measured pKa values in 4 different solvents: water, acetonitrile, dimethyl sulfoxide, and methanol. The electrostatic transform method can be applied to any other solvent if the proton solvation energy is known. It is exclusively based on physicochemical principles, not using any empirical fetch factors or explicit solvent molecules, to obtain agreement with measured pKa values and is therefore ready to be generalized to other solute molecules and solvents. From the computed pKa values, we obtained relative proton solvation energies, which agree very well with the proton solvation energies computed recently by ab initio methods, and used these energies in the present study.

Proton solvation in protic and aprotic solvents.

Protonation pattern strongly affects the properties of molecular systems. To determine protonation equilibria, proton solvation free energy, which is a central quantity in solution chemistry, needs to be known. In this study, proton affinities (PAs), electrostatic energies of solvation, and pKA values were computed in protic and aprotic solvents. The proton solvation energy in acetonitrile (MeCN), methanol (MeOH), water, and dimethyl sulfoxide (DMSO) was determined from computed and measured pKA values for a specially selected set of organic compounds. pKA values were computed with high accuracy using a combination of quantum chemical and electrostatic approaches. Quantum chemical density functional theory computations were performed evaluating PA in the gas-phase. The electrostatic contributions of solvation were computed solving the Poisson equation. The computations yield proton solvation free energies with high accuracy, which are in MeCN, MeOH, water, and DMSO -255.1, -265.9, -266.3, and -266.4 kcal/mol, respectively, where the value for water is close to the consensus value of -265.9 kcal/mol. The pKA values of MeCN, MeOH, and DMSO in water correlates well with the corresponding proton solvation energies in these liquids, indicating that the solvated proton was attached to a single solvent molecule.

Structured water molecules in the binding site of bromodomains can be displaced by cosolvent.

Bromodomains are α-helical bundles of approximately 110 residues that recognize acetylated lysine side chains mainly on histone tails. Bromodomains are known to play an important role in cancer and inflammation, and as such, significant efforts are being made to identify small-molecule inhibitors of these epigenetic reader proteins. Here, explicit solvent molecular dynamics (MD) simulations of two bromodomains (BAZ2B and CREBBP) are used to analyze the water molecules that seem to be conserved at the bottom of the acetyl-lysine binding site in most crystal structures of bromodomains. The MD runs suggest that the occupancy of the structured water molecules is influenced by conformational transitions of the loop that connects helices Z and A. Additional simulations in the presence of 50 molecules of cosolvent (i.e., 440 mM of dimethylsulfoxide, methanol, or ethanol) indicate that some of the structured water molecules can be displaced transiently. The residence time in the acetyl-lysine binding site is calculated to be about 1 ns, 2-5 ns, and 10-30 ns for methanol, ethanol, and dimethylsulfoxide, respectively, while the affinity of the three cosolvents for BAZ2B and CREBBP is in the range of 50-500 mM. The results described have implications for ligand design, suggesting that only structured water molecules that do not exchange with cosolvent should be maintained in crystal structures used for docking campaigns, and that hydroxy substituents should be incorporated in the ligand so as to map the structured water molecules replaced by (m)ethanol.