Directional hydrogen bonding in the MM3 force field: II

Conformational energies of reference organic molecules: benchmarking of common efficient computational methods against coupled cluster theory

We selected 145 reference organic molecules that include model fragments used in computer-aided drug design. We calculated 158 conformational energies and barriers using force fields, with wide applicability in commercial and free softwares and extensive application on the calculation of conformational energies of organic molecules, e.g. the UFF and DREIDING force fields, the Allinger's force fields MM3-96, MM3-00, MM4-8, the MM2-91 clones MMX and MM+, the MMFF94 force field, MM4, ab initio Hartree-Fock (HF) theory with different basis sets, the standard density functional theory B3LYP, the second-order post-HF MP2 theory and the Domain-based Local Pair Natural Orbital Coupled Cluster DLPNO-CCSD(T) theory, with the latter used for accurate reference values. The data set of the organic molecules includes hydrocarbons, haloalkanes, conjugated compounds, and oxygen-, nitrogen-, phosphorus- and sulphur-containing compounds. We reviewed in detail the conformational aspects of these model organic molecules providing the current understanding of the steric and electronic factors that determine the stability of low energy conformers and the literature including previous experimental observations and calculated findings. While progress on the computer hardware allows the calculations of thousands of conformations for later use in drug design projects, this study is an update from previous classical studies that used, as reference values, experimental ones using a variety of methods and different environments. The lowest mean error against the DLPNO-CCSD(T) reference was calculated for MP2 (0.35 kcal mol-1), followed by B3LYP (0.69 kcal mol-1) and the HF theories (0.81-1.0 kcal mol-1). As regards the force fields, the lowest errors were observed for the Allinger's force fields MM3-00 (1.28 kcal mol-1), ΜΜ3-96 (1.40 kcal mol-1) and the Halgren's MMFF94 force field (1.30 kcal mol-1) and then for the MM2-91 clones MMX (1.77 kcal mol-1) and MM+ (2.01 kcal mol-1) and MM4 (2.05 kcal mol-1). The DREIDING (3.63 kcal mol-1) and UFF (3.77 kcal mol-1) force fields have the lowest performance. These model organic molecules we used are often present as fragments in drug-like molecules. The values calculated using DLPNO-CCSD(T) make up a valuable data set for further comparisons and for improved force field parameterization.

Construction of a Gaussian Process Regression Model of Formamide for Use in Molecular Simulations

FFLUX, a novel force field based on quantum chemical topology, can perform molecular dynamics simulations with flexible multipole moments that change with geometry. This is enabled by Gaussian process regression machine learning models, which accurately predict atomic energies and multipole moments up to the hexadecapole. We have constructed a model of the formamide monomer at the B3LYP/aug-cc-pVTZ level of theory capable of sub-kJ mol^-1 accuracy, with the maximum prediction error for the molecule being 0.8 kJ mol^-1. This model was used in FFLUX simulations along with Lennard-Jones parameters to successfully optimize the geometry of formamide dimers with errors smaller than 0.1 Å compared to those obtained with D3-corrected B3LYP/aug-cc-pVTZ. Comparisons were also made to a force field constructed with static multipole moments and Lennard-Jones parameters. FFLUX recovers the expected energy ranking of dimers compared to the literature, and changes in C═O and C-N bond lengths associated with hydrogen bonding were found to be consistent with density functional theory.

Tinker 8: Software Tools for Molecular Design

The Tinker software, currently released as version 8, is a modular molecular mechanics and dynamics package written primarily in a standard, easily portable dialect of Fortran 95 with OpenMP extensions. It supports a wide variety of force fields, including polarizable models such as the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field. The package runs on Linux, macOS, and Windows systems. In addition to canonical Tinker, there are branches, Tinker-HP and Tinker-OpenMM, designed for use on message passing interface (MPI) parallel distributed memory supercomputers and state-of-the-art graphical processing units (GPUs), respectively. The Tinker suite also includes a tightly integrated Java-based graphical user interface called Force Field Explorer (FFE), which provides molecular visualization capabilities as well as the ability to launch and control Tinker calculations.

Molecular self-assembly of substituted terephthalic acids at the liquid/solid interface: investigating the effect of solvent

Self-assembly of three related molecules – terephthalic acid and its hydroxylated analogues – at liquid/solid interfaces (graphite/heptanoic acid and graphite/1-phenyloctane) has been studied using a combination of scanning tunnelling microscopy and molecular mechanics and molecular dynamics calculations. Brickwork-like patterns typical for terephthalic acid self-assembly have been observed for all three molecules. However, several differences became apparent: (i) formation or lack of adsorbed monolayers (self-assembled monolayers formed in all systems, with one notable exception of terephthalic acid at the graphite/1-phenyloctane interface where no adsorption was observed), (ii) the size of adsorbate islands (large islands at the interface with heptanoic acid and smaller ones at the interface with 1-phenyloctane), and (iii) polymorphism of the hydroxylated terephthalic acids’ monolayers, dependent on the molecular structure and/or solvent. To rationalise this behaviour, molecular mechanics and molecular dynamics calculations have been performed, to analyse the three key aspects of the energetics of self-assembly: intermolecular, substrate–adsorbate and solvent–solute interactions. These energetic characteristics of self-assembly were brought together in a Born–Haber cycle, to obtain the overall energy effects of formation of self-assembled monolayers at these liquid/solid interfaces.

Dioxygen binding to Fe-MOF-74: microscopic insights from periodic QM/MM calculations

Accurate MOF-FF parameter sets were determined for the ferrous and ferric forms of an iron-based metal–organic framework (MOF) called Fe-MOF-74. For this purpose, density functional theory (DFT) calculations were applied to truncated cluster models of Fe-MOF-74, and the DFT-calculated geometries and energy derivatives were used for the force-field parameterization. The resultant parameter sets performed remarkably well in reproducing the experimentally determined structure of the MOF. We also performed periodic quantum mechanics (QM) / molecular mechanics (MM) calculations employing a subtractive scheme called ONIOM, with the optimized MOF-FF parameters used for the MM calculations, in an attempt to evaluate the binding energies between O2 and several Fe-MOF-74 variants. The calculated binding energy for Fe-MOF-74 agreed very well with the experimental value, and QM/MM geometry optimization calculations confirmed that the O2-bound complex has a side-on geometry. Our calculations also predicted that, when the two neighboring iron ions around the O2-binding site are replaced with other metal ions (Mg2+, Ni2+, Zn2+, Co2+, or Mn2+), there are noticeable variations in the binding energy, indicating that these substituted metal ions affect the O2 binding indirectly.

Thermodynamics of 4,4'-stilbenedicarboxylic acid monolayer self-assembly at the nonanoic acid-graphite interface

A direct calorimetric measurement of the overall enthalpy change associated with self-assembly of organic monolayers at the liquid-solid interface is for most systems of interest practically impossible. In previous work we proposed an adapted Born-Haber cycle for an indirect assessment of the overall enthalpy change by using terephthalic acid monolayers at the nonanoic acid-graphite interface as a model system. To this end, the sublimation enthalpy, dissolution enthalpy, the monolayer binding enthalpy in vacuum, and a dewetting enthalpy are combined to yield the total enthalpy change. In the present study the Born-Haber cycle is applied to 4,4'-stilbenedicarboxylic acid monolayers. A detailed comparison of these two aromatic dicarboxylic acids is used to evaluate and quantify the contribution of the organic backbone for stabilization of the monolayer at the nonanoic acid-graphite interface.

Combined covalent-electrostatic model of hydrogen bonding improves structure prediction with Rosetta

Interactions between polar atoms are challenging to model because at very short ranges they form hydrogen bonds (H-bonds) that are partially covalent in character and exhibit strong orientation preferences; at longer ranges the orientation preferences are lost, but significant electrostatic interactions between charged and partially charged atoms remain. To simultaneously model these two types of behavior, we refined an orientation dependent model of hydrogen bonds [Kortemme et al. J. Mol. Biol. 2003, 326, 1239] used by the molecular modeling program Rosetta and then combined it with a distance-dependent Coulomb model of electrostatics. The functional form of the H-bond potential is physically motivated and parameters are fit so that H-bond geometries that Rosetta generates closely resemble H-bond geometries in high-resolution crystal structures. The combined potentials improve performance in a variety of scientific benchmarks including decoy discrimination, side chain prediction, and native sequence recovery in protein design simulations and establishes a new standard energy function for Rosetta.

Accurate prediction of polarised high order electrostatic interactions for hydrogen bonded complexes using the machine learning method kriging

As intermolecular interactions such as the hydrogen bond are electrostatic in origin, rigorous treatment of this term within force field methodologies should be mandatory. We present a method able of accurately reproducing such interactions for seven van der Waals complexes. It uses atomic multipole moments up to hexadecupole moment mapped to the positions of the nuclear coordinates by the machine learning method kriging. Models were built at three levels of theory: HF/6-31G(**), B3LYP/aug-cc-pVDZ and M06-2X/aug-cc-pVDZ. The quality of the kriging models was measured by their ability to predict the electrostatic interaction energy between atoms in external test examples for which the true energies are known. At all levels of theory, >90% of test cases for small van der Waals complexes were predicted within 1 kJ mol(-1), decreasing to 60-70% of test cases for larger base pair complexes. Models built on moments obtained at B3LYP and M06-2X level generally outperformed those at HF level. For all systems the individual interactions were predicted with a mean unsigned error of less than 1 kJ mol(-1).

A new pH sensitive fluorescent and white light emissive material through controlled intermolecular charge transfer

A new, pH dependent and water-soluble, conjugated oligomer (amino, trimethylammonium oligophenylene vinylene, ATAOPV) was synthesized with a quaternary ammonium salt and an aromatic amine at the two ends of a π-conjugated oligomer, thus creating a strong dipole across the molecule. A unique white light LED is successfully fabricated from a stimuli responsive organic molecule whose emission properties are dominated by the pH value of the solution through controlled intermolecular charge transfer.

Prediction of enhanced solvent-induced enantioselectivity for a ring opening with a bifurcating reaction path

A fully atomistic molecular dynamics simulation predicts enhance induction of enantiomeric excess in the products of a reaction with a bifurcating reaction coordinate, when run in a chiral solvent.

Shear and friction between carbon nanotubes in bundles and yarns

We perform a detailed density functional theory assessment of the factors that determine shear interactions between carbon nanotubes (CNTs) within bundles and in related CNT and graphene structures including yarns, providing an explanation for the shear force measured in recent experiments (Filleter, T. etal. Nano Lett. 2012, 12, 73). The potential energy barriers separating AB stacked structures are found to be irrelevant to the shear analysis for bundles and yarns due to turbostratic stacking, and as a result, the tube-tube shear strength for pristine CNTs is estimated to be <0.24 MPa, that is, extremely small. Instead, it is pinning due to the presence of defects and functional groups at the tube ends that primarily cause resistance to shear when bundles are fractured in weak vacuum (∼10(-5) Torr). Such defects and groups are estimated to involve 0.55 eV interaction energies on average, which is much larger than single-atom vacancy defects (approximately 0.039 eV). Furthermore, because graphitic materials are stiff they have large coherence lengths, and this means that push-pull effects result in force cancellation for vacancy and other defects that are internal to the CNTs. Another important factor is the softness of cantilever structures relative to the stiff CNTs in the experiments, as this contributes to elastic instability transitions that account for significant dissipation during shear that has been observed. The application of these results to the mechanical behavior of yarns is discussed, providing general guidelines for the manufacture of strong yarns composed of CNTs.

One-way rotation of a molecule-rotor driven by a shot noise

The shot noise of a tunneling current passing through a molecule-motor can sustain a one-way rotation when populating the molecular excited states by tunneling inelastic excitations. We demonstrate that a ratchet-like ground state rotation potential energy curve is not necessary for the rotation to occur. A relative shift in energy difference between the maxima of this ground state and the minima of the excited states is the necessary condition to get to a unidirectional rotation. The rotor speed of rotation and its rotation direction are both controlled by this shift, indicating the necessity of a careful design of both the ground and excited states of the next generation of molecule-motors to be able to generate a motive power at the nanoscale.

Born-Haber cycle for monolayer self-assembly at the liquid-solid interface: assessing the enthalpic driving force

The driving force for self-assembly is the associated gain in free energy with decisive contributions from both enthalpy and entropy differences between final and initial state. For monolayer self-assembly at the liquid-solid interface, solute molecules are initially dissolved in the liquid phase and then become incorporated into an adsorbed monolayer. In this work, we present an adapted Born-Haber cycle for obtaining precise enthalpy values for self-assembly at the liquid-solid interface, a key ingredient for a profound thermodynamic understanding of this process. By choosing terephthalic acid as a model system, it is demonstrated that all required enthalpy differences between well-defined reference states can be independently and consistently assessed by both experimental and theoretical methods, giving in the end a reliable value of the overall enthalpy gain for self-assembly of interfacial monolayers. A quantitative comparison of enthalpy gain and entropy cost reveals essential contributions from solvation and dewetting, which lower the entropic cost and render monolayer self-assembly a thermodynamically favored process.

Experimental-computational study of shear interactions within double-walled carbon nanotube bundles

The mechanical behavior of carbon nanotube (CNT)-based fibers and nanocomposites depends intimately on the shear interactions between adjacent tubes. We have applied an experimental-computational approach to investigate the shear interactions between adjacent CNTs within individual double-walled nanotube (DWNT) bundles. The force required to pull out an inner bundle of DWNTs from an outer shell of DWNTs was measured using in situ scanning electron microscopy methods. The normalized force per CNT-CNT interaction (1.7 ± 1.0 nN) was found to be considerably higher than molecular mechanics (MM)-based predictions for bare CNTs (0.3 nN). This MM result is similar to the force that results from exposure of newly formed CNT surfaces, indicating that the observed pullout force arises from factors beyond what arise from potential energy effects associated with bare CNTs. Through further theoretical considerations we show that the experimentally measured pullout force may include small contributions from carbonyl functional groups terminating the free ends of the CNTs, corrugation of the CNT-CNT interactions, and polygonization of the nanotubes due to their mutual interactions. In addition, surface functional groups, such as hydroxyl groups, that may exist between the nanotubes are found to play an unimportant role. All of these potential energy effects account for less than half of the ~1.7 nN force. However, partially pulled-out inner bundles are found not to pull back into the outer shell after the outer shell is broken, suggesting that dissipation is responsible for more than half of the pullout force. The sum of force contributions from potential energy and dissipation effects are found to agree with the experimental pullout force within the experimental error.

Direct evaluation of individual hydrogen bond energy in situ in intra- and intermolecular multiple hydrogen bonds system

The results of evaluating the individual hydrogen bond (H-bond) strength are expected to be helpful for the rational design of new strategies for molecular recognition or supramolecular assemblies. Unfortunately, there is few obvious and unambiguous means of evaluating the energy of a single H-bond within a multiple H-bonds system. We present a local analytic model, ABEEMσπ H-bond energy (HBE) model based on ab initio calculations (MP2) as benchmark, to directly and rapidly evaluate the individual HBE in situ in inter- and intramolecular multiple H-bonds system. This model describes the HBE as the sum of electrostatic and van der Waals (vdW) interactions which all depend upon the geometry and environment, and the ambient environment of H-bond in the model is accounted fairly. Thus, it can fairly consider the cooperative effect and secondary effect. The application range of ABEEMσπ HBE model is rather wide. This work has discussed the individual H-bond in DNA base pair and protein peptide dimers. The results indicate that the interactions among donor H atom, acceptor atom as well as those atoms connected to them with 1,2 or 1,3 relationships are all important for evaluating the HBE, although the interaction between the donor H atom and the acceptor atom is large. Furthermore, our model quantitatively indicates the polarization ability of N, O, and S in a new style, and gives the percentage of the polarization effect in HBE, which can not be given by fixed partial charge force field.

STRUCTURES AND ENERGETICS OF PORPHYRIN-FULLERENE SUPRAMOLECULAR COMPLEXES

In this paper, the structures and energetics of porphyrin-fullerene supramolecular complexes have been investigated theoretically via molecular mechanics calculations on complexes having different porphyrins. The results of the calculations allow us to delineate the effect of the substituents on the porphyrin ring on the porphyrin-fullerene interaction energies. The calculations also allow us to understand how the van der Waals force affects the structures of the 2:1, 2:2, and 3:1 supramolecular complexes.

STUDY ON COMPLEXES OF TRYPSIN AND ITS INHIBITORS BY MEANS OF ATOM-BOND ELECTRONEGATIVITY EQUALIZATION METHOD FUSED INTO MOLECULAR MECHANICS (ABEEM/MM)

Trypsin is one of the most important enzymes and plays important roles in the regulation of biological processes. The newly developed trypsin–inhibitor interaction potential model based on atom-bond electronegativity equalization method fused into molecular mechanics (ABEEM/MM) was employed to study complexes of trypsin and its inhibitors. Some structural properties, including root-mean-square deviations (RMSD) of bond length, bond angle and key dihedral, and coordinated RMS Shifts of atoms and hydrogen bond, were studied using ABEEM/MM method and compared with OPLS-AA force-field. At the same time, comparative study on the charges of the hydrogen atoms at the tail of the ligand was also discussed in order to investigate the effect of hydrogen bond between trypsin and its ligand. This work demonstrates that ABEEM/MM model can well reproduce good structures for these trypsin–inhibitor complexes with rather small RMSD of bond length, bond angles, key dihedrals, and RMS Shifts of atomic coordinates with respect to the experimental crystal structures, compared with the results from the OPLS-AA method. The charges obtained by ABEEM/MM model are in good accordance with those from an ab initio calculation. Moreover, both the polarization and the salt-bridge effects have been taken into account. It is shown that ABEEM charges can properly describe the electrostatic interactions between the protein trypsin and its inhibitors.

Incorporation dynamics of molecular guests into two-dimensional supramolecular host networks at the liquid-solid interface

The objective of this work is to study both the dynamics and mechanisms of guest incorporation into the pores of 2D supramolecular host networks at the liquid-solid interface. This was accomplished by adding molecular guests to prefabricated self-assembled porous monolayers and the simultaneous acquisition of scanning tunneling microscopy (STM) topographs. The incorporation of the same guest molecule (coronene) into two different host networks was compared, where the pores of the networks either featured a perfect geometric match with the guest (for trimesic acid host networks) or were substantially larger than the guest species (for benzenetribenzoic acid host networks). Even the moderate temporal resolution of standard STM experiments in combination with a novel injection system was sufficient to reveal clear differences in the incorporation dynamics in the two different host networks. Further experiments were aimed at identifying a possible solvent influence. The interpretation of the results is aided by molecular mechanics (MM) and molecular dynamics (MD) simulations.

Polarizable Atomic Multipole-based Molecular Mechanics for Organic Molecules

An empirical potential based on permanent atomic multipoles and atomic induced dipoles is reported for alkanes, alcohols, amines, sulfides, aldehydes, carboxylic acids, amides, aromatics and other small organic molecules. Permanent atomic multipole moments through quadrupole moments have been derived from gas phase ab initio molecular orbital calculations. The van der Waals parameters are obtained by fitting to gas phase homodimer QM energies and structures, as well as experimental densities and heats of vaporization of neat liquids. As a validation, the hydrogen bonding energies and structures of gas phase heterodimers with water are evaluated using the resulting potential. For 32 homo- and heterodimers, the association energy agrees with ab initio results to within 0.4 kcal/mol. The RMS deviation of hydrogen bond distance from QM optimized geometry is less than 0.06 Å. In addition, liquid self-diffusion and static dielectric constants computed from molecular dynamics simulation are consistent with experimental values. The force field is also used to compute the solvation free energy of 27 compounds not included in the parameterization process, with a RMS error of 0.69 kcal/mol. The results obtained in this study suggest the AMOEBA force field performs well across different environments and phases. The key algorithms involved in the electrostatic model and a protocol for developing parameters are detailed to facilitate extension to additional molecular systems.

Novel nonlinear knowledge-based mean force potentials based on machine learning

The prediction of 3D structures of proteins from amino acid sequences is one of the most challenging problems in molecular biology. An essential task for solving this problem with coarse-grained models is to deduce effective interaction potentials. The development and evaluation of new energy functions is critical to accurately modeling the properties of biological macromolecules. Knowledge-based mean force potentials are derived from statistical analysis of proteins of known structures. Current knowledge-based potentials are almost in the form of weighted linear sum of interaction pairs. In this study, a class of novel nonlinear knowledge-based mean force potentials is presented. The potential parameters are obtained by nonlinear classifiers, instead of relative frequencies of interaction pairs against a reference state or linear classifiers. The support vector machine is used to derive the potential parameters on data sets that contain both native structures and decoy structures. Five knowledge-based mean force Boltzmann-based or linear potentials are introduced and their corresponding nonlinear potentials are implemented. They are the DIH potential (single-body residue-level Boltzmann-based potential), the DFIRE-SCM potential (two-body residue-level Boltzmann-based potential), the FS potential (two-body atom-level Boltzmann-based potential), the HR potential (two-body residue-level linear potential), and the T32S3 potential (two-body atom-level linear potential). Experiments are performed on well-established decoy sets, including the LKF data set, the CASP7 data set, and the Decoys “R”Us data set. The evaluation metrics include the energy Z score and the ability of each potential to discriminate native structures from a set of decoy structures. Experimental results show that all nonlinear potentials significantly outperform the corresponding Boltzmann-based or linear potentials, and the proposed discriminative framework is effective in developing knowledge-based mean force potentials. The nonlinear potentials can be widely used for ab initio protein structure prediction, model quality assessment, protein docking, and other challenging problems in computational biology.

Agent-based modeling for the 2D molecular self-organization of realistic molecules

We extend our previously developed agent-based (AB) algorithm to the study of the self-assembly of a fully atomistic model of experimental interest. We study the 2D self-assembly of a rigid organic molecule (1,4-benzene-dicarboxylic acid or TPA), comparing the AB results with Monte Carlo (MC) and MC simulated annealing, a technique traditionally used to solve the global minimization problem. The AB algorithm gives a lower energy configuration in the same simulation time than both of the MC simulation techniques. We also show how the AB algorithm can be used as a part of the protocol to calculate the phase diagram with less computational effort than standard techniques.

Challenges in protein folding simulations: Timescale, representation, and analysis

Experimental studies of protein folding processes are frequently hampered by the fact that only low resolution structural data can be obtained with sufficient temporal resolution. Molecular dynamics simulations offer a complementary approach, providing extremely high resolution spatial and temporal data on folding processes. The effectiveness of such simulations is currently hampered by continuing questions regarding the ability of molecular dynamics force fields to reproduce the true potential energy surfaces of proteins, and ongoing difficulties with obtaining sufficient sampling to meaningfully comment on folding mechanisms. We review recent progress in the simulation of three common model systems for protein folding, and discuss how recent advances in technology and theory are allowing protein folding simulations to address their current shortcomings.

Biomolecular Structure and Modeling: Historical Perspective

physics, chemistry, and biology have been connected by a web of causal explanation organized by induction-based theories that telescope into one another. … Thus, quantum theory underlies atomic physics, which is the foundation of reagent chemistry and its specialized offshoot biochemistry, which interlock with molecular biology — essentially, the chemistry of organic macromolecules — and hence, through successively higher levels of organization, cellular, organismic, and evolutionary biology. … Such is the unifying and highly productive understanding of the world that has evolved in the natural sciences.

Molecular dynamics simulations of N-methylacetamide (NMA) in water by the ABEEM/MM model

The N-methylacetamide (NMA) is a very interesting kind of compound and often serves as a model of the peptide bond. The interaction between NMA and water provides a convenient prototype for the solvation of peptides in aqueous solutions. We have carried out molecular dynamics (MD) simulations of a NMA molecule in water under 1 atm and 298 K. The simulations make use of the newly developed NMA–water fluctuating charge ABEEM/MM potential model ( Yang, Z. Z.; Qian, P. J. Chem. Phys. 2006, 125, 064311 ), which is based on the combination of the atom-bond electronegativity equalization method (ABEEM) and molecular mechanics (MM). This model has been successfully applied to NMA–water gas clusters, NMA(H2O)n (n = 1–6), and accurately reproduced many static properties. For the NMA–water ABEEM/MM potential model, two characters must be emphasized in the simulations. Firstly, the model allows the charges in system to fluctuate, responding to the ambient environment. Secondly, for two major types of intermolecular hydrogen bonds, which are the hydrogen bond forming between the lone-pair electron on amide oxygen and the water hydrogen, and the one forming between the lone-pair electron on water oxygen and the amide hydrogen, we take special treatments in describing the electrostatic interaction by the use of the parameters klpO=,H and klpO–,HN–, respectively, which explicitly describe the short-range interaction of hydrogen bonds in the hydrogen bond interaction region. All sorts of properties have been studied in detail, such as, radial distribution function, energy distribution, ABEEM charge distribution and dipole moment, and so on. These simulation results show that the ABEEM/MM-based NMA–water potential model appears to be robust, giving the solution properties in excellent agreement with other dynamics simulations on similar systems.

VARICELLA: a variable-cell direct space method for structure determination from powder diffraction data

A direct space method for structure determination from powder diffraction data is proposed. Employing a hybrid Monte Carlo algorithm for generating the random conformations of a flexible molecular model, and sampling in a modified multicanonical statistical ensemble, it allows for variable cell parameters during an iterative search process. The acceptance-rejection criterion involves both a disagreement factor between the calculated and the experimental diffraction profiles and a modified crystal energy so that the space of tentative solutions can be widely explored while maintaining some physical meaningfulness of the proposals. Allowing the cell to be variable requires the zero shift to be treated as an optimizing parameter; this, in turn, requiring the disagreement factor to be based on the Fourier transform of the spectrum. The algorithm is presented in both a serial and a parallel version, the latter presenting several advantages, such as the possibility to probe different structures at a time while keeping them far from each other in the space defined by suitable order parameters. The method is built up and carefully tested by using, as a case study, a crystal of 3-ethyl 2,3-exo-disyndiotactic norbornene heptamer recently determined by single crystal x-ray diffraction techniques.

Force field bias in protein folding simulations

Long timescale (>1 micros) molecular dynamics simulations of protein folding offer a powerful tool for understanding the atomic-scale interactions that determine a protein's folding pathway and stabilize its native state. Unfortunately, when the simulated protein fails to fold, it is often unclear whether the failure is due to a deficiency in the underlying force fields or simply a lack of sufficient simulation time. We examine one such case, the human Pin1 WW domain, using the recently developed deactivated morphing method to calculate free energy differences between misfolded and folded states. We find that the force field we used favors the misfolded states, explaining the failure of the folding simulations. Possible further applications of deactivated morphing and implications for force field development are discussed.

Optimized molecular dynamics force fields applied to the helix-coil transition of polypeptides

Obtaining the correct balance of secondary structure propensities is a central priority in protein force-field development. Given that current force fields differ significantly in their alpha-helical propensities, a correction to match experimental results would be highly desirable. We have determined simple backbone energy corrections for two force fields to reproduce the fraction of helix measured in short peptides at 300 K. As validation, we show that the optimized force fields produce results in excellent agreement with nuclear magnetic resonance experiments for folded proteins and short peptides not used in the optimization. However, despite the agreement at ambient conditions, the dependence of the helix content on temperature is too weak, a problem shared with other force fields. A fit of the Lifson-Roig helix-coil theory shows that both the enthalpy and entropy of helix formation are too small: the helix extension parameter w agrees well with experiment, but its entropic and enthalpic components are both only about half the respective experimental estimates. Our structural and thermodynamic analyses point toward the physical origins of these shortcomings in current force fields, and suggest ways to address them in future force-field development.

Synthesis and antimitotic activity of novel 2-methoxyestradiol analogs--Part II

The syntheses and antimitotic activity of several novel 18a-homo-analogs of 2-methoxyestradiol are described. Structural modifications of the parent 2-methoxy-18a-homoestradiol include introduction of unsaturation in the D-ring and methylation of the 17-OH. Of seven analogs synthesized, one has demonstrated superior biological activities compared to 2-methoxyestradiol. The relationship between biological activity and the conformational preference of the 13-ethyl group as determined by computational analysis is discussed.