Theoretical gain optimization studies in 10.6 µm CO2-N2 gasdynamic lasers. IV. Further results of parametric study

Based on a method proposed by Reddy and Shanmugasundaram, similar solutions have been obtained for the steady inviscid quasi-one-dimensional nonreacting flow in the supersonic nozzle of CO2-N2-H2O and CO2-N2-He gasdynamic laser systems. Instead of using the correlations of a nonsimilar function NS for pure N2 gas, as is done in previous publications, the NS correlations are computed here for the actual gas mixtures used in the gasdynamic lasers. Optimum small-signal optical gain and the corresponding optimum values of the operating parameters like reservoir pressure and temperature and nozzle area ratio are computed using these correlations. The present results are compared with the previous results and the main differences are discussed. Journal of Applied Physics is copyrighted by The American Institute of Physics.

Theoretical gain optimization in CO2-N2-H2 gasdynamic lasers with two-dimensional wedge nozzles

Attention is given to the results of optimization studies with a 16-micron CO2-N2-H2 GDL employing two-dimensional wedge nozzles. The optimum value of the achievable gain reaches 12.7 percent/cm on the P(15) line for a 30:50:20 percent respective apportionment of the aforementioned gases. The corresponding optimum values for reservoir pressure and area ratio are computed as functions of reservoir temperature, and presented graphically.

Conformational studies on cyclic dipeptides

An analysis of 11 crystal structures of cyclic dipeptides so far reported in the literature is made, with main reference to the internal parameters of these molecules. Preferred conformations of the side chains of cyclic dipeptides with different α-amino acid residues have been studied by classical energy calculations. The possible conformations of the DKP ring are also studied. The significance of the non-bonded interaction in deciding the pathway for conformational change has also been investigated. The agreement between theoretical results and experimental observations is quite good, both with respect to the conformation of these molecules as well as the enthalpy difference as estimated from n.m.r. studies between different conformers.

Free energy calculations of glycosaminoglycan - Protein interactions

Glycosaminoglycans (GAGs) are complex highly charged linear polysaccharides that have a variety of roles in biological processes. We report the first use of molecular dynamics (MD) free energy calculations using the MM/PBSA method to investigate the binding of GAGs to protein molecules, namely the platelet endothelial cell adhesion molecule 1 (PECAM-1) and annexin A2. Calculations of the free energy of the binding of heparin fragments of different sizes reveal the existence of a region of low GAG-binding affinity in domains 5-6 of PECAM-1 and a region of high affinity in domains 2-3, consistent with experimental data and ligand-protein docking studies. A conformational hinge movement between domains 2 and 3 was observed, which allows the binding of heparin fragments of increasing size (pentasaccharides to octasaccharides) with an increasingly higher binding affinity. Similar simulations of the binding of a heparin fragment to annexin A2 reveal the optimization of electrostatic and hydrogen bonding interactions with the protein and protein-bound calcium ions. In general, these free energy calculations reveal that the binding of heparin to protein surfaces is dominated by strong electrostatic interactions for longer fragments, with equally important contributions from van der Waals interactions and vibrational entropy changes, against a large unfavorable desolvation penalty due to the high charge density of these molecules.

Can current force fields reproduce ring puckering in 2-O-sulfo-α-l-iduronic acid? A molecular dynamics simulation study

The monosaccharide 2-O-sulfo-α-l-iduronic acid (IdoA2S) is one of the major components of glycosaminoglycans. The ability of molecular mechanics force fields to reproduce ring-puckering conformational equilibrium is important for the successful prediction of the free energies of interaction of these carbohydrates with proteins. Here we report unconstrained molecular dynamics simulations of IdoA2S monosaccharide that were carried out to investigate the ability of commonly used force fields to reproduce its ring conformational flexibility in aqueous solution. In particular, the distribution of ring conformer populations of IdoA2S was determined. The GROMOS96 force field with the SPC/E water potential can predict successfully the dominant skew-boat to chair conformational transition of the IdoA2S monosaccharide in aqueous solution. On the other hand, the GLYCAM06 force field with the TIP3P water potential sampled transitional conformations between the boat and chair forms. Simulations using the GROMOS96 force field showed no pseudorotational equilibrium fluctuations and hence no inter-conversion between the boat and twist boat ring conformers. Calculations of theoretical proton NMR coupling constants showed that the GROMOS96 force field can predict the skew-boat to chair conformational ratio in good agreement with the experiment, whereas GLYCAM06 shows worse agreement. The omega rotamer distribution about the C5–C6 bond was predicted by both force fields to have torsions around 10°, 190°, and 360°.

Free energy calculations of the interactions of c-Jun-based synthetic peptides with the c-Fos protein

The c-Fos–c-Jun complex forms the activator protein 1 transcription factor, a therapeutic target in the treatment of cancer. Various synthetic peptides have been designed to try to selectively disrupt the interaction between c-Fos and c-Jun at its leucine zipper domain. To evaluate the binding affinity between these synthetic peptides and c-Fos, polarizable and nonpolarizable molecular dynamics (MD) simulations were conducted, and the resulting conformations were analyzed using the molecular mechanics generalized Born surface area (MM/GBSA) method to compute free energies of binding. In contrast to empirical and semiempirical approaches, the estimation of free energies of binding using a combination of MD simulations and the MM/GBSA approach takes into account dynamical properties such as conformational changes, as well as solvation effects and hydrophobic and hydrophilic interactions. The predicted binding affinities of the series of c-Jun-based peptides targeting the c-Fos peptide show good correlation with experimental melting temperatures. This provides the basis for the rational design of peptides based on internal, van der Waals, and electrostatic interactions.

Calculations of the free energy of interaction of the c-Fos−c-Jun coiled coil: Effects of the solvation model and the inclusion of polarization effects

The leucine zipper region of activator protein-1 (AP-1) comprises the c-Jun and c-Fos proteins and constitutes a well-known coiled coil protein−protein interaction motif. We have used molecular dynamics (MD) simulations in conjunction with the molecular mechanics/Poisson−Boltzmann generalized-Born surface area [MM/PB(GB)SA] methods to predict the free energy of interaction of these proteins. In particular, the influence of the choice of solvation model, protein force field, and water potential on the stability and dynamic properties of the c-Fos−c-Jun complex were investigated. Use of the AMBER polarizable force field ff02 in combination with the polarizable POL3 water potential was found to result in increased stability of the c-Fos−c-Jun complex. MM/PB(GB)SA calculations revealed that MD simulations using the POL3 water potential give the lowest predicted free energies of interaction compared to other nonpolarizable water potentials. In addition, the calculated absolute free energy of binding was predicted to be closest to the experimental value using the MM/GBSA method with independent MD simulation trajectories using the POL3 water potential and the polarizable ff02 force field, while all other binding affinities were overestimated.

Theoretical Studies on Coupled Electron and Proton Transfer in Cytochrome c Oxidase

The complexity of life is based on an effective energy transduction machinery, which has evolved during the last 3.5 billion years. In aerobic life, the utilization of the high oxidizing potential of molecular oxygen powers this machinery. Oxygen is safely reduced by a membrane bound enzyme, cytochrome c oxidase (CcO), to produce an electrochemical proton gradient over the mitochondrial or bacterial membrane. This gradient is used for energy-requiring reactions such as synthesis of ATP by F0F1-ATPase and active transport. In this thesis, the molecular mechanism by which CcO couples the oxygen reduction chemistry to proton-pumping has been studied by theoretical computer simulations. By building both classical and quantum mechanical model systems based on the X-ray structure of CcO from Bos taurus, the dynamics and energetics of the system were studied in different intermediate states of the enzyme. As a result of this work, a mechanism was suggested by which CcO can prevent protons from leaking backwards in proton-pumping. The use and activation of two proton conducting channels were also enlightened together with a mechanism by which CcO sorts the chemical protons from pumped protons. The latter problem is referred to as the gating mechanism of CcO, and has remained a challenge in the bioenergetics field for more than three decades. Furthermore, a new method for deriving charge parameters for classical simulations of complex metalloenzymes was developed.

Food aid and food production: A theoretical analysis

It has been claimed that food aid leads to permanent dependency as it depresses domestic food prices and thus farmers find it profitable to take land out of food production and into more lucrative activities. This paper develops two alternative scenarios under which the hypothesis about the damaging effect of food aid may not be true. Under the first scenario, it is argued that food production in developing countries is often low due to unfavourable trade policies and if food aid is tied to the removal of bias against the agricultural sector, food aid will not have any disincentive effect on food production. The second exercise argues that the revenue raised by the recipient government by selling aid could be used for R and D in agricultural production.

Which One Is Preferred: Myers-Saito Cyclization of Ene-Yne-Allene or Garratt-Braverman Cyclization of Conjugated Bisallenic Sulfone? A Theoretical and Experimental Study

A competitive scenario between Myers-Saito (MS) and Garraff-Braverman (GB) cyclization has been created in a molecule. High-level computations indicate a preference for GB over MS cyclization. The activation energies for the rate-determining steps of the GB and MS cyclizations were found to be the same (24.4 kcal/mol) at the B3LYP/6-31G* level of theory; thus, from the kinetic point of view, both reactions are feasible. However, the main biradical intermediate GB2 of the GB reaction is 6.2 kcal/mol lower in energy than the biradical MS2, which is the main intermediate of MS reaction, so GB cyclization is thermodynamically favored over MS cyclization. To verify the prediction by computational techniques, bisenediynyl sulfones 1-4 and bisenediynyl sulfoxide 17 were synthesized. Under basic conditions, these molecules isomerized to a system possessing both the ene-yne-allene and the bisallenic sulfone. The isolation of only one product, identified as the corresponding naphthalene- or benzene-fused sulfone 8-11, indicated the occurrence of GB cyclization as the sole reaction pathway. No product corresponding to the MS cyclization pathway could be isolated. Though the theoretical prediction showed a preference for the GB pathway over the MS pathway, the exclusive preference for GB over MS cyclization is very striking. Further analysis showed that the intramolecular self-quenching nature of the GB pathway may play an important role in the complete preference for this reaction. Apart from the mechanistic studies, these sulfones showed DNA cleavage activity that had an inverse relation with the reactivity order. Our findings are important for the design of artificial DNA-cleaving agents.

Simple theoretical analysis of the Einstein's photoemission from quantum confined superlattices

In this paper, we study the Einstein's photoemission from III-V, II-VI, IV-VI and HgTe/CdTe quantum well superlattices (QWSLs) with graded interfaces and quantum well effective mass superlattices in the presence of a quantizing magnetic field on the basis of newly formulated dispersion relations in the respective cases. Besides, the same has been studied from the afore-mentioned quantum dot superlattices and it appears that the photoemission oscillates with increasing carrier degeneracy and quantizing magnetic field in different manners. In addition, the photoemission oscillates with film thickness and increasing photon energy in quantum steps together with the fact that the solution of the Boltzmann transport equation will introduce new physical ideas and new experimental findings under different external conditions. The influence of band structure is apparent from all the figures and we have suggested three applications of the analyses of this paper in the fields of superlattices and microstructures.

Cu-II-Azide Polymers of Cu-3 and Cu-6 Building Units: Synthesis,Structures, and Magnetic Exchange Mechanism

Two new neutral copper-azido polymers [Cu-3(N-3)(6)(tmen)(2)](n)(1)and [Cu-6(N-3)(12)(deen)(2)](n) (2) [tmen = N,N,N, N-tetramethylethylenediamine and deen = N,N-diethylethylenediamine] have been synthesized by using lower molar equivalents of the chelating diamine ligands with Cu(NO3)(2)center dot 3H(2)O and an excess of NaN3. The single crystal X-ray structure shows that in the basic unit of the 1D complex 1, the three Cu-II ions are linked by double end-on azido bridges with Cu-N-EO-Cu angles on both sides of the magnetic exchange critical angle of 108 degrees. Complex 2 is a 3D framework of a basic u-6 cluster. Cryomagnetic susceptibility measurements over a wide range of temperature exhibit dominant ferromagnetic behavior in both the complexes. Density functional theory calculations (B3LYP functional) have been performed on the trinuclear unit to provide a qualitative theoretical interpretation of the overall ferromagnetic behavior shown by the complex 1.

Type II beta-turn conformation of pivaloyl-L-prolyl-a-aminoiso-butyryl-N-methylamide: Theoretical, spectroscopic and X-ray studies

Pivaloyl-L-Pro-Aib-N-methylamihdaes been shown to possess one intramolecular hydrogen bond in (CD&SO solution, by 'H-nmr methods, suggesting the existence of p-turns, with Pro-Aib as the corner residues. Theoretical conformational analysis suggests that Type II P-turn conformations are about 2 kcal mol-' more stable than Type 111 structures. A crystallographic study has established the Type I1 /%turn in the solid state. The molecule crystallizes in the space group P21 with a = 5.865 8, b = 11.421 A, c = 12.966 A, /3 = 97.55", and 2 = 2. The structure has been refined to a final R value of 0.061. The Type I1 p-turn conformation is stabilized by an intramolecular 4 - 1 hydrogen bond between the methylamide NH and the pivaloyl CO group. The conformational angles are @pro= -57.8", $pro = 139.3', @Aib = 61.4', and $Ajb = 25.1'. The Type 11 /%turn conformation for Pro-Aib in this peptide is compared with the Type I11 structures observed for the same segment in larger peptides.