Conformational Energy Research Articles

Conformational analysis of acetamide in the ground and lowest excited electronic states

For acetamide molecule (CH3CONH2) in the ground (S0) and lowest excited singlet (S1) and triplet (T1) electronic states calculations of equilibrium geometry parameters, harmonic vibrational frequencies and barriers to conformational transitions (also conformer energy differences in excited states) using following ab initio methods: MP2, CCSD(T), CASSCF, CASPT2 and MRCI were performed. One-, two- and three-dimensional potential energy surface (PES) sections by different large amplitude motions (LAM) coordinates were calculated by means of MP2/aug-cc-pVTZ (S0) and CASPT2/cc-pVTZ (S1, T1). As a result of electronic excitation, both CCON and CNH2 fragments become pyramidal. On 2D PES sections by torsion (CN) and inversion coordinates there are six minima forming three pairs of enantiomers. Using PES sections different anharmonic vibrational problems were solved and the frequencies of large amplitude vibrations were estimated.

Glutamine Amide Flip Elicits Long Distance Allosteric Responses in the LOV Protein Vivid.

Light-oxygen-voltage (LOV) domains sense blue light through the photochemical formation of a cysteinyl-flavin covalent adduct. Concurrent protonation at the flavin N5 position alters the hydrogen bonding interactions of an invariant Gln residue that has been proposed to flip its amide side chain as a critical step in the propagation of conformational change. Traditional molecular dynamics (MD) and replica-exchange MD (REMD) simulations of the well-characterized LOV protein Vivid (VVD) demonstrate that the Gln182 amide indeed reorients by ∼180° in response to either adduct formation or reduction of the isoalloxazine ring to the neutral semiquinone, both of which involve N5 protonation. Free energy simulations reveal that the relative free energies of the flipped Gln conformation and the flipping barrier are significantly lower in the light-adapted state. The Gln182 flip stabilizes an important hinge-bβ region between the PAS β-sheet and the N-terminal cap helix that in turn destabilizes an N-terminal latch region against the PAS core. Release of the latch, observed both experimentally and in the simulations, is known to mediate light-induced VVD dimerization. This computational study of a LOV protein, unprecedented in its agreement with experiment, provides an atomistic view of long-range allosteric coupling in a photoreceptor.

Synthesis, DNA binding and in-vitro cytotoxicity studies on novel bis-pyrazoles

A new series of bis-pyrazoles 6a-t were synthesized from 3,5-dimethyl pyrazole using sequential approach. All these compounds were characterized by IR, 1H NMR, 13C NMR and mass spectral data. The interaction of newly synthesized bis-pyrazoles with DNA was investigated through molecular docking and absorption spectroscopic technique. Among all bis-pyrazoles compounds, the 6h compound showed lower conformational energy through in silico analysis. The interaction of each molecule in this series 6a-t with the various concentrations of DNA was examined through the UV–visible spectroscopic studies. The UV–visible spectroscopy studies on the specific binding of compound 6a, 6b, 6g, 6h, 6d, 6i, 6k, 6n, 6s with DNA have exhibited spectral shifts and the results were discussed. In further the compounds 6a-t were subjected to the in-vitro cytotoxicity studies against human pancreatic adenocarcinoma, human non-small cell lung carcinoma cell lines. Among the screened compounds, N-(3-isopropoxy-1-isopropyl-4-(3,5-dimethyl-1H-pyrazol-1-yl)-1H-pyrazol-5-yl)cyclobutane carboxamide and N-(5′-Isopropoxy-2′-isopropyl-3,5-dimethyl-2′H-[1,4′] bipyrazolyl-3′-yl)-dimethane sulfonamide were found as lead molecules since they have exhibited promising activity against both the cancer cell lines used in this study, whereas the compounds 4-(trifluoromethyl)-N-(3-isopropoxy-1-isopropyl-4-(3,5-dimethyl-2H-pyrrol-2-yl)-1H-pyrazol-5-yl)benzamide and 2,6-difluoro-N-(3-isopropoxy-1-isopropyl-4-(3,5-dimethyl-2H-pyrrol-2-yl)-1H-pyrazol-5-yl) benzamide were found to be active against the pancreatic cell line only. Rest all the other compounds were found to exhibit moderate to good activity towards both the cell lines.

Temperature-Induced Misfolding in Prion Protein: Evidence of Multiple Partially Disordered States Stabilized by Non-Native Hydrogen Bonds.

The structural basis of pathways of misfolding of a cellular prion (PrPC) into the toxic scrapie form (PrPSC) and identification of possible intermediates (e.g., PrP*) still eludes us. In this work, we have used a cumulative ∼65 μs of replica exchange molecular dynamics simulation data to construct the conformational free energy landscapes and capture the structural and thermodynamic characteristics associated with various stages of the thermal denaturation process in human prion protein. The temperature-dependent free energy surfaces consist of multiple metastable states stabilized by non-native contacts and hydrogen bonds, thus rendering the protein prone to misfolding. We have been able to identify metastable conformational states with high β-content (∼30-40%) and low α-content (∼10-20%) that might be precursors of PrPSC oligomer formation. These conformations also involve participation of the unstructured N-terminal domain, and its role in misfolding has been investigated. All the misfolded or partially unfolded states are quite compact in nature despite having large deviations from the native structure. Although the number of native contacts decreases dramatically at higher temperatures, the radius of gyration and number of intraprotein hydrogen bonds and contacts remain relatively unchanged, leading to stabilization of the misfolded conformations by non-native interactions. Our results are in good agreement with the established view that the C-terminal regions of the second and third helices (H2 and H3, respectively) of mammal prions might be the Achilles heels of their stability, while separation of B1-H1-B2 and H2-H3 domains seems to play a key role, as well.

Protein-Protein Complex Structure Prediction using the Solution Theory in the Energy Representation

Proteins conduct their functions through interactions with other molecules. Thus, accurate prediction of protein complex is a key to understand functional mechanisms. We have developed a method to evaluate binding free energy differences of complex models generated by docking prediction using the all-atom molecular dynamics simulation and solution theory in the energy representation. By evaluating binding free energy difference with the method, we previously show that “near-native” models similar to crystal structure are successfully selected as the lowest energy structures in two protein-protein complexes, bovine trypsin with CMTI-1 squash inhibitor and RNase SA with barstar [Takemura, J. Chem. Phys. 2012]. The method requires relatively short MD simulations (2 ns) and can calculate binding free energy difference among several hundreds of complex models. As a protein-protein complex structure prediction method, we combined the evaluation method with CyClus [Omori, Proteins, 2013] that performs fast clustering/reranking using a cylindrical approximation of interface and improves the results of the rigid-body docking.The initial complex models were prepared using the protein-protein rigid body docking program, ZDOCK. The generated complex models were then clustered and reranked using CyClus. Free energy analysis for top 100 or 300 complex models returned by CyClus were conducted using conformational energies, solute entropies, and solvation free energies calculated with the solution theory in the energy representation. Our analysis improved the results of CyClus and the complex models similar to the native structure have the lowest binding free energies, suggesting that this procedure is effective in protein-protein complex structure predictions. We further improved the procedure by placing interface waters into protein-protein interface before free energy evaluations.

Physical Binding of a Tobacco-Specific Carcinogen (NNK) Metabolite to the Human TP53 Gene with Ramifications for DNA Damage and Mutations

The most prevalent and carcinogenic of the tobacco-specific nitrosamines is 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK). NNK is a potent carcinogen that causes numerous types of cancers, most commonly lung adenomas and adenocarcinomas. The diazonium ion formed through the enzymatic metabolism of NNK by cytochrome P450s, can bind chemically to DNA (DNA damage). These adducts, if not repaired, can cause mutations that lead to tumorigenesis. We do not have a complete understanding of the binding of the diazonium ion to DNA which is important for future preventative strategies of NNK-associated cancers. The focus of this research is to computationally model and understand the physical interactions between the NNK diazonium ion and exon 5 of the human TP53 gene. We used a Lamarckian genetic algorithm to minimize a physics-based free energy function that allows us to estimate the free energy differences between the NNK diazonium ion along with TP53 in an unbound state and when they are physically bound. The free energy function employs pair-wise terms that evaluate dispersion, Pauli repulsions, hydrogen bonding, electrostatics and desolvation, along with a term for the lost conformational energy when the molecules bind. We will present the most probable NNK diazonium ion-TP53 physical binding sites and discuss how DNA damage (pyridyloxobutylated adducts in particular) may be occurring when it is exposed to this tobacco-specific carcinogen.

Colossal Strain Release by Conformational Energy Up-Conversion in a Compressed Molecular Crystal

External stimuli trigger conformational up-conversion in the molecules of di-p-tolyl disulfide (CH3–C6H4–S−)2, compensating the stress and absorbing its energy. These mechanochemical transformations explain at the molecular level the lubricating performance of di-p-tolyl disulfide and suggest new applications for storing energy. High pressure induces the conversion either in the solid state, as a phase transition, or at the nucleation stage on the solution-solid interphase as a detour of crystallization preference. The discontinuous transition at 1.6 GPa reduces the symmetry of ambient-pressure monoclinic phase α (space group P21) to triclinic phase β (space group P1). The recrystallizations above 0.45 GPa yield a new polymorph γ (space group P21/c), most spectacular in the conformational up-conversion absorbing 6 kJmol–1.

Revised RNA Dihedral Parameters for the Amber Force Field Improve RNA Molecular Dynamics.

The backbone dihedral parameters of the Amber RNA force field were improved by fitting using multiple linear regression to potential energies determined by quantum chemistry calculations. Five backbone and four glycosidic dihedral parameters were fit simultaneously to reproduce the potential energies determined by a high-level density functional theory calculation (B97D3 functional with the AUG-CC-PVTZ basis set). Umbrella sampling was used to determine conformational free energies along the dihedral angles, and these better agree with the population of conformations observed in the protein data bank for the new parameters than for the conventional parameters. Molecular dynamics simulations performed on a set of hairpin loops, duplexes and tetramers with the new parameter set show improved modeling for the structures of tetramers CCCC, CAAU, and GACC, and an RNA internal loop of noncanonical pairs, as compared to the conventional parameters. For the tetramers, the new parameters largely avoid the incorrect intercalated structures that dominate the conformational samples from the conventional parameters. For the internal loop, the major conformation solved by NMR is stable with the new parameters, but not with the conventional parameters. The new force field performs similarly to the conventional parameters for the UUCG and GCAA hairpin loops and the [U(UA)6A]2 duplex.

Ensemble simulations: folding, unfolding and misfolding of a high-efficiency frameshifting RNA pseudoknot

Massive all-atom molecular dynamics simulations were conducted across a distributed computing network to study the folding, unfolding, misfolding and conformational plasticity of the high-efficiency frameshifting double mutant of the 26 nt potato leaf roll virus RNA pseudoknot. Our robust sampling, which included over 40 starting structures spanning the spectrum from the extended unfolded state to the native fold, yielded nearly 120 μs of cumulative sampling time. Conformational microstate transitions on the 1.0 ns to 10.0 μs timescales were observed, with post-equilibration sampling providing detailed representations of the conformational free energy landscape and the complex folding mechanism inherent to the pseudoknot motif. Herein, we identify and characterize two alternative native structures, three intermediate states, and numerous misfolded states, the latter of which have not previously been characterized via atomistic simulation techniques. While in line with previous thermodynamics-based models of a general RNA folding mechanism, our observations indicate that stem-strand-sequence-separation may serve as an alternative predictor of the order of stem formation during pseudoknot folding. Our results contradict a model of frameshifting based on structural rigidity and resistance to mechanical unfolding, and instead strongly support more recent studies in which conformational plasticity is identified as a determining factor in frameshifting efficiency.

Rigidity and flexibility in the tetrasaccharide linker of proteoglycans from atomic-resolution molecular simulation.

Proteoglycans (PGs) are covalent conjugates between protein and carbohydrate (glycosaminoglycans). Certain classes of glycosaminoglycans such as chondroitin sulfate/dermatan sulfate and heparan sulfate utilize a specific tetrasaccharide linker for attachment to the protein component: GlcAβ1-3Galβ1-3Galβ1-4Xylβ1-O-Ser. Toward understanding the conformational preferences of this linker, the present work used all-atom explicit-solvent molecular dynamics (MD) simulations combined with Adaptive Biasing Force (ABF) sampling to determine high-resolution, high-precision conformational free energy maps ΔG(φ, ψ) for each glycosidic linkage between constituent disaccharides, including the variant where GlcA is substituted with IdoA. These linkages are characterized by single, predominant (> 97% occupancy), and broad (45° × 60° for ΔG(φ, ψ) < 1 kcal/mol) free-energy minima, while the Xyl-Ser linkage has two such minima similar in free-energy, and additional flexibility from the Ser sidechain dihedral. Conformational analysis of microsecond-scale standard MD on the complete tetrasaccharide-O-Ser conjugate is consistent with ABF data, suggesting (φ, ψ) probabilities are independent of the linker context, and that the tetrasaccharide acts as a relatively rigid unit whereas significant conformational heterogeneity exists with respect to rotation about bonds connecting Xyl to Ser. © 2017 Wiley Periodicals, Inc.

Computational Prediction and Rationalization, and Experimental Validation of Handedness Induction in Helical Aromatic Oligoamide Foldamers.

Metadynamics simulations were used to describe the conformational energy landscapes of several helically folded aromatic quinoline carboxamide oligomers bearing a single chiral group at either the C or N terminus. The calculations allowed the prediction of whether a helix handedness bias occurs under the influence of the chiral group and gave insight into the interactions (sterics, electrostatics, hydrogen bonds) responsible for a particular helix sense preference. In the case of camphanyl-based and morpholine-based chiral groups, experimental data confirming the validity of the calculations were already available. New chiral groups with a proline residue were also investigated and were predicted to induce handedness. This prediction was verified experimentally through the synthesis of proline-containing monomers, their incorporation into an oligoamide sequence by solid phase synthesis and the investigation of handedness induction by NMR spectroscopy and circular dichroism.

Computational Design of Experiment Unveils the Conformational Reaction Coordinate of GH125 α-Mannosidases.

Conformational analysis of enzyme-catalyzed mannoside hydrolysis has revealed two predominant conformational itineraries through B2,5 or 3H4 transition-state (TS) conformations. A prominent unassigned catalytic itinerary is that of exo-1,6-α-mannosidases belonging to CAZy family 125. A published complex of Clostridium perfringens GH125 enzyme with a nonhydrolyzable 1,6-α-thiomannoside substrate mimic bound across the active site revealed an undistorted 4C1 conformation and provided no insight into the catalytic pathway of this enzyme. We show through a purely computational approach (QM/MM metadynamics) that sulfur-for-oxygen substitution in the glycosidic linkage fundamentally alters the energetically accessible conformational space of a thiomannoside when bound within the GH125 active site. Modeling of the conformational free energy landscape (FEL) of a thioglycoside strongly favors a mechanistically uninformative 4C1 conformation within the GH125 enzyme active site, but the FEL of corresponding O-glycoside substrate reveals a preference for a Michaelis complex in an OS2 conformation (consistent with catalysis through a B2,5 TS). This prediction was tested experimentally by determination of the 3D X-ray structure of the pseudo-Michaelis complex of an inactive (D220N) variant of C. perfringens GH125 enzyme in complex with 1,6-α-mannobiose. This complex revealed unambiguous distortion of the -1 subsite mannoside to an OS2 conformation, matching that predicted by theory and supporting an OS2 → B2,5 → 1S5 conformational itinerary for GH125 α-mannosidases. This work highlights the power of the QM/MM approach and identified shortcomings in the use of nonhydrolyzable substrate analogues for conformational analysis of enzyme-bound species.

Intramolecular backbone···backbone hydrogen bonds in polypeptide conformations. The other way around: ɛ-turn.

In this study, we performed a detailed literature survey of the ɛ-turn in peptides and proteins. This three-dimensional structural feature is characterized by an eleven-membered pseudo-cycle closed by an intramolecular backbone…backbone H-bond. Interestingly, in this motif the direction of the N-H…O = C H-bond runs opposite to that of the much more popular and extensively investigated α-, β-, and γ-turns. We did not authenticate unequivocally the ɛ-turn main-chain reversal topology in any linear short peptide. However, it is frequently observed in small cyclic peptides formed by four, five, and six amino acid residues with stringent geometric requirements. Rather surprisingly, ɛ-turns do occur in proteins, although to a relatively moderate extent, as an isolated feature or in the turn segment of hairpin motifs based on two antiparallel, pleated β-strands. Moreover, the ɛ-turn may also host not only the seven-membered, intramolecularly H-bonded, pseudo-cycle termed γ-turn, either of the classic or inverse type, but also one (or even two) cis peptide bond(s) or a β-bulge conformation. Based on their ϕ, ψ backbone torsion angles, we were able to classify the protein ɛ-turns in six different families. Conformational energy computations using the DFT methodology were also performed on the ɛ-turns adopted by the amino acid triplet -Gly-Gly-Gly- (Gly is the most commonly found residue at each of the three positions in our analysis of proteins). Again, in this computational study, six families of turns were identified, but only some of them resemble rather closely those extracted from our investigation on proteins.

Structure and conformational analysis of 2-hydroxy-5-isobutyl-1,3,2-dioxaborinane

The structure of 2-hydroxy-5-isobutyl-1,3,2-dioxaborinane, one of the major autooxidation products of 5-isobutyl-2-isopropyl-1,3,2-dioxaborinane, has been studied by 1H, 13C, and 11B NMR spectroscopy together with X-ray analysis. Molecules of the title compounds adopt a sofa conformation with equatorial isobutyl substituent. The ring interconversion path, free conformational energy, and optimal conformation of the isobutyl group (corresponding to the X-ray diffraction data) were determined by DFT quantum chemical calculations at the PBE/3ξ level of theory.

Conformational studies of 3-[(2,2-dimethylhydrazinyl)methylene]-pentane-2,4-dione using vibrational and NMR spectra, X-ray analysis and ab initio calculations

The conformers of push–pull 3-[(2,2-dimethylhydrazinyl)methylene]-pentane-2,4-dione (CH3)2NNHCHC(COCH3)2 (DMHMP) have been studied experimentally by NMR and vibrational spectroscopy and theoretically by ab initio calculations at MP2 and DFT B3LYP levels in various basis sets. The NMR spectra were obtained in chloroform and dimethylsulfoxide and the IR and Raman spectra of DMHMP as a solid and as a solute in various less and more polar solvents at room temperature have been recorded.DMHMP was prepared as a pure solid and the data from X-ray analysis revealed that DMHMP exists in solid state as EZa conformer with an intramolecular hydrogen bond. The geometries and relative energies of possible conformers of DMHMP were evaluated at the both levels of theory in several basis sets and compared with the data from X-ray analysis.According to the NMR spectra the studied compound exists as a single entity. On the other hand vibrational spectra revealed that in less polar DMHMP solutions the presence of the second less polar ZZa conformer is possible, whereas in more polar solvent only one EZa conformer is observed. The influence of the environment polarity on this conformational equilibrium is discussed with respect to the SCRF solvent effect calculations using IEFPCM model. The observed IR and Raman bands were compared with calculated MP2/cc-pVTZ harmonic vibrational frequencies and assigned on the basis of potential energy distribution.

Computing conformational free energy differences in explicit solvent: An efficient thermodynamic cycle using an auxiliary potential and a free energy functional constructed from the end points.

Many biomolecules undergo conformational changes associated with allostery or ligand binding. Observing these changes in computer simulations is difficult if their timescales are long. These calculations can be accelerated by observing the transition on an auxiliary free energy surface with a simpler Hamiltonian and connecting this free energy surface to the target free energy surface with free energy calculations. Here, we show that the free energy legs of the cycle can be replaced with energy representation (ER) density functional approximations. We compute: (1) The conformational free energy changes for alanine dipeptide transitioning from the right-handed free energy basin to the left-handed basin and (2) the free energy difference between the open and closed conformations of β-cyclodextrin, a "host" molecule that serves as a model for molecular recognition in host-guest binding. β-cyclodextrin contains 147 atoms compared to 22 atoms for alanine dipeptide, making β-cyclodextrin a large molecule for which to compute solvation free energies by free energy perturbation or integration methods and the largest system for which the ER method has been compared to exact free energy methods. The ER method replaced the 28 simulations to compute each coupling free energy with two endpoint simulations, reducing the computational time for the alanine dipeptide calculation by about 70% and for the β-cyclodextrin by > 95%. The method works even when the distribution of conformations on the auxiliary free energy surface differs substantially from that on the target free energy surface, although some degree of overlap between the two surfaces is required. © 2016 Wiley Periodicals, Inc.

Biasing Simulations of DNA Base Pair Parameters with Application to Propellor Twisting in AT/AT, AA/TT, and AC/GT Steps and Their Uracil Analogs.

An accurate and efficient implementation of the six DNA base pair parameters as order parameters for enhanced sampling simulations is presented. The parameter definitions are defined by vector algebra operations on a reduced atomic set of the base pair, and correlate very well with standard definitions. Application of the model is illustrated by umbrella sampling simulations of propeller twisting within AT/AT, AA/TT, and AC/GT steps and their uracil analogs. Strong correlations are found between propeller twisting and a number of conformational parameters, including buckle, opening, BI/BII backbone configuration, and sugar puckering. The thymine methyl group is observed to notably alter the local conformational free energy landscape, with effects within and directly upstream of the thymine containing base pair.

En route towards the peptide γ-helix: X-ray diffraction analyses and conformational energy calculations of Adm-rich short peptides.

We performed the solution-phase synthesis of a set of model peptides, including homo-oligomers, based on the 2-aminoadamantane-2-carboxylic acid (Adm) residue, an extremely bulky, highly lipophilic, tricyclic, achiral, Cα -tetrasubstituted α-amino acid. In particular, for the difficult peptide coupling reaction between two Adm residues, we took advantage of the Meldal's α-azidoacyl chloride approach. Most of the synthesized Adm peptides were characterized by single-crystal X-ray diffraction analyses. The results indicate a significant propensity for the Adm residue to adopt γ-turn and γ-turn-like conformations. Interestingly, we found that a -CO-(Adm)2 -NH- sequence is folded in the crystal state into a regular, incipient γ-helix, at variance with the behavior of all of the homo-dipeptides from Cα -tetrasubstituted α-amino acids already investigated, which tend to adopt either the β-turn or the fully extended conformation. Our density functional theory conformational energy calculations on the terminally blocked homo-peptides (n=2-8) fully confirmed the crystal-state data, strongly supporting the view that this rigid Cα -tetrasubstituted α-amino acid residue is largely the most effective building block for γ-helix induction, although to a limited length (anti-cooperative effect). Copyright © 2016 European Peptide Society and John Wiley & Sons, Ltd.

Multiconfiguration Pair-Density Functional Theory: A New Way To Treat Strongly Correlated Systems.

The electronic energy of a system provides the Born-Oppenheimer potential energy for internuclear motion and thus determines molecular structure and spectra, bond energies, conformational energies, reaction barrier heights, and vibrational frequencies. The development of more efficient and more accurate ways to calculate the electronic energy of systems with inherently multiconfigurational electronic structure is essential for many applications, including transition metal and actinide chemistry, systems with partially broken bonds, many transition states, and most electronically excited states. Inherently multiconfigurational systems are called strongly correlated systems or multireference systems, where the latter name refers to the need for using more than one ("multiple") configuration state function to provide a good zero-order reference wave function. This Account describes multiconfiguration pair-density functional theory (MC-PDFT), which was developed as a way to combine the advantages of wave function theory (WFT) and density functional theory (DFT) to provide a better treatment of strongly correlated systems. First we review background material: the widely used Kohn-Sham DFT (which uses only a single Slater determinant as reference wave function), multiconfiguration WFT methods that treat inherently multiconfigurational systems based on an active space, and previous attempts to combine multiconfiguration WFT with DFT. Then we review the formulation of MC-PDFT. It is a generalization of Kohn-Sham DFT in that the electron kinetic energy and classical electrostatic energy are calculated from a reference wave function, while the rest of the energy is obtained from a density functional. However, there are two main differences with respent to Kohn-Sham DFT: (i) The reference wave function is multiconfigurational rather than being a single Slater determinant. (ii) The density functional is a function of the total density and the on-top pair density rather than being a function of the spin-up and spin-down densities. In work carried out so far, the multiconfigurational wave function is a multiconfiguration self-consistent-field wave function. The new formulation has the advantage that the reference wave function has the correct spatial and spin symmetry and can describe bond dissociation (of both single and multiple bonds) and electronic excitations in a formally and physically correct way. We then review the formulation of density functionals in terms of the on-top pair density. Finally we review successful applications of the theory to bond energies and bond dissociation potential energy curves of main-group and transition metal bonds, to barrier heights (including pericyclic reactions), to proton affinities, to the hydrogen bond energy of water dimer, to ground- and excited-state charge transfer, to valence and Rydberg excitations of molecules, and to singlet-triplet splittings of radicals. We find that that MC-PDFT can give accurate results not only with complete-active-space multiconfiguration wave functions but also with generalized-active-space multiconfiguration wave functions, which are practical for larger numbers of active electrons and active orbitals than are complete-active-space wave functions. The separated-pair approximation, which is a special case of generalized active space self-consistent-field theory, is especially promising. MC-PDFT, because it requires much less computer time and storage than pure WFT methods, has the potential to open larger and more complex strongly correlated systems to accurate simulation.

Effect of alcohol on the structure of cytochrome C: FCS and molecular dynamics simulations.

Effect of ethanol on the size and structure of a protein cytochrome C (Cyt C) is investigated using fluorescence correlation spectroscopy (FCS) and molecular dynamics (MD) simulations. For FCS studies, Cyt C is covalently labeled with a fluorescent probe, alexa 488. FCS studies indicate that on addition of ethanol, the size of the protein varies non-monotonically. The size of Cyt C increases (i.e., the protein unfolds) on addition of alcohol (ethanol) up to a mole fraction of 0.2 (44.75% v/v) and decreases at higher alcohol concentration. In order to provide a molecular origin of this structural transition, we explore the conformational free energy landscape of Cyt C as a function of radius of gyration (Rg) at different compositions of water-ethanol binary mixture using MD simulations. Cyt C exhibits a minimum at Rg ∼ 13 Å in bulk water (0% alcohol). Upon increasing ethanol concentration, a second minimum appears in the free energy surface with gradually larger Rg up to χEtOH ∼ 0.2 (44.75% v/v). This suggests gradual unfolding of the protein. At a higher concentration of alcohol (χEtOH > 0.2), the minimum at large Rg vanishes, indicating compaction. Analysis of the contact map and the solvent organization around protein indicates a preferential solvation of the hydrophobic residues by ethanol up to χEtOH = 0.2 (44.75% v/v) and this causes the gradual unfolding of the protein. At high concentration (χEtOH = 0.3 (58% v/v)), due to structural organization in bulk water-ethanol binary mixture, the extent of preferential solvation by ethanol decreases. This causes a structural transition of Cyt C towards a more compact state.