Calculations of biomolecules

Since the discovery of the nuclear Overhauser effect (NOE, see previous section) [4, 5] and scalar coupling constants [36, 37] decades ago, NMR-derived structure calculations of biomolecules largely depended on the measurement of these two parameters [38]. Recently it became possible to use cross-correlated relaxation (CCR) to directly measure angles between bond vectors [39] (see also Chapt 7). In addition, residual dipolar couplings of weakly aligned molecules were discovered to measure the orientation of bond vectors relative to the alignment tensor (see Sect 16.5). Measurement of cross-correlated relaxation was described experimentally earlier for homonuclear cases [40, 41] and is widely used in solid-state NMR [42 14]. 

Solvation energies for other multipoles inside a spherical cavity, including corrections due to salt effects, can be found, for example in Ref. 29. Analytical solutions of the Poisson equation for some other cavities, such as ellipse or cylinder, are also known [2] but are of little use in solvation calculations of biomolecules. For cavities of general shape only numerical solution of the Poisson and Poisson-Boltzmann equations is possible. There are two well-established approaches to the numerical solution of these equations the finite difference and the finite element methods. 

The predictive power of computer simulation relies heavily on the accuracy offeree field and the efficiency of phase space sampling. The broader and comprehensive applications of biomolecular simulation have highlighted the limitations of the existing force fields, and there is an urgent need to develop next generation force field that includes electrostatic polarization for biomolecules. The desired polarizable or polarized force fields should be ideally based on quantum mechanical calculations of biomolecules, which is a challenging task for computational chemists. 

Quantum chemical calculations of biomolecules and their environment... 

Many computational chemistry techniques are extremely computer-intensive. Depending on the type of calculation desired, it could take anywhere from seconds to weeks to do a single calculation. There are many calculations, such as ah initio analysis of biomolecules, that cannot be done on the largest computers in existence. Likewise, calculations can take very large amounts of computer memory and hard disk space. In order to complete work in a reasonable amount of time, it is necessary to understand what factors contribute to the computer resource requirements. Ideally, the user should be able to predict in advance how much computing power will be needed. 

Empirical energy functions can fulfill the demands required by computational studies of biochemical and biophysical systems. The mathematical equations in empirical energy functions include relatively simple terms to describe the physical interactions that dictate the structure and dynamic properties of biological molecules. In addition, empirical force fields use atomistic models, in which atoms are the smallest particles in the system rather than the electrons and nuclei used in quantum mechanics. These two simplifications allow for the computational speed required to perform the required number of energy calculations on biomolecules in their environments to be attained, and, more important, via the use of properly optimized parameters in the mathematical models the required chemical accuracy can be achieved. The use of empirical energy functions was initially applied to small organic molecules, where it was referred to as molecular mechanics [4], and more recently to biological systems [2,3]. 

As discussed in many previous studies of biomolecules, the treatment of electrostatic interactions is an important issue [69, 70, 84], What is less widely appreciated in the QM/MM community, however, is that a balanced treatment of QM-MM electrostatics and MM-MM electrostatics is also an important issue. In many implementations, QM-MM electrostatic interactions are treated without any cut-off, in part because the computational cost is often negligible compared to the QM calculation itself. For MM-MM interactions, however, a cut-off scheme is often used, especially for finite-sphere type of boundary conditions. This imbalanced electrostatic treatment may cause over-polarization of the MM region, as was first discussed in the context of classical simulations with different cut-off values applied to solute-solvent and solvent-solvent interactions [85], For QM/MM simulations with only energy minimizations, the effect of over-polarization may not be large, which is perhaps why the issue has not been emphasized in the past. As MD simulations with QM/MM potential becomes more prevalent, this issue should be emphasized. 

The biomolecule sensitivity is defined as the magnitude of the resonant mode shift for a given biomolecule surface density, i.e. fi/L/ap. Upon determining the biomolecule sensitivity, the detection limit (DL) can be calculated. The DL is the minimum surface density of biomolecules that can possibly be measured by the sensing device ... 

Note The calculation of relative molecular mass, Mr, of organic molecules exceeding 2000 u is significantly influenced by the basis it is performed on. Both the atomic weights of the constituent elements and the natural variations in isotopic abundance contribute to the differences between monoisotopic- and relative atomic mass-based values. In addition, they tend to characteristically differ between major classes of biomolecules. This is primarily because of molar carbon content, e.g., the difference between polypeptides and nucleic acids is about 4 u at Mr = 25,000 u. Considering terrestrial sources alone, variations in the isotopic abundance of carbon lead to differences of about 10-25 ppm in Mr which is significant with respect to mass measurement accuracy in the region up to several 10 u. [41]... 

Dithiins are the only biomolecules found in nature that are formally nonaromatic living organisms tend to avoid synthesizing antiaromatic compounds because of their thermodynamic and kinetic instability. Ab initio calculations of... 

The development of molecular mechanics, which incorporates quanmm mechanical data into a simplified mathematical framework derived from the classical equations of motion to permit reasonable calculations on biomolecules of large size. 

The ground-state wave function of cytosine has been calculated by practically all the semiempirical as well as nonempirical methods. Here, we shall discuss the application of these methods to interpret the experimental quantities that can. be calculated from the molecular orbitals of cytosines and are related to the distribution of electron densities in the molecules. The simplest v-HMO method yielded a great mass of useful information concerning the structure and the properties of biological molecules including cytosines. The reader is referred to the book1 Quantum Biochemistry for the application of this method to interpret the physicochemical properties of biomolecules. Here we will restrict our attention to the results of the v-SCF MO and the all-valence or all-electron treatments of cytosines. 

This volume was developed from the second international symposium on NMR chemical shifts. This meeting was organized by the editors at the 216th National Meeting of the American Chemical Society, Boston, Massachusetts, August 23-26,1998, and was cosponsored by the ACS Divisions of Computers in Chemistry and Physical Chemistry. The symposium included four extended lectures by Jameson, Ando, Oldfield, and Nicholas, which are included in this volume as Chapters 1-4, respectively. These lectures provide a convenient review of the current state of the art in the calculation of the NMR chemical shielding (Jameson, Chapter 1), and its application to different areas of chemistry polymers (Ando), biomolecules (Oldfield), and catalysis (Nicholas). 

On the other hand, there are several clear perspectives for future improvements and extensions of COSMO-RS. One of the most obvious perspectives is the improvement of the underlying quantum chemical methods. While density functional theory appears to have reached its limit regarding the quality of the electrostatics, and hence of the COSMO polarization charge densities, there will be an increase in the availability of higher correlated ab initio methods like coupled cluster calculations at affordable computational cost. Quantum chemical calculation of local polarizability and eventually of suitable descriptors for dispersion forces should provide additional information about the strength of local surface interactions and can be used to improve the various surface interaction functionals. At the other end, the quantum chemical COSMO calculations for larger biomolecules and enzymes, which have just become available at reasonable... 

Computational methods are increasingly applied to calculate pXa values of biomolecules. The data obtained for the gas phase are subsequently tried to be converted into values also relevant for solution conditions. These efforts, frequently made to understand processes occurring in catalytic centers of enzymes (213), have occasionally also included nucleotides (214). 

For MM+ (energy calculations of small biomolecules or ligands) Choose either Bond dipoles or Atomic charges (assigned via Builds Set charge) for use in the calculations of nonbounded Electrostatic interactions. Select None (calculate all nonbonded interactions recommended for small molecules), Switched or Shifted for Cutoffs (for large molecules). 

---